包實秋+許東+黃選平+譚曉頌

摘 要：針對高超聲速飛行器氣動光學(xué)效應(yīng)及其校正研究的技術(shù)需求，設(shè)計開發(fā)了一套氣動光學(xué)效應(yīng)仿真系統(tǒng)，用于對氣動光學(xué)效應(yīng)進(jìn)行全過程仿真。系統(tǒng)基于“模塊化”、“可擴(kuò)展”和“可視化”的設(shè)計思路，采用“仿真平臺+數(shù)據(jù)庫+功能插件”的系統(tǒng)實現(xiàn)方案，在VisualStudio2008開發(fā)環(huán)境下，完成了氣動光學(xué)效應(yīng)仿真系統(tǒng)平臺、數(shù)據(jù)庫及相關(guān)功能插件的開發(fā)工作。系統(tǒng)支持用戶根據(jù)自身需求自主編制仿真任務(wù)，并可在接口規(guī)范的約束下，自主定制功能模塊對系統(tǒng)進(jìn)行擴(kuò)展。本系統(tǒng)可以為解決高超聲速飛行器氣動光學(xué)效應(yīng)及其校正方法的理論研究與工程應(yīng)用問題提供理論依據(jù)和技術(shù)手段。

關(guān)鍵詞：氣動光學(xué);仿真軟件;任務(wù)編制;模塊化;可擴(kuò)展性

中圖分類號：TP311.5 文獻(xiàn)標(biāo)識碼：A 文章編號：1673-5048（2014）05-0041-05

0 引 言

高超聲速飛行器是當(dāng)今世界各國的研究熱點之一。當(dāng)飛行器在大氣層中高速飛行時，會產(chǎn)生劇烈的氣動光學(xué)效應(yīng)，嚴(yán)重影響光學(xué)成像過程，是高超聲速飛行器光學(xué)成像探測必須解決的瓶頸問題。

到目前為止，國內(nèi)外很多研究機(jī)構(gòu)已對氣動光學(xué)效應(yīng)機(jī)理及其校正方法開展了大量的研究工作，也報道了一些研究成果。這些研究包括氣動光學(xué)效應(yīng)機(jī)理[1]、測試試驗[2-3]、數(shù)值模擬[4-5]和退化圖像校正[6-7]等方面。

在此基礎(chǔ)上，以美國為代表的西方國家對氣動光學(xué)效應(yīng)的研究，己開始向工程應(yīng)用方向發(fā)展。如美國TeledyneBrownEngineering公司[8]、HIPEC

評估中心以及Lockheed公司等研究機(jī)構(gòu)分別研制了AOQCode、DNAOS[9]等氣動光學(xué)效應(yīng)仿真軟件，用于對高超聲速飛行器的氣動光學(xué)效應(yīng)進(jìn)行仿真分析和探測性能預(yù)測，并已經(jīng)在相關(guān)型號的武器研制中發(fā)揮了重要作用。

國內(nèi)對氣動光學(xué)效應(yīng)也開展了相關(guān)的研究工作，在氣動光傳輸[10]、氣動熱輻射[11]、氣動加熱[12]、退化質(zhì)量評估[13]等方面取得了一定的研究進(jìn)展。但是到目前為止，還沒有出現(xiàn)實用化的、能夠?qū)鈩庸鈱W(xué)效應(yīng)全過程進(jìn)行完整仿真和分析的系統(tǒng)平臺。在實際的工程應(yīng)用中，用戶難以真正系統(tǒng)、全面地考慮氣動光學(xué)效應(yīng)對成像探測帶來的影響問題。

本文在當(dāng)前理論研究的基礎(chǔ)上，設(shè)計并開發(fā)了一套氣動光學(xué)仿真系統(tǒng)。該系統(tǒng)能夠?qū)Ω叱曀倮@流、流場與頭罩氣動熱耦合、流場和光學(xué)頭罩的氣動光傳輸、流場和光學(xué)頭罩的氣動熱輻射、光學(xué)探測系統(tǒng)成像等進(jìn)行全過程數(shù)值模擬和仿真分析，并能夠?qū)鈩庸鈱W(xué)效應(yīng)引起的像質(zhì)退化進(jìn)行評估和分析，能夠為氣動光學(xué)效應(yīng)仿真理論研究和工程應(yīng)用提供完整的解決方案。

1 系統(tǒng)總體設(shè)計方案

