摘 要：提出了一个逐次逼近的计算方法，扩展了对LED自由曲面透镜的计算能力。首先用光通量线方法计算来自LED的n条光线的角度，使它们携带着相同的能量；再由照度分布要求求得像面上n个目标落点的位置。设定一个尝试平面为透镜的最后一面，并任意给定此平面上的n个尝试折射点位置，从而得到n个总的光线偏折角。按给定的透镜各面偏折光线能力的权重把总偏折角分配到各面，得到每个表面前后的光线角度，再用折射定律求得构成各面的n个连续的小面。同时可得到光线的修正方向，用于下一次循环计算。这样逐次逼近计算，直至误差满足预先给定的任意小量的要求。新方法可以用于多种计算，包括像面照度分布事先给定、一个或多个透镜的形状的同步计算、远场或近场照明、照明面是平面或曲面，或者以上这些情况的结合。该方法有很强的计算能力，简单易行，且以一个很快的收敛速度达到相应的精度。

关键词：逐次逼近法；近场；光通量线；二次光学设计；LED；偏折能力

引言

计算LED透镜的形状有多种非成像光学的设计方法[1-5]。其中光通量线(以下称LFR)方法[4]可以计算给定任意照度分布的透镜，而偏折能力分配法[5]则可以同步计算多表面透镜。同时用这两种方法是一个计算透镜的优秀方法，其计算步骤如下：

1) 在极坐标中把LED光强的空间分布曲线分为n份，但各份的光通量必须相等，分割线的方向就是透镜的n条入射线的方向，构成一个n维的入射线角度数列。

2) 因为每条光线光通量相同，像面的照度分布要求可以转换为光线在像面上的分布，从而得到像面上n条光线落点的位置。假定透镜尺度相对很小(远场条件)，得到n个出射线方向。也构成一个n维出射线方向数列。

实际作业中，根据线路排列方式(三角形、水平)移动稳固滑块，调整横杆与立杆顶部的垂直距离满足中相导线与两边相导线的垂直距离后，旋紧第二锁紧件固定住滑块。

3) 透镜出射线和入射线的方向之差称为总偏折角，由求得的上述两个序列后可以求得总偏折角序列。设透镜表面数为m，把此n个总偏折角每个都按各面的权重再分到m个透镜的表面上，这就得到了(m+1)×n个角度，也就是得到了每个面前后的入射角和出射角。

4) 用折射定律对每个表面依次计算各得n个小线段连成一条折线，共m条折线，构成透镜的m个表面。当n很大时就是一个平滑的透镜表面。图1为双面LED透镜m=60、n=50的示例。

其实李世豫一直为潘际銮的成就骄傲。2008年年初，潘际銮完成了全国第一条高铁京津城际铁路的焊接工程，高铁验收时，潘际銮特意请夫人去体验。回来后，夫人的脸上神采飞扬：“验收的时候坐到高铁司机旁边，我心里高兴极了。高铁的速度真快，我们从来没坐过那么快的火车。”

图1 LFR方法计算单透镜

Fig.1 Single lens calculated with LFR method

上述LFR方法假定了LED尺度很小，也就是满足远场条件。但如果像面比较近，即LED的尺寸相比照明尺度不是非常小，不能满足远场条件，计算就会有误差。为了解决这一问题，我们提出了用于透镜设计的逐步逼近的方法，本文着重叙述逐步逼近方法的概念，具体的公式和程序不在本文中出现。

1 逐步逼近法

1)动机。上述LFR方法在处理一般问题时都符合“远场”条件，此时透镜上折射点的位置的不同对像方光分布基本无影响，我们只需关心光通量线的方向而不是折射点的位置。但是在“近场”情况下，光线在透镜上的位置会影响照明的结果。近场问题的处理应该用光线在表面上的位置来表征。但这对透镜形状的计算却带来很大的困难。其原因在于折射定律擅长的只是求解介质表面的方向和折射线的方向之间的关系。在许多领域的实际工作中，求解的方程十分复杂或者是无法直接求解时往往可以使用逐次逼近的方法[6]，此法在大量实际问题中都得到了广泛的应用。由此我们想到，如果能把这种逼近法的思想用到LED透镜的计算中，也许可以解决原来解决不了的求解问题。本文的结果说明这是可以做到的。

