| 查看: 541 | 回复: 5 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】小世界网络的模拟程序
|
|||
|
自己用fortran写了一个程序，但运行速度太慢， 如果哪位有，请发到wqfscience@yahoo.com.cn 不胜感激。 [ Last edited by lvzhu2007 on 2008-9-18 at 15:42 ] |
回复此楼» 猜你喜欢
329求调剂
已经有13人回复
-
复试调剂
已经有4人回复
-
求调剂
已经有3人回复
-
材料调剂
已经有10人回复
-
292求调剂
已经有11人回复
-
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有7人回复
-
363求调剂
已经有5人回复
-
0703总分331求调剂
已经有4人回复
-
材料专硕(0856) 339分求调剂
已经有3人回复
-
求调剂
已经有3人回复
2楼2008-09-18 15:41:26
3楼2008-09-18 18:39:13
4楼2008-09-18 19:15:58
5楼2008-10-03 09:35:43
★
wuli8(金币+1,VIP+0):鼓励新虫！(*^__^*) ……
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 |
6楼2008-11-12 10:55:48
