董菁菁，郭震寧，林介本，黃智煒，黃明波，李建功

（華僑大學(xué) 信息科學(xué)與工程學(xué)院，福建 廈門 361021）

發(fā)光二極管輻射模型的快速建立與驗(yàn)證

董菁菁，郭震寧，林介本，黃智煒，黃明波，李建功

（華僑大學(xué) 信息科學(xué)與工程學(xué)院，福建 廈門 361021）

基于Source Property Generator軟件，根據(jù)配光曲線圖快速獲得相應(yīng)的光源數(shù)據(jù)，提出運(yùn)用非線性回歸分析阻尼最小二乘（L－M）法，由模擬公式得到的光源數(shù)據(jù)與由配光曲線（或測(cè)量）獲得的光源數(shù)據(jù)之間的誤差，以尋求模擬公式參數(shù)的變化步長(zhǎng)；然后，由誤差均方根與歸一化互相關(guān)兩個(gè)判斷標(biāo)準(zhǔn)來(lái)重建輻射模型的精確性.結(jié)果顯示：誤差均方根低于1%，歸一化互相關(guān)接近100%.該方法只需提供光源的配光曲線圖即可以獲得相應(yīng)的光強(qiáng)分布函數(shù)和輻射模型，如果前面的函數(shù)設(shè)置沒(méi)有出現(xiàn)錯(cuò)誤，一般只要一次就能通過(guò)驗(yàn)證；而且只要改變某些要求參數(shù)，即可便捷地設(shè)計(jì)出一系列同類型透鏡的輻射模型函數(shù).

光學(xué)設(shè)計(jì)；發(fā)光二極管；輻射模型；阻尼最小二乘法

目前，對(duì)發(fā)光二極管（LED）真實(shí)光源設(shè)計(jì)的需求不斷加強(qiáng)，不僅體現(xiàn)在近場(chǎng)、中場(chǎng)、遠(yuǎn)場(chǎng)方面［1－4］，而且在制作發(fā)光二極管透鏡（包括一次光學(xué)設(shè)計(jì)與二次光學(xué)設(shè)計(jì)）或燈具時(shí)，需根據(jù)不同光源的輻射模型來(lái)設(shè)計(jì)透鏡或燈具以達(dá)到所需的照明效果［5－7］.然而，由于發(fā)光二極管芯片類型和發(fā)光二極管透鏡形狀的多樣性，導(dǎo)致建立其精確的輻射模型要比建立傳統(tǒng)光源輻射模型更加困難.發(fā)光二極管光源廠家提供的一般是發(fā)光二極管的配光曲線，較少提供相應(yīng)的光源數(shù)據(jù)和輻射模型.發(fā)光二極管光源數(shù)據(jù)可以通過(guò)CCD對(duì)發(fā)光二極管光源進(jìn)行檢測(cè)獲得［8－10］，然后，將光源數(shù)據(jù)導(dǎo)入光學(xué)仿真軟件以獲得發(fā)光二極管輻射模型 .但這不能確切地獲取其輻射模型函數(shù)（也稱強(qiáng)度分布函數(shù)），給發(fā)光二極管光源模擬、發(fā)光二極管用作光纖傳輸?shù)男盘?hào)光源［11］、透鏡設(shè)計(jì)、燈具制作，以及發(fā)光二極管陣列設(shè)計(jì)等帶來(lái)困難［2－4］.因此，快速獲得發(fā)光二極管的光源數(shù)據(jù)并建立其輻射模型及相應(yīng)函數(shù)是十分必要的.MORENO等［12］提出根據(jù)廠商提供的光源數(shù)據(jù)，通過(guò)逐個(gè)改變函數(shù)參數(shù)的方法獲得發(fā)光二極管輻射模型函數(shù)，其誤差均方根（RMS）在1%～5%之間，歸一化互相關(guān)（NCC）在99.84%～99.97%之間.本文為解決發(fā)光二極管輻射模型難建立的問(wèn)題，提出一種快速、有效地建立發(fā)光二極管輻射模型的方法.