2)可行性。处理近场问题的困难在于已知LED发出的入射光线的方向和像面的落点无法解得透镜的形状。但如果能够尝试性地给出出射方向就可以解得初始的形状。于是我们提出先给定一个尝试平面作为透镜的最后一面，并给定该面上各光线的尝试折射点，再按照度分布要求给定LFR在像面的落点，这就可以得到出射光线的尝试方向，这样就可以用原有的LFR方法计算得到初始透镜形状。然而再逐步修正透镜使之逐渐满足“落点”的要求。经过多次逼近计算，可以得到符合折射定律和落点要求的透镜形状。本文的结果说明这是可行的方法。

3)难点。这个方法的难点是如何给出下一次逼近的光线方向。对此我们用最为简单的方法很好地解决了这个问题。就是简单地把尝试计算得到的新的折射点到目标落点的方向作为下一次计算的方向。

2 逐步逼近法计算LED透镜的步骤

2.1 计算入射LFR角度序列

考虑一个旋转对称单透镜，平面像面、均匀照度分布、极少的光线(n=8)。此例并非某个实际设计的例子，只是为了形象地说明逼近计算方法而设定的。

设LED位于坐标原点，透镜前表面中心在A点(x=10 mm)，后表面在B点(x=30 mm)，弧形像面中心在C点(x=70 mm)，见图2。

图2 LED的LFR及光线目标落点

Fig.2 LED’s LFR and target light drop points

按照LRF方法，首先需要求得LED发出的光线的方向序列，这些光线携带相同的能量。由于照度正比于单位面积上的光通量，而光通量正比于光强与立体角的乘积，因此光强分布曲线对总光束角的积分的值将正比于总光通量。根据这一规律，可以把LED发出的光能按照携带的光通量相等而分成8份(以n=8为例)，即把光强分布曲线对总光束角的积分的值分成8等份。它们的角度间隔就是9条具有相同能量的LFR的方向，这便是透镜入射线方向序列。示意性的结果见图2左部原点出发的9条光线簇，大致可以看到这里是大光强处光线密度大，小光强处光线密度小,符合光强空间分布规律。LFR计算的具体方法[4]在此不详述。

2.2 设定像面和光线落点

由于LFR是带有相同能量的光线，因此像面的光线落点的密度分布决定了像面的照度分布，这就不难根据事先给定的照度分布求得LFR的落点沿像面的分布。

在照度均匀分布的情况下，只要将位于D处的被照明平面均匀分成8份即可得到9个落点，得到图2所示的9个小圆形分割点。

对非均匀照度分布时，可以事先给定照度分布函数，再按照落点的密度正比于照度的原则来划分。

2.3 设置尝试平面及尝试折射点

这是逐步逼近法的重点。在透镜的后表面位置C点即x=30处给定一个任意长度例如70 mm的垂线，代表尝试平面直径。在这个垂线上任意给定9个小圆形试探折射点。

图3 入射和出射的LFR角度序列

Fig.3 Angles sequence of incident and exit LRF

把图3中C平面和D平面上的9×2个小圆一一相连，连线(图3中的点划线)的方向代表光线从LED最后一个表面出发方向，构成了透镜出射线方向序列。这样就得到了透镜的入射线和出射线两个尝试角度序列。

有了这两个角度序列，我们就可以用LFR的方法进行计算求得透镜的形状了。对于单透镜的情况，设两面的折射能力相等，则可求得上两个角度的平均值，就得到第三个角度序列，那就是在透镜内部的光线角度序列(图3中没有画出这一个序列)。有了这9×3个方向，我们就可以顺序用折射定律计算得到透镜前后表面共8×2个小面的角度和位置，结果得到了一个弯月形的透镜，见图4的粗实线。

图4 第一次计算结果和修正角度

Fig.4 Calculation results after first cycle and correction angles for next cycle

2.4 逐次逼近法计算透镜形状

