不確定性傳播的混沌多項式方法研究進(jìn)展

熊芬芬，陳江濤，任成坤，張立，李澤賢

1 北京理工大學(xué) 宇航學(xué)院，北京 100081

2 中國空氣動力研究與發(fā)展中心，四川 綿陽 621000

0 引 言

復(fù)雜工程產(chǎn)品（例如飛行器、船舶、汽車等）高水平、高效率的開發(fā)設(shè)計對國民經(jīng)濟(jì)及國防事業(yè)的發(fā)展有著舉足輕重的作用。為了縮短設(shè)計周期、降低開發(fā)成本、滿足產(chǎn)品不斷提升的性能需求和更新?lián)Q代需要，20 世紀(jì)60 年代中期以來，有限元分析（FEA）和計算流體力學(xué)（CFD）等計算機(jī)仿真技術(shù)和優(yōu)化方法被廣泛應(yīng)用于復(fù)雜系統(tǒng)的設(shè)計。模型確認(rèn)、不確定性分析和優(yōu)化設(shè)計是復(fù)雜工程產(chǎn)品基于仿真的開發(fā)設(shè)計中涉及的兩大關(guān)鍵內(nèi)容。前者主要在于保證仿真模型的高保真度，從而為可靠的產(chǎn)品性能分析和設(shè)計提供保障；后者旨在通過設(shè)計提高產(chǎn)品性能在各種不確定性下的可靠性和穩(wěn)健性。對于二者而言，需要解決的首要問題皆為不確定性傳播（uncertainty propagation, UP），這也一直是工程優(yōu)化領(lǐng)域重要的理論課題之一。然而，隨著工程系統(tǒng)設(shè)計的復(fù)雜化、多學(xué)科化，仿真分析模型計算規(guī)模和計算量增長，響應(yīng)函數(shù)高維、強(qiáng)非線性，給不確定性傳播帶來“維數(shù)災(zāi)難”、精度低、可靠性差等諸多難題。

作為一種高效的概率不確定性傳播和量化理論，混沌多項式（polynomial chaos, PC）方法[1-2]由于其堅實的數(shù)學(xué)基礎(chǔ)和良好性能，近些年在學(xué)術(shù)界和工業(yè)界受到廣泛關(guān)注?；煦缍囗検椒椒▽嵸|(zhì)上相當(dāng)于將隨機(jī)變量表示為一組正交多項式的加權(quán)和，構(gòu)建一個隨機(jī)代理模型，不確定性傳播就直接在這個代理模型上進(jìn)行。它能夠?qū)哂腥我夥植碱愋偷碾S機(jī)變量實現(xiàn)較為精確的近似，且理論上當(dāng)條件滿足，即可獲得指數(shù)收斂速度。而且，一旦混沌多項式模型構(gòu)建完成，輸出響應(yīng)的統(tǒng)計矩、失效概率以及概率密度函數(shù)都能非常方便地得到。相比于傳統(tǒng)的蒙特卡羅仿真（Monte Carlo simulation, MCS）方法，混沌多項式方法在保證精度的前提下，可大幅降低計算量；相比于傳統(tǒng)的一階[3]、二階[4-5]可靠性分析方法，對非線性函數(shù)具有更高的不確定性傳播精度，且無需函數(shù)的導(dǎo)數(shù)信息，應(yīng)用起來更加靈活。目前，混沌多項式方法在機(jī)械[6]、土木[7-8]、材料[9]、電力電子[10-11]、汽車[12]、航空航天[13-19]、船舶[20-23]、控制[14,24-25]等領(lǐng)域都得到了廣泛研究和應(yīng)用，在航空航天和船舶領(lǐng)域應(yīng)用尤為廣泛。

由于構(gòu)建混沌多項式模型的計算量隨著不確定性輸入維數(shù)的增加呈指數(shù)增長，高維下面臨嚴(yán)重的“維數(shù)災(zāi)難”問題，這也是目前阻礙混沌多項式方法在工程問題中廣泛應(yīng)用的最大障礙之一。因此，如何解決或緩解混沌多項式的“維數(shù)災(zāi)難”難題，降低計算量，一直是學(xué)術(shù)界的研究熱點。另一方面，混沌多項式方法基于概率理論，僅能處理隨機(jī)不確定性，然而實際問題中存在大量由于知識或數(shù)據(jù)不足而導(dǎo)致的認(rèn)知不確定性，例如在CFD 建模中由于認(rèn)知不夠，湍流模型的建立存在模型不確定性和相關(guān)參數(shù)選取的不確定性。因此，如何充分利用混沌多項式的優(yōu)勢，將其擴(kuò)展為能處理認(rèn)知不確定性，是關(guān)注的另一熱點。本文將首先給出不確定性傳播的基本概念及其在模型確認(rèn)和不確定性優(yōu)化中的作用，然后系統(tǒng)地介紹混沌多項式方法的構(gòu)建原理及其各類變種，接著針對當(dāng)前混沌多項式應(yīng)用中面臨的“維數(shù)災(zāi)難”、計算量大、無法直接處理認(rèn)知不確定性等問題介紹相應(yīng)的解決策略，最后對混沌多項式的研究進(jìn)行展望。

1 基本概念

1.1 不確定性傳播

不確定性傳播也稱不確定性分析，是研究各種不確定性對產(chǎn)品系統(tǒng)性能（泛指系統(tǒng)輸出）的影響規(guī)律的方法。簡而言之，不確定性傳播就是在給定輸入的不確定性信息下，如何估算輸出響應(yīng)的不確定性信息。不確定性主要分為兩大類：隨機(jī)不確定性和認(rèn)知不確定性，前者無法控制或減少，主要基于概率理論進(jìn)行研究，而后者主要由認(rèn)知或數(shù)據(jù)不足導(dǎo)致，可以減少。目前有多種研究認(rèn)知不確定性的方法，例如模糊[26]、區(qū)間[27]、證據(jù)[28-29]理論等。由于概率方法發(fā)展成熟，具有堅實的數(shù)學(xué)基礎(chǔ)，很多情況下結(jié)合貝葉斯定理也可以處理認(rèn)知不確定性，因此目前概率方法應(yīng)用最為廣泛。本文介紹的混沌多項式方法建立在隨機(jī)概率空間，主要用于處理隨機(jī)不確定性。圖1展示了概率不確定性傳播的基本概念，從數(shù)學(xué)上具體描述為：在隨機(jī)輸入x=[x1, ···,xd]存在不確定性的情況下（此時x1, ···,xd的不確定性可以用其概率密度函數(shù)、累積分布函數(shù)、或均值和方差描述），計算輸出響應(yīng)y的不確定性信息，包括均值、方差、失效概率、概率密度函數(shù)（PDF）等。

圖1 概率不確定性傳播示意圖Fig. 1 Probabilistic uncertainty propagation

1.2 模型確認(rèn)

