本文提出 了一種優(yōu)化非序列光學系統(tǒng)的方法。 推薦的方法是使用像素插值(Pixel Interpolation)、探測器數(shù)據(jù)合集(光照時刻數(shù)據(jù))和正交下降優(yōu)化器。 例如，優(yōu)化一個自由曲面反射鏡，使 LED 的亮度從23 Cd 到大于250 Cd只需幾步。 (聯(lián)系我們獲取文章附件)

簡介

OpticStudio 的優(yōu)化功能允許用戶通過將系統(tǒng)參數(shù)設為變量，在評價函數(shù)編輯器中定義性能標準來改進設計。這個過程會對設計產(chǎn)生巨大的影響，所以選擇合適的變量和標準非常重要。序列模式和非序列模式中可用的標準類型有所不同。本文為非序列系統(tǒng)的優(yōu)化提供了一種建議方式。

例如，通過優(yōu)化自由曲面反射鏡，最大限度地將 LED 的亮度從23 Cd 提高到大于250 Cd，只需幾分鐘。

阻尼最小二乘法與正交下降法對比

OpticStudio 中有兩種局部優(yōu)化算法 :阻尼最小二乘法(DLS)和 正交下降法(OD)。DLS 運用數(shù)值微分計算，在一個較小的評價函數(shù)設計的解空間里確定優(yōu)化方向。這種梯度方法是為光學系統(tǒng)設計專門開發(fā)的，被推薦用于所有成像和經(jīng)典光學優(yōu)化問題。然而，在純非序列系統(tǒng)優(yōu)化中，由于采用像素探測器進行探測，DLS 的優(yōu)化效果較差。并且評價函數(shù)本身是不連續(xù)的，這也可能導致梯度搜尋方法失敗。

下面是當評價函數(shù)只有一個變量時，對非序列系統(tǒng)的評價函數(shù)進行查看。

可以看出，很長一段區(qū)間內評價函數(shù)根本沒有變化，發(fā)生的變化是突然且不連續(xù)的。這使得通過梯度搜尋方法進行優(yōu)化變得困難。

正交下降優(yōu)化利用變量的正交化和解空間的離散采樣來降低評價函數(shù)值。OD 算法不計算評價函數(shù)的數(shù)值微分。對于評價函數(shù)存在原本噪聲的系統(tǒng)而言，例如非序列系統(tǒng)，OD 通常比 DLS 算法要好。它在照度最大化、亮度增強和均勻性優(yōu)化等優(yōu)化問題中非常有用。

像素插值和非相干強度數(shù)據(jù)（NSDD）

除了使用的特定算法外，OpticStudio 中還有幾個能夠顯著提升非序列系統(tǒng)優(yōu)化效率的功能。

如上所述，由于探測器的像素特性，非序列的解空間往往是不連續(xù)的。如果一條給定光線中的能量只分配給一個像素，那么當系統(tǒng)更改導致光線在該像素內的任何位置移動時，就沒有數(shù)量上的差異。因此，當一條光線跨越邊界進入一個新的像素時就需要對不連續(xù)的評價函數(shù)求微分，優(yōu)化非常困難。

這可以通過觀察入射到探測器上的一條光線來說明。下面的通用繪圖顯示了探測器上的輻照度質心如何隨著光線位置的變化而變化。

避免由像素探測器引起的量化效應的一種方法是使用像素插值。不是將100%的光線能量分配給單個被擊中的像素，而是根據(jù)光線在像素內的截距位置將一部分能量分配給相鄰像素。因此，當系統(tǒng)變化導致光線在像素上移動時，評價函數(shù)會發(fā)生明顯的變化?？梢栽谔綔y器的 物體屬性( Object Properties )> 類型(Type )中啟用像素插值。

如果我們在啟用了像素插值的探測器上掃描一條光線，會發(fā)現(xiàn)輻照度質心，以及大多數(shù)標準值的變化都會是連續(xù)的，并且可以很容易地使用 DLS。

可以利用優(yōu)化操作數(shù) NSDD 計算評價函數(shù)中所述的輻照度質心。NSDD 代表非序列非相干強度數(shù)據(jù)，是讀取非相干探測器數(shù)據(jù)最有效的操作數(shù)。NSDC 等效于相干計算。NSDD 操作數(shù)的語法如下:

