Fortran 95中如何正确设置随机数种子生成可复现的高斯白噪声序列？

问题根源 #

你之前把随机数种子初始化放在box_muller_mod子程序内，导致每次调用该子程序时都会重置随机数生成器的状态，全程生成完全相同的随机数序列——因为每次调用都从同一个初始种子开始生成。

实现步骤 #

1. 把种子初始化移到主程序开头，仅执行一次 #

随机数生成器的状态是全局的，只需在程序启动时初始化一次种子，后续所有随机数调用都会基于这个初始状态持续推进，不会重复。

2. Box-Muller子程序仅负责高斯噪声生成，不处理种子 #

子程序只需要调用Fortran内置的随机数生成函数（random_number）获取均匀分布随机数，再通过Box-Muller公式转换为高斯分布即可。

3. 支持手动指定种子，方便更换序列 #

通过手动设置种子数组，可以固定随机数序列；更换种子值就能得到不同的噪声序列。

示例代码 #

主程序（初始化种子+调用Box-Muller） #

program brownian_dynamics
    use box_muller_mod
    implicit none
    integer :: nsize, i, n_steps
    integer, allocatable :: seed(:)
    real :: noise

    ! 设置随机数种子（仅执行一次）
    call random_seed(size=nsize)
    allocate(seed(nsize))
    seed = 12345  ! 更换此值可得到不同随机序列
    call random_seed(put=seed)
    deallocate(seed)

    n_steps = 1000  ! 时间步数量
    do i = 1, n_steps
        ! 调用Box-Muller生成高斯白噪声
        noise = generate_gaussian_noise()
        ! 后续布朗动力学计算逻辑...
        print *, "Step ", i, " Noise: ", noise
    end do

end program brownian_dynamics

Box-Muller模块（仅负责噪声生成） #

module box_muller_mod
    implicit none
    private
    public :: generate_gaussian_noise
    real :: z1, z2
    logical :: has_spare = .false.  ! 缓存一个生成的高斯数，提高效率

contains

    function generate_gaussian_noise() result(gaussian)
        real :: gaussian, u1, u2, r
        if (has_spare) then
            ! 使用上一次生成的缓存值
            gaussian = z2
            has_spare = .false.
        else
            ! 生成两个(0,1)均匀分布随机数
            do
                call random_number(u1)
                call random_number(u2)
                u1 = 2.0*u1 - 1.0
                u2 = 2.0*u2 - 1.0
                r = u1**2 + u2**2
                if (r > 0.0 .and. r < 1.0) exit  ! 确保在单位圆内
            end do
            r = sqrt(-2.0*log(r)/r)
            z1 = u1*r
            z2 = u2*r
            gaussian = z1
            has_spare = .true.
        end if
    end function generate_gaussian_noise

end module box_muller_mod

关键说明 #

• random_seed是Fortran内置的随机数种子控制子程序，通过put参数传入自定义种子数组，可固定随机序列；若不指定put，则系统会自动生成随机种子（不同运行可能不同）。
• Box-Muller算法会一次性生成两个高斯分布随机数，用has_spare变量缓存其中一个，避免重复计算，提升效率。
• 确保主程序中种子初始化仅执行一次，所有随机数调用（包括子程序内的）都会共享同一个全局随机状态。

内容的提问来源于stack exchange，提问作者QuestionTheAnswer