Matlab：Fox_Li数值迭代法求激光器谐振腔本征模式

最新推荐文章于 2026-06-06 08:43:32 发布

原创

最新推荐文章于 2026-06-06 08:43:32 发布 · 1.1w 阅读

本文探讨了谐振腔经典理论的局限性及新型结构谐振腔的模式求解方法，介绍了Fox-Li迭代法等数值模拟方法在平行平面腔、正支共焦非稳腔和方形镜对称共焦腔中的应用。通过实例展示了不同腔型的本征模式计算过程及结果，强调激光谐振腔中光场分布仅与腔结构相关。

谐振腔的经典理论仅给出了部分简单腔型的模式解析解。对于平行平面腔以及在激光器的不断发展过程中所涌现的许多新型结构谐振腔通常是没有解析结果的 ,必须采用各种数值模拟方法进行求解 ,如 Fox-Li 迭代法、快速傅立叶变换法(FFT)、等效透镜波导法、特征向量法、有限元法(FEM)和有限差分法(FDM)等。
代码：

%% Fox-Li数值迭代法求谐振腔的本征模式   该程序单位统一为：m
function main ()
clc;    clear;      close all;
%% 参数设置
step = 300;             %采样点数
%   /* 平行平面腔参数 */
lambda1 = 1064e-9;      %波长1064nm
k1 = 2*pi/lambda1;
L1 = 100*lambda1;       %腔长100λ
a1 = 25*lambda1;        %反射镜半径 25λ
F1 = a1^2/lambda1/L1;   %菲涅耳数 6.25
x1 = linspace(-a1,a1,step);
%   /* 正支共焦非稳腔参数 */
lambda2 = 10.6e-6;      %波长10.6μm
k2 = 2*pi/lambda2;
L2 = 15;                %腔长15m
Mx = 2; My = 2; M = 2;  %光腔的几何放大率
Nc = 0.83;
a2_1 = 2e-2;            %镜一尺寸
a2_2 = a2_1;            %镜二尺寸
x2 = linspace(-a2_1,a2_1,step);
%   /* 方形镜对称共焦腔参数 */
lambda3 = 1064e-9;      %波长1064nm
k3 = 2*pi/lambda3;
L3 = 100*lambda3;       %腔长100λ
a3 = 25*lambda3;        %反射镜半径 25λ
F3 = a3^2/lambda3/L3;   %菲涅耳数 6.25
x3 = linspace(-a3,a3,step);
%% 迭代一次
E0 = ones(1,step);  %均匀平面波
E0_1 = E0;
E0_2 = E0;
E0_3 = E0;
E1 = zeros(1,step);%先分配内存
E2 = zeros(1,step);%先分配内存
E3 = zeros(1,step);%先分配内存
for j = 1:step
    E1(j) = pxpmq(L1,a1,k1,lambda1,x1(j),E0_1,step);        %E1代表平行平面腔
    E2(j) = zzgjfwq(L2,a2_2,lambda2,Mx,x2(j),E0_2,step);    %E2代表正支共焦非稳腔
    E3(j) = dcgjq(L3,a3,k3,lambda3,x3(j),E0_3,step);        %E3代表方形镜对称共焦腔
end
phase_E1_1 = angle(E1);
I1_1 = E1.*conj(E1);    I1_1 = I1_1/max(max(I1_1));
phase_E2_1 = angle(E2);
I2_1 = E2.*conj(E2);    I2_1 = I2_1/max(max(I2_1));
phase_E3_1 = angle(E3);
I3_1 = E3.*conj(E3);    I3_1 = I3_1/max(max(I3_1));
%正支共焦非稳腔三维计算
% E0 = ones(1,step);  %均匀平面波
% E2