隨機(jī)行走之計(jì)算機(jī)模擬——王慶(中國(guó)科學(xué)技術(shù)大學(xué))

1、RANDOM WALK王慶PB11204037物理化學(xué)研討第三小組12013/12/24 Tuesday一維晶格隨機(jī)行走模擬（C）#include #include #include #include #includevoid main() FILE*outfile; if(outfile=fopen(outfile.dat,w)=NULL) printf(cannot open this filen); exit(0); float a=0.0; float b=0.0; int x=0; int i=1; srand(int)time(NULL);while(i0.5) x=x+1; el

2、se x=x-1; printf(%d %dn,i,x); fprintf(outfile,%d %dn,i,x); i=i+1; fclose(outfile);22013/12/24 Tuesday一維隨機(jī)行走模擬1000個(gè)點(diǎn)32013/12/24 Tuesday4一維隨機(jī)行走模擬10000個(gè)點(diǎn)2013/12/24 Tuesday5一維隨機(jī)行走模擬100000個(gè)點(diǎn)2013/12/24 Tuesday#include#includeusing namespace std;int main() FILE*outfile; if(outfile=fopen(outfile.dat,w)=NULL

3、) printf(cannot open this filen); exit(0); int x=0,y=0,i=0; int step; cout起點(diǎn):(0,0)endl; fprintf(outfile,%d %d %dn,i,x,y); srand(time(NULL); for(i=1;i100;i+) step=rand()%4; if(step=0) x+; couti向前endl; cout位置(x,y)endl; fprintf(outfile,%d %d %dn,i,x,y); else if(step=1) x-; couti向后endl; cout位置(x,y)endl;

4、 fprintf(outfile,%d %d %dn,i,x,y); else if(step=2) y-; couti向左endl; cout位置(x,y)endl; fprintf(outfile,%d %d %dn,i,x,y); else if(step=3) y+; couti向右endl; cout位置(x,y)endl; fprintf(outfile,%d %d %dn,i,x,y); printf(*n); printf(以上數(shù)據(jù)是二維隨機(jī)行走%d步的結(jié)果！n,i); fclose(outfile);6二維晶格隨機(jī)行走模擬（C+）2013/12/24 Tuesday二維隨機(jī)行

5、走100步位置 (X,Y)與時(shí)間T關(guān)系圖72013/12/24 Tuesday8二維隨機(jī)行走100步距離R與時(shí)間T關(guān)系圖2013/12/24 Tuesday9二維隨機(jī)行走10000步位置 (X,Y)與時(shí)間T關(guān)系圖2013/12/24 Tuesday10二維隨機(jī)行走10000步距離R與時(shí)間T關(guān)系圖2013/12/24 Tuesday#include#include#includeusing namespace std;int main() FILE*outfile; if(outfile=fopen(outfile.dat,w)=NULL) printf(cannot open this fil

6、en); exit(0); int x=0,y=0,i=0; int step; long double r=0,max=0; cout起點(diǎn):(0,0)endl; fprintf(outfile,%d %d %dn,i,x,y); srand(time(NULL); for(i=1;i10E8;i+)/此處設(shè)置循環(huán)次數(shù) step=rand()%4; if(step=0) x+; couti向前endl; cout位置(x,y)endl; fprintf(outfile,%d %d %dn,i,x,y);11二維晶格沙漠中醉漢隨機(jī)行走模型 else if(step=1) x-; couti向后e

7、ndl; cout位置(x,y)endl; fprintf(outfile,%d %d %dn,i,x,y); else if(step=2) y-; couti向左endl; cout位置(x,y)endl; fprintf(outfile,%d %d %dn,i,x,y); else if(step=3) y+; couti向右endl; cout位置(x,y)max) max=r; if(max=40) break;/此處設(shè)置沙漠半徑 printf(最遠(yuǎn)距離是：%.2fn,max); if(max=40) printf(上天憐憫啊，我終于走出沙漠了！n);/此處設(shè)置沙漠半徑 else p

8、rintf(我是故意不走出沙漠的，哼！n); printf(* * * * * * * * * * * * * * * * * * *n); printf(以上數(shù)據(jù)是二維隨機(jī)行走%d步的結(jié)果！n,i); fclose(outfile);2013/12/24 Tuesday12醉漢沙漠隨機(jī)行走模型某次運(yùn)行結(jié)果預(yù)設(shè)（速度：1m/s，沙漠半徑是10000m，完全隨機(jī)晶格線路）走1000步，計(jì)算機(jī)運(yùn)行約4秒，最遠(yuǎn)距離23.09；走10000步，計(jì)算機(jī)運(yùn)行約28秒，最遠(yuǎn)距離127.88；走100000步，計(jì)算機(jī)運(yùn)行約270秒，最遠(yuǎn)距離345.04；走1000000步。2013/12/24 Tuesday

9、13二維晶格自規(guī)避隨機(jī)行走GSAW2013/12/24 Tuesday以(0,0)為起點(diǎn)的GSAW模擬情形?？梢?jiàn)該模型很容易就會(huì)“進(jìn)入死胡同”從而終止程序。14二維晶格自規(guī)避隨機(jī)行走GSAW2013/12/24 Tuesday15蜂窩狀晶格隨機(jī)行走示意圖2013/12/24 Tuesday#include #include #include #include #include#define PI 3.141592654main()FILE*fengwo; if(fengwo=fopen(fengwo.dat,w)=NULL) printf(cannot open this filen);exi

10、t(0); double x=0.0,y=0.0; int i=1;int step; 16srand(int)time(NULL); for(i=1;i100;i+)if(i%2=1) step=rand()%3; if(step=0) y=y+1.0; else if(step=1) x=x+0.5;y=y-cos(PI/6); else if(step=2) x=x-0.5;y=y-cos(PI/6); elsestep=rand()%3; if(step=0) y=y-1.0; else if(step=1) x=x+0.5;y=y+cos(PI/6); else if(step=2) x=x-0.5;y=y+cos(PI/6); printf(%f %fn,x,y); fprintf(fengwo,%f %