| 查看: 520 | 回复: 5 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】小世界网络的模拟程序
|
|||
|
自己用fortran写了一个程序，但运行速度太慢， 如果哪位有，请发到wqfscience@yahoo.com.cn 不胜感激。 [ Last edited by lvzhu2007 on 2008-9-18 at 15:42 ] |
回复此楼» 猜你喜欢
过年走亲戚时感受到了所开私家车的鄙视链
已经有10人回复
-
今年春晚有几个节目很不错，点赞！
已经有12人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有13人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有12人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复
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