clear
clc
close all
N=256; %---采样点数
L=0.5*10^-3; %---采样大小
wo=0.1*10^-3; %---输入高斯光束腰
dx=L/N;
row=N;
wave_length = 0.6328e-6; %---待测激光波长
f=0.4; %---透镜焦距
l=1; %---拓扑荷数
% figure;imagesc(theta);
length=L; %---采样尺寸
dx = length/(N); %----x空间取样间隔
dy = length/(N); %----y空间取样间隔
sx_range = -(length/2-dx/2):dx:(length/2); %---x空间取样范围
sy_range = -(length/2-dy/2):dy:(length/2); %---y空间取样范围
% aa=-1.2800e-04:dx:1.2800e-04;
[sx,sy] = meshgrid(sx_range,sy_range); %---空间取样-分割为二维场
h=i*wave_length*f*exp(i*atan2(sx,sy))./(2*pi*(sx.^2+sy.^2));
% h=(2*pi*(sx.^2+sy.^2));
% h=LBwindow(h,0,0,fix(0.0025/dx));
% hh=i*wave_length*f*exp(i*theta);
% figure;imagesc(abs(h));
% figure;imagesc(angle(h));
% % figure;plot(abs(h(ceil(N/2),:)));
% figure;plot(real(h(:,ceil(N/2))));
% figure;plot(real(h(ceil(N/2),:)));
% figure;plot(imag(h(:,ceil(N/2))));
% figure;plot(imag(h(ceil(N/2),:)));
% E0=ones(N,N);
% E0=LBwindow(E0,0,0,80);
% figure;imagesc(E0);colormap(gray);
% E1=conv2(E0,h,'same');
% % figure;imagesc(abs(E1));colormap(gray);
% figure;imagesc(E1.*conj(E1));colormap(gray);
figure;
for winp=fix(N)
h=i*wave_length*f*exp(i*l*atan2(sx,sy))./(2*pi*(sx.^2+sy.^2));
h=LBwindow(h,0,0,winp);
ii=20;
E00=ones(N,N);
E0=LBwindow(E00,0,0,ii);
% E_out=LBwindow(E00,0,0,ii-fix(ii/1.5));
% E0=E0-E_out;
% E0=E0.*exp(i*E0);
subplot(1,2,1);imagesc(sx_range*1000,sy_range*1000,E0.*conj(E0));colormap(gray);
set(gca,'FontSize',20);
xlabel('x/mm','FontSize',20);
ylabel('y/mm','FontSize',20);
E1=conv2(E0,h,'same');
subplot(1,2,2);imagesc(sx_range*1000,sy_range*1000,E1.*conj(E1));colormap(gray);
set(gca,'FontSize',20);
xlabel('x/mm','FontSize',20);
ylabel('y/mm','FontSize',20);
pause(0.01);
end
评论1