2楼: Originally posted by snoopyzhao at 2011-09-04 14:56:30:
每次都一样是什么意思？同一个程序，每次运行结果当然应该是一样的啊……

3楼: Originally posted by 贺仪 at 2011-09-04 15:07:58:
但我是个随机过程 就不应该一样的 。我用matlab得到的是不同的，但是用fortran得到的就是相同的。在取种子的时候seed取时间能不能精确到毫秒呢？要能到就可以了

4楼: Originally posted by snoopyzhao at 2011-09-04 15:29:13:
你的意思是，你在程序中取了很多随机数，然后，随机数是一样的吗？

另外，种子与时间有啥关系？你以时间作种吗？

给段代码看看吧……

6楼: Originally posted by asaka at 2011-09-05 02:44:38:
你把
call random_seed ()
放在do循环外头！

7楼: Originally posted by 贺仪 at 2011-09-05 09:24:07:
也试过了，最前面也放过作用不大，后来只有让程序等一段时间再运行 有点时间差也就好了。不过好像与变量声明也有关系
real：：rd
与
real rd
有什么区别呢？
后者是不可以，利用时间产生随机数的 加了call  ...

program test
   real :: r
   
!  call random_seed()
   call init_random_seed()
   
   do i=1,10
      call random_number(r)
      write(*,*) r
   end do
   
end program test

subroutine init_random_seed()
   integer :: i, n, clock
   integer, dimension(:), allocatable :: seed
   
   call random_seed(size = n)
   allocate(seed(n))
   
   call system_clock(count=clock)
   
   seed = clock + 37 * (/ (i - 1, i = 1, n) /)
   call random_seed(put = seed)
   
   deallocate(seed)
end subroutine init_random_seed