1 設(shè)計(jì)原理

為了精確地演示發(fā)光二極管（LED）的光強(qiáng)分布和空間輻射模型，必須考慮光發(fā)射面（芯片、芯片陣列、熒光面）經(jīng)反射杯反射和封裝透鏡折射的光線分布.在數(shù)學(xué)上，這種模型可用若干個(gè)高斯函數(shù)或余弦函數(shù)來(lái)描述［12］.

當(dāng)光源輻射模型旋轉(zhuǎn)對(duì)稱時(shí)，發(fā)光二極管的強(qiáng)度模型為若干個(gè)高斯函數(shù)的總和，有

或者為若干個(gè)余弦函數(shù)的疊加，有

式（1），（2）中：I（θ）為強(qiáng)度，一給定方向上單位立體角的光通量；θ為以發(fā)光二極管為中心的坐標(biāo)系的極角；g1i，g2i，g3i或c1i，c2i，c3i為函數(shù)的參數(shù).

對(duì)于旋轉(zhuǎn)非對(duì)稱封裝導(dǎo)致在兩垂直方位角方向上具有不同輻射模型的發(fā)光二極管，其強(qiáng)度角變化模型可為

或余弦函數(shù)為

式（3），（4）中：φ為方位角；θ為以發(fā)光二極管為中心的坐標(biāo)系的極角；g1i，g2i，g3i，g4i或c1i，c2i，c3i，c4i為函數(shù)的參數(shù).

式（1）～（4）對(duì)大多發(fā)光二極管是相當(dāng)精確的，可運(yùn)用于建立普通發(fā)光二極管光源模型.高斯函數(shù)或余弦函數(shù)的個(gè)數(shù)主要取決于發(fā)光二極管輻射模型的形狀.因此，采用誤差均方根（RMS）和歸一化互相關(guān)（NCC）［10］兩個(gè)標(biāo)準(zhǔn)來(lái)判斷所建模型的精確性.即

式（5）～（6）中：I（θj，φk）m，I（θj，φk）c分別為在（θj，φk）點(diǎn)處模擬與實(shí)驗(yàn)測(cè)量得到的光強(qiáng)值；ˉIm，ˉIe分別為模擬和實(shí)驗(yàn)測(cè)量得到的平均光強(qiáng)值；j為以發(fā)光二極管為中心的坐標(biāo)系的極角θ的取樣個(gè)數(shù)；k為發(fā)光二極管方位角的取樣個(gè)數(shù)；N 為總的（θj，φk）點(diǎn)個(gè)數(shù).

通常RMS用來(lái)度量曲線擬合的好壞，即RMS越小，說(shuō)明曲線擬合得越好，參數(shù)選得越合理，誤差RMS一般小于5%［13］.NCC表示模擬值與測(cè)量值的相似性，對(duì)于大多數(shù)運(yùn)用，其值高于99%［3］.

2 設(shè)計(jì)方法

發(fā)光二極管有很多種不同的輻射模型，以致其很難找到某一公式去確切定義.因此，基于Source Property Generator軟件，根據(jù)配光曲線獲得相應(yīng)的光源數(shù)據(jù)，運(yùn)用相應(yīng)的理論對(duì)其非線性回歸分析，最后將所得到的光源模型進(jìn)行實(shí)驗(yàn)驗(yàn)證.文中以對(duì)稱型高斯函數(shù)的模型為例.

2.1 函數(shù)個(gè)數(shù)的決定與參數(shù)的選取

運(yùn)用方程（1），令i＝1，g＝［g11，g21，g31］＝［1，0°，30°］，且保持其他參數(shù)不變，只改變某一參數(shù)，其光強(qiáng)分布如圖1所示.

圖1 不同參數(shù)的光強(qiáng)分布示意圖Fig.1 Intensity distribution schematic of different parameters