對于基于仿真的復(fù)雜產(chǎn)品性能分析和優(yōu)化設(shè)計，要想保證精度，構(gòu)建高可信度的仿真模型是關(guān)鍵。然而，由于物理過程的復(fù)雜性及人們的認(rèn)知偏差，物理建模與數(shù)值模擬始終存在不確定性。以CFD 計算為例，其存在幾何模型、模型假設(shè)、模型參數(shù)（例如湍流模型封閉系數(shù)）、迭代方法等各類不確定性，從而嚴(yán)重影響數(shù)值模擬結(jié)果的可信度, 使得決策需承擔(dān)很大風(fēng)險。以飛行器結(jié)構(gòu)有限元分析數(shù)值模擬為例，傳統(tǒng)的確定性有限元分析方法并未考慮有限元建模與試驗測量中普遍存在的不確定性，在實際應(yīng)用中僅以變量的均值對問題進(jìn)行分析與描述?；蛘咭試@均值附近的“安全系數(shù)”這類概念粗略估計問題的隨機(jī)性，基于工程師個人的工程經(jīng)驗得到“安全”或“不安全”這類結(jié)論。若以此作為飛行器結(jié)構(gòu)設(shè)計的關(guān)鍵性決策依據(jù)，顯然是不盡合理的。若不確定性無法得到有效量化，使用與真實結(jié)果存在較大差異的數(shù)值模擬進(jìn)行設(shè)計，極有可能導(dǎo)致真實系統(tǒng)達(dá)不到預(yù)期的性能要求，引入潛在風(fēng)險。對于數(shù)值模擬模型的確認(rèn)，由于存在大量模型不確定性和參數(shù)不確定性，所以必須首先分析不確定性對數(shù)值模擬結(jié)果的影響，實現(xiàn)有效的不確定性傳播，對仿真結(jié)果進(jìn)行不確定性量化，評估其可信度，進(jìn)而進(jìn)行模型修正或重選，確保仿真結(jié)果的高可信度。事實上，模型修正通常涉及不確定性傳播的反問題求解，不確定性傳播也是其中的核心之一。不確定性量化分析在模型確認(rèn)總體技術(shù)路線中占據(jù)重要地位，該問題的研究是近年來國內(nèi)外研究人員關(guān)注的熱點。事實上，美國機(jī)械工程師協(xié)會在其計算固體力學(xué)驗證和確認(rèn)（V&V）項目指南中，已將不確定性量化分析列為模型確認(rèn)中必不可少的關(guān)鍵研究內(nèi)容[30]。圖2 以CFD 模型確認(rèn)為例，展示了不確定性傳播在其中的作用，在模型確認(rèn)過程中要反復(fù)調(diào)用不確定性傳播模塊進(jìn)行不確定性量化，這也是其中最耗時的部分。實際上由于不確定性高維、數(shù)值模擬較為耗時，“維數(shù)災(zāi)難”是不確定性傳播面臨的最大難題，因此提高其計算效率尤為重要。

圖2 不確定性傳播在模型確認(rèn)中的作用Fig. 2 UP in model validation

由于混沌多項式的高效性，其在數(shù)值模擬模型確認(rèn)中應(yīng)用非常廣泛，尤其對于CFD 數(shù)值模擬等通常較為耗時的問題，混沌多項式方法能大幅降低計算量，因此其在CFD 模型確認(rèn)的不確定性量化中應(yīng)用尤為廣泛。Meldi 等[31]將混沌多項式方法用于空間演化混合層流大渦模擬（LES）的不確定性量化，并分析了不同模擬參數(shù)的敏感性。王瑞利等[32]將非嵌入式混沌多項式方法應(yīng)用于爆轟模型，對平面爆轟問題和散心爆轟問題進(jìn)行了不確定性分析；趙輝等[33]發(fā)展了非嵌入式混沌多項式方法，研究了湍流模型系數(shù)的不確定性對RAE2822 翼型跨聲速繞流模擬的影響和材料物性參數(shù)的不確定性對燒蝕熱響應(yīng)預(yù)測的影響，為工程中多變量不確定性量化問題提供了有效的解決方案；劉全等[34]將混沌多項式方法應(yīng)用于拉式流體模型，對Sod 激波管問題進(jìn)行了不確定度量化；Enderle 等[35]將非嵌入式混沌多項式方法應(yīng)用于湍流噴霧燃燒模擬，通過對所用的雷諾平均Navier-Stokes 模型進(jìn)行不確定性分析實現(xiàn)了對變量降維，進(jìn)而實現(xiàn)了模型優(yōu)化。Schaefer 等[36-37]將混沌多項式方法用于跨聲速壁面束縛流湍流模型封閉系數(shù)的不確定性量化和敏感度分析，提出了一種基于混沌多項式方法的用于工業(yè)級規(guī)模氣動分析的不確定性量化框架，考慮了湍流模型封閉系數(shù)、來流工況、網(wǎng)格收斂誤差等不確定性。

1.3 不確定性優(yōu)化設(shè)計

確定性條件下的優(yōu)化技術(shù)已經(jīng)成功運用到諸多工程設(shè)計問題中，但是設(shè)計條件的變化，例如載荷、材料特性和操作環(huán)境的變化，往往使得工程系統(tǒng)存在很大的不確定性。以往由于數(shù)學(xué)處理和計算速度等方面的原因，通常將這些不確性量作為確定量處理，導(dǎo)致產(chǎn)品的設(shè)計性能對這些不確定性因素非常敏感，穩(wěn)健性低，或在不確定性下產(chǎn)品性能無法滿足約束，設(shè)計失效，極有可能導(dǎo)致災(zāi)難性的后果。因此，產(chǎn)生了不確定性下的優(yōu)化設(shè)計。根據(jù)設(shè)計理念的不同，可將其分為穩(wěn)健優(yōu) 化 設(shè)計（robust design optimization, RDO）[38-39]（見圖3）和基于可靠性的優(yōu)化設(shè)計（reliabilitybased design optimization, RBDO）[40-41]（見圖4），也有將二者結(jié)合的，這就是所謂的穩(wěn)健可靠性優(yōu)化。通過在設(shè)計過程中考慮設(shè)計變量、設(shè)計參數(shù)、設(shè)計決策和系統(tǒng)分析模型等不確定性因素的影響，來降低不確定性對系統(tǒng)性能的影響。如圖3 所示，穩(wěn)健優(yōu)化所得最優(yōu)解，雖然相對確定性優(yōu)化損失了一定的產(chǎn)品性能，但是卻對不確定性的敏感程度更低，具有更強(qiáng)的穩(wěn)健性。如圖4 所示，基于可靠性的優(yōu)化所得最優(yōu)解離約束邊界相對較遠(yuǎn)，雖然一定程度上可能會導(dǎo)致設(shè)計的保守，但是能保證設(shè)計始終落在可行性區(qū)域內(nèi)，系統(tǒng)是安全的；而確定性最優(yōu)解由于設(shè)計變量或參數(shù)存在不確定性，最優(yōu)解很大可能違反約束。

圖3 確定性最優(yōu)和穩(wěn)健性最優(yōu)Fig. 3 Deterministic and robust optimums

圖4 確定性最優(yōu)和可靠性最優(yōu)Fig. 4 Deterministic and reliability-based optimums

圖5 展示了不確定性優(yōu)化的流程圖。由圖可知，不確定性傳播是其中必不可少的模塊，用于計算目標(biāo)和約束函數(shù)的均值、標(biāo)準(zhǔn)差或失效概率。不確定性傳播的精度和效率幾乎決定了整個優(yōu)化的精度和效率[42]。如上所述，對于復(fù)雜工程系統(tǒng)設(shè)計，仿真分析模型通常較為耗時、非線性程度高且高維，因此計算量大是其不確定性傳播面臨的主要難題。

圖5 不確定性傳播在不確定性優(yōu)化中的地位Fig. 5 UP in design optimization under uncertainty

