光强分布MATLAB.docx
《光强分布MATLAB.docx》由会员分享，可在线阅读，更多相关《光强分布MATLAB.docx（10页珍藏版）》请在冰豆网上搜索。
光强分布MATLAB
、高斯光束模拟的焦点附近光场分布
forn=1:
10
1=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,'-'）;
holdon
end
右半边：
forn=1:
10
1=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,'-'）;
holdon
plot（z,r2,'-'）;
holdon
end
全图：
clearall
forn=1:
10
1=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,'-'）;axisequal
holdon
plot（z,r2,'-'）;axisequal
holdon
end
xlabel（'Z（R/m'）;
ylabel（'r（R/m'）;
title（'焦点附近光场分布'）
带空间坐标:
clearall
forn=1:
10
1=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））.*z
A2））））;
r2=-w0*sqrt（0.5.*（1+（wl/（pi*w0A2））.*z.A2）.*log（1./（l*（1+（wl/（pi*w0A2））.*
z.A2））））;
plot（z,r1,'-'）;axisequal
holdon
plot（z,r2,'-'）;axisequal
holdon
end
xlabel（'Z（R/m'）;
ylabel（'r（R/m'）;
title（'焦点附近光场分布'）
Z'AI
近光塌为舟
012
cm—COXAY
-2-10123
XAxisTitle
、焦点附近的三维光场分布
三、高斯光束在传播过程中，焦点处的光强分布
1、高斯光束在光腰处的三维光强分布图
clearall
N=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）;%将极坐标网格转化为直角坐标网格
I=exp（-2*rho.A2./w0.A2）;%高斯光束在归一化的光强值
surf（x,y,I）;
colormaphot;%colormapgray线性灰度图colormap
hot热色颜色图去掉这句命令为彩图
shadinginterp;
xlabel（'x/m'）;
ylabel（'y/m'）;
zlabel（'intensity/a.u.'）;
axis（[-3*w03*w0-3*w03*w001]）
title（'焦点附近三维光场分布'）
colorbar;
2、高斯光束在光腰处的二维光强分布图
clearall
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.A2+Y.A2）;
I=exp（-2*R.A2./w0.A2）;
imagesc（x,y,l,[O1]）;axisequal
colormaphot;%colormapgray线性灰度图colormap
hot热色颜色图去掉这句命令为彩图
colorbar;
xlabel（'x/m'）;
ylabel（'y/m'）;
axis（[-3*w03*w0-3*w03*w0]）
title（'焦点附近三维光场分布'）
击』*202AE
3、高斯光束在传播过程中，焦点处的三维光强分布
%高斯光束的光腰半径
%高斯光束的光波长
%高斯光束的瑞利长度
%高斯光束的传输范围
clearall
w0=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）.A2）;
l=w0A2./wz.A2.*exp（-2*pr.A2./wz.A2）;
surf（pz,pr,I）;
shadinginterp;
xlabel（'Z（为/m'）;
ylabel（'r（R/m'）;
zlabel（'intensity/a.u.'）;
axis（[-3*ZR3*ZR-Inflnf-InfInf]）
title（'焦点附近三维光场分布'）
colorbar;
colormaphot;%colormapgray线性灰度图colormaphot热色颜
色图
view
（2）
%显示二维视图
oyoao?
0&
05
□i.
0J
0.2
01
近三维也馬労帘
-4-3-2-101234
Z（A>rtn*1（/
I1O議点削近三维也场労帝
-4-3-2-101234
Z（A}rtn
升级会员 