光強分布MATLAB

1、光強分布MATLAB作者:日期:、高斯光束模擬的焦點附近光場分布for n=1:10I=0.1*n;a=sqrt(1./I-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:10I=0.1*n;a=sqrt(1./I-1);z=0:0.0001:a;r1=sqrt(0.5.*(1+z.A2).*log(1./(I.*(1+z.A2);r2=-sqrt(0.5.*(1+z.A2).*log(1./(I.*(1+z.A2);plot(z,r1, -);hold o

2、nplot(z,r2, -);hold on end全圖：clear allfor n=1:10I=0.1*n;a=sqrt(1./I-1);z=-a:0.00001*a:a;r1=sqrt(0.5.*(1+z.A2).*log(1./(I.*(1+z.A2);r2=-sqrt(0.5.*(1+z.A2).*log(1./(I.*(1+z.A2);plot(z,r1, -); axis equal hold onplot(z,r2, -); axis equal hold on endxlabel(Z( /rd)ylabel(r( /m詠)title(焦點附近光場分布)帶空間坐標:clear

3、allfor n=1:10I=0.1*n;w0=200e-9;wl=780e-9;a=pi*w0A2*sqrt(1./I-1)/wl;z=-a:0.0001*a:a;r1=w0*sqrt(0.5.*(1+(wl/(pi*w0A2).*z.A2).*log(1./(I.*(1+(wl/(pi*w0A2).*z.A2);r2=-w0*sqrt(0.5.*(1+(wl/(pi*w0A2).*z.A2).*log(1./(I.*(1+(wl/(pi*w0A2).*z.A2); plot(z,r1, -); axis equal hold onplot(z,r2, -); axis equal hold

4、 on endxlabel(Z(入)/m);ylabel(r( 入)/m);title(焦點附近光場分布,)定4M城無場方和保點附近光場分布Z(A)-2 -3-2-101X Axis Title、焦點附近的三維光場分布三、高斯光束在傳播過程中，焦點處的光強分布1、高斯光束在光腰處的三維光強分布圖clear allN=200;w0=200e-9;%高斯光束的光腰半徑r=linspace(0,3*w0,N);eta=linspace(0,2*pi,N);rho,theta=meshgrid(r,eta);%生成極坐標網格x,y=pol2cart(theta,rho);%將極坐標網格轉化為直角坐標網

5、格I=exp(-2*rh0 A2./w0 A2);%高斯光束在歸一化的光強值surf(x,y,I);colormap hot;%colormap gray 線性灰度圖colormap hot 熱色顏色圖去掉這句命令為彩圖shading interp;xlabel(x /m);ylabel(y /m);zlabel(intensity /a.u.);axis(-3*w0 3*w0 -3*w0 3*w0 0 1)title(焦點附近三維光場分布)colorbar;n tsuvc-2、高斯光束在光腰處的二維光強分布圖clear allN=200;w0=200e-9;%高斯光束的光腰半徑x=linsp

6、ace(-3*w0,3*w0,N);y=linspace(-3*w0,3*w0,N);X=meshgrid(x,y);Y=meshgrid(y,x);Y=Y;R=sqrt(X,A2+Y.A2);I=exp(-2*RA2./w0 A2);imagesc(x,y,I,0 1); axis equalcolormap hot;%colormap gray 線性灰度圖 colormap hot 熱色顏色圖去掉這句命令為彩圖colorbar;xlabel(x /m);ylabel(y /m);axis(-3*w0 3*w0 -3*w0 3*w0)title(焦點附近三維光場分布)4點.網上_爆均分書6咽

7、Qm4x 10-3、高斯光束在傳播過程中，焦點處的三維光強分布clear allw0=200e-9;%高斯光束的光腰半徑wavelength=780e-9;ZR=pi*w0A2/wavelength;L=3*ZR;%高斯光束的光波長%高斯光束的瑞利長度%高斯光束的傳輸范圍N=200;z=linspace(-L,L,N);r=linspace(-4*w0,4*w0,N);pr,pz=meshgrid(r,z);wz=w0*sqrt(1+(wavelength*pz/pi/w0A2)-2);I=w0A2./wz.A2.*exp(-2*pr.A2./wz.A2);surf(pz,pr,I);shading interp;xlabel(Z(/m);ylabel(r(/rrK)zlabel(intensity /a.u.); axis(-3*ZR 3*ZR -Inf Inf -Inf Inf) title(焦點附近三維光場分布,)colorbar;colorm