對氣動光學(xué)效應(yīng)的研究結(jié)果表明，氣動光學(xué)效應(yīng)主要包括：高超聲速繞流引起的流場光傳輸效應(yīng)、流場高溫氣體分子的熱輻射效應(yīng)、光學(xué)頭罩溫度的升高和應(yīng)力應(yīng)變引起的頭罩光傳輸效應(yīng)、光學(xué)頭罩熱輻射效應(yīng)等[14]。因此，氣動光學(xué)效應(yīng)仿真的研究內(nèi)容可劃分為六部分，如圖1所示。

這六部分涵蓋了氣動光學(xué)效應(yīng)影響的全過程。通過仿真，可以得到氣動光學(xué)效應(yīng)給紅外圖像所帶來的“抖動”、“模糊”、以及強(qiáng)烈的“熱背景干擾”等影響，并以相關(guān)指標(biāo)對上述影響進(jìn)行定量的分析。

針對上述需求，圖2給出了系統(tǒng)的總體設(shè)計方案。

氣動光學(xué)仿真系統(tǒng)被劃分為“仿真參數(shù)”，“仿真計算”，“數(shù)據(jù)庫”三部分。

“數(shù)據(jù)庫”部分負(fù)責(zé)仿真系統(tǒng)的數(shù)據(jù)存儲與管理，其存儲了仿真過程所需的全部參數(shù)，計算模塊的配置信息，系統(tǒng)平臺的配置信息以及仿真過程中所生成的中間結(jié)果。

“仿真參數(shù)”部分為用戶提供了配置參數(shù)的接口，用戶可以通過該部分的模塊完成特定參數(shù)的選擇、修改、添加以及刪除。用戶所進(jìn)行的所有設(shè)定都會同步至數(shù)據(jù)庫進(jìn)行保存。仿真參數(shù)主要包括如下內(nèi)容：

（1）目標(biāo)圖像：用于仿真的原始圖像，最后會將氣動光學(xué)效應(yīng)作用在此圖像上，進(jìn)而得到退化結(jié)果。

（2）頭罩材料參數(shù)：飛行器光學(xué)頭罩材料相關(guān)參數(shù)，如比熱容、折射率等。

（3）飛行彈道參數(shù)：描述飛行器飛行的軌跡。

（4）大氣物理模型：提供數(shù)值模擬時所采用的大氣模型，根據(jù)相關(guān)模型計算大氣的溫度、壓強(qiáng)等參數(shù)。

（5）頭罩結(jié)構(gòu)參數(shù)：頭罩結(jié)構(gòu)的約束參數(shù)，通過這些參數(shù)可以完成頭罩幾何結(jié)構(gòu)的建模。

“仿真計算”部分負(fù)責(zé)完成氣動光學(xué)仿真系統(tǒng)的仿真計算任務(wù)，采用基于模塊化的思想進(jìn)行設(shè)計，主要包括以下模塊：

（1）網(wǎng)格建模：根據(jù)具體的頭罩結(jié)構(gòu)參數(shù)以及網(wǎng)格建模要求，實現(xiàn)對仿真對象的自動網(wǎng)格建模，并將建模結(jié)果數(shù)據(jù)保存在相應(yīng)文件中。

（2）流固耦合數(shù)值模擬：以“頭罩材料參數(shù)”、“飛行彈道參數(shù)”、“大氣物理模型”和網(wǎng)格建模結(jié)果為輸入?yún)?shù)，完成流場CFD數(shù)值模擬與流場-頭罩氣動熱耦合分析，模擬飛行時頭罩本身以及周圍流場的溫度、壓力、密度等狀態(tài)量的分布并保

（3）流場光線追跡：根據(jù)流場密度的數(shù)值模擬仿真結(jié)果，計算流場密度與折射率分布，并對光線在流場中的傳輸過程進(jìn)行仿真計算，獲得流場出射光線的波前畸變，完成對流場光傳輸效應(yīng)分析。

（4）流場熱輻射計算：根據(jù)流場的數(shù)值模擬結(jié)果、大氣物理模型以及氣體分子光譜輻射的相關(guān)理論，對高溫繞流場中的氣體輻射進(jìn)行仿真計算，獲得流場的熱輻射分布及輻射出射，完成對流場熱輻射效應(yīng)分析。

