用MATLAB复现高斯光束通过双透镜系统：从ABCD矩阵到可视化光斑演变

用MATLAB复现高斯光束通过双透镜系统从ABCD矩阵到可视化光斑演变在光学工程和激光技术领域理解高斯光束在复杂光学系统中的传输特性至关重要。本文将带您一步步实现高斯光束通过双透镜系统的完整MATLAB仿真从ABCD矩阵理论推导到动态光斑演变可视化让抽象的光学公式转化为直观的图形输出。1. 高斯光束与ABCD定律基础高斯光束是激光光学中最基本的传输模式其特性由束腰半径$w_0$和波前曲率半径$R(z)$描述。在自由空间传播时高斯光束的复参数q定义为q(z) z i*z_R其中$z_R\pi w_0^2/\lambda$是瑞利长度。当光束通过光学元件时ABCD定律告诉我们q (A*q B)/(C*q D)薄透镜的传输矩阵为矩阵元素物理意义表达式A横向放大率1B等效传播距离0C光学功率-1/fD角度变换因子1对于自由空间传播距离d其传输矩阵为function M space(d) M [1, d; 0, 1]; end2. MATLAB仿真环境搭建2.1 核心函数实现我们需要三个关键函数来构建仿真系统自由空间传播矩阵生成器function M space(d) % 生成自由空间传播矩阵 % d: 传播距离(mm) M [1, d; 0, 1]; end薄透镜变换矩阵生成器function M lens(f) % 生成薄透镜变换矩阵 % f: 透镜焦距(mm)正值为凸透镜负值为凹透镜 M [1, 0; -1/f, 1]; end高斯光束状态计算器function [w, R] beamPropagate(M, q, lambda, M2) % 计算通过光学元件后的光束参数 % M: ABCD矩阵 % q: 输入复参数 % lambda: 波长(mm) % M2: 光束质量因子 q_new (M(1,1)*q M(1,2))/(M(2,1)*q M(2,2)); w sqrt(lambda*M2/pi/abs(imag(1/q_new))); % 光斑半径 R 1/real(1/q_new); % 波前曲率半径 end2.2 参数初始化设置建议使用结构体组织仿真参数便于管理和修改params.lambda 1.064e-3; % 波长(mm) params.M2 1.5; % 光束质量因子 params.w0 0.1; % 初始束腰半径(mm) params.f1 -50; % 透镜1焦距(mm) params.f2 30; % 透镜2焦距(mm) params.d1 60; % 初始束腰到透镜1距离(mm) params.d2 50; % 透镜间距(mm) params.d3 40; % 透镜2到观察面距离(mm)3. 完整传输路径仿真实现3.1 分阶段传播模拟将光束传输分为五个阶段进行模拟初始束腰到透镜1z linspace(0, params.d1, 1000); w zeros(size(z)); for i 1:length(z) M space(z(i)); [w(i), ~] beamPropagate(M, q0, params.lambda, params.M2); end通过透镜1M_lens1 lens(params.f1); q (M_lens1(1,1)*q M_lens1(1,2))/(M_lens1(2,1)*q M_lens1(2,2));透镜1到透镜2z linspace(0, params.d2, 1000); w_lens1_to_lens2 zeros(size(z)); for i 1:length(z) M space(z(i)); [w_lens1_to_lens2(i), ~] beamPropagate(M, q, params.lambda, params.M2); end3.2 可视化结果输出使用MATLAB绘图函数展示完整传输路径figure(Position, [100, 100, 800, 600]) hold on; % 绘制各段传播曲线 plot(z1, w1, b, LineWidth, 2); plot(z2, w2, r, LineWidth, 2); plot(z3, w3, g, LineWidth, 2); % 标记透镜位置 xline(params.d1, --k, Lens 1); xline(params.d1params.d2, --k, Lens 2); % 图形修饰 xlabel(传播距离 (mm)); ylabel(光束半径 (mm)); title(高斯光束通过双透镜系统传输); legend(初始传播, 透镜1后传播, 透镜2后传播); grid on;4. 参数影响分析与优化4.1 关键参数敏感性分析通过修改以下参数观察系统响应参数典型范围对光束影响初始束腰0.05-0.5 mm决定光束发散角和瑞利长度透镜1焦距-100至100 mm控制第一次光束变换的强度透镜间距10-200 mm影响两透镜间的光束演化光束质量1.0-3.0决定光束发散特性和最终光斑大小4.2 交互式参数探索创建交互式界面方便参数调整function interactiveSimulator() % 创建图形界面 fig uifigure(Name, 高斯光束仿真); % 添加控制组件 w0Slider uislider(fig, Position, [100 350 200 3],... Limits, [0.05 0.5], Value, 0.1); f1Slider uislider(fig, Position, [100 300 200 3],... Limits, [-100 100], Value, -50); % 添加回调函数 w0Slider.ValueChangedFcn (src,event) updatePlot(); f1Slider.ValueChangedFcn (src,event) updatePlot(); function updatePlot() % 获取当前参数值 current_w0 w0Slider.Value; current_f1 f1Slider.Value; % 更新仿真并重绘 % ...仿真代码... end end5. 高级应用与扩展5.1 多透镜系统建模通过串联多个透镜矩阵可以模拟更复杂的光学系统% 三透镜系统示例 lens_system {... lens(100), space(50),... lens(-75), space(80),... lens(50)}; % 计算整体传输矩阵 M_total eye(2); for i 1:length(lens_system) M_total M_total * lens_system{i}; end5.2 光束质量评估指标除了光斑尺寸还可以计算以下指标光束参数乘积BP w * theta; % w为光斑半径theta为发散角M²因子验证M2_measured (pi * w^2 * theta) / lambda;聚焦特性分析[~, focus_idx] min(w); focus_position z(focus_idx);在实际项目中我发现合理设置采样点数量对结果精度影响很大。对于快速变化的区域如束腰附近需要增加局部采样密度% 非均匀采样示例 z_near_waist linspace(waist_pos-10, waist_pos10, 500); z_other [linspace(0, waist_pos-10, 200), linspace(waist_pos10, total_len, 200)]; z unique(sort([z_near_waist, z_other]));