简介这份资源面向光学工程、光电技术及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 一般取 50200太少面型会有棱角太多计算量上升且拟合噪声变大。目标面半径 R、光源到目标面距离 H、透镜折射率 n、透镜中心厚度 t 这几个参数决定了整个面型的走向改一个就要重跑全套。2.3 代码包里的模块划分与数据流拿到这套实现包先别急着跑主脚本花五分钟看一下目录结构会省很多事。典型的结构是文件/目录作用关键参数main_design.m主入口串起建模、追迹、评价目标半径、距离、折射率source_model.m光源发光强度模型半强角、朗伯阶数surface_solver.m自由曲面逐点求解环带数 M、迭代容差ray_trace.m蒙特卡洛光线追迹光线数、接收面网格evaluate.m照度均匀度与光斑评价均匀度阈值、评价区域results/计算结果与配光数据照度矩阵、均匀度数值数据流是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.59PMMA 约 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和追迹结果叠在一张图上评价项设计值实测值偏差来源中心照度计算值实测值光源实际光通量均匀度0.850.78扩展光源装配误差光斑半径0.5 m0.53 m折射率批次差异边缘角度60°63°透镜装配偏心偏差在 10% 以内一般可以接受超过就要回头查光源模型和装配间隙。我自己的习惯是每次设计完都强制走一遍“导出点云 → 仿真复核 → 实测对标”三步哪怕项目再急也不跳过因为跳过的那次往往就是后面返工的那次。希望这套 LED 透镜二次光学设计 Matlab 实现包能帮你把配光这件事从玄学变成可复现的流程少熬几个改透镜的夜。本文还有配套的精品资源点击获取