蒙特卡洛方法模拟光声成像图像重建

该程序通过蒙特卡洛模拟光声成像的过程 ，能够反应光子的运动。

function output=iterat（inputstep） 
input=[000000];
x=input（1）;y=input（2）;z=input（3）;ux=input（4）;uy=input（5）;uz=input（6）; 
g=0.8; 
step=0.49383;
step=step*log10（rand（））/log10（2.71828）; 
if （step>5） 
    step=5; 
end 
l=step;%步长 
 
e1=rand（1）; 
e2=rand（1）; 
if（g==0） 
    costheta=2*e1-1; 
else 
    costheta=（1+g^2-（（1-g^2）/（1-g+2*g*e1））^2）/2/g; 
end 
fai=2*pi*e2; 
sinfai=sin（fai）; 
cosfai=cos（fai）; 
sintheta=（1-（costheta）^2）^0.5; 
if （abs（uz）<=0.99999） 
    uxx=（sintheta/（1-uz^2）^0.5）*（ux*uz*cosfai-uy*sinfai）+ux*costheta; 
    uyy=（sintheta/（1-uz^2）^0.5）*（uy*uz*cosfai+ux*sinfai）+uy*costheta; 
    uzz=-（1-uz^2）^0.5*sintheta*cosfai+uz*costheta; 
else 
    uxx=sintheta*sinfai 
    uyy=sintheta*sinfai 
    uzz=sign（uz）*costheta 
end 
uuxx=uxx/（uxx^2+uyy^2+uzz^2）^0.5; 
uuyy=uyy/（uxx^2+uyy^2+uzz^2）^0.5; 
uuzz=uzz/（uxx^2+uyy^2+uzz^2）^0.5; 
xx=x+l*uuxx; 
yy=y+l*uuyy; 
zz=z+l*uuzz; 
output=[xx yy zz uuxx uuyy uuzz];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        965  2010-11-30 13:29  mc模拟光声成像\iterat.m

     文件       5495  2010-11-30 14:52  mc模拟光声成像\mc.m

     目录          0  2001-05-14 16:32  mc模拟光声成像

----------- ---------  ---------- -----  ----

                 6460                    3