代码:

clc
clear all
close all
%%  L-G光束双缝干涉
N = 300;            %取样点数
lambda = 632e-9;    %波长632nm
k = 2*pi/lambda;    %波数
x = linspace(-2e-5,2e-5,N);
y = linspace(-2e-5,2e-5,N);
[X,Y] = meshgrid(x,y);
[theta,r] = cart2pol(X,Y);
w0 = 3e-1;          %束腰
d = 2e-4;           %双缝间距200μm
D = 9e-4;           %双缝与观察屏之间的距离900μm
p = 1;
Z_R = pi*w0^2/lambda;      %瑞利长度
z = 0;
w_z = w0*sqrt(1+(z/Z_R)^2);%光束在z位置的半径
figure;
for m = -4:4E = sqrt(2*factorial(p)/pi/(p+factorial(abs(m))))*(1/w_z)*(sqrt(2)*r/w_z).^abs(m)....*exp(-r.^2/w_z^2).*laguerre(p,abs(m),2*r.^2/w_z^2).*exp(-1i*m*theta).*exp(-1i*k*z)....*exp(-1i*k*r.^2*z/2/(z^2+Z_R^2))*exp(-1i*(2*p+abs(m)+1)*atan(z/Z_R));I = E.*conj(E);I_1 = 4*I.*cos(pi*X*d/lambda/D+delta_phi(m,Y)/2);subplot(3,3,m+5)h1 = pcolor(X,Y,I_1);colorbar;set(h1,'edgecolor','none','facecolor','interp');title(['m = ',num2str(m)]);colormap(gray);axis square;
end
suptitle('拉盖尔-高斯光束双缝干涉')   %为图一添加总标题
%% delta_phi
function result = delta_phi(m,y)%delphs=@(L,y)L.*(2.*pi-2.*atan(a./y)).*(y>0&y<8)+L.*(2.*atan(-a./y)).*(y>-8&y<0);%m为拓扑荷数,改变拓扑荷数会使图形中的明暗条纹分布发生极大的改变result = m*2*(0.5*pi+atan(1e7*y));
end
%% 拉盖尔多项式(文献5中的公式)
function result = laguerre(p,l,x)
result = 0;
if p == 0result = 1;
elseif p == 1result = 1+abs(l)-x;
elseresult = (1/p)*((2*p+l-1-x).*laguerre(p-1,abs(l),x)-(p+l-1)*laguerre(p-2,abs(l),x));
end
end

运行结果:

△phi(y)那里我看的参考文献中只有phase-y图像,没有具体的关系式,这个△phi(y)公式是参考一位大佬的代码
博客连接
涡旋光及其部分干涉图样MATLAB模拟
参考文献:
[1]徐丽娟. 涡旋光束的产生及特性研究[D]. 浙江大学, 2014.
[2]Sztul H I , Alfano R R . Double-slit interference with Laguerre-Gaussian beams[J]. Optics Letters, 2006, 31(7):999-1001.
[3]陈志婷. 涡旋光束的特性研究[D]. 燕山大学, 2013.
[4]赵麒, 白忠臣, 周骅,等. 拉盖尔-高斯光束作用下熔石英温度及应力研究%Research of temperature and thermal stress of fused silica irradiated by Laguerre-Gaussian beam[J]. 激光技术, 2018, 042(001):121-126.
[5]石业娇. 面向Fredholm微分方程的广义拉盖尔多项式求解方法[J]. 湘潭大学自然科学学报, 2018, 040(001):31-35.

