%?本程序計算圓形高斯光束入射到位于?束腰處的?圓盤產生的衍射及伯松亮斑
close?all;clear?all;
w?=?0.25;????????????????????%?入射光束腰，束腰半徑為0.25mm
lamda?=?0.632e-3;???????????%?入射激光波長
k?=?2*pi/lamda;?????????????%?波數
q?=?1i*pi*w^2/lamda;????????%?束腰處的q參數
N?=?500;????????????????????%?光強采樣精度
x?=?linspace（-0.50.5N）;???%?觀察光斑的坐標范圍
z?=?10;?????????????????????%?衍射傳播的距離
E0?=?exp（-1i*k*x.^2/2/q）;???%?束腰處的光場復振幅分布

[X?Y]?=?meshgrid（xx）;??????%?坐標網格化
[Ex?Ey]?=meshgrid（E0E0）;???%?光場二維化
E1?=?Ex.*Ey;????????????????%?束腰處的二維光場
figure;imagesc（abs（E1）.^2）;?%?圓形高斯光束
colormap?hot;axis?image;

Disk?=?one