| 查看: 628 | 回复: 5 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】小世界网络的模拟程序
|
|||
|
自己用fortran写了一个程序，但运行速度太慢， 如果哪位有，请发到wqfscience@yahoo.com.cn 不胜感激。 [ Last edited by lvzhu2007 on 2008-9-18 at 15:42 ] |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有25人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有8人回复
-
分享一下我之前已中青C的计划书的filecode
已经有4人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有3人回复
-
为什么网上很多人说本周 12号出结果
已经有6人回复
-
听说今天filecode变了
已经有49人回复
-
帮忙看看fileCode
已经有12人回复
-
FileCode能看出啥？
已经有17人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有6人回复
-
filecode
已经有9人回复
★
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
2楼2008-09-18 15:41:26
3楼2008-09-18 18:39:13
4楼2008-09-18 19:15:58