由于混沌多項式的高效性，其在不確定性優(yōu)化設(shè)計中也得到了大量應(yīng)用。蔡宇桐等[19]針對壓氣機(jī)葉片制造中的加工誤差，應(yīng)用非嵌入式混沌多項式進(jìn)行不確定性量化，得到了增強(qiáng)葉片氣動性能穩(wěn)健性的優(yōu)化方向；胡晚亭等[8]將混沌多項式方法應(yīng)用于Winkler 地基沉降計算，得到參數(shù)變化對地基沉降可靠性的影響，進(jìn)而明確了設(shè)計方向；李冬琴等[22]利用混沌多項式法分析多維隨機(jī)不確定性因素對船舶優(yōu)化方案的影響，完成了船舶多學(xué)科穩(wěn)健設(shè)計優(yōu)化研究，有效減少和避免了船舶設(shè)計優(yōu)化方案失效的可能性。Lie 等[43]將混沌多項式方法與禁忌搜索算法結(jié)合起來對天線陣列反問題進(jìn)行穩(wěn)健優(yōu)化，相較于拉丁超立方抽樣方法極大提高了計算效率并且保證了穩(wěn)健性能的計算精度；Mandur 等[44]在模型參數(shù)不確定性傳播的穩(wěn)健優(yōu)化問題中使用混沌多項式方法，展現(xiàn)出了比蒙特卡羅法更強(qiáng)大的計算精度和計算效率優(yōu)勢；Huang 等[45]針對探測器火星大氣進(jìn)入的軌跡穩(wěn)健優(yōu)化問題，采用混沌多項式方法獲得認(rèn)知不確定條件下軌跡性能的近似解析函數(shù)，建立了一個雙環(huán)嵌套的穩(wěn)健優(yōu)化模型。

2 混沌多項式理論基礎(chǔ)

混沌多項式方法分為嵌入式（intrusive）和非嵌入式（non-intrusive）[46]。嵌入式混沌多項式方法主要應(yīng)用于動力學(xué)系統(tǒng)的不確定性傳播，其主要思想是將不確定性源進(jìn)行混沌多項式模型表達(dá)，然后將其代入到動力學(xué)微分方程，利用Galerkin投影，將隨機(jī)微分方程轉(zhuǎn)換為一組更高維的以混沌多項式模型系數(shù)為狀態(tài)量的確定性微分方程組，通過求解該微分方程組得到混沌多項式系數(shù)，從而完成不確定性傳播[47]。由于嵌入式混沌多項式方法需要進(jìn)入系統(tǒng)方程內(nèi)部，對方程進(jìn)行擴(kuò)維修改，不可避免地會引入數(shù)值誤差。另一方面，當(dāng)前的系統(tǒng)分析模型通常都經(jīng)過了反復(fù)標(biāo)定驗證和確認(rèn)，對其進(jìn)行模型內(nèi)部調(diào)整可能性非常小。此外，隨著仿真軟件的大量應(yīng)用，復(fù)雜工程產(chǎn)品設(shè)計分析基本都涉及黑箱型響應(yīng)函數(shù)。因此，嵌入式混沌多項式方法應(yīng)用非常有限，非嵌入式混沌多項式方法將響應(yīng)函數(shù)看作一個黑箱，僅關(guān)注輸入和輸出的映射關(guān)系，構(gòu)建輸出響應(yīng)的混沌多項式代理模型，從而實現(xiàn)不確定性傳播。由于簡單易于操作，相比嵌入式混沌多項式，其研究應(yīng)用更加廣泛。本文主要介紹非嵌入式混沌多項式方法。

2.1 基本原理

以響應(yīng)函數(shù)y=g(x) (x=[x1,···,xi,···,xd])為例，混沌多項式理論下輸出響應(yīng)表示為p階截斷的混沌多項式模型。

式中：bi為 待求的混沌多項式系數(shù)； Φi為正交多項式；正交多項式的總項數(shù)為P+1=(d+p)!/d!p!；ξ=[ξ1,···,ξi,···,ξd]， 為標(biāo)準(zhǔn)隨機(jī)向量，其中 ξi與原隨機(jī)變量xi存 在一定的轉(zhuǎn)換關(guān)系，具體與xi和 ξi的分布形式有關(guān)。

一旦混沌多項式系數(shù)bi求出，就可直接在混沌多項式模型上運行MCS，得到輸出響應(yīng)y的隨機(jī)概率特性，包括均值 μy、 方差 σy、概率分布等。在式(1) 的基礎(chǔ)上，也可得到低階統(tǒng)計矩關(guān)于混沌多項式系數(shù)的如下解析表達(dá)：

式中，E 表示求期望，是一個運算符。

對于混沌多項式方法，主要涉及兩方面內(nèi)容：正交多項式 Φi的構(gòu)建和混沌多項式系數(shù)bi的求取，其中后者涉及響應(yīng)函數(shù)的調(diào)用，是計算量消耗的主要來源。

2.2 正交多項式構(gòu)建

混沌多項式方法最早由Wiener[48]提出，以埃爾米特（Hermite）正交多項式作為基函數(shù)，被稱為Wiener 混沌多項式。根據(jù)Cameron-Martin 理論，它對正態(tài)分布型輸入具有指數(shù)收斂速度，但是對于其它的分布類型，則收斂速度明顯降低。這是因為Hermite 正交多項式的權(quán)函數(shù)剛好與正態(tài)分布的概率密度函數(shù)具有相同形式。利用上述特點，Xiu 等[2]根據(jù)Askey 方案中各隨機(jī)分布類型概率密度函數(shù)與正交多項式權(quán)函數(shù)一一對應(yīng)的關(guān)系，針對不同隨機(jī)輸入類型，采用相應(yīng)類型的一元正交多項式作為基函數(shù)，對其進(jìn)行直接張量積操作（tensor product）得到由多元正交多項式Φi(ξ)(i=0,1,···,P)構(gòu)建的混沌多項式模型，將混沌多項式理論擴(kuò)展到了廣義混沌多項式（general PC,gPC）方法，使其對Askey 方案中所包含的正態(tài)、均勻、指數(shù)、貝塔（Beta）、伽馬（Gamma）等多種隨機(jī)分布類型均具有指數(shù)收斂速度。當(dāng)輸入變量不滿足Askey 方案中所列的5 種分布時，則通常將其首先轉(zhuǎn)換為標(biāo)準(zhǔn)正態(tài)分布，然后選取Hermite正交多項式作為基函數(shù)。gPC 方法極大地拓展了混沌多項式理論的應(yīng)用范圍，使其迅速被大量應(yīng)用。

對于一般的常見分布類型，gPC 方法能夠得到較好的不確定性傳播結(jié)果。但是在實際工程應(yīng)用中，存在諸多不屬于甚至與Askey 方案相差較遠(yuǎn)的分布類型，例如質(zhì)量、擴(kuò)散系數(shù)、剛度系數(shù)等物理參數(shù)均服從對數(shù)正態(tài)分布。此時，通常需要利用不確定性變換將其轉(zhuǎn)化為Askey 方案中的分布類型，即：將原隨機(jī)變量xi(i=1,···,d)轉(zhuǎn)換為Askey 方案中所示的某類標(biāo)準(zhǔn)隨機(jī)變量 ξi，而變換自然會引入一定的誤差，降低不確定性分析的精度。針對隨機(jī)分布類型的復(fù)雜性和多樣性特點，不少研究提出針對任意的隨機(jī)分布類型，自行構(gòu)建正交多項式。Witteveen 等[49]提出了基于Gram-Schmidt 正交分解的混沌多項式方法；Zhang 等[50]和Xu 等[51]提出了基于斯蒂爾吉斯（Stieltjes）過程的混沌多項式方法。這些研究分別利用Gram-Schmidt 正交分解和Stieltjes 過程，構(gòu)建相對于各維隨機(jī)輸入分布的最優(yōu)一元正交多項式基函數(shù)。在此基礎(chǔ)上，類似于gPC 通過張量積操作，構(gòu)建p階 截 斷的 多 元正 交 多項 式 Φi(x)(i=0,1,···,P)。相比于gPC，這些混沌多項式方法適用范圍更加廣泛，可應(yīng)對任意隨機(jī)輸入分布類型，無需利用變換，提高了收斂速度。

