伪随机数抽取随机数序列——随机数产生的改进

言深深

众所周知fortran内部内置一张随机数表，我们选取的时候就是从随机数表的某一项开始，往下顺次抽取需要的个数，但是往往系统内置的随机数相邻两个数产生的时候内在动力学相关性比较大，为此我们需要重新挑选随机数，于是我用随机数作为选取状态变量，随机从已经产生的m个伪随机数中选择n个数：
parameter(m=10000,n=1000,ns=100,k=1)
!m：产生伪随机数的序列长度
!n：产生新的随机数个数
!ns：构建随机位置序列的初始选择，可以随便写，保证每次产生的随机数不一样
!x数组存放的是原始伪随机数
!y数组存放的是新随机数
!r数组里面存放的是抽取随机数的标签
!k是随机数振幅

real x(m),y(n),r(n)
iseed=1
do i=1,m
x(i)=ran(iseed)
enddo

do i=ns,ns+n-1
r(i-ns+1)=mod(int(x(i)*m),n)
        if(r(i-ns+1).eq.0..and.i.ne.1)then
        r(i-ns+1)=int((r(i-ns+1)+r(1))/2)        !此处判断抽取的随机序号不为0
        endif
enddo

do i=1,n
        y(i)=x(r(i))*k
enddo

open(1,file="ran rand.txt")
do i=1,n
write(1,"(f20.5)") y(i)
enddo
close(1)
end

言深深

什么叫七天憋六个字，深深体会到了！没有课题，没有目的，没有方向的悲哀啊~~~~

传说中的谁

鄙视自己抢沙发的。下午笨笨也抽个时间出一个帖子了

尽头的尽头

支持深深啊

言深深

传说中的谁 发表于 2011-10-8 12:21

                               
登录/注册后可看大图

鄙视自己抢沙发的。下午笨笨也抽个时间出一个帖子了

不要威胁我啊，压力已经很大了！

言深深

尽头的尽头 发表于 2011-10-8 13:09

                               
登录/注册后可看大图

支持深深啊

友好合作

mofangbao

想法不错哟，原来fortran自带随机数表的啊，之前都不知道。貌似判断为0之后直接使用r(1)不可靠的 万一第一个数字就等于0...，呵呵,不过概率比较小啦，用另外生成的随机数做标签也可以试试

言深深

mofangbao 发表于 2011-10-8 14:16

                               
登录/注册后可看大图

想法不错哟，原来fortran自带随机数表的啊，之前都不知道。貌似判断为0之后直接使用r(1)不可靠的 万一第一个 ...

呵呵，我今天上午想半天，没解决掉这个bug，等我想到好办法的时候会修正一下的啦~~~~没想到还是被你看出来了哈哈

757084134

真的不错啊