由圖1（a）可知：g11決定模型的幅值，隨著g11的增加，幅值隨之增大；由圖1（b）可知：當(dāng)g21為正時(shí)，g21決定模型的最大強(qiáng)度值處的輻射角（φ）；隨著｜g21｜的增加，兩個(gè)最大峰分得越開；當(dāng)g21為負(fù)時(shí)，其模型為g21為正時(shí)的補(bǔ)集.由圖1（c）可知：g31決定模型的半強(qiáng)度角，隨著g31的增加，輻射角的有效寬度增大.綜合圖1并根據(jù)已知的配光曲線來(lái)決定高斯函數(shù)或余弦函數(shù)的個(gè)數(shù)，同時(shí)得到參數(shù)g11，g21，g31的初始值.

2.2 光源數(shù)據(jù)的獲得

將配光曲線圖復(fù)制到Source Property Generator軟件中，根據(jù)配光曲線圖選取相應(yīng)的坐標(biāo)系；然后，開始描點(diǎn)，最終獲得輻射角及其相應(yīng)的光強(qiáng)值.

2.3 模型的建立

根據(jù)函數(shù)個(gè)數(shù)與參數(shù)的初始設(shè)定，運(yùn)用非線性回歸分析阻尼最小二乘（Levenberg－Marquart，簡(jiǎn)稱L－M）法［14］，由模擬公式的光源數(shù)據(jù)與根據(jù)配光曲線圖獲得的光源數(shù)據(jù)（又稱測(cè)量數(shù)據(jù)）之間的誤差來(lái)尋求模擬公式參數(shù)變化的步長(zhǎng)；同時(shí)，由RMS與NCC兩個(gè)判斷標(biāo)準(zhǔn)來(lái)重建輻射模型的精確性.

3 算法流程

阻尼最小二乘法可用下式求搜索方向，有

式（7）中：J＝J（x），f＝f（x），而f（x）為模擬數(shù)據(jù)與測(cè)量數(shù)據(jù)的差值；f（x＋h）＝f（x）＋J（x）h＋O（‖h‖）2；J∈Rm×n雅克比函數(shù)，（J（x））i，j＝（x）；hl，m為步長(zhǎng)；μ 為阻尼因子，它可以控制hl，m的大小方向.當(dāng)μ＝0時(shí)，即為Gauss－Newton法；而當(dāng)μ趨于無(wú)窮時(shí)，I趨于零矢量，即為最速下降法.因此，只要給一個(gè)足夠大的μ，f（x＋hl，m）＜f（x）就始終為真，而即使遇到影響Gauss－Newton法有效性的病態(tài)二次項(xiàng)，也可通過(guò)阻尼因子μ來(lái)控制.LM法給出的是介于Gauss－Newton法和最速下降法之間的搜索方向.對(duì)于尋找最佳的函數(shù)系數(shù)，該法相對(duì)靈活、簡(jiǎn)單、快捷.

根據(jù)得到所需的函數(shù)、函數(shù)個(gè)數(shù)與函數(shù)參數(shù)g的初始值，運(yùn)用優(yōu)化標(biāo)準(zhǔn)RMS，NCC與非線性回歸分析L－M法獲得最佳參數(shù)g，其流程圖如圖2所示.

圖2 算法流程圖Fig.2 Arithmetic flow chart

4 實(shí)驗(yàn)分析與驗(yàn)證

4.1 基本模型

選擇3個(gè)基本的旋轉(zhuǎn)對(duì)稱輻射模型，120°的近似Lambertian型采用方程（2）計(jì)算（i＝1），Batwing型采用方程（2）計(jì)算（i＝2），Side emitter型采用方程（1）計(jì)算（i＝3），相關(guān)的參數(shù)設(shè)置如表1所示.

將輻射模型模擬光強(qiáng)分布公式的數(shù)據(jù)與測(cè)量數(shù)據(jù)進(jìn)行比較，結(jié)果如圖3所示.圖3（a），（b），（c）中的右圖均為空間立體輻射模型，觀看默認(rèn)方位角為37.5°，俯角為30°.