可以看到，第一次循环后试探平面自动变成了曲面，从其上的折射点发出的光线(实线)的确和试探方向(点划线)精确平行(这正是LRF方法的结果)，但是却没有射向希望的落点。实际落点与目标落点之差最大达13.6 mm。逐步逼近法可以迅速缩小这一偏差，但如何选择下一次计算的修正量是逼近法成功与否的关键。

由图1可知，S8与S13集中在第二象限，它们的豆味(f)最为突出；S1，S2，S3，S4，S14，S15和S16聚集在第一、四象限，它们的胡椒(a)、胡椒(f)、花椒(a)、辣椒(a)、烧烤(a)、甜(f)和辣(f)等风味特征比较突出；而S5，S6，S7，S9，S10，S11和S12的鲜(f)比较突出。由图1可知，16种怪味胡豆样品彼此位置相距较远，说明它们的风味特征存在差异，这与方差分析的结论一致。

我们的解决方法是简单地用第一次计算得到的透镜后表面的光线落点和希望的像面落点(图3中D面的小圆形)的连线(图3中虚线)作为第2次计算的像方光线的方向，进行第2次计算，第2次计算结果如图5所示。这时，实际的光线落点和希望落点之差Δmax=0.14 mm，比第一次的结果缩小了约100倍。

图5 第二次计算结果

Fig.5 Calculation results for second cycle

继续上述计算，在本例中经过4次逼近计算后就得到了理想的结果，如图6所示，最大误差仅为Δmax = 0.003 mm。

对照组运用红霉素治疗方式，红霉素（唐山红星药业有限责任公司，国药准字H13022717）口服治疗，3次/d，连续使用10 d作为一个疗程，连续使用3个疗程

图6 最终结果

Fig.6 Final result

增加LFR的数目，如n=40，5次循环后的结果如图7所示。这时Δmax=0.005 mm，曲线也较为平滑了。实际设计中一般n的值为200左右就可以得到可以用于加工的平滑的表面。

图7 n=40时的计算结果

Fig.7 Calculation result of n=40

以上用直观的图形较清晰地给出了本方法的内容，方法中完全不涉及复杂的微分方程等计算，这就不难编写计算机程序。

The non-monotonic behavior of ημd during periods I and IV indicates changes of the predominance of the two ion sources. From Fig. 3 we can state that this phenomenon is important whenever the conductive filament exhibits most of its largest extension because oxygen ions are lacking throughout it.

实践表明，试探面的形状、大小和其上的落点的位置分布的宽容度很大，即使初始给出的是一个非平面的尝试表面和非均匀的尝试点分布，用逐步逼近法都将趋向于同一个结果。

采用HPLC方法来定量分析底物咖啡碱、可可碱、茶碱在发酵过程中的消耗情况，反应液原液稀释至合适倍数后经HPLC进样分析，三种生物碱的定量检测方法采用色谱条件[14,15]：C18柱（4.6×200 mm，5 µm，welchrom），流动相：A：去离子水，B：N, N-二甲基甲酰胺：甲醇：乙酸=40：2：1.5，流速：1 mL/min。进样量：10 μL等条件在278 nm下检测，试验重复六次（n=6），分析并计算结果。

3)经液氮冷浸煤岩损伤机制分析，在含水量较低的情况下，水的冻结速度较缓慢，影响机制以静水压理论与分凝冰理论为主，在煤岩饱水情况下，液氮对煤岩的影响以水冰相变造成的体积膨胀为主，造成对煤基质的损伤。

山东煤机装备集团经过多年理论研究和实验总结出一套先进的粉煤灰湿法选炭工艺，如图1所示。该工艺主要采用了本公司生产的FWX系列浮选柱、BXN系列高频振动斜板浓缩机、DU系列水平带式过滤机、矿浆预处理器、搅拌桶和滚筒筛等设备。

3 非均匀照明的逼近法设计

以上的例子是均匀照度分布。对于非均匀照度分布，只要给出像面照度分布函数也可以进行计算。例如，给定一个中心照度小而边缘照度大的正弦函数分布a+b·sin(x+c),其中a=31、b=30、c=-90°，其照度分布函数如图8所示。