NSDD Surf Det# Pix# Data #Ignored
Surf 定義非序列組的面序號（在純非序列模式中為1）；
Det# 定義所需的探測器的編號（也可以用于清除一個或所有探測器數(shù)據(jù)）；
Pix# 定義要返回的指定像素上的數(shù)值，Data 定義是否返回入射通量、輻照度或強度數(shù)據(jù)；
#Ignored 定義要忽略的邊緣像素數(shù)量。

這些參數(shù)允許優(yōu)化許多標準值：最小光斑大?。ㄗ钚【礁臻g寬度）、最大能量（總光通量)、空間均勻性（所有像素的標準差）、準直度（最小均方根角寬度）等等。有關 NSDD 功能的詳細描述，請參見 OpticStudio 幫助文件部分: 優(yōu)化選項卡( The Optimize Tab )(序列模式界面)> 自動優(yōu)化組( Automatic Optimization Group )> 評價函數(shù)編輯器(Merit Function Editor )(自動優(yōu)化組)> 按類別優(yōu)化操作數(shù)( Optimization Operands by Category )> 非序列光線追跡和探測器操作數(shù)(Non-Sequential Ray Tracing and Detector Operands )。

系統(tǒng)設置

發(fā)光二極管 (LEDs) 是應用廣泛的重要光源。在汽車照明和顯示照明等領域，經(jīng)常需要通過添加輔助光學元件來改變 LED 光源的發(fā)光強度，從而提高 LED 的亮度。

我們將從實際 LED 光源的測量數(shù)據(jù)開始?，F(xiàn)在，我們所需要知道的是，徑向光源(Source Radial)被用來輸入測量到的功率，作為角度的函數(shù)。測量到的單色光源總輸出功率為27 lm，峰值波長為627 nm。光源采用 Sobol 采樣 ，用最少的光線獲得最佳信噪比。

在 系統(tǒng)資源管理器( System Explorer )下，我們將系統(tǒng) 單位( Units )設置如下:

LED 的光通量是以流明為單位來測量的，所以我們選擇該單位進行模擬。因此光照度單位為lm/m 2或 Lux。發(fā)光強度(“亮度”)是用 lm/sr（流明/立體角）或 Cd（坎德拉）來計量的。亮度用 lm/m 2/sr或Cd/m 2 來計量，有時也稱為 nit（尼特）。系統(tǒng)初始設置如下：

LED 光源將光線射向一面平面鏡，然后照亮探測器表面。這個文件可以從本文頂部的示例文件鏈接處下載?？梢钥吹教綔y器上的空間分布和角度分布如下圖:

從圖中可以看出，由于 LED 發(fā)出的光線在反射鏡上有輕微的溢出，所以空間分布和角度分布略有不對稱。這是為了給設計增加一點復雜性而有意為之。

從發(fā)光強度圖中可以看出，41 Cd左右的峰值亮度發(fā)生在極角為27°左右時。近似垂直于探測器表面的光線只有23 Cd 的發(fā)光強度(將在下面討論如何得到這個數(shù)字)。這樣的輪廓對于車前燈照明系統(tǒng)或投影機照明系統(tǒng)來說并不好。一般情況下，我們希望低角度時的光亮度盡可能高，這樣光源就可以投射到遠處。

現(xiàn)在，我們將優(yōu)化反射鏡的形狀，使軸上的亮度最高。為此，我們必須執(zhí)行以下步驟:

1 定義一個可表達我們想要達到的目標的評價函數(shù)。
2 定義反射鏡表面將如何變化。
3 運行優(yōu)化。

評價函數(shù)

根據(jù)設計是否能很好地滿足其指定的工作性能，使用評價函數(shù)來定義光學設計的“質量”。在這種情況下，我們想要在零度角下得到最大的光亮度(發(fā)光強度)，這可以通過操作數(shù) NSDD 和 NSTR 很容易地實現(xiàn)。在本設計中，探測器為物體 3，我們想要得到在0°下的發(fā)光強度。 探測器查看器 ( Detector Viewer )如下所示。

