高斯光束的matlab仿真教学内容

高斯光束的matlab仿真精品文档题目：根据高斯光束数学模型，模拟仿真高斯光束在谐振腔中某一位置处的归一化强度分布并给出其二维、三维强度分布仿真图；用 Matlab读取实际激光光斑照片中所记录的强度数据（读取照片中光斑的一个直径所记录的强度数据即可，Matlab读取照片数据命令为 imread）,用该数据画出图片中激光光斑的强度二维分布图，与之前数学模型仿真图对比。（如同时考虑高斯光束光斑有效截面半径和等相位面特点，仿真高斯光束光强、光斑有效截面半径以及等相位面同时随传播距离 z的变化并给出整体仿真图可酌情加分。）原始光斑如图1所示，用imread命令读入matlab后直接用imshow命令读取即可，CCD采集的高斯光束光强分布图1CCD采集的高斯光束强度分布读入的数据是一个 224X244的矩阵，矩阵中的数值代表光强分布。用读入的数据取中间一行（122行）画出强度分布如图 2所示。收集于网络，如有侵权请联系管理员删除精品文档实验测量高斯曲线18016014012010080604020010015020050图2实验测量高斯曲线用理论上的高斯曲线公式画出理论高斯曲线如图3所示。收集于网络，如有侵权请联系管理员删除精品文档理论高斯曲线10.80.60.40.20-30-20-10010203040-40图3理论高斯曲线文件如下：A=imread(作业激光原理与应用高斯.bmp');A1=A(:,122);x1=1:1:224;x2=-100:1:100;a2=exp(-x2.^2/10);figureimshow(A);axisofftitle('\fontsize{12}CCD采集的高斯光束光强分布');figureplot(x2,a2,'linewidth',1,'color','b');axis([-404001.2])title('\fontsize{12}实验测量高斯曲线')figureplot(x1,A1,'linewidth',1,'color','r')title('\fontsize{12}理论高斯曲线')axis([502000180])画三维强度分布。取图片矩阵的中间层，用 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)));gridonxlabel('x'),ylabel('y'),zlabel('z');title('三维强度分布');再用matlab仿真理论上传播过程中高斯光束的变化这次先给出M文件：%Gaussian_propagation.m%SimulationofdiffractionofGaussianBeamclear;%GaussianBeam%N:samplingnumberN=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;for n=1:N+1for m=1:N+1%Spaceaxisx(m)=(m-1)*dx-L/2;y(n)=(n-1)*dx-L/2;%GaussianBeaminspacedomainGau(n,m)=exp(-(x(m)^2+y(n)^2)/(wo^2)); %FrequencyaxisKx(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;%FreespacetransferfunctionH(n,m)=exp(j/(2*ko)*z*(Kx(m)^2+Ky(n)^2));endend%GaussianBeaminFrequencydomainFGau=fft2(Gau);FGau=fftshift(FGau);%PropagatedGaussianbeaminFrequencydomainFGau_pro=FGau.*H;%PeakamplitudeoftheinitialGaussianbeamPeak_ini=max(max(abs(Gau)));sprintf( 'Initialpeakamplitudeis%f[mm]' ,Peak_ini) %PropagatedGaussianbeaminspacedomainGau_pro=ifft2(FGau_pro);Gau_pro=Gau_pro;%PeakamplitudeofthepropagatedGaussianbeamPeak_pro=max(max(abs(Gau_pro)));sprintf( 'Propagatedpeakamplitudeis%f[mm]' ,Peak_pro) %CalculatedBeamWidth[NM]=min(abs(x));Gau_pro1=Gau_pro(:,M);[N1M1]=min(abs(abs(Gau_pro1)-abs(exp(-1)*Peak_pro)));Bw=dx*abs(M1-M)*10^3;sprintf( 'Beamwidth(numerical)is%f[mm]' ,Bw) %TheoreticalBeamWidthW=(2*z_ray)/ko*(1+(z/z_ray)^2);收集于网络，如有侵权请联系管理员删除精品文档W=(W^0.5)*10^3;sprintf('Beamwidth(theoretical)is%f[mm]',W)%axisinmmscalex=x*10^3;y=y*10^3;figure(1);mesh(x,y,abs(Gau))title('InitialGaussianBeam')xlabel('x[mm]')ylabel('y[mm]')axis([min(x)max(x)min(y)max(y)01])axissquarefigure(2);mesh(x,y,abs(Gau_pro))title('propagatedGaussianBeam')xlabel('x[mm]')ylabel('y[mm]')axis([min(x)max(x)min(y)max(y)01])axissquare程序主要根据高斯光束的传播规律计算传播过程中任意 z处的高斯光强分布。运行结果：Numberofsamples(enterfrom100to500)=500wavelengthoflightin[micrometers]=0.568WaistofGaussianBeamin[mm]=1ans=Rayleighrangeis5530.972982[mm]Propagationlength(z)in[mm]100000ans=Initialpeakamplitudeis1.000000[mm]ans=Propagatedpeakamplitudeis0.210252[mm]收集于网络，如有侵权请联系管理员删除精品文档ans=Beam?width(numerical)is1.940000[mm]ans=Beam?width(theoretical)is18.107635[mm]>>束腰半径处的理想高斯光强分布收集于网络，如有侵权请联系管理员删除精品文档传播1m处的理想高斯光强分布收集于网络，如有侵权请联系管理员删除精品文档传播10m处的理想高斯光强分布传播20m处的理想高斯光强分布收集于网络，如有侵权请联系管理员删除精品文档传播30m处的理想高斯光强分布传播50m