24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

[交流] 【求助】小世界网络的模拟程序

自己用fortran写了一个程序，但运行速度太慢，
如果哪位有，请发到wqfscience@yahoo.com.cn
不胜感激。

[ Last edited by lvzhu2007 on 2008-9-18 at 15:42 ]

回复此楼

» 猜你喜欢

1楼 2008-09-18 15:33:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

huiy_chen

应助: (幼儿园)
在线:
虫号: 631108

★
wuli8(金币+1,VIP+0):鼓励新虫！(*^__^*) ……

我也正在编写关于小世界模型的fortran程序，可是得出来的结总是不对，我的方向是非线性动力学，能否把你的程序给我发下，huiy_chen@stu.snnu.edu.cn，不胜感激！我的代码如下：
program main
implicit none
c ******************************************************
c    creat a small-world network
c    ns is the number of edges
c    n is the number of vertex
c    kh is the nearest edges that connected to vertex i
c ******************************************************
   integer i,j,ns,n,jj,kh,k
real*8 p
   integer G(1:10,1:10)
real*8 R(1:100)
ns=0
write(*,*)"n="
read(*,*) n
write(*,*)"kh="
read(*,*) kh
write(*,*)"p="
read(*,*)p
do j=1,n
   do i=1,n
   G(i,j)=0
   end do
end do

do i=1,n
   do jj=1,kh
      j=i+jj
      if(j.gt.n)j=j-n
      G(i,j)=1
G(j,i)=1
ns=ns+1

end do
   end do
c******************************************************************
c    add edges randomly,ns denotes the degree of vertex
c******************************************************************
do i=1,n-kh-1
   do j=i+kh+1,n
      if(G(i,j).eq.1.d0)  cycle
         call random_seed()
         do k=1,1000
         call random_number(R(k))
         if(R(k).le.p)then
            G(i,j)=1
            G(j,i)=1
            ns=ns+1
         write(*,*) ns
         end if

      end do
   end do
end do
end