上述所有混沌多項式方法都建立在已知隨機(jī)輸入完整的概率分布函數(shù)基礎(chǔ)之上，而在實際工程應(yīng)用中隨機(jī)參數(shù)的信息可能以各種形式存在，如離散的原始數(shù)據(jù)樣本，尤其對于復(fù)雜系統(tǒng)，往往由于價格高昂、耗時太長等，難以得到其完整的概率密度函數(shù)。此時，上述基于各維隨機(jī)輸入完整概率分布函數(shù)的混沌多項式方法不再適用。為此，Oladyshkin 等[52]提出了一種數(shù)據(jù)驅(qū)動混沌多項式（data-driven PC）方法，該方法可應(yīng)對各種類型的隨機(jī)分布以及離散的原始數(shù)據(jù)（無需概率分布函數(shù)），并具有良好的收斂性和精度。Wang 等[53]在此基礎(chǔ)上進(jìn)一步將投影法（Galerkin projection）引入到數(shù)據(jù)驅(qū)動混沌多項式中計算混沌多項式系數(shù)，提出了相應(yīng)的高斯節(jié)點和權(quán)值的計算方法，進(jìn)一步提高了其靈活性。數(shù)據(jù)驅(qū)動混沌多項式方法的主要思想為，根據(jù)各維隨機(jī)輸入變量的離散數(shù)據(jù)或概率分布函數(shù)，計算其統(tǒng)計矩，利用正交多項式的正交性，推導(dǎo)矩匹配（moment-matching）方程，進(jìn)而通過匹配隨機(jī)輸入變量一定階次的統(tǒng)計矩，完成最優(yōu)一元正交多項式的構(gòu)建。需要注意的是，當(dāng)離散數(shù)據(jù)不足時，一元正交多項式的構(gòu)建會引入一定的誤差。

考慮到在實際工程應(yīng)用中，常有隨機(jī)輸入變量相關(guān)的情況。例如，結(jié)構(gòu)中的材料屬性與疲勞屬性[54]、飛行器氣動噪聲與外表面分布的隨機(jī)載荷[55]等。針對上述問題，多數(shù)研究學(xué)者采用變換方法，如正交變換[56]、Rosenblatt 變換[57]和Nataf變換[58]等，將相關(guān)隨機(jī)輸入變量轉(zhuǎn)換為相互獨立的標(biāo)準(zhǔn)正態(tài)隨機(jī)變量。但是這些變換方法均屬于非線性變換，會導(dǎo)致變換后的響應(yīng)函數(shù)呈現(xiàn)強(qiáng)非線性特征，特別是當(dāng)相關(guān)隨機(jī)輸入變量服從復(fù)雜的非正態(tài)概率分布或響應(yīng)函數(shù)非線性較強(qiáng)時，在不確定性分析中會引入非常大的計算誤差。此外，這些轉(zhuǎn)換方法均要依賴隨機(jī)變量完整的分布函數(shù)，當(dāng)相關(guān)隨機(jī)輸入變量以離散數(shù)據(jù)形式存在時，上述變換方法顯然均不再適用。為此，在數(shù)據(jù)驅(qū)動混沌多項式方法基礎(chǔ)上，Lin 等[59]提出了一種能夠直接處理輸入變量相關(guān)性的混沌多項式方法，其能避免變換帶來的誤差，且能應(yīng)對離散數(shù)據(jù)形式的相關(guān)輸入，擴(kuò)展了混沌多項式方法的適用范圍。Paulson 等[60]針對嵌入式混沌多項式方法，利用Gram-Schmidt 正交變換，提出了一種考慮輸入變量相關(guān)性的方法。隨后，Wang 等[61]將該方法應(yīng)用到了隨機(jī)潮流。這類考慮相關(guān)性的數(shù)據(jù)驅(qū)動混沌多項式方法依然采用矩匹配方法建立一元正交多項式基，只不過在求解矩匹配方程的時候引入了相關(guān)隨機(jī)輸入變量混合矩的概念，進(jìn)而考慮變量的相關(guān)性。

這類自行構(gòu)建正交多項式的方法顯然提高了混沌多項式方法的適用范圍，具有更加廣闊的應(yīng)用前景，但是構(gòu)建正交多項式基的過程難免引入近似或數(shù)值等誤差。尤其是針對隨機(jī)變量為復(fù)雜分布（如雙峰或多峰），數(shù)據(jù)驅(qū)動混沌多項式方法可能需要匹配很高階次的統(tǒng)計矩，矩匹配方程的求解極易出現(xiàn)奇異，導(dǎo)致正交多項式基的構(gòu)建精度難以保證。gPC 方法由于實現(xiàn)方便，穩(wěn)健性相對更高，且理論上增加混沌多項式階次到一定程度即可應(yīng)對任意形式的輸入分布，因此，目前依然是應(yīng)用最為廣泛的混沌多項式方法。

2.3 階次確定

混沌多項式模型的構(gòu)建首先必須給定混沌多項式模型的階數(shù)。事實上，混沌多項式模型階數(shù)p對不確定性傳播的計算量和精度具有重要影響。常用的做法是進(jìn)行階數(shù)的收斂性測試，先從較低混沌多項式模型階次出發(fā)進(jìn)行不確定性傳播，增加階次繼續(xù)進(jìn)行不確定性傳播，同時考慮先前樣本的重復(fù)利用，直到相鄰兩次的結(jié)果（如輸出響應(yīng)y的均值和方差）變化不大，則認(rèn)為當(dāng)前混沌多項式階數(shù)滿足要求。但是對于某些特殊問題，非線性程度高，分布函數(shù)甚至呈現(xiàn)非規(guī)則雙峰的情況，則需要高達(dá)p=15 的階數(shù)才能獲得滿意的精度[62]。對于工程不確定性分析，往往涉及非線性黑箱型響應(yīng)函數(shù)，對于混沌多項式模型階次的確定并非易事。Hu 和Youn[63]利用留一法（leaveone-out, Loo）評估混沌多項式的精度，構(gòu)建混沌多項式階次的循環(huán)，不斷增加混沌多項式階次，直至滿足精度要求為止，從而決定最終的混沌多項式階數(shù)。目前這方面的研究并不多，主要原因是對于一般的實際工程系統(tǒng)，通常系統(tǒng)輸出響應(yīng)主要受各維輸入變量以及其低階交叉項的影響，高階交叉項的影響較小[64]，因此，混沌多項式模型階數(shù)p=2 或3 就能滿足工程實際需求，出現(xiàn)高階混沌多項式模型的情況較少。

2.4 系數(shù)求取

混沌多項式系數(shù)求取是混沌多項式進(jìn)行不確定性傳播的關(guān)鍵，目前主要有投影法和回歸法兩類方法。

2.4.1 投影法

利用Galerkin 投影方法，將式(1)兩邊同時依次投影到各正交多項式 Φj(ξ)上，得

根據(jù)內(nèi)積的定義，并利用正交多項式的正交性，整理式(4)，得

2.4.2 回歸法