Matlab:拉盖尔-高斯光束杨氏双缝干涉相关推荐

  1. 厄米高斯光束 matlab,拉盖尔高斯光束_厄米高斯光束MATLAB仿真

    <拉盖尔高斯光束_厄米高斯光束MATLAB仿真>由会员分享,可在线阅读,更多相关<拉盖尔高斯光束_厄米高斯光束MATLAB仿真(4页珍藏版)>请在人人文库网上搜索. 1.激光原 ...

  2. Matlab:实现杨氏双缝干涉仿真

    本文基于光学原理,用Matlab实现杨氏双缝干涉现象仿真. Project Code % 作者:ZQJ % 日期:2021.1.4 星期一 %***********************模拟杨氏干涉 ...

  3. MATLAB:杨氏双缝干涉实验(Young‘s Interference)

    目录 1.概念 2.MATLAB代码实现 1.概念 提出时间:1801年 提出者:托马斯·杨​(Thomas Young) 双缝干涉实验起源于光学研究,以它的发现者​托马斯·杨​(Thomas You ...

  4. 运用ChitGPT AI大模型生成杨氏双缝干涉实验MATLAB代码

    由于杨氏双缝干涉实验涉及到光的波动性和干涉现象,因此需要用到复杂的光学理论和数学分析.MATLAB是一款强大的数学软件,可以方便地进行数值计算和可视化分析,因此可以用MATLAB来模拟和分析杨氏双缝干 ...

  5. 用matlab电磁波双缝干涉实验报告,杨氏双缝干涉实验的MatLab模拟.pdf

    杨氏双缝干涉实验的MatLab模拟 第31卷第6期 光学仪器 V01.31,No.6 2009年12月 OPrICALINSTRlrM[ENTS 文章编号:1005-5630(2009)06-0025 ...

  6. 杨氏双缝干涉 matlab仿真,基于MATLAB的杨氏双缝干涉实验仿真.pdf

    基于MATLAB的杨氏双缝干涉实验仿真.pdf ( ) 第 27卷 湖北师范学院学报 自然科学版 Vol27 第 1期 Jou rnal of Hubei Norm al U n iversity ( ...

  7. 用matlab电磁波双缝干涉实验报告,基于MATLAB的杨氏双缝干涉实验的研究

    1杨氏双缝干涉 如图1所示,设单缝光阑S 到狭缝S 1和S 2距离相等,狭缝S 1和S 2的宽度均为b ,双缝之间遮光部分的宽为a ,双缝间的距离为d (d =a +b ),光屏离双缝的距离为D . ...

  8. 杨氏双缝干涉 matlab仿真,基于MATLAB的杨氏双缝干涉实验仿真

    第 27卷 湖北师范学院学报 (自然科学版 ) Vol27 第 1期 Journal of HubeiNormalUniversity (Natural Science) No11, 2007 基于 ...

  9. 【光学】基于matlab GUI杨氏双缝干涉【含Matlab源码 001期】

    一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[光学]基于matlab GUI杨氏双缝干涉[含Matlab源码 001期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2: 付费 ...

最新文章

  1. Android ActionBar完全解析,使用官方推荐的最佳导航栏(下) .
  2. 数学大神攻克猜字游戏Wordle,求解算法成绩逼近理论极限,连信息论都用上了...
  3. Linux主机WordPress伪静态设置方法
  4. 我国有线电视信源编码现有体制和对高清的制约
  5. 开发人员学Linux(1):VirtualBox中安装CentOS7过程记录
  6. Jafka源码粗略解读之二--关于JMX
  7. 在数组中找重复数、只出现一次的数或丢失数的题目(Leetcode题解-Python语言)
  8. wifi×××破解的密码查看方式
  9. 海康摄像头如何查看IP,重置密码
  10. 【关于测试开发工程师】
  11. ZEMAX | 如何使用 ISO 元件制图工具
  12. CCIE培训分解实验内容
  13. 自己动手写CPU(6)简单算术操作指令
  14. Android 获取毫秒时间戳
  15. android 滤镜 原理,android openglse实现滤镜九宫格
  16. Python作业1-高级数据类型
  17. mysql组复制(MGR)——背景
  18. 识别字体软件测试,2行代码帮你搞定自动化测试的文字识别
  19. 2011年恒生电子校园招聘笔试题目
  20. 微信小程序调用相册和相机

热门文章

  1. 快速了解必要的网络知识
  2. sql 语句left join关联超过2张表时,随着left join一路向左,永远只能是最左边的关联表与新表关联吗?还是母表可以一直on下去?
  3. Juniper training之Junos as a second language
  4. 收藏(网站和软件下载地址)
  5. 记录一下如何购买云服务器
  6. TC Games 电脑玩手机游戏助手全新版本常见问题解决方法汇总
  7. 路由详解(九阳真经)
  8. SQL Server 2005利用分区对海量数据的处理 [转自13590--北极燕鸥]
  9. 产能节拍php,工厂生产线的产能如何提升(实例讲解)
  10. 视频号算法推荐机制,视频号好友点赞立即推送吗?:国仁网络资讯