代碼：

clc clear all close all %% L-G光束雙縫干涉 N = 300; %取樣點數 lambda = 632e-9; %波長632nm k = 2*pi/lambda; %波數 x = linspace(-2e-5,2e-5,N); y = linspace(-2e-5,2e-5,N); [X,Y] = meshgrid(x,y); [theta,r] = cart2pol(X,Y); w0 = 3e-1; %束腰 d = 2e-4; %雙縫間距200μm D = 9e-4; %雙縫與觀察屏之間的距離900μm p = 1; Z_R = pi*w0^2/lambda; %瑞利長度 z = 0; w_z = w0*sqrt(1+(z/Z_R)^2);%光束在z位置的半徑 figure; for m = -4:4E = sqrt(2*factorial(p)/pi/(p+factorial(abs(m))))*(1/w_z)*(sqrt(2)*r/w_z).^abs(m)....*exp(-r.^2/w_z^2).*laguerre(p,abs(m),2*r.^2/w_z^2).*exp(-1i*m*theta).*exp(-1i*k*z)....*exp(-1i*k*r.^2*z/2/(z^2+Z_R^2))*exp(-1i*(2*p+abs(m)+1)*atan(z/Z_R));I = E.*conj(E);I_1 = 4*I.*cos(pi*X*d/lambda/D+delta_phi(m,Y)/2);subplot(3,3,m+5)h1 = pcolor(X,Y,I_1);colorbar;set(h1,'edgecolor','none','facecolor','interp');title(['m = ',num2str(m)]);colormap(gray);axis square; end suptitle('拉蓋爾-高斯光束雙縫干涉') %為圖一添加總標題 %% delta_phi function result = delta_phi(m,y)%delphs=@(L,y)L.*(2.*pi-2.*atan(a./y)).*(y>0&y<8)+L.*(2.*atan(-a./y)).*(y>-8&y<0);%m為拓撲荷數，改變拓撲荷數會使圖形中的明暗條紋分布發生極大的改變result = m*2*(0.5*pi+atan(1e7*y)); end %% 拉蓋爾多項式(文獻5中的公式) function result = laguerre(p,l,x) result = 0; if p == 0result = 1; elseif p == 1result = 1+abs(l)-x; elseresult = (1/p)*((2*p+l-1-x).*laguerre(p-1,abs(l),x)-(p+l-1)*laguerre(p-2,abs(l),x)); end end

運行結果：

△phi(y)那里我看的參考文獻中只有phase-y圖像，沒有具體的關系式，這個△phi(y)公式是參考一位大佬的代碼
博客連接
渦旋光及其部分干涉圖樣MATLAB模擬
參考文獻：
[1]徐麗娟. 渦旋光束的產生及特性研究[D]. 浙江大學, 2014.
[2]Sztul H I , Alfano R R . Double-slit interference with Laguerre-Gaussian beams[J]. Optics Letters, 2006, 31(7):999-1001.
[3]陳志婷. 渦旋光束的特性研究[D]. 燕山大學, 2013.
[4]趙麒, 白忠臣, 周驊,等. 拉蓋爾-高斯光束作用下熔石英溫度及應力研究%Research of temperature and thermal stress of fused silica irradiated by Laguerre-Gaussian beam[J]. 激光技術, 2018, 042(001):121-126.
[5]石業嬌. 面向Fredholm微分方程的廣義拉蓋爾多項式求解方法[J]. 湘潭大學自然科學學報, 2018, 040(001):31-35.

總結

以上是生活随笔為你收集整理的Matlab：拉盖尔-高斯光束杨氏双缝干涉的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。