【光学】涡旋光与球面波的干涉【含Matlab源码 597期】

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【光学】基于matlab涡旋光与球面波的干涉【含Matlab源码 597期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab物理应用（初级版）

备注：
点击上面蓝色字体付费专栏Matlab物理应用（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab物理应用（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、部分源代码

clear
clc
row=1080; %像素点设置成1080pix
g1=zeros(row); %生成1080阶零矩阵，初始化g1
g2=zeros(row); %生成1080阶零矩阵，初始化g1
w0=100; %束腰半径
l=3; %拓扑荷值
for m1=1:row %涡旋光函数
for n1=1:row
if (m1-row/2).^{2+(n1-row/2).}2<row/2row/2
g1(n1,m1)=sqrt((m1-row/2)^2+(n1-row/2)2)^{abs(l)*exp(-((m1-row/2)}2+(n1-row/2)^2)/w02)exp(1ilatan2((n1-row/2),(m1-row/2)));
end
end
end
A0=max(max(abs(g1))); %球面波函数
for m1=1:row
for n1=1:row
if (m1-row/2).^{2+(n1-row/2).}2<row/2*row/2
g2(n1,m1)=A0.exp(-((m1-row/2)^2+(n1-row/2)2)/(row0.2)^{2).*exp(1i*pi*-2*((m1-row/2)}2+(n1-row/2)^2)/(row/6)2);
end
end
end

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 门云阁.MATLAB物理计算与可视化[M].清华大学出版社，2013.

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除