LED透镜二次光学设计Matlab实现：自由曲面建模与光线追迹

简介：这份资源面向光学工程、光电技术及LED封装应用方向的本科生与硕士生，聚焦LED照明系统中的透镜二次光学设计问题，提供一套可直接运行的Matlab实现方案。代码基于光线追迹或非成像光学原理，完成自由曲面或旋转对称透镜的光学参数建模与光强分布优化，输出照度分布图与光效评估数据等典型结果，适配matlab2014a、2019a、2021a多个版本。压缩包共7个文件，约5KB，以m主程序、txt计算结果与说明文档为主，另含py脚本及少量工程配置文件，分别承担核心计算、结果输出与使用指引等用途，下载后无需额外配置即可运行查看设计效果。目前已有63人学习关注。对于开展课程设计、毕业设计或基础科研验证的读者，可借助主程序理解二次光学建模流程，对照输出文件核对光强与照度分布，并参考说明文档快速定位参数含义，从而完成从建模到仿真评估的完整闭环，适合作为LED配光需求下的入门实践参考。

1. 从一次配光翻车说起：这套 LED 透镜二次光学设计 Matlab 实现包到底能干什么

做 LED 灯具结构或光学设计的同行大概率都经历过这种场面：透镜供应商给了一条配光曲线，装到整灯上实测却完全不是那么回事，光斑要么中间发暗、要么边缘出现一圈黄环，客户拿着照度计在测试台上盯着你看。问题往往不在透镜本身，而在二次光学设计阶段没有把光源模型、透镜面型和目标配光三者之间的映射关系算清楚。这套 LED 透镜二次光学设计 Matlab 实现包，就是冲着这个环节来的——它把自由曲面透镜的建模、光线追迹、配光评价串成了一条可复现的 Matlab 流程，附带代码和计算结果，适合做 LED 照明模组、路灯透镜、投光灯配光以及背光匀光的工程师直接拿来改参数跑。

它解决的不是“画一个透镜”的问题，而是“给定目标配光，反推透镜面型并验证”的问题。适合两类人：一类是刚接手二次光学任务、手里只有光源手册和照度目标的新手，能照着代码把流程走通；另一类是做过多款透镜、想把手算经验固化成脚本的老手，可以拿它当参数扫描和方案对比的底座。下面按“原理怎么立住 → 代码怎么跑 → 坑在哪 → 怎么进阶验证”的顺序拆开讲。

2. 二次光学设计的数学底子：从朗伯光源到自由曲面映射

2.1 为什么二次光学不能靠试错堆透镜

LED 芯片本身大多是朗伯或近朗伯发光，半强角在 120° 左右，直接出光会形成一个又宽又散的圆斑，照度均匀度差、眩光高。一次光学（封装透镜）只能把光收一收，真正决定整灯配光形状的是二次光学元件——也就是我们常说的 TIR 透镜、反光杯加透镜、或者自由曲面透镜。二次光学的核心任务，是把光源的发光强度分布 I(θ,φ) 通过透镜两个面的折射和全反射，映射成目标面上的照度分布 E(x,y)。

如果靠“改一个面型参数、跑一次仿真、看一眼光斑”的试错法，一个方案迭代十几轮很正常，而且改到后面参数之间互相耦合，根本说不清哪个面在起作用。所以合格的做法是先建立能量映射关系，再解面型。常见做法是采用非成像光学的边缘光线理论和网格划分法：把光源按角度分成 N 份，把目标面按面积分成 N 份，让每一份光源能量对应一份目标面积，从而得到一组入射-出射光线对，再逐点构建自由曲面。

2.2 能量映射与 Snell 面型构建的基本方程

设光源为点光源近似，发光强度按朗伯分布 I(θ)=I0·cosθ。目标面是一个半径为 R 的圆形区域，要求照度均匀。把光源半角 θmax 内分成 M 个环带，每个环带能量相等：

E_ring = ∫ I(θ) sinθ dθ dφ

目标面同样分成 M 个等面积环带。第 i 个环带对应的出射角 θi 由能量守恒确定。得到 θi 后，用 Snell 定律在透镜前表面逐点求解面型坐标。前表面一般取平面或球面，后表面为自由曲面，光线在后表面折射后射向目标环带。Matlab 里实现时，通常用离散点云加样条拟合的方式生成面型，再导入光学仿真软件做验证。

这套包里代码的核心就是两个循环：外层遍历环带，内层用 Snell 方程解出每个采样点的法向量，再积分得到面型坐标。参数上，M 一般取 50～200，太少面型会有棱角，太多计算量上升且拟合噪声变大。目标面半径 R、光源到目标面距离 H、透镜折射率 n、透镜中心厚度 t 这几个参数决定了整个面型的走向，改一个就要重跑全套。

2.3 代码包里的模块划分与数据流

