光强分布MATLAB.docx

1、光强分布MATLAB、高斯光束模拟的焦点附近光场分布for n=1:101=0.1* n;a=sqrt(1./l-1);z=0:0.0001:a;r=sqrt(0.5.*(1+z.A2).*log(1./(I.*(1+ z. A2);plot(z,r, -);hold onend右半边：for n=1:101=0.1* n;a=sqrt(1./l-1);z=0:0.0001:a;r仁 sqrt(0.5.*(1+z.A2).*log(1./(I.*(1+ z.八2); r2=-sqrt(0.5.*(1+z.A2).*log(1./(I.*(1+ z.八2); plot(z,r1, -);hold

2、 onplot(z,r2, -);hold onend全图：clear allfor n=1:101=0.1* n;a=sqrt(1./l-1);z=-a:0.00001*a:a;r仁 sqrt(0.5.*(1+z.A2).*log(1./(I.*(1+ z. A2);r2=-sqrt(0.5.*(1+z.A2).*log(1./(I.*(1+ z.八2);plot(z,r1, -); axis equalhold onplot(z,r2, -); axis equalhold onendxlabel(Z( R/m);ylabel(r( R/m);title(焦点附近光场分布)带空间坐标:cl

3、ear allfor n=1:101=0.1* n;w0=200e-9;wl=780e-9;a=pi*w0A2*sqrt(1./I-1)/wl;z=-a:0.0001*a:a;r仁w0*sqrt(0.5.*(1+(wl/(pi*w0A2).*z.A2).*log(1./(l*(1+(wl/(pi*w0A2).*zA2);r2=-w0*sqrt(0.5.*(1+(wl/(pi*w0A2).*z.A2).*log(1./(l*(1+(wl/(pi*w0A2).*z.A2);plot(z,r1, -); axis equalhold onplot(z,r2, -); axis equalhold o

4、nendxlabel(Z( R/m);ylabel(r( R/m);title(焦点附近光场分布)ZAI近光塌为舟0 1 2cm COXA Y-2-10123X Axis Title、焦点附近的三维光场分布三、高斯光束在传播过程中，焦点处的光强分布1、高斯光束在光腰处的三维光强分布图clear allN=200;w0=200e-9; %高斯光束的光腰半径r=li nspace(0,3*w0,N);eta=li nspace(0,2*pi,N);rho,theta=meshgrid(r,eta); % 生成极坐标网格x,y=pol2cart(theta,rho); %将极坐标网格转化为直角坐标网

5、格I=exp(-2*rho.A2./w0.A2); %高斯光束在归一化的光强值surf(x,y,I);colormap hot; %colormap gray 线性灰度图 colormaphot热色颜色图 去掉这句命令为彩图shad ing in terp;xlabel(x /m);ylabel(y /m);zlabel(i nten sity /a.u.);axis(-3*w0 3*w0 -3*w0 3*w0 0 1)title(焦点附近三维光场分布)colorbar;2、高斯光束在光腰处的二维光强分布图clear allx=li nspace(-3*w0,3*w0,N);y=li nspa

6、ce(-3*w0,3*w0,N);X=meshgrid(x,y);Y=meshgrid(y,x);Y=Y;R=sqrt(X.A2+Y.A2);I=exp(-2*R.A2./w0.A2);imagesc(x,y,l,O 1); axis equalcolormap hot; %colormap gray 线性灰度图 colormaphot热色颜色图 去掉这句命令为彩图colorbar;xlabel(x /m);ylabel(y /m);axis(-3*w0 3*w0 -3*w0 3*w0)title(焦点附近三维光场分布)击 *2 0 2 A E3、高斯光束在传播过程中，焦点处的三维光强分布%高

7、斯光束的光腰半径%高斯光束的光波长%高斯光束的瑞利长度%高斯光束的传输范围clear allw0=200e-9;wavele ngth=780e-9;ZR=pi*w0A2/wavele ngth;L=3*ZR;N=200;z=li nspace(-L,L,N);r=li nspace(-4*w0,4*w0,N);pr,pz=meshgrid(r,z);wz=w0*sqrt(1+(wavele ngth*pz/pi/w0A2).A2);l=w0A2./wz.A2.*exp(-2*pr.A2./wz.A2);surf(pz,pr,I);shad ing in terp;xlabel(Z(为 /m);ylabel(r( R /m);zlabel(i nten sity /a.u.);axis(-3*ZR 3*ZR -Inflnf -Inf Inf)title（焦点附近三维光场分布）colorbar;colormap hot; %colormap gray 线性灰度图 colormap hot 热色颜色图view(2)%显示二维视图oy oa o?0&05 i.0 J0.20 1近三维也馬労帘-4 -3 -2 -1 0 1 2 3 4Z(Artn *1(/I1O 議点削近三维也场労帝-4 -3 -2 -1 0 1 2 3 4Z(Artn