高斯光束地matlab仿真.docxVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实用标准 题目:根据高斯光束数学模型, 模拟仿真高斯光束在谐振腔中某一位置处的归一化强度分布 并给出其二维、三维强度分布仿真图;用 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 理论高斯曲线 文案大全 实用标准 M文件如下: A=imread( 作业 激光原理与应用 高斯 .bmp ); 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 plot(x2,a2, linewidth ,1, color , b ); axis([-40 40 0 1.2]) 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文件: %Gaussian_propagation.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

文档评论(0)

153****3726 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档