拿到这套实现包，先别急着跑主脚本，花五分钟看一下目录结构会省很多事。典型的结构是：

数据流是：source_model 生成光源采样 → surface_solver 解出面型点云 → ray_trace 把面型当黑匣子做追迹 → evaluate 输出均匀度和光斑图。理解这条链路，后面改参数才知道改哪一层。

3. 把代码跑起来：Matlab 环境配置与主流程实操

3.1 环境准备与依赖检查

这套代码对 Matlab 版本要求不算苛刻，R2018b 及以上基本都能跑，但要用到曲线拟合和插值函数，建议装 Curve Fitting Toolbox 和 Optimization Toolbox。如果只做基础面型求解，Optimization Toolbox 可以不用，但参数扫描那部分会受限。常见做法是先在命令行检查：

% 检查关键工具箱是否可用 ver('curvefit') % 曲线拟合，面型样条需要 ver('optim') % 优化，参数扫描需要 ver('parallel') % 并行，光线追迹加速可选

如果某个工具箱缺失，代码里对应函数会报 undefined，这时候要么装工具箱，要么把相关调用换成基础函数实现。注意不要直接注释掉拟合步骤，面型点云不拟合直接追迹，光斑会出现明显锯齿。

3.2 主脚本参数设置与运行

打开 main_design.m，最先看到的是参数区。以下是我一般会先确认的几个：

% ===== 设计输入参数 ===== R_target = 0.5; % 目标面半径，单位 m H_distance = 3.0; % 光源到目标面距离，单位 m n_lens = 1.59; % 透镜材料折射率，PC 约 1.59，PMMA 约 1.49 theta_max = pi/3; % 光源半角，60 度 M_rings = 120; % 能量环带数 N_rays = 200000; % 追迹光线数 center_thick = 0.008; % 透镜中心厚度，单位 m

参数含义和调整方向：R_target 和 H_distance 决定了配光角度，R/H 越大光束越宽；n_lens 必须和实际材料一致，填错会导致面型整体偏移；theta_max 要覆盖光源实际能量分布，设小了边缘能量被截断，光斑边缘会突然变暗；M_rings 从 80 起调，看面型是否光滑；N_rays 影响追迹噪声，20 万条在普通工作站上跑几分钟，要更快可以降到 5 万先看趋势。

设置完直接运行：

% 运行主设计流程 run('main_design.m')

运行后会在 results 目录下生成照度矩阵和配光曲线数据。如果中途报错，先看是哪个模块，常见的是 surface_solver 里迭代不收敛，这时候把 M_rings 降下来或者放宽容差再试。

3.3 面型求解与光线追迹的关键代码段

面型求解的核心循环大致是这样：

% 逐环带求解自由曲面点坐标 for i = 1:M_rings % 当前环带对应的目标半径 r_i = R_target * sqrt(i / M_rings); % 由能量守恒反解出射角 theta_out = acos(1 - (i / M_rings) * (1 - cos(theta_max))); % 用 Snell 定律解后表面法向量 [normal_vec, point_new] = solve_snell( ... point_old, theta_out, n_lens); % 存储点云 surface_points(i, :) = point_new; point_old = point_new; end % 样条拟合生成光滑面型 surface_spline = fit(surface_points(:,1), surface_points(:,2), 'smoothingspline');

逻辑说明：每个环带先算目标半径 r_i，注意这里用 sqrt 是因为面积和半径平方成正比，等面积划分必须开方；theta_out 由能量守恒反解，保证每个环带能量相等；solve_snell 是自定义函数，输入上一点、出射角和折射率，输出新点坐标和法向量。参数上，point_old 的初值决定面型起点，一般取透镜中心厚度对应的点。拟合那一步的 smoothingspline 平滑参数可以调，太平滑会丢失面型细节，太紧会引入噪声。

追迹部分用蒙特卡洛：

% 蒙特卡洛光线追迹 illuminance = zeros(grid_size, grid_size); for k = 1:N_rays % 按朗伯分布采样光源方向 theta = acos(1 - rand * (1 - cos(theta_max))); phi = 2 * pi * rand; % 追迹到目标面 [hit_x, hit_y] = trace_ray(theta, phi, surface_spline, H_distance); % 累加到照度矩阵 illuminance = accumulate(illuminance, hit_x, hit_y); end % 归一化并计算均匀度 uniformity = calc_uniformity(illuminance);

这段代码里 rand 采样保证光源方向符合朗伯分布，trace_ray 内部做折射和全反射判断，accumulate 把落点映射到网格。N_rays 越大均匀度数值越稳，但超过 50 万条后提升有限，性价比下降。

4. 避坑与排查：二次光学 Matlab 实现里最容易翻车的五件事

4.1 光斑出现同心圆环