（5）頭罩結(jié)構(gòu)及折射率變化分析：根據(jù)流固耦合數(shù)值模擬的結(jié)果，導(dǎo)出頭罩的溫度分布數(shù)據(jù)，求出頭罩內(nèi)部的應(yīng)力應(yīng)變分布數(shù)據(jù)。

（6）頭罩光線追跡：首先根據(jù)頭罩溫度和應(yīng)力應(yīng)變的數(shù)值模擬結(jié)果，計算頭罩折射率分布數(shù)據(jù)，然后計算由于頭罩折射率分布不均勻帶來的光波前畸變，完成對頭罩光傳輸效應(yīng)的分析。

（7）頭罩熱輻射計算：根據(jù)頭罩溫度和應(yīng)力應(yīng)變的數(shù)值模擬結(jié)果，對被嚴(yán)重加熱頭罩的熱輻射情況進(jìn)行計算，通過求解輻射傳輸方程獲得頭罩內(nèi)表面的熱輻射出射分布，完成對頭罩熱輻射效應(yīng)的分析。

（8）光學(xué)成像仿真：綜合前述流場光傳輸、流場熱輻射、頭罩光傳輸和頭罩熱輻射的仿真結(jié)果，根據(jù)成像理論，對目標(biāo)圖像進(jìn)行成像仿真，得出退化的圖像結(jié)果。

（9）圖像退化評估：通過對退化圖像與原始目標(biāo)圖像的對比分析，對圖像信噪比、對比度、制導(dǎo)精度、失真程度等指標(biāo)進(jìn)行評估，為后續(xù)的目標(biāo)識別、跟蹤與制導(dǎo)提供依據(jù)。

2 系統(tǒng)的軟件實現(xiàn)

2.1 系統(tǒng)實現(xiàn)技術(shù)方案

系統(tǒng)實現(xiàn)方案如圖3所示。系統(tǒng)全部模塊插件以DLL形式封裝，并提供標(biāo)準(zhǔn)規(guī)范接口嵌入氣動光學(xué)仿真系統(tǒng)平臺中。系統(tǒng)內(nèi)嵌了ANSYS平臺，ANSYS提供給用戶豐富的二次開發(fā)接口，系統(tǒng)能夠根據(jù)仿真相關(guān)參數(shù)自動生成符合其規(guī)范的命令流，然后調(diào)用ANSYS自動讀取命令流來實現(xiàn)所需功能。

2.2 系統(tǒng)平臺類建模

本文采用UML語言對軟件平臺進(jìn)行建模，模型符合UML標(biāo)準(zhǔn)規(guī)范[15]，系統(tǒng)平臺類模型設(shè)計如圖4所示。

主要類的功能如下：

（1）CApp/CView/CMainFrame/CDoc：MFC基本框架類，分別管理程序運行實例、程序視窗、程序數(shù)據(jù)。

（2）DataBase類：程序與數(shù)據(jù)庫的交互模塊，其中定義了相關(guān)狀態(tài)變量以及數(shù)據(jù)庫指針，方便數(shù)據(jù)庫隨時訪問;對程序所需的數(shù)據(jù)庫操作進(jìn)行了封裝，規(guī)范了數(shù)據(jù)庫接口;同時也方便該接口在模塊中進(jìn)行復(fù)用。

（3）SetUp類：該類完成軟件關(guān)于配置的相關(guān)操作以及維護(hù)與其相關(guān)的變量，通過該類的實例可以完成對軟件的配置操作。

（4）Save&Read類：軟件可以對當(dāng)前執(zhí)行的仿真工程進(jìn)行保存和對已保存的工程進(jìn)行讀取，該類封裝了此功能相應(yīng)操作及變量。

（5）QdList類：軟件左側(cè)模塊名稱列表窗口，該列表列出當(dāng)前可用功能插件。

（6）DataWnd/InfoWnd/OutputWnd類：軟件數(shù)據(jù)顯示窗口、信息顯示窗口、狀態(tài)顯示窗口，此三個窗口對軟件的狀態(tài)進(jìn)行實時顯示。

（7）RunSim類：仿真操作窗口，控制仿真任務(wù)的起停，顯示仿真任務(wù)狀態(tài)。

（8）DLLOperation類：接口類，封裝了軟件平臺對DLL的所有操作。

（9）SelfDefine類：抽象類，代表軟件中自定義的相關(guān)數(shù)據(jù)結(jié)構(gòu)。

（10）DataManage類：維護(hù)軟件正常運行所需的相關(guān)變量，封裝了訪問、修改這些變量的操作。

