获取代码方式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