现象：追迹出来的照度图有一圈一圈的明暗环，像水波纹。原因：环带数 M_rings 太少，或者面型拟合时平滑过度，导致相邻环带之间出现台阶。解决：把 M_rings 提到 150 以上，拟合改用 'thinplateinterp' 或降低平滑参数，同时检查 solve_snell 里法向量是否做了归一化。

4.2 均匀度始终上不去

现象：反复调参数，均匀度卡在 0.6 左右，目标面边缘照度明显低于中心。原因：目标面划分用了等半径而不是等面积，能量映射本身就错了。解决：确认 r_i 计算里有没有 sqrt，等面积划分必须开方；另外检查 theta_max 是否覆盖了光源 95% 以上的能量，截断太多边缘必然暗。

4.3 面型出现自交或尖点

现象：拟合出的透镜面型在某些位置交叉，或者出现尖锐凸起，加工不出来。原因：Snell 迭代步长太大，或者折射率填错导致出射角解算异常。解决：把环带数加大、步长减小；核对 n_lens 和实际材料一致；在 solve_snell 里加一个角度合法性判断，出射角超过全反射临界角时报警而不是硬算。

4.4 追迹结果和仿真软件对不上

现象：Matlab 算出来均匀度 0.85，导入光学仿真软件跑出来只有 0.7。原因：Matlab 里用的是点光源近似，实际 LED 有发光面尺寸；另外透镜厚度和装配间隙没算进去。解决：把光源模型从点光源改成扩展光源，采样时加一个面内随机偏移；把透镜中心厚度和实际装配空气隙补进追迹距离。常见做法是留 10% 的均匀度余量给仿真和实测差异。

4.5 运行速度慢到无法迭代

现象：改一个参数跑一次要十几分钟，参数扫描根本跑不动。原因：N_rays 设太大，且循环里用了大量临时变量和未预分配的数组。解决：先降 N_rays 到 5 万看趋势，确认方向后再用 20 万精算；把 illuminance 矩阵预分配；追迹循环用 parfor 并行，前提是装了 Parallel Computing Toolbox。注意并行时随机数种子要设好，否则每次结果略有波动。

5. 进阶用法：参数扫描、面型导出与实测对标

5.1 用参数扫描找最优折射率和厚度组合

单跑一组参数只能得到一个方案，实际项目里往往要在 PC 和 PMMA 之间选，还要权衡中心厚度和模具成本。我一般会写一个外层扫描脚本，把 n_lens 和 center_thick 做成网格，每组跑一次追迹，记录均匀度和光斑半径：

% 参数扫描：折射率 x 中心厚度 n_list = [1.49, 1.59, 1.65]; t_list = [0.006, 0.008, 0.010]; results = zeros(length(n_list), length(t_list)); for a = 1:length(n_list) for b = 1:length(t_list) params.n_lens = n_list(a); params.center_thick = t_list(b); [uniformity, spot_radius] = run_design(params); results(a, b) = uniformity; end end % 输出最优组合 [max_u, idx] = max(results(:)); [best_a, best_b] = ind2sub(size(results), idx); fprintf('最优折射率 %.2f，中心厚度 %.3f m，均匀度 %.3f\n', ... n_list(best_a), t_list(best_b), max_u);

这段脚本的价值在于把“凭经验选材料”变成“看数据选材料”。参数上，n_list 按实际可选材料填，t_list 按模具和注塑能力填，不要填理论上最优但做不出来的值。扫描完记得把每组的面型点云存下来，后面做对比图用。

5.2 面型点云导出给加工和仿真

Matlab 算出的面型点云要能落到实际加工，必须导出成通用格式。常见做法是导出为 CSV 或 STL：

% 导出面型点云为 CSV csvwrite('results/surface_points.csv', surface_points); % 导出为 STL 供 3D 查看 tri = delaunay(surface_points(:,1), surface_points(:,2)); stlwrite('results/lens_surface.stl', tri, surface_points);

导出后注意两点：一是坐标单位统一成毫米，加工端一般用毫米；二是点云密度要够，太稀 STL 会有明显三角面。导出前把点云按弧长重采样一遍，分布更均匀。

5.3 用实测配光反推验证设计偏差

设计做完不算完，实测对标才是闭环。拿照度计或配光机测出实际光斑，把实测数据读进 Matlab，和追迹结果叠在一张图上：

偏差在 10% 以内一般可以接受，超过就要回头查光源模型和装配间隙。我自己的习惯是每次设计完都强制走一遍“导出点云 → 仿真复核 → 实测对标”三步，哪怕项目再急也不跳过，因为跳过的那次往往就是后面返工的那次。希望这套 LED 透镜二次光学设计 Matlab 实现包能帮你把配光这件事从玄学变成可复现的流程，少熬几个改透镜的夜。

本文还有配套的精品资源，点击获取