用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)

薄透镜的传输矩阵为：

对于自由空间传播距离d，其传输矩阵为：

function M = space(d) M = [1, d; 0, 1]; end

2. MATLAB仿真环境搭建

2.1 核心函数实现

我们需要三个关键函数来构建仿真系统：

1. 自由空间传播矩阵生成器：

function M = space(d) % 生成自由空间传播矩阵 % d: 传播距离(mm) M = [1, d; 0, 1]; end

2. 薄透镜变换矩阵生成器：

function M = lens(f) % 生成薄透镜变换矩阵 % f: 透镜焦距(mm)，正值为凸透镜，负值为凹透镜 M = [1, 0; -1/f, 1]; end

3. 高斯光束状态计算器：

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); % 波前曲率半径 end

2.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 分阶段传播模拟

将光束传输分为五个阶段进行模拟：

1. 初始束腰到透镜1：

z = 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

2. 通过透镜1：

M_lens1 = lens(params.f1); q = (M_lens1(1,1)*q + M_lens1(1,2))/(M_lens1(2,1)*q + M_lens1(2,2));

3. 透镜1到透镜2：

z = 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); end

3.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.d1+params.d2, '--k', 'Lens 2'); % 图形修饰 xlabel('传播距离 (mm)'); ylabel('光束半径 (mm)'); title('高斯光束通过双透镜系统传输'); legend('初始传播', '透镜1后传播', '透镜2后传播'); grid on;

4. 参数影响分析与优化

4.1 关键参数敏感性分析

通过修改以下参数观察系统响应：

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 end

5. 高级应用与扩展

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}; end

5.2 光束质量评估指标

除了光斑尺寸，还可以计算以下指标：

1. 光束参数乘积：

BP = w * theta; % w为光斑半径，theta为发散角

2. M²因子验证：

M2_measured = (pi * w^2 * theta) / lambda;

3. 聚焦特性分析：

[~, focus_idx] = min(w); focus_position = z(focus_idx);

在实际项目中，我发现合理设置采样点数量对结果精度影响很大。对于快速变化的区域（如束腰附近），需要增加局部采样密度：

% 非均匀采样示例 z_near_waist = linspace(waist_pos-10, waist_pos+10, 500); z_other = [linspace(0, waist_pos-10, 200), linspace(waist_pos+10, total_len, 200)]; z = unique(sort([z_near_waist, z_other]));