（11）Module類：模塊管理類，封裝了添加模塊、刪除模塊相關(guān)操作，以及模塊數(shù)量等相關(guān)狀態(tài)變量，用于實現(xiàn)模塊的可擴(kuò)展，支持用戶自主構(gòu)建模塊。

2.3 功能插件類建模

功能插件類模型設(shè)計如圖5所示。

StandModel接口類定義了與插件交互所需要的接口函數(shù)，平臺通過這些函數(shù)接口來使用封裝在插件內(nèi)的仿真運算功能、UI接口、圖標(biāo)等資源。為了使插件功能完整而又獨立于系統(tǒng)平臺，插件還具備以下重要組成部分：

（1）SetUpDlg類：該類實現(xiàn)了插件的參數(shù)配置窗口，通過此窗口對插件所需要的參數(shù)進(jìn)行配置;配置完成后，該類調(diào)用DataBase類將配置信息存入數(shù)據(jù)庫。

（2）RunDlg類：該類封裝了插件的仿真運算功能，并且實現(xiàn)了一個運行窗口來顯示仿真運行的狀態(tài)。仿真運行之前，該類會通過DataBase類讀取數(shù)據(jù)庫中的配置參數(shù)以及相關(guān)輸入?yún)?shù)，完成初始化工作后，進(jìn)行仿真運算。

（3）DataPreProcess類為數(shù)據(jù)預(yù)處理類，在利用有限元數(shù)據(jù)進(jìn)行數(shù)值仿真的功能插件中，該類完成數(shù)據(jù)預(yù)處理功能，提供了對有限元數(shù)據(jù)的快速訪問機(jī)制，和對不同格式有限元數(shù)據(jù)的訪問功能，保證了運算準(zhǔn)確性以及運算效率。Element為實現(xiàn)該類所需要的數(shù)據(jù)結(jié)構(gòu)。

圖6給出了系統(tǒng)所采用的“平臺+插件”架構(gòu)的組件圖。

StandModel為虛類，所有模塊組件必須繼承此類并實現(xiàn)該類，才能加入到平臺中被平臺調(diào)用，而每個組件的接口實現(xiàn)不同，因此其功能相互獨立。

系統(tǒng)架構(gòu)中的接口規(guī)范對外開放，用戶只需要繼承StandModel接口，并在新加入的組件中實現(xiàn)接口中的函數(shù)，然后將模塊添加至平臺即可實現(xiàn)對系統(tǒng)功能模塊的自主定制。 3 仿真示例

基于以上設(shè)計方案開發(fā)的氣動光學(xué)效應(yīng)仿真系統(tǒng)，采用通用的用戶操作界面，符合用戶使用習(xí)慣，易于操作。圖7給出了一個已經(jīng)完成氣動光學(xué)效應(yīng)仿真任務(wù)編排的界面圖。

圖中左側(cè)列表對已經(jīng)添加的模塊進(jìn)行實時顯示，工具欄上具有模塊配置按鈕，可以隨時添加，卸載模塊;創(chuàng)建仿真任務(wù)時，只需從左側(cè)列表中將相關(guān)模塊拖入到中間窗口即可創(chuàng)建相應(yīng)模塊，模塊之間采用鼠標(biāo)連線的操作方式進(jìn)行編排。系統(tǒng)提供數(shù)據(jù)源模塊和數(shù)據(jù)中間存儲模塊，使得仿真任務(wù)可以從任意階段開始，提高了系統(tǒng)使用的靈活性。此外系統(tǒng)還提供任務(wù)狀態(tài)、數(shù)據(jù)實時顯示等輔助功能。

圖8和圖9分別給出了兩組不同條件下的圖像仿真結(jié)果。相應(yīng)的退化質(zhì)量評估結(jié)果如表1、表2所示。

圖8為飛行速度Ma=3，飛行高度20km，攻角0度下的仿真結(jié)果，（a）為未經(jīng)氣動光學(xué)效應(yīng)作用的原始圖像，（b）右圖為氣動光學(xué)效應(yīng)作用下的結(jié)果圖像。圖9為飛行速度Ma=5，飛行高度10 km，攻角53°下的仿真結(jié)果。圖8～9均為球形頭罩成像結(jié)果，在成像過程中，首先利用流場和頭罩光線追跡的計算結(jié)果，得出與光傳輸效應(yīng)等效的PSF，并與圖像做卷積，得出在光傳輸效應(yīng)影響下的圖像，再將流場與頭罩的熱輻射計算結(jié)果進(jìn)行量化并轉(zhuǎn)換成灰度分布，疊加至圖像上。從圖中可以看到，飛行條件越苛刻，氣動光學(xué)效應(yīng)越嚴(yán)重。