圖3 簡(jiǎn)單輻射模型的模擬數(shù)據(jù)與測(cè)量數(shù)據(jù)比較Fig.3 Comparison of simulated data with measured data of simple radiation patterns

表1 簡(jiǎn)單輻射模型的優(yōu)化前后的函數(shù)參數(shù)Tab.1 Function parameters of simple radiation patterns before and after optimization

4.2 相對(duì)復(fù)雜模型

Nichia的NSPW345CS發(fā)光二極管光強(qiáng)分布模型是非旋轉(zhuǎn)對(duì)稱的，所以用方程（3）計(jì)算（i＝2）；LUXEON Rebel DS56暖白、中性白與冷白的發(fā)光二極管光強(qiáng)分布模型是旋轉(zhuǎn)對(duì)稱的，所以用方程（1）計(jì)算（i＝2），相關(guān)參數(shù)設(shè)置如表2所示.

表2 相對(duì)復(fù)雜輻射模型的優(yōu)化前后的函數(shù)參數(shù)Tab.2 Function parameters of complex radiation patterns before and after optimization

兩個(gè)輻射模型的模擬光強(qiáng)分布公式的數(shù)據(jù)與測(cè)量數(shù)據(jù)比較，如圖4所示.圖4（a），（b）中的右圖均為空間立體輻射模型.

圖4 相對(duì)復(fù)雜輻射模型的模擬數(shù)據(jù)與測(cè)量數(shù)據(jù)比較Fig.4 Comparison of simulated data with measured data of complex radiation model

綜上圖3，4可得，所模擬的光強(qiáng)分布的精確度很高，RMS與NCC遠(yuǎn)超過(guò)設(shè)計(jì)的所需精度，而且實(shí)際操作靈活、簡(jiǎn)單、方便.

5 結(jié)論

基于Source Property Generator軟件，根據(jù)配光曲線圖獲得相應(yīng)的發(fā)光二極管光源數(shù)據(jù).提出運(yùn)用非線性回歸分析（L－M）法，由模擬公式的光源數(shù)據(jù)與根據(jù)配光曲線圖（或測(cè)量）獲得的光源數(shù)據(jù)（又稱測(cè)量數(shù)據(jù)）之間的誤差來(lái)尋求模擬公式參數(shù)變化的步長(zhǎng)，同時(shí)由RMS與NCC兩個(gè)判斷標(biāo)準(zhǔn)來(lái)重建輻射模型的精確性.結(jié)果表明：RMS低于1%，NCC接近100%.

采用這種方法只需要光源的配光曲線圖就可以獲得相應(yīng)的光強(qiáng)分布函數(shù)和輻射模型，如果前面的函數(shù)個(gè)數(shù)與參數(shù)設(shè)置沒(méi)有出現(xiàn)錯(cuò)誤，一般只要一次就能通過(guò)驗(yàn)證，這極大地節(jié)省了設(shè)計(jì)的周期；而且更重要的一點(diǎn)，只要改變某些要求參數(shù)，便可以輕而易舉地設(shè)計(jì)出一系列的同類型透鏡的輻射模型函數(shù).該法具有相對(duì)靈活、簡(jiǎn)單、方便和精確性較高等特點(diǎn).

［1］MORENO I，TZONCHEV R I.Effects on illumination uniformity due to dilution on arrays of LEDs［J］.SPIE，2004，5529：268－275.

［2］MORENO I，AVENDA?O－ALEJO M，TZONCHEV R I.Designing light－emitting diode arrays for uniform nearfield irradiance［J］.Applied Optics，2006，45（10）：2265－2272.

［3］SUN Ching－cherng，LEE Tsung－xian，MA Shih－h(huán)sin，et al.Precise optical modeling for LED lightiing verified by cross correlation in the midfield region［J］.Optical Letters，2006，31（14）：2193－2195.