Study on the variation of direct shearing test of pile-soil interface on excess pore water pressure in clayey soil

图8 正弦函数给定为像面照度分布

Fig.8 Sine function is given as the illumination distribution on image surface

计算得其透镜形状和LFR分布为图9所示。其光线落点的密度沿弧面的分布是按图8的正弦分布由中心向边缘变大。

图9 按照事先给定照度分布的计算结果

Fig.9 Calculation results with illumination distribution given in advance

对于其他类型的照度分布，只要给出分布函数，可以套用上面的方法。好在实际的工程中都不需要用复杂的照度分布，因此给出函数也不是难事。

4 被照面为非平面的逼近法设计

考虑到在近距离照明中有时会遇到像面不是平面的情况，用平面代替弧面必然会带来误差，而用我们的逼近方法计算却可以容易地处理非平面像面的问题。

被照明面是一个曲面并没有给我们的计算带来任何困难。只要在这个曲面上按照照度分布的要求求得一个个光线的落点的坐标就可以了。唯一的困难是在数学上描述这个曲面，实际工作中这也不是难事。

以下用一个球表面均匀照度分布为例。假定其圆弧可以用下列参数方程表示：

x(x0,R,t)=x0+R·cos(t)，y(y0,R,t)=y0+R·sin(t)

其中圆心的坐标为(x0,y0)，R为半径，t为辐角。

在本例中，被照明面是中心在C点(x=70 mm)，半径为130 mm的凸弧形表面，圆心在x=200 mm，y=0。

将此弧面均分42份，得到图 10所示的43个分割点。若希望非均匀照度分布，则需要按照落点密度正比于照度的原则求得落点位置。

“风雨说”以黄进为首提出的，他们认为广东丹霞地貌及北国山山顶的上凹坑形成主要是由风化、雨水的滴蚀、溶蚀作用下产生的，称之为“风化、雨滴溶蚀说”，简称“风雨说”。

图10 非平面的被照明面

Fig.10 Non-Plane Image Surface

对一个可以用参数方程表示的二次曲面，逼近方法同上。其他形状也是用参数方程表示方便。

5 多透镜的逼近法设计

图11 LFR方法计算双透镜

Fig.11 Double lens calculated with LFR method

LRF方法中计算多个透镜的方法同样可以用在本方法中。某些LED透镜若做成双透镜会比单透镜更为合适[5]。SSWYD-800型无影灯[7]有一款透镜就做成了双透镜。其光线较少时的计算结果见图11。为了显示清楚些，图形用了n=20的计算结果。图11为无影灯的工作距离为1 000 mm时，用文献[7]的总体设计方法，用逐步逼近法同步计算的LED无影灯的双透镜四个面的计算结果。我们也用文献[4]的方法按远场条件进行了计算，得到结果一致，这是因为此时满足远场条件了。

6 结论

我们用逼近法计算LED自由曲面透镜，从建立一个虚拟的尝试表面并给定尝试点出发，给出了修正计算的方法，进行多次循环计算，避开了求解方程等复杂的问题，不但可以得到很高的计算精度，还具有快速的收敛速度。

从计算的角度看，本文的方法实际上解决了一个从“方向”转为计算“落点”的问题，以前是由计算光线方向来求透镜形状[4，5]，而现在近场问题则是要由计算光线落点来求透镜形状。

逼近法把计算透镜的LFR方法扩展到既可以用于近场计算，又可以用于远场计算，既可以计算平面像面又可以计算非平面像面。同时逼近法保留了原方法的几个重要的优点，如可以用于非均匀的照度分布或事先给定的任意像面照度分布的计算，以及同步计算双透镜系统。这给LED透镜二次光学设计带来了方便。

