Matlab:拉盖尔高斯光束的干涉
代码:
clc
clear all
close all
%% 拉盖尔-高斯光束与平面光波干涉
N = 300; %取样点数
lambda = 632e-9; %波长632nm
k = 2*pi/lambda; %波数
w0 = 3;
x = linspace(-1e-4,1e-4,N); y = x;
[X,Y] = meshgrid(x,y);
[theta,r] = cart2pol(X,Y);p = 0;
Z_R = pi*w0^2/lambda; %瑞利长度
z = 0;
w_z = w0*sqrt(1+(z/Z_R)^2);%光束在z位置的半径
figure;
for m = -4:4E1 = exp(-1i*k*X); %平面波E2 = 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));c1 = E1+E2;E_1 = c1.*conj(c1);subplot(3,3,m+5)h1 = pcolor(X,Y,E_1);colorbar;set(h1,'edgecolor','none','facecolor','interp');title(['m = ',num2str(m)]);colormap(gray);axis square;
end
suptitle('拉盖尔-高斯光束与平面光波干涉') %为图一添加总标题
%% 涡旋光束与球面光波干涉
N = 200; %取样点数
lambda = 632e-9; %波长632nm
k = 2*pi/lambda; %波数
x = linspace(-2e-3,2e-3,N); y = x;
[X,Y] = meshgrid(x,y);
[theta,r] = cart2pol(X,Y);
Z = 1;
Z_R = pi*w0^2/lambda; %瑞利长度
z = 0;
w_z = w0*sqrt(1+(z/Z_R)^2);%光束在z位置的半径
figure;
for m = -4:4E3 = exp(-1i*k*Z*(1+0.5*X.^2/Z^2+0.5*Y.^2/Z^2)); %球面波E4 = 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));c2 = E3+E4;E_2 = c2.*conj(c2);subplot(3,3,m+5)h2 = pcolor(X,Y,E_2);colorbar;set(h2,'edgecolor','none','facecolor','interp');title(['m = ',num2str(m)]);colormap(gray);axis square;
end
suptitle('拉盖尔高斯光束与球面光波干涉') %为图二添加总标题
%% 拓扑荷数相反的拉盖尔高斯光束互相干涉
N = 200; %取样点数
lambda = 632e-9; %波长632nm
k = 2*pi/lambda; %波数
x = linspace(-8,8,N); y = x;
[X,Y] = meshgrid(x,y);
[theta,r] = cart2pol(X,Y);
Z_R = pi*w0^2/lambda; %瑞利长度
z = 0;
w_z = w0*sqrt(1+(z/Z_R)^2);%光束在z位置的半径
figure;
for m = 1:9E5 = 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));E6 = 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));c2 = E5+E6;E_3 = c2.*conj(c2);subplot(3,3,m)h3 = pcolor(X,Y,E_3);colorbar;set(h3,'edgecolor','none','facecolor','interp');title(['m = ',num2str(m)]);colormap(gray);axis square;
end
suptitle('拓扑荷数相反的拉盖尔高斯光束互相干涉') %为图三添加总标题
%% 拉盖尔多项式(文献3中的公式)
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
运行结果:
参考文献:
[1]徐丽娟. 涡旋光束的产生及特性研究[D]. 浙江大学, 2014.
[2]赵麒, 白忠臣, 周骅,等. 拉盖尔-高斯光束作用下熔石英温度及应力研究%Research of temperature and thermal stress of fused silica irradiated by Laguerre-Gaussian beam[J]. 激光技术, 2018, 042(001):121-126.
[3]石业娇. 面向Fredholm微分方程的广义拉盖尔多项式求解方法[J]. 湘潭大学自然科学学报, 2018, 040(001):31-35.
Matlab:拉盖尔高斯光束的干涉相关推荐
- 厄米高斯光束 matlab,拉盖尔高斯光束_厄米高斯光束MATLAB仿真
<拉盖尔高斯光束_厄米高斯光束MATLAB仿真>由会员分享,可在线阅读,更多相关<拉盖尔高斯光束_厄米高斯光束MATLAB仿真(4页珍藏版)>请在人人文库网上搜索. 1.激光原 ...
- Matlab:拉盖尔-高斯光束杨氏双缝干涉
代码: clc clear all close all %% L-G光束双缝干涉 N = 300; %取样点数 lambda = 632e-9; %波长632nm k = 2*pi/lambda; % ...
- 横向剪切干涉matlab仿真,横向剪切干涉检测非球面波前数值模拟
引言随着先进光学制造技术和干涉检测技术的不断发展,中.大型光学非球面在天文光学.空间光学.激光大气传输.惯性约束聚变(ICF)以及国防.军工等领域得到了越来越多的应用,并且已成为该领域起支撑作用的关键 ...
- 【光学】基于matlab GUI双孔干涉【含Matlab源码 2119期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[光学]基于matlab GUI双孔干涉[含Matlab源码 2119期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2: 付费 ...
- 【湍流】基于matlab模拟拉盖尔高斯光束传播的光强
✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信.
- Matlab:拉盖尔-高斯光束光阑衍射
代码: function main() clc; clear; close all; %% /* 参数设置 */ Nxy = 200; %x,y坐标采样点 lambda = 632e-9; %波长为6 ...
- matlab牛顿环gif,牛顿环干涉实验的 Matlab模拟
<牛顿环干涉实验的 Matlab模拟>由会员分享,可在线阅读,更多相关<牛顿环干涉实验的 Matlab模拟(5页珍藏版)>请在人人文库网上搜索. 1.牛顿环干涉牛顿环干涉实验的 ...
- 【光学】基于matlab模拟双孔干涉附matlab代码
1 内容介绍 通过Matlab软件编程,实现光学双缝干涉的计算机仿真.仿真结果对学生理解光学原理的基本概念很有帮助,提高了教学效果. 2 部分代码 %REDME!!!该仿真模拟了双孔干事实验中光屏逐 ...
- matlab模拟高斯光束波前相位分布,光束波前校正与信标研究
光束波前校正与信标研究 [摘要]: 大气湍流引起光波波前的随机扰动,是限制光波大气传输和天文成像分辨率的主 要因素.20世纪70年代发展起来的自适应光学技术能主动校正光波相位畸变,提高 光束远场能量集 ...
- 【湍流】基于Matlab模拟高斯光束在湍流大气中传输仿真,得到大气湍流相位屏、以及光斑强度变化
✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信.
最新文章
- 实现apply call 方法
- 关于级数∑(x n-x n-1)一致收敛性的一点儿理解
- php 导出csv字符串,PHP CSV字符串到数组
- 021_jdbc-mysql入门
- Net中的Request和Response对象的理解
- 1022词法分析实验总结
- Fibonacci数列整除性质的组合证明
- Java继承中成员变量和成员函数的覆盖
- mysql查球队胜场和败场sql_NBA球队_nba有哪些球队_NBA史上实力最弱的球队是哪个?用Python+SQL我们找到了答案...
- 剖析HotSpot的Launcher
- 每天一点正则表达式复习(二)
- 测试抑郁症软件,App Store 上的“抑郁症测试 - 心理测试”
- 易语言 内存调用html文件夹,将文件移到指定文件夹 易语言移动文件到指定文件夹内...
- Office server webs app 集成(JAVA)
- linux打开chm文件,安装kchmviewer,在linux下查看chm文件
- 【计算机网络】计算机网络的体系结构
- 关于Android中RSA数字签名的理解及使用
- 做自己喜欢做的事,需要勇气,也需要能力
- Go语言系列——31-自定义错误、32-panic和recover、33-函数是一等公民(头等函数)、34-反射、35-读取文件、36-写入文件
- 内部乙方部门的转型之困
热门文章
- 阿里云python面试题_一看“左程云:200道算法与数据结构”,二刷“阿里云:70+算法题、30种大厂笔试高频知识点”,3月过去终于挺进我梦中的字节!...
- 2015061906 - firebug安装和下载(2)
- BMP(DIB)图片格式
- 汇编指令大全(带注释)
- 程序员如何写工作日志?
- 如何在Internet Explorer中使用VARCHART XGantt控件?
- 【技术贴】从51下载的网站代码asp源码怎么运行?怎么打开?
- 绿色清爽win7电脑主题 +炫酷紫色win7主题+美女写真电脑主题
- windows 7 RC(7106.0.090408)下载另附windows7驱动收集整合(5月31日更新)
- eclipse安卓插件ADT下载地址