基于Matlab模拟涡旋光波和平面光波干涉仿真.zip

clc clear all close all lambda=1.064e-6; k=2*pi/lambda; x=(-1:0.01:1)*1e-4; y=x; [X,Y]=meshgrid(x,y); phi=atan2(Y,X); for l=1:5 E1=exp(-i*k*X); E2=exp(i*l*phi); c=E1+E2; E=c.*conj(c); E=(E-min(min(E)))/(max(max(E))-min(min(E))); subplot(2,3,l) pcolor(x,y,E) shading interp colormap gray axis off axis square; title(num2str(l)) % waitforbuttonpress end