现实的就医环境中，由于很多人缺乏对发烧这一症状的常识性认识，这就给滥用激素提供了市场。人们在就医时常常对医生说：“我没时间养病，我不能休息，我得把烧快点压下去。”鉴于当前紧张的医患关系，不少医生为迎合病人的这一需求，就会给病人输注激素强行退烧。于是地塞米松被当作退烧药加入输液中便成了十分普遍的现象。实际上，地塞米松不是退烧药，它属于糖皮质激素的一种，也就是老百姓常说的激素。激素能起到抑制免疫的作用，免疫系统不工作了，烧也就退了，但它同时阻碍了人体自身对抗疾病的能力，使得病菌趁机生长繁殖引发新的感染，也就是说，容易使人的体质变差，可能会导致严重感染。

7 讨论

1)设想被照明面位于x=1 000 m远处，则其上1 mm的落点位置误差其角度差别只有0.0001°，可以让程序使用12位有效数字，角度精度可以非常高，则即使多次运算也不会带来有影响的误差。图11说明本方法在非近场的情况同样适用。

2)新方法的必要条件还在于这种逼近能否收敛。实际工作证明此法可以收敛，而且收敛速度快，我们计算的多个例子都是循环次数在10以下时计算误差小于0.001 mm。

3)逐步逼近方法可以推广到其他难以一次性直接计算的问题中去，如一个面既有反射又有透射的LED照明系统[8，9]。此时，一般方法难以求解得到曲面的形状，但逼近法可以计算得到。

参考文献

[1] WANG Lin, QIAN Keyuan,LUO Yi. Discontinuous free-form lens design for prescribed irradiance [J]. Applied Optics, 2007, 46(18)：3716-3723.

[2] Harald Ries.Tailored freeform optical surface [J]. Journal of the Optical Society of America， 2002, 19(3)： 590-595.

[3] 丁毅，顾培夫.实现均匀照明的自由曲面反射器 [J]. 光学学报，27(3)：540-544.

[4] 周士康，陈春根，许礼，等. 光通量线方法用于LED二次光学设计[J]. 照明工程学报，2016, 27(1)：101-111.

[5] 周士康，陈春根.LED多表面LED二次光学透镜设计的偏折力分配法[J]. 照明工程学报，2018, 29(2)：47-51.

[6] 赵根榕.逐次逼近法[M]. 北京：科学普及出版社，1965.

[7] ZHOU Shikang, CHEN Chungen. Overall Optical Design of LED Surgical Luminaires of Variable Light Spot.Optics，2019，8(1):1-6.

[8] 周士康，陈春根. 具有反射透射两用表面的LED自由曲面透镜的逼近法计算[J]，照明工程学报，2017, 28(3)：97-101.

[9] 周士康，崔小红，程德诗. 大功率LED白光准直光源透镜的设计与应用[C]. 第十三届全国LED产业发展与技术研讨会，2012.

Successive Approximation Calculations of LED Lens

ZHOU Shikang, DU Jin, CHEN Chungen

(Shanghai Sansi Electronic Engineering Co. Ltd., Shanghai 201199, China)

Abstract：A successive approximation method is proposed in this paper, which greatly extends the computational power of the LED lens shape with free form surface. In this method, the angles of the n rays from the LED are calculated firstly using the LFR (Light Flux Ray) method, so that they carry the same energy. The positions of the target n drop points of LFR on the image surface can be obtained according to the requirements of the illumination distribution. An attempt plane is set as the last surface of the lens and the positions of n attempt refraction points on this plane are given arbitrarily, so as to get the n total deflection angles of each ray. These angles are distributed to each surface according to the weight of the deflection light ability given previously, then the ray angles before and after each surface are obtained. Then, by using the law of refraction, the consecutive n facets are obtained which constitute the lens surface. At the same time, the corrected ray directions can be obtained for next cycle calculation. So the successive approximation is computed until the error satisfies any small amount specified in advance. This new method can be used for a wide variety of calculation, including the illumination distribution was given beforehand, the shape with one or more lens are calculated synchronously, the lighting surface location is near or far field, the shape of the lighting surface is flat or curved etc.， or their combination of the above. Therefore it has very strong calculating ability and very wide scope of applications. This method is simple and easy to be achieved, and any high computational accuracy can be achieved with a fast convergence speed.

Key words:successive approximation method; near field; light flux ray; secondary optical design; light emitting diode; deflection light ability