光強(qiáng)分布MATLAB資料

一、高斯光束模擬的焦點(diǎn)附近光場(chǎng)分布forn=1:10I=0.1*n;a=sqrt(1./I-1);z=0:0.0001:a;r=sqrt(0.5.*(1+z.^2).*log(1./(I.*(1+z.^2))));plot(z,r,'-');holdonend右半邊：forn=1:10I=0.1*n;a=sqrt(1./I-1);z=0:0.0001:a;r1=sqrt(0.5.*(1+z.^2).*log(1./(I.*(1+z.^2))));r2=-sqrt(0.5.*(1+z.^2).*log(1./(I.*(1+z.^2))));plot(z,r1,'-');holdonplot(z,r2,'-');holdonend全圖：clearallforn=1:10I=0.1*n;a=sqrt(1./I-1);z=-a:0.00001*a:a;r1=sqrt(0.5.*(1+z.^2).*log(1./(I.*(1+z.^2))));r2=-sqrt(0.5.*(1+z.^2).*log(1./(I.*(1+z.^2))));plot(z,r1,'-');axisequalholdonplot(z,r2,'-');axisequalholdonendxlabel('Z(λ)/m');ylabel('r(λ)/m');title('焦點(diǎn)附近光場(chǎng)分布')二、焦點(diǎn)附近的三維光場(chǎng)分布三、高斯光束在傳播過(guò)程中，焦點(diǎn)處的光強(qiáng)分布1、高斯光束在光腰處的三維光強(qiáng)分布圖clearallN=200;w0=200e-9;%高斯光束的光腰半徑r=linspace(0,3*w0,N);eta=linspace(0,2*pi,N);[rho,theta]=meshgrid(r,eta);%生成極坐標(biāo)網(wǎng)格[x,y]=pol2cart(theta,rho);%將極坐標(biāo)網(wǎng)格轉(zhuǎn)化為直角坐標(biāo)網(wǎng)格I=exp(-2*rho.^2./w0.^2);%高斯光束在歸一化的光強(qiáng)值surf(x,y,I);colormaphot;%colormapgray線性灰度圖colormaphot熱色顏色圖去掉這句命令為彩圖shadinginterp;xlabel('x/m');ylabel('y/m');zlabel('intensity/a.u.');axis([-3*w03*w0-3*w03*w001])title('焦點(diǎn)附近三維光場(chǎng)分布')colorbar;2、高斯光束在光腰處的二維光強(qiáng)分布圖clearallN=200;w0=200e-9;%高斯光束的光腰半徑x=linspace(-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.^2+Y.^2);I=exp(-2*R.^2./w0.^2);imagesc(x,y,I,[01]);axisequalcolormaphot;%colormapgray線性灰度圖colormaphot熱色顏色圖去掉這句命令為彩圖colorbar;xlabel('x/m');ylabel('y/m');axis([-3*w03*w0-3*w03*w0])title('焦點(diǎn)附近三維光場(chǎng)分布')3、高斯光束在傳播過(guò)程中，焦點(diǎn)處的三維光強(qiáng)分布clearallw0=200e-9;%高斯光束的光腰半徑wavelength=780e-9;%高斯光束的光波長(zhǎng)ZR=pi*w0^2/wavelength;%高斯光束的瑞利長(zhǎng)度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/w0^2).^2);I=w0^2./wz.^2.*exp(-2*pr.^2./wz.^2);surf(pz,pr,I);shadinginterp;xlabel('Z(λ)/m');ylabel('r(λ)/m');zlabel('intensity/a.u.');axis([-3*ZR3*ZR-InfInf-InfInf])title