计算机模拟第21课.ppt免费

6.8 雨滴运动的计算机模拟; ntrial=10000; hmin=10; hmax=200; for ix=hmin:10:hmax stepsm=0; sigsq=0; x=0; y=0; z=ix; for jx=1:ntrial [x,y,nsteps]=rainrw(x,y,z); stepsm=stepsm+nsteps; sigsq=sigsq+x*x+y*y; end stepsm=stepsm/ntrial; sigma=sqrt(sigsq/ntrial); result=[ix,stepsm,sigma] end ;function [x,y,nsteps]=rainrw(x,y,z) for jx=1:10000 p=rand; if p=0.55 z=z-1; elseif p=0.6 z=z+1; elseif p=0.7 y=y+1; elseif p=0.8 y=y-1; elseif p=0.9 x=x+1;