Xiong 等[62]在上述SRSM 基礎(chǔ)上，進(jìn)一步提出了加權(quán)隨機(jī)響應(yīng)面方法（Weighted SRSM），通過引入樣本權(quán)值的概率，考慮其在概率空間的分布特性，在相同樣本的情況下，提高了不確定性傳播的精度。與確定性下的響應(yīng)面相似，回歸中樣本的性能很大程度上決定了響應(yīng)面的好壞，因此決定了SRSM 不確定性傳播的精度。關(guān)于抽樣方法，Hosder 等[69]從精度和收斂性方面，對隨機(jī)抽樣、拉丁超立方設(shè)計（Latin hypercube design, LHD）和Hemmasy 抽樣方法進(jìn)行了綜合的分析比較，并推薦用兩倍于未知混沌多項式系數(shù)個數(shù)的樣本（即N=2（P+1））可以得到比較滿意的結(jié)果。Isukapalli[67]提出在Wiener 混沌多項式方法中運用Hermite 積分節(jié)點，也就是Hermite 正交多項式的根，作為回歸樣本來求解系數(shù)。在該方法的基礎(chǔ)上，可從Askey 方案中所列正交多項式出發(fā)，分別求解p+1 階正交多項式的根，對這些根（一維空間）進(jìn)行直接張量積操作，得到多維空間的全因子設(shè)計樣本，然后在這些樣本點上進(jìn)行最小二次回歸，得到混沌多項式系數(shù)。這種抽樣策略使得樣本點大多數(shù)集中在概率空間中的高頻率區(qū)域，類似于重要性抽樣的原理，因此提高了混沌多項式系數(shù)的估算精度。由于這種抽樣方法對一維樣本采用直接張量積，得到的多維樣本的組合數(shù)為(p+1)d，呈指數(shù)增長。出于計算量的考慮，目前仍局限于啟發(fā)式地選取其中部分樣本進(jìn)行回歸，通常選取概率空間出現(xiàn)頻率較大的樣本點。也有研究提出利用單項求容積法則（monomial cubature rules, MCR）來產(chǎn)生回歸樣本，構(gòu)建混沌多項式模型，這就是所謂的PC-MCR 方法。該方法運用MCR 產(chǎn)生樣本點，由于樣本點數(shù)目少，可以全部用來估算混沌多項式系數(shù)，不僅保證了不確定性傳播的精度，而且與前面提到的高斯積分點方法相比大大減少了所需的樣本數(shù)。

回歸法中由于涉及矩陣求逆等運算，針對高維問題可能會較為繁瑣耗時，因此回歸法通常適用于低維問題。隨機(jī)抽樣、拉丁超立方設(shè)計、高斯積分點、MCR 抽樣方法在混沌多項式系數(shù)計算方面具有相似的精度，MCR 方法雖然計算效率方面最為可觀，但由于在目前的科學(xué)計算軟件中并未見到成熟算法，而且可選取的樣本形式有多種，具體選何種較為主觀，因此其應(yīng)用相對較少。

3 緩解“維數(shù)災(zāi)難”的方法

混沌多項式方法具有較高的精度和效率，但是其計算量通常隨著隨機(jī)輸入的維數(shù)呈指數(shù)增長，高維下存在“維數(shù)災(zāi)難”問題，這也是目前混沌多項式在實際應(yīng)用中面臨的最大難題。關(guān)于該問題，目前產(chǎn)生了諸多方法去緩解“維數(shù)災(zāi)難”。

3.1 多項式基截斷

為了解決維數(shù)災(zāi)難，最常用的方法是減少全階混沌多項式模型中的正交多項式項數(shù)。其中，雙曲線截斷（hyperbolic truncation）策略[70]是最常用的一種，其主要思想為：一般的實際工程系統(tǒng)的輸出響應(yīng)主要受各維輸入變量以及其低階交叉項的影響，受高階交叉項的影響較小[64]，從而可以直接去除某些多項式項，達(dá)到降低計算量的目的。

傳統(tǒng)的全階混沌多項式模型從一元正交多項式基出發(fā)，利用直接張量積構(gòu)建多元正交多項式，當(dāng)混沌多項式模型階數(shù)為p時，需要滿足

式中：ai為 多元正交多項式 Φi中第i維變量對應(yīng)的一元正交多項式的階數(shù)；p為混沌多項式模型的階數(shù)，也就是式(1)中的正交多項式Φi(i=0,1,···,P)的最高階次。

在直接張量積下，式(1)中正交多項式的總項數(shù)為P+1=(d+p)!/d!p!，可見總項數(shù)隨著d的增加呈指數(shù)增長。為了降低計算量，考慮以下雙曲線截斷策略，去除部分高階交叉多項式項，以達(dá)到減少混沌多項式系數(shù)個數(shù)、降低計算量的目的。

式中：q為自定義稀疏因子；I為滿足不等式的a=(a1,···,ad)的集合。

以p={3,4,5,6}和q={1,0.75,0.5}為例，在二維空間中上述截斷方法產(chǎn)生的階數(shù)組合如圖6 所示，圖中藍(lán)色的圓圈表示滿足式（10）中不等式要求的a=(a1,···,ad)的集合。顯然，通過引入雙曲線截斷策略，確實去除了一部分高階交叉項，對降低計算量具有重要作用。而且，相同的階數(shù)p下，q越小，去除的高階交叉項越多，計算量降低越明顯。當(dāng)q=1（圖6中第一行）時，式（10）退化為常規(guī)的全階混沌多項式方法。

圖6 不同稀疏因子下一元正交多項式的階數(shù)組合Fig. 6 Order combinations of orthogonal polynomial with different sparse factors

除了這類直接去除高階交叉項多項式項的截斷策略，Hampton 和Doostan[71]提出了各向異性多項式階數(shù)（anisotropic polynomial order）的概念，在構(gòu)建混沌多項式模型的過程中，認(rèn)為 Φi中對應(yīng)于各維變量的一元多項式的最高階次并非相同，通過引入誤差指數(shù)，自適應(yīng)地搜尋最優(yōu)的各維正交多項式階數(shù)a=(a1,···,ad)，使之滿足

可見，當(dāng)pi=p(i=1,···,d)時，該方法退化為全階混沌多項式模型。實際中，輸出響應(yīng)關(guān)于各維隨機(jī)變量的非線性是不同的，因此上述考慮各向異性的多項式階數(shù)確定方法能夠在計算量一定的情況下，將更多資源配置于非線性較強(qiáng)的維度上，從而提高精度。

多項式基截斷策略能夠在一定程度上降低計算量，但是在進(jìn)行截斷的時候稀疏因子的選取較為主觀。如果選取的值過小，雖然可以很大程度降低計算量，但是由于截除了過多的高階交叉項，精度的損失較為嚴(yán)重。因此，該方法適合于高階交叉項影響較小的問題。

3.2 稀疏重構(gòu)

稀疏重構(gòu)的主要思想是在全階混沌多項式模型的基礎(chǔ)上，通過去除對輸出響應(yīng)影響不大的正交多項式項 Φi(i=0,1,···,P)，減少混沌多項式系數(shù)的個數(shù)，從而降低計算量，是當(dāng)前應(yīng)對混沌多項式“維數(shù)災(zāi)難”難題最有效的途徑之一。稀疏重構(gòu)下構(gòu)建的稀疏混沌多項式模型可表示為

其中，混沌多項式系數(shù)的個數(shù)為P1+1， 顯然P1+1 ?P+1。

