2021-10-31-MATLAB 调试代码

%在这里记录一下代码
保存了自己的代码，有很大问题啊，设置的长和宽分别是960和1200，但是生成的图不是这么大像素的分辨率，还没有想通，这里调用了一个函数，不能公开。这个代码也只有我自己能用哈哈哈。

clear all
cm = 0.01;
um = 1e-6;
M_a=1200;
N_a=960;
radius=M_a/2;
nm = 1e-9;
lambda = 638.8*nm;
Pitch_spatial= 5*um;
 Pitch = 20*um;


xa=ones(M_a,1)*(-M_a/2:M_a/2-1);%--生成256行 （-128,127）的数  xa如果不乘以Pitch，在测试平面波光斑偏移时，才有偏移
ya=(-M_a/2:M_a/2-1)'*ones(1,M_a);%--生成256列（-128,127）的数
%--坐标转换--%



 %涡旋光探针函数
[theta,rou]=cart2pol(xa,ya);%将笛卡尔坐标转换为极坐标
 l=2;
ll=-10;


%  vox4=exp(j*llll*theta);
% 
%+vox2+vox3;%+vox4;


%--------------
L_rang = [-5,5];

%OAM谱的L维上的采样点的间隔
L_deta = N_a/(L_rang(2)-L_rang(1)+1);

%OAM谱的R维上的采样间隔
R_deta =20;

m = 1;

 for r= radius-1:-R_deta:3
    
      n = 1;
%---环形窄带滤波---%
      narrow_r_band = grating(M_a, N_a, Pitch_spatial, 1, r)