[交流] matlab的蒙特卡洛模拟提高运算速率 已有1人参与

%蒙特卡洛模拟光子在水雾中的 传播 kext=2.4056; ks=2.4053; ka=0.0003; r=250*10^(-6); n0=19646104; c=pi*(r^2)*n0*ks; for m=10:10:250 L=m; s1=0; s2=0; z=zeros(1,10000); x=zeros(1,10000); y=zeros(1,10000); n=1; while(n>=1&&n<=10000)     if((z(1,n)>=L)&&((x(1,n))^2+(y(1,n))^2<2.88))                 s1=s1+1;         n=n+1;     else if((z(1,n)>=L)&&((x(1,n))^2+(y(1,n))^2>=2.88))         n=n+1;         else if(rand(1)<=(ka/kext))             n=n+1;             s2=s2+1;         else syms t si;      % t=solve(int(4.40189*10^(-8)+1.39435*10^(-8)*sin((pi/1.58425)*(si-2.33743)),si,0,t)-rand(1)*1.3792*10^(-7));         y1=4.40189*10^(-8)+1.39435*10^(-8)*sin(pi/1.58425*(si-2.33743));         R=rand(1)*1.379*10^(-7);         y2=int(y1,si,0,t)-R;         t=vpasolve(y2);                         rand1=-log(rand(1))/c;         rand2=rand(1)*2*pi;           rand3=t-pi/2;         z(1,n)=z(1,n)+rand1*cos(rand3);         x(1,n)=x(1,n)+rand1*cos(rand2)*sin(rand3);         y(1,n)=y(1,n)+rand1*sin(rand2)*sin(rand3);             n=n+0;         end     end     end end fprintf('Distance:%i\n',L); fprintf('The number of photons arriving at the probe plane:%i\n',s1); end 本程序需要对10000个光子进行模拟， 若先对1000个进行模拟，出来一个L对应的是s1耗时5个小时。 那改为10000个光子，出来所有L=10:10:250对应的s1就可能得几个月了 求问如何提高运行效率

回复此楼