對于稀疏混沌多項式方法，關(guān)鍵在于如何去發(fā)掘非重要的正交多項式項。Blatman 和Sudret[72]提出了一種自適應(yīng)算法來自動檢測重要的正交多項式項，認(rèn)為那些能夠最大程度降低預(yù)測誤差的正交多項式項為重要的項并予以保留，從而減少混沌多項式系數(shù)的個數(shù)。隨后，他們又提出利用最小角回歸（least angle regression, LAR）方法去發(fā)掘非重要的正交多項式項，進(jìn)而將其從全階混沌多項式模型中去除，以降低計算混沌多項式系數(shù)的計算量[73]。王豐剛[74]對基于LAR 的稀疏混沌多項式方法展開了研究，通過諸多算例測試發(fā)現(xiàn)，隨著維數(shù)和混沌多項式階數(shù)的增長，全階混沌多項式所需的樣本數(shù)量明顯增加，導(dǎo)致混沌多項式求解系數(shù)的回歸矩陣出現(xiàn)“病態(tài)”，難以求解，然而稀疏混沌多項式則幾乎不受此影響，同時還能保持較高的計算精度。Hu 和Youn[63]針對工程可靠性分析和設(shè)計，提出在全階混沌多項式模型的正交多項式集合中，通過引入誤差指數(shù)，自適應(yīng)地篩選誤差指數(shù)最大的雙變元正交多項式項來構(gòu)建稀疏混沌多項式模型。陳光宋等[75]采用最小絕對收縮和選擇算子（LASSO）回歸自動選擇混沌多項式的重要項及其展開系數(shù)，進(jìn)一步由混沌多項式系數(shù)解析獲得全局靈敏度系數(shù)。Cheng等[76]提出了一種可觀測響應(yīng)保持同倫下的自同態(tài)調(diào)制（D-MORPH）算法，用于構(gòu)建稀疏混沌多項式模型。Diaz 等[77]、Tsilifis 等[78]和陳江濤等[79]分別研究了利用壓縮感知（compressed sensing）方法構(gòu)建稀疏混沌多項式模型，壓縮感知方法是圖像和信號處理領(lǐng)域興起的新方法，能夠高效地重構(gòu)稀疏信號，需要的采樣點數(shù)目小于自由度的個數(shù)。若混沌多項式展開是稀疏的，可通過求解以下優(yōu)化問題得到：

式中，0 范數(shù)表示b中非零元素的個數(shù)，式中所涉及的變量符號參見式(7)?？紤]到式(13) 為NPhard 問題，求解很難，且實際應(yīng)用中考慮到測量噪聲的情況，不要求 ψb=Y精確滿足，則式(13)一般變?yōu)?/p>

在求解過程中，需要指定截斷誤差ε。除了以上通過發(fā)掘全階混沌多項式模型中非重要的正交多項式來構(gòu)建稀疏混沌多項式的方法，也有研究將混沌多項式表示為高維模型表達(dá)（high dimensional model representation, HDMR）的架構(gòu)，認(rèn)為各變量單獨作用和雙變量共同作用的HDMR函數(shù)就已經(jīng)可以較為精確地描述系統(tǒng)的輸出，忽略交叉高階項，從而構(gòu)建所謂的稀疏混沌多項式模型，以達(dá)到降低混沌多項式模型中正交多項式項數(shù)的目的[80]。

LAR、壓縮感知等稀疏重構(gòu)方法在構(gòu)建稀疏模型的過程中，需要通過迭代發(fā)掘重要的正交多項式項，尤其對于非線性較強(qiáng)且維數(shù)較高（d＞10）的問題，收斂過程可能非常慢，且收斂過程受所選取的樣本點和問題的非線性程度影響非常大。

3.3 稀疏網(wǎng)格數(shù)值積分

3.1 和3.2 節(jié)旨在通過減小全階混沌多項式模型中正交多項式的項數(shù)來降低計算量。稀疏網(wǎng)格數(shù)值積分則通過在投影法計算混沌多項式系數(shù)中采取Smoyak 算法，生成積分節(jié)點，相對于全因子節(jié)點大為減小積分節(jié)點的個數(shù)，從而緩解“維數(shù)災(zāi)難”難題。關(guān)于這方面的研究目前非常多。Winokur[81]針對混沌多項式研究了一種自適應(yīng)稀疏網(wǎng)格數(shù)值積分方法，大大降低了計算量，并將其成功應(yīng)用于2004 年9 月伊萬颶風(fēng)穿過墨西哥灣時的海洋環(huán)流模型中。Wu 等[82]將稀疏網(wǎng)格數(shù)值積分用于混沌多項式，求解了不確定性下的翼型氣動優(yōu)化。Xiong 等[83]將基于稀疏網(wǎng)格數(shù)值積分的混沌多項式方法應(yīng)用于火箭彈穩(wěn)健優(yōu)化。稀疏網(wǎng)格數(shù)值積分確實能在一定程度上降低中、低維（d＜10）不確定性傳播問題的混沌多項式系數(shù)的計算量，根據(jù)維數(shù)可以預(yù)估函數(shù)調(diào)用次數(shù)，而且具有較強(qiáng)的穩(wěn)健性，是當(dāng)前緩解維數(shù)災(zāi)難較為可靠的方法，但對于高維問題其應(yīng)對能力非常有限，計算量依然非常大。

4 多可信度混沌多項式

在基于仿真的工程設(shè)計中，分析模型（例如CFD和FEA）往往具有高度非線性、計算耗時的特點。采用混沌多項式方法直接基于高精度仿真模型進(jìn)行不確定性傳播，也同樣面臨計算量大的問題。對于工程系統(tǒng)設(shè)計，由于學(xué)科分工愈來愈細(xì)，分析方法和仿真建模手段逐漸多元化。隨著設(shè)計進(jìn)程的推進(jìn)，往往伴隨有多種不同精度和計算量的分析模型產(chǎn)生。例如，某小型飛機(jī)涉及多個固定翼和螺旋槳的氣動耦合分析，可利用多種多可信度氣動仿真工具實現(xiàn)，包括簡單低階模型（二維渦模型和葉素理論螺旋槳模型）、中精度模型（渦格法和Euler CFD，其中螺旋槳模型均采用激勵盤）、高精度模型（基于雷諾平均Navier-Stokes方程的 CFD）。此外，還有地面試驗數(shù)據(jù)和飛行試驗數(shù)據(jù)。通常認(rèn)為模型的精度越高，耗費代價越大，其可生成的樣本數(shù)量也越少。為降低計算量，充分利用多個分析模型，產(chǎn)生了多可信度建模（multi-fidelity modeling）的思想，也稱多模型融合（model fusion），在大量廉價的低精度樣本基礎(chǔ)上，利用少量高精度樣本點為導(dǎo)引或修正，通過融合不同精度和計算量的樣本數(shù)據(jù)，建立多可信度代理模型（multi-fidelity metamodel），能夠在保證代理模型精度的同時，盡可能地降低計算量。在該方面具有代表性的研究有Kennedy 與O'Hagan提出的基于高斯隨機(jī)過程及上述差值（或比值）模型的多層級co-kriging 方法[84]等。

由于多可信度建模理論在降低計算量方面效果顯著，有學(xué)者基于混沌多項式方法開展了多可信度不確定性傳播的研究，通過以少量高精度樣本為引導(dǎo)，融合大量低精度樣本數(shù)據(jù)并建立不同精度模型之間的修正混沌多項式模型，來提高低精度混沌多項式模型的可信度，從而達(dá)到降低計算量的目的。Ng 和Eldred[85]提出了基于混沌多項式和稀疏網(wǎng)格數(shù)值積分的多可信度不確定性傳播方法；Palar 等[86]提出了基于最小二次回歸的多可信度混沌多項式方法。上述多可信度混沌多項式方法采用加法修正的多可信度建模策略，即在低精度混沌多項式模型上進(jìn)行加法項混沌多項式修正，雖有效解決了特定的問題，但依然存在不足。例如：Ng 和Eldred 提出的方法需要采集稀疏網(wǎng)格數(shù)值積分點處的樣本，其數(shù)量和位置都不是任意的；Matteo[87]提出的方法因利用線性回歸法計算混沌多項式系數(shù)，雖然在采樣策略上相對靈活，但需要高/低精度樣本點嵌套。為此，Berchier 提出用低精度混沌多項式模型來預(yù)測非嵌套樣本點處的響應(yīng)值，但這必然會引入一定的預(yù)測誤差，降低了多可信度混沌多項式模型的精度。Yan 和Zhou[88]針對貝葉斯推理反問題，提出了一種基于自適應(yīng)抽樣的加法修正多可信度混沌多項式方法。Cheng 等[89]提出了一種基于高斯過程回歸的多層級多可信度稀疏混沌多項式方法。Wang等[90]基于高斯隨機(jī)過程，將多層級co-kriging 方法從確定性多可信度建模領(lǐng)域擴(kuò)展到不確定性量化，構(gòu)建一種了多可信度混沌多項式方法，研究表明，相較于常用的基于加法修正的多可信度混沌多項式方法精度大幅提高，且該方法能夠應(yīng)對模型的精度水平為非層次型的情況。除了以上兩類基于加/乘法修正和高斯隨機(jī)過程的多可信度混沌多項式方法，還有研究提出利用輸出空間映射（output space mapping）技術(shù)實現(xiàn)多可信度混沌多項式模型構(gòu)建[16]，其主要思想為建立低精度混沌多項式模型到高精度混沌多項式模型的映射關(guān)系，映射關(guān)系最常見的為線性映射。