［4］MORENO I，MU?OZ J，LVANOV R.Uniform illumination of distant targets using a spherical light－emitting diode array［J］.Optical Engineering，2007，46（3）：033001－1－033001－7.

［5］劉旭，李海峰.現(xiàn)代投影顯示技術(shù)［M］.浙江：浙江大學(xué)出版社，2009：84－113.

［6］WANG Kai，CHEN Fei，LIU Zong－yuan，et al.Design of compact freeform lens for application specific light－emitting diode packaging［J］.Optics Express，2010，18（2）：413－425.

［7］丁毅，鄭臻榮，顧培夫.實(shí)現(xiàn)LED照明的自由曲面透鏡設(shè)計(jì)［J］.光子學(xué)報(bào)，2009，38（6）：1486－1489.

［8］CASSARLY W J.LED modeling：Pros and cons of common methods［J］.NASA Tech Briefs，2002（4）：Ⅱa－2a.

［9］MORENO I，SUN Ching－cherng.Three－dimensional measurement of light－emitting diode radiation pattern：A rapid estimation［J］.Measurement Science and Technology，2009，20（7）：1－6.

［10］吳寶寧，李宏光，俞兵，等.LED光學(xué)參數(shù)測(cè)試方法研究［J］.應(yīng)用光學(xué)，2007，28（4）：513－516.

［11］李玉權(quán)，崔敏.光波導(dǎo)理論與技術(shù)［M］.北京：人民郵電出版社，2002：55－61.

［12］MORENO I，SUN Ching－cherng.Modeling the radiation pattern of LEDs［J］.Optics Express，2008，16（3）：1808－1819

［13］DANIEL C，WOOD F S.Fitting equation to data：Computer analysis of multifactor data［M］.2nd ed.New York：Wiley，1999.

［14］MADSEN K，NIELSEN H B，TINGLEFF O.Methods for non－linear least squares problems［R］.2nd ed.Copenhagen：Technical University of Denmark，2004.

Quickly Modeling of Light Emitting Diode Radiation Pattern and Certifying

DONG Jing－jing，GUO Zhen－ning，LIN Jie－ben，HUANG Zhi－wei，HUANG Ming－bo，LI Jian－gong

（College of Information Science and Engineering，Huaqiao University，Xiamen 361021，China）

According to the angular intensity distribution，light emitting diode（LED）intensity datasheets are obtained by using Source Property Generator software.This paper firstly offers the way using nonlinear regression Levenberg－Marquart method which seeks corresponding mutative step length of the parameter of simulating equation，according to the error between the intensity value of simulating equation and Measurement，and changes it.Meanwhile，using two estimation criterion，namely，root mean square（RMS）and normalized cross correlation（NCC），rebuild the accuracy of LED radiation pattern.Results show that RMS is less than 1%，NCC is near 100%.Only providing the map of LED source′s the angular intensity distribution can gain corresponding intensity distribution equation and radiation pattern immediately，if there is not any wrong in function setting.Most important point is that a series of similar types of LENs′Radiation Pattern equations can be obtained by changing some parameters which need to be solved.

optical design；light emitting diode；radiation pattern；levenberg－marquart method

TN 312.802

A

1000－5013（2012）02－0143－06

2011－04－19

郭震寧（1958－），男，教授，主要從事半導(dǎo)體發(fā)光器件及其光學(xué)設(shè)計(jì)與應(yīng)用研究.E－mail：znguo2003＠yahoo.com.cn.

福建省科技計(jì)劃重點(diǎn)項(xiàng)目（2009H0034）；福建省自然科學(xué)基金資助項(xiàng)目（2010J01338）；福建省廈門市科技計(jì)劃重點(diǎn)項(xiàng)目（3502Z20093033）；福建省泉州市科技計(jì)劃重點(diǎn)項(xiàng)目（2008G7）

（責(zé)任編輯：錢筠 英文審校：吳逢鐵）