這張圖顯示了光線從 x 和 y 軸的-90°到+90°入射到探測器上光線的角度范圍。沒有光線的角度超過35°，因為 LED 沒有超過這個角度發(fā)射。峰值強度出現(xiàn)在27°左右。我們感興趣的是接近0°的光功率。針對這樣的分布有兩個標準 :RMS 角寬度和質心發(fā)光強度。設置 RMS 角寬度是要獲得準直光線（即具有相同入射角的光線），質心目標是使入射角為零。下面的評價函數(shù)可看到這個角度下的發(fā)光強度。

第一個 NSDD 操作數(shù)讀取探測器物體0，該物體不存在；沒有物體0可以存在。這是此操作數(shù)的一個特殊用法: 在 OpticStudio 中用于清除所有探測器。另外，可以通過定義一個負數(shù)(例如：Det# = -3，只清除探測器3)來清除單一探測器。這在定義了多個探測器的系統(tǒng)中非常有用。

接下來，操作數(shù) NSTR 告訴 OpticStudio 去完成光線追跡。第二和第三個操作數(shù) NSDD 讀出探測器物體3，質心 x 和 y (Pix# = -6， -7)，數(shù)據(jù)項2，即功率/單位立體角。注意，我們的目標是質心發(fā)光強度(角度)，而不是質心照度(空間)。第四個操作數(shù) NSDD 讀取所有像素數(shù)據(jù)的 RMS 角寬度。此外，最后一個操作數(shù) NSDD 讀取了中心像素(5101)的強度，以便進行比較；注意，沒有分配權重，操作數(shù)就不會對評價函數(shù)有貢獻。這個值大約是23 Cd。

最后一個 NSDD 操作數(shù)與一個 OPGT 操作數(shù)相結合，以保持探測器上的最小光通量。我們把這個設為25，因為這是探測器上的初始光通量。如果這個操作數(shù)不存在，那么完全移除反射鏡就有可能使評價函數(shù)值為0！如果沒有能量落在探測器上，質心強度和 RMS 半徑為零，這就是我們的目標。這個“解”強調了定義良好的評價函數(shù)的重要性。在優(yōu)化過程中，OpticStudio 將嘗試將評價函數(shù)變?yōu)榱?，而不管這在系統(tǒng)中的物理意義。

自由曲面反射鏡

自由曲面通常被描述為多個低階多項式，如樣條曲線或貝塞爾曲線。它們通常用來描述渦輪葉片、車身和船殼等形狀。

然而，在光學系統(tǒng)設計時，保留基本圓錐曲線部分的概念，并使該部分向自由曲面形式偏離，是非常有助于設計的。原因很快就會揭曉。我們將使用 擴展多項式面( Extended Polynomial Surface )物體，這個表面可由以下形式的方程描述:

第一項是光學設計中所鐘愛的標準圓錐非球面，用于設計球面、橢圓、拋物面、雙曲面、反射面等等。第二項表示從該表面出發(fā)的一系列越來越高階的多項式形變。多項式是 x 和 y 的冪級數(shù)，第一項是 x，然后是 y，接著是 x*x, x*y, y*y 等。有2個1階項，3 個2階項，4個3階項，以此類推。最高階為20，由此可創(chuàng)建最多230個多項式非球面系數(shù)。位置值 x 和 y 除以歸一化半徑得到的多項式系數(shù)是無量綱的。在本設計中，多項式的最大階數(shù)限制在20項以內，因此最大的自由曲面偏差為 x 0y 5和 x 5y 0。在實際設計中這既不是必要的也不是建議的，只是在設計過程中可供選擇而已。

現(xiàn)在，如果我們使用全局繪圖，可以在掃描反射鏡曲率半徑時顯示中心像素強度。

這張圖說明了優(yōu)化非序列系統(tǒng)的困難，以及定義良好的評價函數(shù)的必要性。如果我們看一下整體的評價函數(shù)值與曲率半徑的關系，就能明白為什么質心和均方根半徑是更好的優(yōu)化目標。

現(xiàn)在我們的評價函數(shù)正確地定義了設計標準，我們將使用局部和全局優(yōu)化算法進行優(yōu)化，對比 DLS 和 OD 的優(yōu)化結果。

優(yōu)化

OpticStudio 包含兩個“全局”優(yōu)化例程，用于搜索更大的解空間區(qū)域。全局優(yōu)化算法采用基因演算法優(yōu)化、多起點計算和局部優(yōu)化相結合的方法，能夠有效地搜索到低評價函數(shù)區(qū)域的多維參數(shù)空間。 錘形(Hammer)優(yōu)化器也使用基因演算法和局部優(yōu)化器，一旦全局搜索找到一個有希望的參數(shù)空間區(qū)域，就會對設計進行全面的優(yōu)化。

評價函數(shù)的初始值為14.8，0°時的亮度為 23Cd。

我們先使用 DLS 使用局部搜索例程(Optimize > Optimize!)，并將其與使用 OD 算法的優(yōu)化結果進行比較。最終的錘形優(yōu)化將在這兩種情況下執(zhí)行。

我們的評價函數(shù)和起始系統(tǒng)已經(jīng)定義好了，剩下的唯一步驟就是分配變量。我們有22個變量：曲率半徑、圓錐系數(shù)和20個多項式系數(shù)。分配這些參數(shù)變量狀態(tài)，并使用帶有自動循環(huán)數(shù)的DLS算法啟動局部優(yōu)化器。

優(yōu)化了一段時間后(約4.4分鐘)，OpticStudio 找到了解。評價函數(shù)降到了大約6.69，中心像素亮度為238 Cd。此次優(yōu)化說明使用像素插值和定義良好的評價函數(shù)，即使是DLS算法也可以在非序列解空間中有效地工作。

用一個新名稱保存結果文件以便進行比較，然后再次打開起始文件。這一次，分配所有22個變量并使用 OD 算法進行優(yōu)化。根據(jù)前面對兩種局部優(yōu)化算法的比較，我們可以期望這種優(yōu)化能夠更快地進行，并得到更好的解。在這種情況下，此優(yōu)化比 DLS 運行時間稍長(約6.5分鐘)，但得到了更小的評價函數(shù)值(6.68)。

為了測試我們得到的是一個可信的最優(yōu)解，而不是陷入一個較高的局部最小值，我們將運行錘形優(yōu)化器來優(yōu)化這兩個系統(tǒng)。下表整理的結果表明，DLS 優(yōu)化得到了一定的改善，但中心像素亮度仍低于 OD 算法優(yōu)化得到的像素亮度；進一步說明了這兩個局部優(yōu)化例程的強大功能。在這些測試中使用的機器擁有英特爾四核CPU (2.90 GHz)和 16GB 運行內存:

初始輻照度和輻射強度分布，以及使用 DLS 和 OD 算法優(yōu)化的系統(tǒng)如下圖所示。

起始系統(tǒng):

DLS 優(yōu)化:

OD 優(yōu)化:

結果表明，兩種優(yōu)化方法得到的解相似。但值得注意的是，雖然 OD 算法的優(yōu)化時間比 DLS 稍長，但所獲得的性能優(yōu)于 DLS 優(yōu)化和30分鐘錘形優(yōu)化。事實上，優(yōu)化一個包含22個變量和幾個不同照度目標的系統(tǒng)只需要幾分鐘！我們已經(jīng)在短時間內從一個平面鏡得到一個完全優(yōu)化的解決方案，之后將需要我們進行公差分析和生成系統(tǒng)的圖紙。對具有相同數(shù)量變量的序列系統(tǒng)均方根半徑進行優(yōu)化需要同樣的時間。這個例子清楚地說明了 OD 算法在優(yōu)化非序列系統(tǒng)時的效率。在本文頂部的示例文件鏈接中，您可以找到 OD 優(yōu)化文件。

請記住，我們從未直接針對中心像素亮度進行優(yōu)化，而是優(yōu)化均方根半徑和質心位置。我們總是可以對這個標準施加有限的權重，但是均方根半徑很可能因此而增加。您會發(fā)現(xiàn)，如果您的目標是中心像素強度而不是照明二階矩，那么系統(tǒng)性能會明顯變差。