高斯光束的matlab仿真.docVIP

高斯光束的matlab仿真 高斯光束的matlab仿真 PAGE / NUMPAGES 高斯光束的matlab仿真 题目：根据高斯光束数学模型， 模拟仿真高斯光束在谐振腔中某一位置处的归一化强度分布 并给出其二维、三维强度分布仿真图；用 Matlab 读取实际激光光斑照片中所记录的强度数 据（读取照片中光斑的一个直径所记录的强度数据即可， Matlab 读取照片数据命令为 imread ）, 用该数据画出图片中激光光斑的强度二维分布图，与之前数学模型仿真图对比。 （如同时考 虑高斯光束光斑有效截面半径和等相位面特点， 仿真高斯光束光强、 光斑有效截面半径以及 等相位面同时随传播距离 z 的变化并给出整体仿真图可酌情加分。 ） 原始光斑如图 1 所示，用 imread 命令读入 matlab 后直接用 imshow 命令读取即可， CCD 采集的高斯光束光强分布 图 1 CCD采集的高斯光束强度分布 读入的数据是一个 224 X 244 的矩阵，矩阵中的数值代表光强分布。用读入的数据取中间一行（ 122 行）画出强度分布如图 2 所示。 实验测量高斯曲线 180 160 140 120 100 80 60 40 20 0 100 150 200 50 图 2 实验测量高斯曲线 用理论上的高斯曲线公式画出理论高斯曲线如图 3 所示。 理论高斯曲线 1 0.8 0.6 0.4 0.2 0 -30 -20 -10 0 10 20 30 40 -40 图 3 理论高斯曲线 文件如下： A=imread( 作业 激光原理与应用 A1=A(:,122); x1=1:1:224; x2=-100:1:100; a2=exp(-x2.^2/10); figure imshow(A); axis off title( \fontsize{12}CCD 采集的高斯光束光强分布 figure  高斯 );  .bmp ); plot(x2,a2, linewidth ,1,color ,b ); axis([-40 40 0 ]) title( \fontsize{12} 实验测量高斯曲线 ) figure plot(x1,A1, linewidth ,1,color ,r ) title( \fontsize{12} 理论高斯曲线 ) axis([50 200 0 180]) 画三维强度分布。取图片矩阵的中间层，用 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 on xlabel(x),ylabel(y),zlabel(z); title( 三维强度分布 ); 再用 matlab 仿真理论上传播过程中高斯光束的变化 这次先给出 M 文件： % %Simulation of diffraction of Gaussian Beam clear; %Gaussian Beam %N:sampling number N=input( Number of samples(enter from 100 to 500)= ); L=10*10^-3; Ld=input( wavelength of light in [micrometers]= ); Ld=Ld*10^-6; ko=(2*pi)/Ld; wo=input( Waist of Gaussian Beam in [mm]= ); wo=wo*10^-3; z_ray=(ko*wo^2)/2*10^3; sprintf( Rayleigh range is %f [mm] ,z_ray) z_ray=z_ray*10^-3; z=input( Propagation length (z) in [mm] ); z=z*10^-3; %dx:step size dx=L/N; for n=1:N+1 for m=1:N+1 %Space axis x(m)=(m-1)*dx-L/2; y(n)=(n-1)*dx-L/2; %Gaussian Beam in space domain Gau(n,m)=exp(-(x(m)^2+y(n)^2)/(wo^2)); %Frequency axis Kx(m)=(2*pi*(m-1))/(N*dx)-((2*pi*(N))/(N*dx))/2; Ky(n)=(2*pi*(n-1))/(N*dx)