由于多可信度建模在降低計算量方面具有巨大潛能，而且切合目前產(chǎn)品設(shè)計過程中存在多種類型分析模型或數(shù)據(jù)的現(xiàn)狀，多可信度混沌多項式研究較為活躍。目前，基于加法修正的多可信度混沌多項式方法由于形式和實施簡單，應(yīng)用最多。基于高斯隨機(jī)過程的多可信度混沌多項式方法精度和靈活性高，但由于利用高斯隨機(jī)過程，建模中參數(shù)估計存在計算繁瑣和穩(wěn)健性不足的缺點，這也是目前基于高斯隨機(jī)過程的建模方法存在的共性問題。整體而言，多可信度混沌多項式方法在降低計算量方面具有巨大潛力，但是與確定性領(lǐng)域的多可信度建模方法相似，需要多可信度分析模型預(yù)測趨勢一致，存在泛化能力不足的問題，而且如何配置各個多可信度分析模型的樣本個數(shù)以保證不確定性傳播的精度，目前也缺乏科學(xué)系統(tǒng)的方法。

5 靈敏度分析

5.1 局部方法

穩(wěn)健優(yōu)化在工程設(shè)計得到廣泛應(yīng)用，混沌多項式由于精度高被大量用于穩(wěn)健優(yōu)化?；谔荻鹊膶?yōu)算法由于效率高在穩(wěn)健優(yōu)化中應(yīng)用很多，而在基于梯度尋優(yōu)的穩(wěn)健優(yōu)化中，需要計算目標(biāo)和約束函數(shù)相對于設(shè)計變量的局部敏度信息，也稱設(shè)計靈敏度（design sensitivity）。此外，該設(shè)計靈敏度還能為設(shè)計決策者提供參考，從而在實際中對產(chǎn)品質(zhì)量波動進(jìn)行管控。常規(guī)的做法是直接利用差分法計算這些設(shè)計靈敏度，但這必然會消耗一定的計算量，尤其當(dāng)設(shè)計變量高維時，計算量較大。

混沌多項式用于局部靈敏度分析的研究不多。Ren 等[91]針對基于數(shù)值積分計算混沌多項式系數(shù)的穩(wěn)健優(yōu)化，推導(dǎo)了設(shè)計靈敏度，無需調(diào)用任何額外的響應(yīng)函數(shù)。該方法被用于基于混沌多項式的翼型氣動穩(wěn)健優(yōu)化，相比于直接依賴差分法的穩(wěn)健優(yōu)化，可降低高達(dá)30%的計算量?？紤]到系數(shù)的求取基于Galerkin 投影和高斯數(shù)值積分，見式(5)，其中分子的期望值 E[yΦi(ξ)]計算都將表示為高斯節(jié)點上的函數(shù)響應(yīng)值的加權(quán)和，最終設(shè)計靈敏度的計算歸結(jié)到期望值關(guān)于節(jié)點的偏導(dǎo)?E[yΦi(ξ)]/?xj.ij，其中xj.ij為高斯積分節(jié)點。分別針對FFNI 和SGNI，針對不同的節(jié)點分布情況，結(jié)合拉格朗日插值近似技術(shù)，利用節(jié)點上的函數(shù)值，在不調(diào)用任何響應(yīng)函數(shù)的情況下，推導(dǎo)了半解析形式的設(shè)計靈敏度。

5.2 全局方法

全局靈敏度分析（global sensitivity analysis, GSA）能夠量化各不確定性因素對系統(tǒng)響應(yīng)不確定性的影響程度，進(jìn)而可適當(dāng)忽略那些影響小的不確定性因素，在不影響精度的前提下降低不確定性傳播的計算量，目前已成為提高不確定性傳播效率的有效途徑之一。全局靈敏度分析方法大體可以分為基于回歸的方法和基于方差的方法[92]，其中基于Sobol'靈敏度指數(shù)的方差分析法因其簡單有效的特點得到了廣泛應(yīng)用。

混沌多項式方法將隨機(jī)輸出響應(yīng)表示為一組正交多項式的加權(quán)組合，基于混沌多項式模型可解析得到響應(yīng)方差，從而非常方便地計算Sobol'靈敏度指數(shù)。最終的靈敏度指數(shù)是關(guān)于混沌多項式系數(shù)的解析表達(dá)，相當(dāng)于混沌多項式系數(shù)計算的副產(chǎn)品。該方法相比于MCS 計算量大為降低，因此基于混沌多項式的Sobol'方差分析法成為一種較為常用的全局靈敏度分析方法。

基于混沌多項式方法的Sobol'靈敏度指數(shù)表示為

式中，DPC為總方差，可通過式(3) 計算得到。由式(15) 可知，一旦混沌多項式模型構(gòu)建好，就可在混沌多項式模型基礎(chǔ)上進(jìn)行全局靈敏度分析，得到關(guān)于混沌多項式系數(shù)的靈敏度指數(shù)的表達(dá)，無需任何額外的函數(shù)調(diào)用。

Sudret[92]將廣義混沌多項式應(yīng)用到全局靈敏度分析中，得到了關(guān)于混沌多項式系數(shù)的解析式的Sobol'靈敏度指數(shù)。Palar 等[93]基于加法修正多可信度混沌多項式方法，進(jìn)行了基于Sobol'靈敏度指數(shù)的全局靈敏度分析。Cheng 等[94]基于支持向量機(jī)回歸提出了用于全局靈敏度分析的自適應(yīng)混沌多項式方法。王晗等[95]提出了基于稀疏多項式混沌展開的孤島微電網(wǎng)全局靈敏度分析方法，并將其用于準(zhǔn)確、快速地辨識影響系統(tǒng)運行狀態(tài)的關(guān)鍵輸入隨機(jī)變量。卜令澤[96]對基于混沌多項式的全局靈敏度分析方法展開深入研究，為大型復(fù)雜結(jié)構(gòu)的靈敏度與可靠度分析提供了新思路。為識別影響自動裝填機(jī)構(gòu)剛度的核心關(guān)鍵參數(shù)，孫佳等[97]采用基于混沌多項式展開的全局靈敏度分析方法，從32 個自動裝填機(jī)構(gòu)參數(shù)中提取出了6 個影響剛度的核心關(guān)鍵參數(shù)。王娟[98]開展了基于混沌多項式的Sobol'全局靈敏度分析，考慮了輸入變量線性相關(guān)性的情況。也有研究首先在全階混沌多項式模型上進(jìn)行雙曲線截斷，構(gòu)建稀疏混沌多項式模型，在此基礎(chǔ)上進(jìn)行全局靈敏度分析，進(jìn)而發(fā)掘Sobol'靈敏度指數(shù)較小的非重要變量，在不確定性傳播中不予以考慮，從而實現(xiàn)降維，然后在降維后的變量空間構(gòu)建階次更高的稀疏混沌多項式，達(dá)到了提高精度并降低計算量的目的[85]。

6 混合不確定性傳播

