Matlab：涡旋光束与平面光波、球面光波干涉

matlab代码：

clc
clear all
close all
%%  涡旋光束与平面光波干涉
N = 300;            %取样点数
lambda = 632e-9;    %波长632nm
k = 2*pi/lambda;    %波数
x = linspace(-1e-2,1e-2,N);
y = linspace(-1e-2,1e-2,N);
[X,Y] = meshgrid(x,y);
[theta,r] = cart2pol(X,Y);

figure;
for m = -4:4
    E1