高斯光束的matlab仿真.docx

题目：根据高斯光束数学模型，模拟仿真高斯光束在谐振腔中某一位置处的归一化强度分布并给出其二维、三维强度分布仿真图；用Matlab读取实际激光光斑照片中所记录的强度数据（读取照片中光斑的一个直径所记录的强度数据即可，Matlab读取照片数据命令为imread），用该数据画出图片中激光光斑的强度二维分布图，与之前数学模型仿真图对比。

（如同时考虑高斯光束光斑有效截面半径和等相位面特点，仿真高斯光束光强、光斑有效截面半径以及等相位面同时随传播距离z的变化并给出整体仿真图可酌情加分。）

原始光斑如图1所示，用imread命令读入matlab后直接用imshow命令读取即可，

CCD采集的高斯光束光强分布

图1CCD采集的高斯光束强度分布

读入的数据是一个224X244的矩阵，矩阵中的数值代表光强分布。用读入的数据取中间一行（122行）画出强度分布如图2所示。

180

160

140

120

100

80

60

40

20

0

理论高斯曲线50 100 150 200

理论高斯曲线

图2实验测量高斯曲线

用理论上的高斯曲线公式画出理论高斯曲线如图3所示。

1

0.8

0.6

0.4

0.2

0

-40 -30 -20 -10 0 10 20 30 40

图3理论高斯曲线

M文件如下:

A=imread('D:\documents、作业、激光原理与应用、高斯.bmp');

A1=A(:,122);

x1=1:1:224;

x2=-100:1:100;

a2=exp(-x2."2/10);

figure

imshow(A);

axisoff

title('\fontsize{12}CCD采集的高斯光束光强分布’)；

figure

plot(x2,a2,Tinewidth',1,'color','b');

axis([-404001.2])

title('\fontsize{12}实验测量高斯曲线')

figure

plot(x1,A1,'linewidth',1,'color','r')

title('\fontsize{12}理论高斯曲线')

axis([502000180])

画三维强度分布。取图片矩阵的中间层，用mesh命令画出三维图如图4所示。

图4三维强度分布

由于读入的图片有一行白边，需要手动去除掉，否则三维图会有一边整体竖起来，影响观察。最终的M文件如下。

A=imread('D:\documents\作业、激光原理与应用、高斯.bmp');

[high,width,color]=size(A);

x=1:width;

y=1:high-1;

mesh(x',y',double(A(2:224,:,1)));

gridon

xlabel('x'),ylabel('y'),zlabel('z');

title('三维强度分布')；

再用matlab仿真理论上传播过程中高斯光束的变化

这次先给出M文件：

%Gaussian_propagation.m

%SimulationofdiffractionofGaussianBeamclear;

%GaussianBeam

%N:samplingnumber

N=input('Numberofsamples(enterfrom100to500)=');

L=10*10"-3;

Ld=input('wavelengthoflightin[micrometers]」)；

Ld=Ld*10"-6;

ko=(2*pi)/Ld;

wo=input('WaistofGaussianBeamin[mm]=');

wo=wo*10“-3;

z_ray=(ko*wo“2)/2*10“3;

sprintf('Rayleighrangeis%f[mm]',z_ray)z_ray=z_ray*10、3;

z=input('Propagationlength(z)in[mm]');

z=z*10、3;%dx:stepsizedx=L/N;

forn=1:N+1

form=1:N+1

%Spaceaxis

x(m)=(m-1)*dx-L/2;

y(n)=(n-1)*dx-L/2;

%GaussianBeaminspacedomain

Gau(n,m)=exp(-(x(m)“2+y(n)“2)/(wo“2));%Frequencyaxis

Kx(m)=(2*pi*(m-1))/(N*dx)-((2*pi*(N))/(N*dx))/2;

Ky(n)=(2*pi*(n-1))/(N*dx)-((2*pi*(N))/(N*dx))/2;

%Freespacetransferf