當(dāng)考慮隨機(jī)（aleatory）和認(rèn)知（epistemic）混合不確定性時，常規(guī)的不確定性傳播過程變?yōu)橐粋€典型的雙層循環(huán)，外層考慮認(rèn)知不確定性，內(nèi)層實質(zhì)為認(rèn)知不確定性固定于某值情況下的隨機(jī)不確定性傳播。由于混沌多項式方法的高效性，近年來出現(xiàn)了不少將混沌多項式方法用于處理隨機(jī)和認(rèn)知混合不確定性傳播的研究工作。其基本思想為：外層利用概率盒（P-box）、證據(jù)、模糊、區(qū)間以及似然估計等方法考慮認(rèn)知不確定性，內(nèi)層進(jìn)行基于混沌多項式的隨機(jī)不確定性傳播，從而避免大量調(diào)用耗時的真實響應(yīng)函數(shù)，降低計算量。

6.1 概率盒方法

圖7 基于概率盒和混沌多項式的混合不確定性傳播Fig. 7 Mixed UP with PC and P-box

6.2 證據(jù)理論

證 據(jù) 理 論（evidence theory）是 由Dempster 率先提出的，后經(jīng)Shafer 系統(tǒng)完善，故又稱為Dempster-Shafer 理論[102]，它是對經(jīng)典概率理論的一種擴(kuò)展，使用概率邊界反映所有可能結(jié)果集合冪集的信任度。證據(jù)理論通過識別框架、基本可信度分配（basic probability assignment, BPA）、可信度函數(shù)Bel()和 似真度函數(shù)Pl()這3 個基本概念構(gòu)成了不確定性建模框架。讀者可參考文獻(xiàn)[29, 94, 103]獲取關(guān)于證據(jù)理論的具體介紹。利用證據(jù)理論做混合不確定性傳播時，基于混沌多項式求取輸出響應(yīng)的累積信度函數(shù)（cumulative belief function,CBF）和累積似真度函數(shù)（cumulative plausibility function, CPF）來表征識別框架內(nèi)系統(tǒng)響應(yīng)的不確定特性，如圖8 所示。在得到輸出響應(yīng)y的CBF 和CPF 之后，可進(jìn)一步建立適合不確定性優(yōu)化設(shè)計的類似于均值和方差的目標(biāo)準(zhǔn)則[104]。

圖8 證據(jù)理論下輸出響應(yīng)y 的CBF 和CPF 曲線Fig. 8 CBF and CPF of output response with Dempster-Shafer theory

6.3 模糊數(shù)

模糊數(shù)（fuzzy number）是常規(guī)實數(shù)的一般化，其含義是它不引用一個值，而是引用一組可能的值，其中每個可能的值都有自己的權(quán)重（范圍0～1），稱為隸屬函數(shù)。在進(jìn)行混合不確定性傳播時，可以將模糊數(shù)看作一種靈活的P-box 的形式，即對每個 αi-cut 的區(qū)間進(jìn)行基于混沌多項式的不確定性傳播[105-106]，最終得到系統(tǒng)輸出響應(yīng)y的隸屬度函數(shù)。圖9 為大致過程的示意圖，其中縱軸表示隸屬度函數(shù)p(x)值 ， α ∈[0,1]是 α-cut 水平。

圖9 基于混沌多項式和模糊理論的混合不確定性傳播Fig. 9 Mixed uncertainty propagation with PC and fuzzy theory

6.4 區(qū)間理論

區(qū)間模型一般定義如下：

式中，上標(biāo)I，L，U 分別表示區(qū)間、區(qū)間下界和區(qū)間上界?；趨^(qū)間數(shù)的混合不確定性傳播可轉(zhuǎn)換為在區(qū)間變量的范圍內(nèi)求解優(yōu)化問題[107-108]。對于任意的區(qū)間變量所在區(qū)間上的某個值，固定認(rèn)知不確定性變量于該值，構(gòu)建隨機(jī)變量的混沌多項式模型，從而可方便地得到輸出響應(yīng)的均值、方差或失效概率。進(jìn)一步通過尋優(yōu)，得到輸出響應(yīng)均值、方差或失效概率的最大值和最小值。

6.5 似然方法

不同于以上4 種方法，似然方法（likelihood）[109]最終依然采用概率理論來進(jìn)行不確定性傳播。似然估計允許同時處理點數(shù)據(jù)和區(qū)間數(shù)據(jù)，它將以稀疏點數(shù)據(jù)和區(qū)間數(shù)據(jù)存在的認(rèn)知不確定性輸入變量表示為概率形式的不確定性，通過最大化似然函數(shù)，得到認(rèn)知不確定性變量的概率密度函數(shù)（PDF），該P(yáng)DF 相當(dāng)于是在當(dāng)前已知數(shù)據(jù)下，該變量的平均PDF。概率表征使得該方法可非常方便地應(yīng)用于各類基于概率的不確定性分析及優(yōu)化設(shè)計理論和方法，而基于概率的方法具有嚴(yán)格的理論基礎(chǔ)，發(fā)展較為成熟，因此應(yīng)用起來非常方便。

7 總結(jié)與展望

混沌多項式方法作為一種較為成熟高效的不確定性傳播方法，近幾年得到了長足發(fā)展，已經(jīng)在各個領(lǐng)域得到了廣泛應(yīng)用。在穩(wěn)健優(yōu)化設(shè)計領(lǐng)域，由于需要在優(yōu)化的每個迭代點構(gòu)建混沌多項式模型，對于高維且非線性較強(qiáng)的優(yōu)化問題，往往涉及大量迭代，計算量顯著上升，而數(shù)值模擬模型確認(rèn)不涉及該迭代，相比于傳統(tǒng)的不確定性傳播方法，混沌多項式的優(yōu)勢凸顯，因此混沌多項式方法在數(shù)值模擬模型確認(rèn)方面最具應(yīng)用潛力。相比其它混沌多項式的變種，廣義混沌多項式方法由于具有相對較強(qiáng)的穩(wěn)健性，針對不同的輸入分類和函數(shù)形式通常能得到較為滿意的結(jié)果，是目前應(yīng)用最廣泛的混沌多項式方法。但其大部分應(yīng)用基本停留于較簡單的工程問題，問題涉及的規(guī)模及維數(shù)都不太高（d＜10），高維問題鮮有報道。在應(yīng)對“維數(shù)災(zāi)難”方面，稀疏網(wǎng)格方法穩(wěn)健性好、泛化能力相對較強(qiáng)，但僅適用于中、低維問題（d＜10）；稀疏重構(gòu)策略計算量降低最為顯著，是最具應(yīng)用潛力的途徑之一，但在適用范圍和泛化能力上還需改進(jìn)。多可信度混沌多項式是降低計算量的有效途徑之一，但是對于高、低精度模型樣本數(shù)目的確定還缺乏較為科學(xué)有效的方法，目前基本停留在不斷試湊的模式，且泛化能力較差。此外，實際問題大都涉及黑箱型響應(yīng)函數(shù)，圍繞混沌多項式開發(fā)較為通用化的程序軟件，根據(jù)實際計算資源預(yù)算自動確定有效的混沌多項式階次和樣本點個數(shù)，以方便實際工程應(yīng)用，建立混沌多項式在各類典型問題應(yīng)用的演示實例，這也是需要進(jìn)一步探索的問題。實際工程中可能涉及高達(dá)上百維的不確定性傳播問題，目前的混沌多項式方法依然會面臨較為突出的“維數(shù)災(zāi)難”問題，這方面還需挖掘新的理論和方法，例如嘗試采用小樣本深度學(xué)習(xí)技術(shù)構(gòu)建輸入輸出的深度神經(jīng)網(wǎng)絡(luò)模型，取代原系統(tǒng)響應(yīng)分析，以降低計算量。