參考文獻(xiàn)：

[1]殷興良.現(xiàn)代光學(xué)新分支學(xué)科———氣動光學(xué)[J].中國工程科學(xué)，2006，7（12）：1-6.

[2]SuttonGW.AeroOpticalFoundationsandApplications[J].AIAAJournal，1985，23（10）：1525-1537.

[3]JumperEJ，F(xiàn)itzgeraldEJ.RecentAdvancesinAeroOptics[J].ProgressinAerospaceSciences，2001，37（3）：299-339.

[4]GordeyevS，JumperE.FluidDynamicsandAeroOptics ofTurrets[J].ProgressinAerospaceSciences，2010，46（8）：388-400.

[5]TromeurE，GarnierE，SagautP，etal.LargeEddySimulationsofAeroOpticalEffectsinaTurbulentBoundary Layer[J].JournalofTurbulence，2003，4（5）：1-22.

[6]洪漢玉，張?zhí)煨?基于多分辨率盲目去卷積的氣動光學(xué)效應(yīng)退化圖像復(fù)原算法[J].計算機(jī)學(xué)報，2004，27（7）：952-963.

[7]LuX，ZhangT，HongH.ImageCorrectionMethodwith PixelDeviationCausedbyAeroOpticsEffects[J].InfraredandLaserEngineering，2007，36（5）：758.

[8]KathmanAD，BrooksLC，KalinDA，etal.ATimeIntegratedImageModelforAeroOpticAnalysis[C]∥AIAA andSDIO，AunualInterceptorTechnocogyConference，Huntsville，AL，1992.

[9]SuttonGW，PondJE，SnowR，etal.HypersonicInterceptorAeroOpticsPerformancePredictions[J].Journalof SpacecraftandRockets，1994，31（4）：592-599.

[10]陳澄，費錦東.側(cè)窗頭罩高速層流流場光學(xué)傳輸效應(yīng)數(shù)值模擬[J].紅外與激光工程，2006，34（5）：548-552.

[11]安永泉，趙剡.紅外成像制導(dǎo)中的氣動熱輻射效應(yīng)機(jī)理[J].紅外與激光工程，2011，40（7）：1199-1204.

[12]夏剛，劉新建，程文科，等.鈍體高超聲速氣動加熱與結(jié)構(gòu)熱傳遞耦合的數(shù)值計算[J].國防科技大學(xué)學(xué)報，2003，25（1）：35-39.

[13]黎明，楊杰.基于支持向量機(jī)的湍流退化圖像加速復(fù)原算法[J].紅外與毫米波學(xué)報，2009，28（6）：472-475.

[14]楊文霞，蔡超，丁明躍，等.超音速/高超音速飛行器湍流流場氣動光學(xué)效應(yīng)分析[J].光電工程，2009，36（1）：88-92.

[15]RumbaughJ，JacobsonI，BoochG.UnifiedModeling LanguageReferenceManual[M].2nded.ThePearson HigherEducation，2004.

[2]SuttonGW.AeroOpticalFoundationsandApplications[J].AIAAJournal，1985，23（10）：1525-1537.

[3]JumperEJ，F(xiàn)itzgeraldEJ.RecentAdvancesinAeroOptics[J].ProgressinAerospaceSciences，2001，37（3）：299-339.

[4]GordeyevS，JumperE.FluidDynamicsandAeroOptics ofTurrets[J].ProgressinAerospaceSciences，2010，46（8）：388-400.

[5]TromeurE，GarnierE，SagautP，etal.LargeEddySimulationsofAeroOpticalEffectsinaTurbulentBoundary Layer[J].JournalofTurbulence，2003，4（5）：1-22.

[6]洪漢玉，張?zhí)煨?基于多分辨率盲目去卷積的氣動光學(xué)效應(yīng)退化圖像復(fù)原算法[J].計算機(jī)學(xué)報，2004，27（7）：952-963.

[7]LuX，ZhangT，HongH.ImageCorrectionMethodwith PixelDeviationCausedbyAeroOpticsEffects[J].InfraredandLaserEngineering，2007，36（5）：758.

