Z?=?exp(-X.^2-Y.^2);
Z?=?mat2gray(abs(Z).^2);
figure
imshow(Z,[]);
%
[X,?Y]?=?meshgrid((-5:.1:5)*1e-3);
R?=?sqrt(X.^2+Y.^2);
T?=?atan(Y./X);
A?=?0.003;
w?=?1e-3;
l?=?1;
u?=?A.*(2^(0.5)*R/w).^abs(l).*exp(-R.^2/w^2).*exp(1i*l*T);
%?看妳計算出來的u圖像中心點值是nan(求atan(0)導致的),因此這裏把它設置為0
u(isnan(u))?=?0;
u?=?mat2gray(abs(u).^2);
figure
imshow(u,[]);
%?拼接圖像
im?=?[u?Z?u];
figure;?imshow(im,?[])