高斯光束地matlab仿真

高斯光束地matlab仿真

ID:36730637

大小:424.42 KB

页数:15页

时间:2019-05-14

高斯光束地matlab仿真_第1页
高斯光束地matlab仿真_第2页
高斯光束地matlab仿真_第3页
高斯光束地matlab仿真_第4页
高斯光束地matlab仿真_第5页
资源描述:

《高斯光束地matlab仿真》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

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

2、束强度分布读入的数据是一个224X244的矩阵,矩阵中的数值代表光强分布。用读入的数据取中间一行(122行)画出强度分布如图2所示。文案大全实用标准图2实验测量高斯曲线用理论上的高斯曲线公式画出理论高斯曲线如图3所示。文案大全实用标准图3理论高斯曲线M文件如下:A=imread('D:documents作业激光原理与应用高斯.bmp');A1=A(:,122);x1=1:1:224;x2=-100:1:100;a2=exp(-x2.^2/10);figureimshow(A);axisofftitle('fontsize{12}CCD采集的高斯光束光强分布');figureplot

3、(x2,a2,'linewidth',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作业激光

4、原理与应用高斯.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=inpu

5、t('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^-

6、3;%dx:stepsizedx=L/N;forn=1:N+1form=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;%Freespacetransferfunction文案大全实用标准H(n,

7、m)=exp(j/(2*ko)*z*(Kx(m)^2+Ky(n)^2));endend%GaussianBeaminFrequencydomainFGau=fft2(Gau);FGau=fftshift(FGau);%PropagatedGaussianbeaminFrequencydomainFGau_pro=FGau.*H;%PeakamplitudeoftheinitialGaussianbeamPea

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。