光强分布MATLABWord格式文档下载.docx

光强分布MATLABWord格式文档下载.docx

《光强分布MATLABWord格式文档下载.docx》由会员分享，可在线阅读，更多相关《光强分布MATLABWord格式文档下载.docx（9页珍藏版）》请在冰豆网上搜索。

r=sqrt（0.5.*（1+z.^2）.*log（1./（I.*（1+z.^2））））;

plot（z,r,'

-'

）;

holdon

end

右半边：

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,'

plot（z,r2,'

全图：

clearall

z=-a:

0.00001*a:

axisequal

xlabel（'

Z（λ）/m'

ylabel（'

r（λ）/m'

title（'

焦点附近光场分布'

）

带空间坐标：

w0=200e-9;

wl=780e-9;

a=pi*w0^2*sqrt（1./I-1）/wl;

0.0001*a:

r1=w0*sqrt（0.5.*（1+（wl/（pi*w0^2））.*z.^2）.*log（1./（I.*（1+（wl/（pi*w0^2））.*z.^2））））;

r2=-w0*sqrt（0.5.*（1+（wl/（pi*w0^2））.*z.^2）.*log（1./（I.*（1+（wl/（pi*w0^2））.*z.^2））））;

二、焦点附近的三维光场分布

三、高斯光束在传播过程中，焦点处的光强分布

1、高斯光束在光腰处的三维光强分布图

N=200;

%高斯光束的光腰半径

r=linspace（0,3*w0,N）;

eta=linspace（0,2*pi,N）;

[rho,theta]=meshgrid（r,eta）;

%生成极坐标网格

[x,y]=pol2cart（theta,rho）;

%将极坐标网格转化为直角坐标网格

I=exp（-2*rho.^2./w0.^2）;

%高斯光束在归一化的光强值

surf（x,y,I）;

colormaphot;

%colormapgray线性灰度图colormaphot热色颜色图去掉这句命令为彩图

shadinginterp;

x/m'

y/m'

zlabel（'

intensity/a.u.'

axis（[-3*w03*w0-3*w03*w001]）

焦点附近三维光场分布'

colorbar;

2、高斯光束在光腰处的二维光强分布图

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]）;

axis（[-3*w03*w0-3*w03*w0]）

3、高斯光束在传播过程中，焦点处的三维光强分布

wavelength=780e-9;

%高斯光束的光波长

ZR=pi*w0^2/wavelength;

%高斯光束的瑞利长度

L=3*ZR;

%高斯光束的传输范围

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）;

Z（λ）/m'

r（λ）/m'

axis（[-3*ZR3*ZR-InfInf-InfInf]）

%colormapgray线性灰度图colormaphot热色颜色图

view

（2）%显示二维视图