[8]KathmanAD，BrooksLC，KalinDA，etal.ATimeIntegratedImageModelforAeroOpticAnalysis[C]∥AIAA andSDIO，AunualInterceptorTechnocogyConference，Huntsville，AL，1992.

[9]SuttonGW，PondJE，SnowR，etal.HypersonicInterceptorAeroOpticsPerformancePredictions[J].Journalof SpacecraftandRockets，1994，31（4）：592-599.

[10]陳澄，費錦東.側(cè)窗頭罩高速層流流場光學(xué)傳輸效應(yīng)數(shù)值模擬[J].紅外與激光工程，2006，34（5）：548-552.

[11]安永泉，趙剡.紅外成像制導(dǎo)中的氣動熱輻射效應(yīng)機(jī)理[J].紅外與激光工程，2011，40（7）：1199-1204.

[12]夏剛，劉新建，程文科，等.鈍體高超聲速氣動加熱與結(jié)構(gòu)熱傳遞耦合的數(shù)值計算[J].國防科技大學(xué)學(xué)報，2003，25（1）：35-39.

[13]黎明，楊杰.基于支持向量機(jī)的湍流退化圖像加速復(fù)原算法[J].紅外與毫米波學(xué)報，2009，28（6）：472-475.

[14]楊文霞，蔡超，丁明躍，等.超音速/高超音速飛行器湍流流場氣動光學(xué)效應(yīng)分析[J].光電工程，2009，36（1）：88-92.

[15]RumbaughJ，JacobsonI，BoochG.UnifiedModeling LanguageReferenceManual[M].2nded.ThePearson HigherEducation，2004.

[2]SuttonGW.AeroOpticalFoundationsandApplications[J].AIAAJournal，1985，23（10）：1525-1537.

[3]JumperEJ，F(xiàn)itzgeraldEJ.RecentAdvancesinAeroOptics[J].ProgressinAerospaceSciences，2001，37（3）：299-339.

[4]GordeyevS，JumperE.FluidDynamicsandAeroOptics ofTurrets[J].ProgressinAerospaceSciences，2010，46（8）：388-400.

[5]TromeurE，GarnierE，SagautP，etal.LargeEddySimulationsofAeroOpticalEffectsinaTurbulentBoundary Layer[J].JournalofTurbulence，2003，4（5）：1-22.

[6]洪漢玉，張?zhí)煨?基于多分辨率盲目去卷積的氣動光學(xué)效應(yīng)退化圖像復(fù)原算法[J].計算機(jī)學(xué)報，2004，27（7）：952-963.

[7]LuX，ZhangT，HongH.ImageCorrectionMethodwith PixelDeviationCausedbyAeroOpticsEffects[J].InfraredandLaserEngineering，2007，36（5）：758.

[8]KathmanAD，BrooksLC，KalinDA，etal.ATimeIntegratedImageModelforAeroOpticAnalysis[C]∥AIAA andSDIO，AunualInterceptorTechnocogyConference，Huntsville，AL，1992.

[9]SuttonGW，PondJE，SnowR，etal.HypersonicInterceptorAeroOpticsPerformancePredictions[J].Journalof SpacecraftandRockets，1994，31（4）：592-599.

[10]陳澄，費錦東.側(cè)窗頭罩高速層流流場光學(xué)傳輸效應(yīng)數(shù)值模擬[J].紅外與激光工程，2006，34（5）：548-552.

[11]安永泉，趙剡.紅外成像制導(dǎo)中的氣動熱輻射效應(yīng)機(jī)理[J].紅外與激光工程，2011，40（7）：1199-1204.

[12]夏剛，劉新建，程文科，等.鈍體高超聲速氣動加熱與結(jié)構(gòu)熱傳遞耦合的數(shù)值計算[J].國防科技大學(xué)學(xué)報，2003，25（1）：35-39.

[13]黎明，楊杰.基于支持向量機(jī)的湍流退化圖像加速復(fù)原算法[J].紅外與毫米波學(xué)報，2009，28（6）：472-475.

[14]楊文霞，蔡超，丁明躍，等.超音速/高超音速飛行器湍流流場氣動光學(xué)效應(yīng)分析[J].光電工程，2009，36（1）：88-92.

[15]RumbaughJ，JacobsonI，BoochG.UnifiedModeling LanguageReferenceManual[M].2nded.ThePearson HigherEducation，2004.