高斯光束的matlab仿真

1、题目：根据高斯光束数学模型， 模拟仿真高斯光束在谐振腔中某一位置处的归一化强度分布并给出其二维、三维强度分布仿真图；用Matlab读取实际激光光斑照片中所记录的强度数Matlab读取照片数据命令为imread）,据（读取照片中光斑的一个直径所记录的强度数据即可,（如同时考用该数据画出图片中激光光斑的强度二维分布图，与之前数学模型仿真图对比。虑高斯光束光斑有效截面半径和等相位面特点,仿真高斯光束光强、光斑有效截面半径以及16 / 12等相位面同时随传播距离 z的变化并给出整体仿真图可酌情加分。） 原始光斑如图1所示，用imread命令读入 matlab后直接用imshow命令读取即可,CCD采集

2、的高斯光束光强分布图1 CCD采集的高斯光束强度分布读入的数据是一个 224 X 244的矩阵，矩阵中的数值代表光强分布。用读入的数据取中间一 行（122行）画出强度分布如图 2所示。实验测量高斯曲线18005016014012010080604020100150200图2实验测量高斯曲线用理论上的高斯曲线公式画出理论高斯曲线如图3所示。理论高斯曲线0.201020300-40-30-20-1040图3理论高斯曲线M文件如下:作业激光原理与应用高斯.bmp');A=imread(A1=A(:,122); x1=1:1:224; x2=-100:1:100;a2=ex

3、 p(-X2A2/10); figure imshow(A);axis offti tle('fon tsize12CCDfigurep lot(x2,a2,采集的高斯光束光强分布);'li newidth',1, 'color' , 'b');axis(-40 40 0 1.2)ti tle('fon tsize12figure plot(x1,A1,实验测量高斯曲线')'li newidth'ti tle('fon tsize12axis(50 200 0 180),1, 'color&

4、#39;理论高斯曲线，丫 )画三维强度分布。取图片矩阵的中间层，用mesh命令画出三维图如图 4所示。图4三维强度分布由于读入的图片有一行白边， 需要手动去除掉， 最终的M文件如下。否则三维图会有一边整体竖起来，影响观察。高斯.bmP');作业激光原理与应用high, width, color = size(A);x=1:width;y=1:high-1;mesh(x', y', double(A(2:224,:,1);grid onxlabel('x'),ylabel('y'),zlabel('z');title(

5、9;三维强度分布');再用 matlab 仿真理论上传播过程中高斯光束的变化 这次先给出 M 文件：%Gaussian_propagation.m%Simulation of diffraction of Gaussian Beamclear;%Gaussian Beam%N:sampling numberN=input('Number of samples(enter from 100 to 500)=');L=10T0A-3;Ld=input( 'wavelength of light in micrometers=');Ld=Ld*10A-6; k

6、o=(2*pi)/Ld; wo=input( 'Waist of Gaussian Beam in mm=' wo=wo*10A-3;z_ray=(ko*wo2)/2*103;);sprintf( 'Rayleigh range is %f mm' z_ray=z_rayT0A-3;z=input( 'Propagation length (z) in mm',z_ray);z=z*10A-3; %dx:ste p size dx=L/N;for n=1:N+1for m=1:N+1%Space axisx(m)=(m-1)*dx-L/2;y(n

7、)=(n-1)*dx-L/2;%Gaussian Beam in space domainGau( n,m)=ex p(-(x(m)A2+y( n) 2)/(wo2);Kx(m)=(2*pi*(m-1)/(N*dx)-(2*pi*(N)/(N*dx)/2;%Frequency axisendKy(n)=(2*pi*(n-1)/(N*dx)-(2*pi*(N)/(N*dx)/2;%Free space transfer functionH(n ,m)=ex pa/(2*ko)*z*(Kx(m)A2+Ky(门)人2); end%Gaussian Beam in Frequency domain F

8、Gau=fft2(Gau);FGau=fftshift(FGau);%Propagated Gaussian beam in Frequency domain FGau_pro=FGau.*H;%Peak amplitude of the initial Gaussian beam Peak_ini=max(max(abs(Gau);sprintf( 'Initial peak amplitude is %f mm',Peak_ini)%PropagatedGaussian beam in space domain Gau_pro=ifft2(FGau_pro);Gau_pro

9、=Gau_pro;%Peak amplitude of the propagated Gaussian beam Peak_pro=max(max(abs(Gau_pro);mm' ,Peak_pro) %Calculated Beamsprintf( 'Propagated peak amplitude is %f WidthN M=min(abs(x);Gau_pro1=Gau_pro(:,M);N1 M1=min(abs(abs(Gau_pro1)-abs(exp(-1)*Peak_pro);Bw=dx*abs(M1-M)T0A3;sprintf( 'Beam w

10、idth(numerical) is %fmm'W=(2*z_ray)/ko*(1+(z/z_ray)A2);W=(W0.5)*103;sprintf( 'Beam width(theoretical) is %fmm'x=x*103;y=y*10A3； figure(1);mesh(x,y,abs(Gau)title(xlabel(ylabel(,Bw) %Theoretical Beam Width,W) %axis in mm scale'Initial Gaussian Beam''x mm')'y mm')axi

11、s(min(x) max(x) min(y) max(y) 0 1) axis square figure(2);mesh(x,y,abs(Gau_pro)'propagated Gaussian Beam''x mm')'y mm')title(xlabel(ylabel(axis(min(x) max(x) min(y) max(y) 0 1) axis square程序主要根据高斯光束的传播规律计算传播过程中任意z处的高斯光强分布。运行结果：Number of samples(enter from 100 to 500)=500 wavel

12、ength of light in micrometers=0.568 Waist of Gaussian Beam in mm=1 ans =Rayleigh range is 5530.972982 mmProp agation len gth in mm100000 ans =In itial peak amp litude is 1.000000 mm ans =Prop agated p eak amp litude is 0.210252 mm ans =Beam?width( nu merical) is 1.940000mm ans =Beam?width(theoretical) is 18.107635mm>>束腰半径处的理想高斯光强分布0.3 XO.S.propagted Gaussian Bewm0-2、50y |mmX mnn0