1. 光学系统像差仿真概述
作为一名光学工程师，我经常需要评估和优化光学系统的成像质量。像差分析是这一过程中不可或缺的环节。传统的光学设计软件虽然功能强大，但往往价格昂贵且不够灵活。经过多年实践，我发现MATLAB是一个极佳的选择，它不仅能实现专业光学软件的绝大部分功能，还能根据特定需求进行深度定制。
光学像差主要分为两大类：单色像差和色差。单色像差包括球差、彗差、像散、场曲和畸变等，它们与光学系统的几何形状有关；色差则是由不同波长的光在介质中折射率不同引起的。理解这些像差的成因和特性，对于设计高质量光学系统至关重要。
在MATLAB中实现像差仿真有几个显著优势：首先，可以完全控制仿真的每个环节；其次，能够方便地与其他分析工具集成；最重要的是，可以建立完整的仿真流程，实现从设计到分析的全自动化。下面我将详细介绍如何利用MATLAB构建这样一个仿真系统。
2. Zernike多项式建模
2.1 Zernike多项式理论基础
Zernike多项式是描述光学波前畸变的标准数学工具，它定义在单位圆上，具有正交性，这使得各项系数相互独立。这种特性使得Zernike多项式特别适合用于像差分析和波前重构。
一个典型的Zernike多项式由径向和角向两部分组成：
code复制Z_n^m(ρ,θ) = R_n^m(ρ)cos(mθ) (m≥0)
Z_n^{-m}(ρ,θ) = R_n^m(ρ)sin(mθ) (m<0)
其中，ρ是归一化径向坐标(0≤ρ≤1)，θ是方位角，n是径向阶数，m是角向频率。径向多项式R_n^m(ρ)可以通过递归关系计算得到。
2.2 MATLAB实现细节
在MATLAB中实现Zernike多项式时，递归计算是最稳定的方法。下面是我优化过的实现代码：
matlab复制function Z = zernike(n, m, rho, theta)
% 参数验证
if n < 0 || abs(m) > n || mod(n-abs(m),2) ~= 0
error('Invalid Zernike indices (n,m)');
end
% 计算径向多项式
R = zeros(size(rho));
for s = 0:(n-abs(m))/2
coeff = (-1)^s * factorial(n-s) / ...
(factorial(s) * factorial((n+abs(m))/2-s) * ...
factorial((n-abs(m))/2-s));
R = R + coeff * rho.^(n-2*s);
end
% 组合角向部分
if m >= 0
Z = R .* cos(m*theta);
else
Z = R .* sin(abs(m)*theta);
end
% 归一化处理
if m == 0
norm = sqrt(n+1);
else
norm = sqrt(2*(n+1));
end
Z = Z * norm;
end
这段代码有几个关键优化点：
- 增加了参数有效性检查
- 使用循环代替递归提高计算效率
- 实现了标准归一化
- 支持复数形式的输出
2.3 像差模式对应关系
不同的Zernike项对应不同的像差模式：
| Zernike项 | 像差类型 | 物理意义 |
|---|---|---|
| Z_1^0 | 平移 | 整体波前平移 |
| Z_2^ | 倾斜 | 波前倾斜 |
| Z_2^0 | 离焦 | 焦点位置偏移 |
| Z_3^ | 彗差 | 非对称像差 |
| Z_3^ | 三叶草像差 | 高阶非对称像差 |
| Z_4^0 | 球差 | 径向对称高阶像差 |
| Z_4^ | 像散 | 正交方向不同焦距 |
理解这种对应关系对于像差分析和校正至关重要。在实际应用中，我们通常需要组合多个Zernike项来准确描述复杂的光学像差。
3. 光线追迹与波前计算
3.1 光线追迹基本原理
光线追迹是光学系统仿真的核心，它通过追踪光线通过光学系统的路径来计算像差。最基本的追迹包括：
- 折射：遵循斯涅尔定律
- 反射：遵循反射定律
- 传播：直线传播
对于旋转对称系统，我们可以简化计算，只考虑子午面内的光线。下面是一个优化的光线追迹函数：
matlab复制function [x_out, y_out, theta_out] = trace_ray(x_in, y_in, theta_in, surface)
% surface结构体包含:
% .R - 曲率半径
% .n1 - 入射侧折射率
% .n2 - 出射侧折射率
% .d - 到下一面的距离
% 计算入射点坐标
if isinf(surface.R) % 平面
x_hit = x_in;
y_hit = y_in + (surface.position - x_in) * tan(theta_in);
else % 球面
% 解光线与球面的交点方程
% 详细计算过程省略...
end
% 计算法线方向
if isinf(surface.R)
normal = [0; 1];
else
normal = [x_hit; y_hit - surface.R] / surface.R;
end
% 折射计算
n_ratio = surface.n1 / surface.n2;
cosI = dot([sin(theta_in); cos(theta_in)], normal);
sinI2 = n_ratio^2 * (1 - cosI^2);
if sinI2 > 1 % 全反射
% 反射处理...
else
% 折射方向计算...
end
% 传播到下一面
x_out = x_hit + surface.d * tan(theta_out);
y_out = y_hit + surface.d;
end
3.2 波前计算与可视化
通过光线追迹得到的光线交点数据可以重建波前。常用的方法是：
- 计算理想波前与实际波前的光程差(OPD)
- 对OPD进行Zernike多项式拟合
- 分析各项系数
matlab复制% 生成采样网格
[x,y] = meshgrid(linspace(-1,1,256));
rho = sqrt(x.^2 + y.^2);
theta = atan2(y,x);
valid = rho <= 1;
% 计算波前
wavefront = zeros(size(rho));
wavefront(valid) = 0.1*zernike(4,0,rho(valid),theta(valid)) + ... % 球差
0.05*zernike(3,1,rho(valid),theta(valid)); % 彗差
% 可视化
figure;
surf(x,y,wavefront,'EdgeColor','none');
xlabel('X'); ylabel('Y'); zlabel('波前误差(λ)');
title('波前像差分布');
colormap jet; colorbar;
3.3 性能优化技巧
光线追迹计算量很大，特别是需要追迹大量光线时。以下是一些优化建议：
- 向量化计算：避免使用循环，尽量用矩阵运算
- 并行计算：使用parfor或GPU加速
- 光线排序：按入射高度排序，提高缓存命中率
- 近似计算：对远轴光线使用近似公式
matlab复制% 并行光线追迹示例
num_rays = 10000;
rays = rand(num_rays, 2); % 随机光线
% 使用并行计算
parfor i = 1:num_rays
[x(i), y(i)] = trace_ray(rays(i,1), rays(i,2), ...);
end
4. 单色像差仿真实现
4.1 球差仿真与分析
球差是最常见的轴上像差，表现为边缘光线的焦点比近轴光线更靠近透镜。在MATLAB中，我们可以精确模拟球差：
matlab复制% 定义光学系统
f = 100; % 焦距(mm)
D = 25; % 孔径(mm)
NA = D/(2*f); % 数值孔径
% 生成光线
num_rays = 500;
h = linspace(0, D/2, num_rays); % 入射高度
% 追迹光线
for i = 1:num_rays
% 近轴近似计算理想像点
x_ideal(i) = h(i) * f / (D/2);
% 实际光线追迹
[x_real(i), ~] = trace_spherical_lens(h(i), f, D);
end
% 计算横向球差
TA = x_real - x_ideal;
% 绘制结果
figure;
plot(h, TA*1000, 'LineWidth', 2); % 转换为微米
xlabel('入射高度(mm)');
ylabel('横向像差(μm)');
title('球差曲线');
grid on;
球差的校正方法包括：
- 使用非球面
- 采用双胶合透镜
- 优化透镜形状因子
4.2 彗差仿真与分析
彗差是轴外点的主要像差，表现为不对称的弥散斑。仿真彗差需要离轴光线追迹：
matlab复制% 定义离轴角度
field_angle = 5; % 度
% 生成光线网格
[x,y] = meshgrid(linspace(-D/2,D/2,20));
valid = sqrt(x.^2 + y.^2) <= D/2;
x = x(valid); y = y(valid);
% 追迹光线
for i = 1:length(x)
% 计算入射方向
dir = [tand(field_angle); y(i)/f; x(i)/f];
dir = dir/norm(dir);
% 追迹光线
[x_img(i), y_img(i)] = trace_off_axis(x(i), y(i), dir, lens_system);
end
% 绘制点列图
figure;
plot(x_img, y_img, 'o');
axis equal; grid on;
xlabel('X方向(μm)'); ylabel('Y方向(μm)');
title(['彗差点列图(视场角=' num2str(field_angle) '°)']);
彗差的校正方法包括：
- 满足正弦条件
- 使用对称结构
- 控制透镜弯曲
4.3 像散与场曲仿真
像散和场曲通常一起分析。像散表现为子午和弧矢焦点的分离，场曲则是像面弯曲。
matlab复制% 定义视场点
field_points = linspace(0, 10, 5); % mm
% 初始化存储
T_spot = zeros(length(field_points), 2); % 子午焦点
S_spot = zeros(length(field_points), 2); % 弧矢焦点
% 对每个视场点追迹光线
for f = 1:length(field_points)
% 子午光线追迹
[T_spot(f,1), T_spot(f,2)] = trace_tangential(field_points(f));
% 弧矢光线追迹
[S_spot(f,1), S_spot(f,2)] = trace_sagittal(field_points(f));
end
% 计算场曲
field_curve_T = T_spot(:,2); % 子午场曲
field_curve_S = S_spot(:,2); % 弧矢场曲
% 绘制结果
figure;
plot(field_points, field_curve_T, 'r-', 'LineWidth', 2);
hold on;
plot(field_points, field_curve_S, 'b--', 'LineWidth', 2);
xlabel('视场位置(mm)');
ylabel('焦点位置(mm)');
legend('子午焦点', '弧矢焦点');
title('场曲与像散');
grid on;
5. 色差仿真与分析
5.1 轴向色差仿真
轴向色差表现为不同波长光线的焦点位置不同。仿真时需要追迹多个波长：
matlab复制% 定义波长(μm)
wavelengths = [0.486, 0.587, 0.656]; % F,d,C线
colors = {'b', 'g', 'r'};
% 追迹每个波长
figure; hold on;
for w = 1:length(wavelengths)
% 追迹轴上光线
[spot, ~] = trace_wavelength(wavelengths(w), 0);
% 绘制点列图
plot(spot(:,1), spot(:,2), [colors{w} '.'], 'MarkerSize', 10);
end
axis equal; grid on;
title('轴向色差弥散斑');
legend('486nm', '587nm', '656nm');
5.2 倍率色差仿真
倍率色差表现为不同波长成像大小不同：
matlab复制% 定义离轴视场点
field_point = 5; % mm
% 初始化存储
image_heights = zeros(length(wavelengths), 1);
% 追迹每个波长
for w = 1:length(wavelengths)
% 追迹主光线
[~, image_heights(w)] = trace_chief_ray(wavelengths(w), field_point);
end
% 计算倍率色差
chromatic_diff = image_heights - image_heights(2); % 相对于d线
% 绘制结果
figure;
bar(wavelengths, chromatic_diff*1000); % 转换为微米
xlabel('波长(μm)');
ylabel('倍率色差(μm)');
title('倍率色差分析');
grid on;
5.3 色差校正方法
在MATLAB中模拟色差校正通常涉及：
- 双胶合透镜设计
- 特殊材料选择
- 衍射光学元件
matlab复制% 双胶合透镜设计示例
% 定义材料
material1 = 'N-BK7';
material2 = 'F2';
% 计算阿贝数
V1 = get_abbe_number(material1);
V2 = get_abbe_number(material2);
% 优化曲率半径
R = optimize_doublet(V1, V2, wavelengths);
% 验证校正效果
[axial_CA, lateral_CA] = evaluate_chromatic_aberrations(R, wavelengths);
disp(['轴向色差: ' num2str(axial_CA) ' mm']);
disp(['倍率色差: ' num2str(lateral_CA) ' mm']);
6. 高级仿真功能实现
6.1 相位屏模拟大气湍流
大气湍流是影响天文观测的主要因素，可以用随机相位屏模拟：
matlab复制function phase = atmospheric_phase_screen(N, D, r0, L0)
% N: 网格点数
% D: 屏幕直径(m)
% r0: 相干长度(m)
% L0: 外尺度(m)
% 生成频率网格
delta_f = 1/D;
fx = (-N/2:N/2-1)*delta_f;
[Fx, Fy] = meshgrid(fx);
% 计算功率谱
f = sqrt(Fx.^2 + Fy.^2);
PSD = 0.023*r0^(-5/3) * f.^(-11/3) .* exp(-(f*L0).^2);
PSD(f==0) = 0;
% 生成随机相位
rand_phase = randn(N) + 1i*randn(N);
H = sqrt(PSD) .* fftshift(rand_phase);
phase = real(ifft2(ifftshift(H))) * N*delta_f;
end
6.2 自适应光学校真
自适应光学系统通过变形镜校正波前畸变：
matlab复制% 生成畸变波前
wavefront = generate_distorted_wavefront();
% 波前传感器测量
[slopes_x, slopes_y] = shack_hartmann(wavefront);
% 重构波前
z_coeff = reconstruct_wavefront(slopes_x, slopes_y);
% 控制变形镜
dm_commands = z_coeff_to_dm(z_coeff);
apply_dm_commands(dm_commands);
% 验证校正效果
corrected_wavefront = measure_wavefront();
residual_error = rms(corrected_wavefront(:));
disp(['残余波前误差: ' num2str(residual_error) ' λ']);
6.3 衍射效应模拟
考虑衍射效应可以提高仿真精度：
matlab复制% 定义孔径
lambda = 0.55e-3; % mm
D = 10; % 孔径直径(mm)
f = 100; % 焦距(mm)
% 计算点扩散函数(PSF)
x = linspace(-1,1,512)*lambda*f/D*10;
[X,Y] = meshgrid(x);
r = sqrt(X.^2 + Y.^2)/(lambda*f/D);
% 圆形孔径衍射
PSF = (2*besselj(1,pi*r)./(pi*r)).^2;
PSF(r==0) = 1;
% 显示PSF
figure;
imagesc(x,x,PSF);
colormap hot; axis image;
xlabel('X (mm)'); ylabel('Y (mm)');
title('衍射受限PSF');
7. 仿真结果分析与可视化
7.1 点列图分析
点列图直观显示成像质量：
matlab复制% 生成光线追迹数据
[spots_x, spots_y] = ray_trace_system(optical_system);
% 绘制点列图
figure;
scatter(spots_x, spots_y, 10, 'filled');
axis equal; grid on;
xlabel('X方向(μm)'); ylabel('Y方向(μm)');
title('点列图分析');
% 计算RMS半径
rms_radius = sqrt(mean(spots_x.^2 + spots_y.^2));
disp(['RMS半径: ' num2str(rms_radius) ' μm']);
7.2 MTF计算与显示
调制传递函数(MTF)评价系统分辨率：
matlab复制% 计算PSF
PSF = calculate_psf(optical_system);
% 计算MTF
MTF = abs(fftshift(fft2(PSF)));
MTF = MTF/max(MTF(:));
% 显示MTF
figure;
freq = linspace(-1,1,size(MTF,1)) * (1/(2*lambda*f/D));
plot(freq, MTF(:,size(MTF,2)/2), 'LineWidth', 2);
xlabel('空间频率(cycles/mm)');
ylabel('对比度');
title('衍射MTF');
grid on;
7.3 波前Zernike分析
分解波前到Zernike多项式：
matlab复制function coeff = zernike_fit(wavefront, max_order)
% wavefront: 输入波前
% max_order: 最大Zernike阶数
% 生成Zernike基
[x,y] = meshgrid(linspace(-1,1,size(wavefront,2)));
rho = sqrt(x.^2 + y.^2);
theta = atan2(y,x);
valid = rho <= 1;
% 构建设计矩阵
A = [];
for n = 0:max_order
for m = -n:2:n
Z = zernike(n, m, rho(valid), theta(valid));
A = [A Z(:)];
end
end
% 最小二乘拟合
coeff = A \ wavefront(valid);
end
8. 工程应用案例
8.1 手机镜头像差优化
现代手机镜头需要平衡体积和成像质量：
matlab复制% 定义优化目标
target_MTF = 0.3; % 在Nyquist频率处
weight = [1, 0.5, 0.2]; % 中心、0.7视场、边缘权重
% 优化变量: 曲率半径、厚度、材料
vars = [10, -8, 2, 1.5, ...]; % 初始猜测
% 优化
options = optimset('Display','iter');
opt_vars = fmincon(@(x) lens_cost(x, target_MTF, weight), ...
vars, [], [], [], [], lb, ub, [], options);
% 评估结果
[MTF, spots] = evaluate_lens(opt_vars);
disp(['优化后MTF: ' num2str(MTF)]);
8.2 望远镜系统仿真
天文望远镜需要考虑大气湍流和跟踪误差：
matlab复制% 定义望远镜参数
D = 1; % 口径(m)
f = 10; % 焦距(m)
seeing = 1; % 视宁度(arcsec)
% 生成湍流相位屏
r0 = (0.98*lambda/seeing/4.85e-6)^(6/5); % 相干长度
phase = atmospheric_phase_screen(512, D, r0, 30);
% 模拟星像
star_image = simulate_star_image(phase, D, f);
% 显示结果
figure;
imagesc(star_image);
axis image; colorbar;
title(['模拟星像(seeing=' num2str(seeing) ' arcsec)']);
8.3 显微物镜设计评估
高NA物镜需要精确模拟矢量效应：
matlab复制% 定义物镜参数
NA = 0.95;
M = 40; % 放大率
tube_length = 200; % 筒长(mm)
% 矢量PSF计算
[PSF, OTF] = vector_psf(NA, lambda, M, tube_length);
% 评估分辨率
cutoff_freq = 2*NA/lambda;
disp(['理论截止频率: ' num2str(cutoff_freq) ' cycles/um']);
% 显示PSF
figure;
subplot(1,2,1);
imagesc(PSF); title('矢量PSF');
subplot(1,2,2);
plot(OTF); title('矢量OTF');
9. 性能优化与加速技巧
9.1 并行计算实现
利用MATLAB并行计算工具箱加速光线追迹：
matlab复制% 初始化并行池
if isempty(gcp('nocreate'))
parpool('local', 4); % 使用4个核心
end
% 分割光线数据
num_rays = 1e6;
ray_batches = ceil(num_rays / 1e4);
% 并行追迹
parfor b = 1:ray_batches
start_idx = (b-1)*1e4 + 1;
end_idx = min(b*1e4, num_rays);
batch_results{b} = trace_ray_batch(rays(start_idx:end_idx,:));
end
% 合并结果
results = cat(1, batch_results{:});
9.2 GPU加速
对于大规模计算，GPU可以显著提升速度：
matlab复制% 检查GPU可用性
if gpuDeviceCount > 0
disp('GPU加速可用');
device = gpuDevice;
disp(device);
else
error('未检测到兼容的GPU');
end
% 将数据传输到GPU
gpu_rays = gpuArray(rays);
gpu_lens = gpuArray(lens_params);
% 在GPU上执行计算
gpu_results = arrayfun(@trace_ray_gpu, gpu_rays, gpu_lens);
% 将结果传回CPU
results = gather(gpu_results);
9.3 内存优化
处理大型波前数据时的内存管理技巧：
matlab复制% 使用内存映射文件处理大数据
filename = 'wavefront.dat';
fileID = fopen(filename, 'w');
fwrite(fileID, zeros(4096^2, 'single'), 'single');
fclose(fileID);
m = memmapfile(filename, 'Format', 'single', 'Writable', true);
% 分块处理
block_size = 512;
for i = 1:block_size:4096
for j = 1:block_size:4096
block = m.Data(i:i+block_size-1, j:j+block_size-1);
% 处理数据块...
m.Data(i:i+block_size-1, j:j+block_size-1) = processed_block;
end
end
10. 实用技巧与常见问题
10.1 提高仿真精度的方法
-
光线采样策略：
- 使用准蒙特卡洛(QMC)采样代替随机采样
- 在边缘区域增加采样密度
- 自适应采样：根据前一轮结果调整采样
-
数值稳定性处理：
- 使用高精度算术(vpa)
- 对临界角进行特殊处理
- 添加微小扰动避免奇点
matlab复制% 高精度计算示例
R = vpa('100.1234567890123456789');
n1 = vpa('1.5168');
n2 = vpa('1.0000');
theta = vpa('45')*pi/vpa('180');
% 精确折射计算
theta2 = asin(n1/n2*sin(theta));
disp(['精确折射角: ' char(theta2*180/pi) '°']);
10.2 常见错误与调试
-
光线丢失问题：
- 检查孔径限制
- 验证光线-表面交点计算
- 添加光线追踪可视化
-
能量不守恒：
- 检查反射/透射系数
- 验证吸收处理
- 添加能量平衡检查
-
数值不稳定：
- 减小步长
- 使用更稳定的算法
- 添加数值阻尼
matlab复制% 光线追迹调试示例
function [x_out, debug_info] = trace_ray_debug(x_in, params)
debug_info = struct();
% 记录初始状态
debug_info.initial = x_in;
% 第一步: 表面交点
[x_hit, valid] = find_intersection(x_in, params);
debug_info.hit_point = x_hit;
debug_info.valid = valid;
if ~valid
x_out = NaN;
return;
end
% 第二步: 法线计算
normal = compute_normal(x_hit, params);
debug_info.normal = normal;
% 第三步: 折射/反射
[x_out, is_reflected] = refract_reflect(x_in, normal, params);
debug_info.final = x_out;
debug_info.reflected = is_reflected;
end
10.3 与其他工具集成
-
与Zemax交互：
- 通过DDE或API连接
- 导入/导出ZMX文件
- 批量分析多个配置
-
Python混合编程：
- 使用MATLAB Engine API
- 调用Python光学库(如PyOptics)
- 数据交换通过HDF5或NPY格式
matlab复制% 调用Python示例
pe = pyenv;
if pe.Status == "NotLoaded"
pyenv('Version','C:\Python39\python.exe');
end
% 调用Python光学计算
result = py.scipy.fft.fft2(double(PSF));
% 将结果转回MATLAB
MTF = abs(double(result));
11. 扩展应用与前沿方向
11.1 自由曲面光学仿真
自由曲面光学元件需要特殊处理方法：
matlab复制% 定义自由曲面
[x,y] = meshgrid(linspace(-1,1,100));
z = 0.1*sin(2*pi*x).*cos(2*pi*y) + 0.05*(x.^3 + y.^3);
% 计算法线
[dx,dy] = gradient(z);
normals = -dx + -dy + 1;
normals = normals ./ sqrt(dx.^2 + dy.^2 + 1);
% 光线-曲面交点
function [hit, t] = ray_freeform_intersection(ray, surface)
% 使用牛顿迭代法求交点
tol = 1e-6;
max_iter = 100;
% 初始猜测
t = estimate_initial_t(ray, surface);
for i = 1:max_iter
% 计算残差和雅可比...
% 迭代更新...
if norm(residual) < tol
break;
end
end
hit = ray.origin + t * ray.direction;
end
11.2 超表面光学仿真
超表面的相位调制需要特殊处理：
matlab复制% 定义超表面单元
unit_size = 0.5; % μm
phase_profile = @(x,y) 2*pi*mod(sqrt(x.^2 + y.^2), unit_size)/unit_size;
% 计算远场衍射
lambda = 0.55; % μm
f = 1000; % μm
U = exp(1i*phase_profile(X,Y));
h = fft2(U);
I = abs(h).^2;
% 显示衍射图案
figure;
imagesc(I);
title('超表面衍射图案');
11.3 计算成像系统仿真
计算成像将光学与图像处理结合：
matlab复制% 光学前向模型
object = im2double(imread('target.png'));
PSF = calculate_psf(optical_system);
image = conv2(object, PSF, 'same');
% 重建算法
estimated_PSF = estimate_psf(image);
reconstructed = deconvlucy(image, estimated_PSF, 10);
% 评估重建质量
psnr_val = psnr(reconstructed, object);
ssim_val = ssim(reconstructed, object);
disp(['PSNR: ' num2str(psnr_val) ', SSIM: ' num2str(ssim_val)]);
12. 完整项目实例
12.1 双高斯物镜仿真
双高斯物镜是经典的中等视场设计：
matlab复制% 定义镜头参数
lens = struct();
lens(1).R = 50; lens(1).d = 5; lens(1).n = 1.5168;
lens(2).R = -50; lens(2).d = 2; lens(2).n = 1.0000;
lens(3).R = -30; lens(3).d = 5; lens(3).n = 1.5168;
lens(4).R = 100; lens(4).d = 30; lens(4).n = 1.0000;
lens(5).R = 100; lens(5).d = 5; lens(5).n = 1.5168;
lens(6).R = -30; lens(6).d = 2; lens(6).n = 1.0000;
lens(7).R = -50; lens(7).d = 5; lens(7).n = 1.5168;
lens(8).R = 50; lens(8).d = 0; lens(8).n = 1.0000;
% 光线追迹分析
[spots, wavefront] = analyze_lens(lens, 10); % 10°视场
% 优化设计
opt_vars = optimize_lens(lens, @cost_function);
% 显示优化结果
display_lens_performance(opt_vars);
12.2 激光扩束系统仿真
激光扩束系统需要严格控制波前质量：
matlab复制% 定义扩束系统
M = 5; % 扩束倍数
lambda = 1.064; % μm
input_beam = gaussian_beam(lambda, 1); % 1mm束腰
% 设计望远镜系统
[telescope, beam_path] = design_beam_expander(M, input_beam);
% 分析波前畸变
input_wavefront = measure_wavefront(input_beam);
output_wavefront = propagate_through_system(input_wavefront, telescope);
% 计算Strehl比
strehl = calculate_strehl_ratio(output_wavefront);
disp(['系统Strehl比: ' num2str(strehl)]);
12.3 全景成像系统仿真
全景镜头需要特殊的光学设计：
matlab复制% 定义鱼眼镜头
fov = 180; % 度
f = 1; % 焦距(mm)
image_height = f * tand(fov/2);
% 光线追迹
angles = linspace(0, fov, 100);
for i = 1:length(angles)
ray = [0, 0, sind(angles(i)), cosd(angles(i))];
[x(i), y(i)] = trace_fisheye(ray, lens_system);
end
% 绘制成像图
figure;
plot(angles, sqrt(x.^2 + y.^2));
xlabel('入射角(度)');
ylabel('像高(mm)');
title('鱼眼镜头投影特性');
grid on;
13. 项目部署与自动化
13.1 创建GUI界面
使用MATLAB App Designer创建用户友好界面：
matlab复制% 创建主界面
classdef OpticalSimulatorApp < matlab.apps.AppBase
properties (Access = public)
UIFigure matlab.ui.Figure
LensTable matlab.ui.control.Table
RayTraceButton matlab.ui.control.Button
ResultAxes matlab.ui.control.UIAxes
end
methods (Access = private)
function rayTraceButtonPushed(app, ~)
% 从表格获取镜头数据
lens_data = app.LensTable.Data;
% 执行光线追迹
[spots, wavefront] = simulate_optical_system(lens_data);
% 显示结果
plot(app.ResultAxes, spots(:,1), spots(:,2), 'o');
axis(app.ResultAxes, 'equal');
grid(app.ResultAxes, 'on');
end
end
end
13.2 生成可执行文件
使用MATLAB Compiler打包应用程序：
matlab复制% 创建编译配置
cfg = coder.config('exe');
cfg.CustomInclude = {'optical_simulator.h'};
cfg.CustomSource = {'optical_simulator.c'};
% 定义主函数
main_file = 'optical_simulator_main.m';
% 执行编译
codegen('-config', cfg, main_file);
13.3 自动化测试框架
建立自动化测试确保仿真可靠性：
matlab复制classdef OpticalSystemTest < matlab.unittest.TestCase