胡政文，張保強(qiáng)，鄧振鴻

廈門大學(xué) 航空航天學(xué)院，廈門 361000

由于制造、測量、計算及模型本身中存在的各種誤差，航空航天系統(tǒng)參數(shù)化建模和多學(xué)科優(yōu)化需要考慮不確定性的影響[1-4]。在不確定性分析時，不確定性參數(shù)按照類別大致可分為隨機(jī)不確定性參數(shù)和認(rèn)知不確定性參數(shù)[5]。不同類型不確定性參數(shù)混合作用與單種不確定性參數(shù)作用導(dǎo)致的模型輸出不確定性存在明顯的差異，因此模型分析過程中須同時考慮隨機(jī)和認(rèn)知不確定性的影響[6-8]。Williamson和Downs[9]在累積分布函數(shù)的基礎(chǔ)上引入概率盒來處理隨機(jī)和認(rèn)知不確定性混合量化問題。概率盒可以讓隨機(jī)和認(rèn)知不確定性在模型中分別傳播[10]，并保留2種不確定性各自的特性，因此在混合不確定性量化分析中得到了廣泛應(yīng)用[11-13]。

同時由于涉及多層次和多學(xué)科，航空航天系統(tǒng)參數(shù)的維數(shù)通常較多[14]。學(xué)科間耦合作用下系統(tǒng)通常呈現(xiàn)高度非線性，難以直接判斷輸入和輸出之間的關(guān)系，高維參數(shù)給系統(tǒng)研究分析帶來了挑戰(zhàn)[15]。因此系統(tǒng)建模前通常須考慮減少輸入?yún)?shù)的維數(shù)，以往研究中主要有兩種思路處理維數(shù)簡化問題：一種是參數(shù)選擇，這種思路認(rèn)為在所有的輸入?yún)?shù)中只有少數(shù)真正與模型輸出相關(guān)，其余參數(shù)對模型輸出不確定性影響較小可以忽略，參數(shù)選擇通常采用靈敏度分析的方法[16]；第二種是降維方法，與參數(shù)選擇相比，降維方法假設(shè)所有輸入?yún)?shù)都影響模型的輸出結(jié)果，但只在參數(shù)的某些特定組合下才能最大程度體現(xiàn)出來，而降維方法的目標(biāo)就是構(gòu)造這類組合[17]。

靈敏度分析主要研究輸入不確定性對輸出不確定性的貢獻(xiàn)大小。隨著模型復(fù)雜性的增加，輸入?yún)?shù)的靈敏度信息逐漸成為建模時需要考慮的因素之一[18-20]。靈敏度分析方法主要分為局部靈敏度分析和全局靈敏度分析，全局靈敏度分析方法主要有Sobol’指標(biāo)、矩靈敏度指標(biāo)等[21-22]。傳統(tǒng)的靈敏度分析方法大都基于概率框架，然而實際中模型通常包含認(rèn)知不確定，須同時考慮區(qū)間的靈敏度分析[21, 23]。Oberguggenberger等[24]闡述了概率和非概率框架下的靈敏度分析方法。近年來，隨著概率盒在混合不確定性量化中的應(yīng)用，概率盒下的全局靈敏度研究越來越多。Ferson和Troy Tucker[25]研究了在混合不確定性概率邊界分析時通過縮減輸入?yún)?shù)的不確定性進(jìn)行靈敏度分析。Song等[26]提出一種擴(kuò)展蒙特卡羅模擬方法，計算概率盒下各輸入的基于方差的全局靈敏度指標(biāo)。Bi等[27]使用Bhattacharyya距離定量描述概率盒下的不確定性參數(shù)靈敏度。不確定性縮減法分析概率盒全局靈敏度具有易于理解、工程中實現(xiàn)較為容易的優(yōu)勢，逐漸得到研究人員的關(guān)注[28]。

另一方面，降維方法在航空發(fā)動機(jī)、衛(wèi)星、汽車、芯片、圖像處理等許多領(lǐng)域中得到了廣泛應(yīng)用[14]。Fodor[29]回顧了常用的降維方法，包括主成分分析、因子分析、隨機(jī)投影法、連續(xù)拓?fù)溆成洹⒒貧w等。主成分分析是普遍使用的線性降維方法之一[30-31]，該方法基于參數(shù)的協(xié)方差矩陣，可以有效地降低均方誤差，也被稱為本征正交分解、Karhunen-Loève (KL) 變換、霍特林變換和經(jīng)驗正交函數(shù)法等[32-33]。Constantine等[34]通過構(gòu)造輸出的梯度協(xié)方差矩陣來判斷輸入?yún)?shù)空間中的主要方向，提出了活躍子空間降維方法?；钴S子空間法推廣了傳統(tǒng)的主成分分析思想，對處理高維復(fù)雜系統(tǒng)具有一定的優(yōu)勢[35]。Lukaczyk等[36]將活躍子空間法應(yīng)用于ONERA-M6跨聲速機(jī)翼的形狀設(shè)計優(yōu)化中，在設(shè)計空間區(qū)域內(nèi)找到了最佳設(shè)計變量。Jefferson等[37]將活躍子空間法應(yīng)用在綜合水文模型的靈敏度分析。Hu等[38]使用活躍子空間法對衛(wèi)星系統(tǒng)進(jìn)行多學(xué)科優(yōu)化設(shè)計。

以往處理高維參數(shù)的簡化問題的研究中，一般使用靈敏度分析或降維方法兩者之一，對于高維混合不確定性系統(tǒng)單種方法降維效果不明顯。本文將靈敏度分析和降維方法結(jié)合在一個框架中，提出一種結(jié)合概率盒全局靈敏度和活躍子空間的跨層降維方法。案例分析中將跨層降維方法應(yīng)用到NASA多學(xué)科不確定性量化挑戰(zhàn)問題中，驗證提出的跨層降維方法的有效性。

1 跨層降維方法基本理論

1.1 概率盒全局靈敏度分析

概率盒定義為一組包含變量所有可能累積分布函數(shù)(Cumulative Distribution Function, CDF)曲線的邊界，可以有效表示混合不確定性影響下的參數(shù)特征，圖1為典型概率盒示意圖。概率盒不確定性分析提供了一種較為直觀的方法來描述受到隨機(jī)和認(rèn)知不確定性綜合作用的輸出特性。

圖1 典型概率盒示意圖Fig.1 Schematic of typical probability box

不確定性縮減法可以用于概率盒框架下不確定性參數(shù)的全局靈敏度評估。輸入?yún)?shù)的靈敏度取決于參數(shù)自身不確定性的大小以及對模型輸出不確定性的影響，概率盒全局靈敏度可以表示為[25]

(1)

式中：B為原始基準(zhǔn)輸入?yún)?shù)；T為縮減不確定性后的輸入?yún)?shù)；unc(·)表示輸出不確定性的度量;S為全局靈敏度指標(biāo)。若unc(·)為標(biāo)量，則靈敏度指標(biāo)S也是標(biāo)量，由此可以對輸入?yún)?shù)進(jìn)行靈敏度排序。輸入?yún)?shù)不確定性縮減后，式(1)中分子unc(T)的值會下降。與基于方差的Sobol’指標(biāo)不同，由于不確定性的減少，對于所有輸入?yún)?shù)，unc(T)均小于unc(B)，因此靈敏度指標(biāo)S范圍是[0,1]。對各輸入?yún)?shù)依次應(yīng)用不確定性縮減法，由式(1)可以得到各輸入?yún)?shù)的S值，將所有的S按照大小進(jìn)行排序即得到對應(yīng)參數(shù)的靈敏度排序。

輸入?yún)?shù)不確定性的縮減有多種方法，不同縮減方式得到的不確定性度量unc(·)的結(jié)果不同。通常有以下2種方式：① 使用特定的概率分布替代不確定性輸入?yún)?shù)，特定的概率分布法只消除參數(shù)的認(rèn)知不確定性，不影響其隨機(jī)不確定性；② 使用固定值替代不確定性輸入?yún)?shù)，固定值法同時消除參數(shù)的隨機(jī)不確定性和認(rèn)知不確定性。

不確定性度量unc(·)也有多種不同的定義方法(如方差、矩等)，不同的定義方式可以用于解決不同的靈敏度分析問題[25]。概率盒是參數(shù)所有可能累計分布函數(shù)的包絡(luò)線，沿著概率盒的上界和下界進(jìn)行積分可以得到概率盒的面積。本文結(jié)合概率盒的特征，采用概率盒的面積作為輸出的不確定性度量unc(·)，此時不確定性縮減法的流程見圖2，式(1)轉(zhuǎn)化為

圖2 概率盒靈敏度分析的不確定性縮減法流程圖Fig.2 Flow chart of pinching method for probability box sensitivity analysis

(2)

式中：A(·)表示概率盒面積。

不確定性縮減法主要比較縮減輸入?yún)?shù)不確定性后輸出概率盒的面積占比下降值，進(jìn)而判斷參數(shù)的靈敏度信息。

1.2 活躍子空間法

假設(shè)原始輸入?yún)?shù)的維數(shù)為n，表示為x=[x1,x2,…,xn]T，降維方法的目的是在滿足損失數(shù)據(jù)信息最少的條件下，通過某種方法將原始的n維參數(shù)x用m維的s=[s1,s2,…,sm]T(m

(3)

對于參數(shù)x，其標(biāo)準(zhǔn)參數(shù)空間可以視為邊界為1的超立方體。非標(biāo)準(zhǔn)的參數(shù)空間可以通過正則化轉(zhuǎn)化為標(biāo)準(zhǔn)參數(shù)空間。假設(shè)在x的參數(shù)空間中隨機(jī)抽樣，樣本數(shù)為N，則輸入?yún)?shù)x可以表示為矩陣X={xi,j:1≤i≤n,1≤j≤N}。對于任意的xi,j，可正則化為

(4)

目標(biāo)函數(shù)f梯度的協(xié)方差矩陣定義為

(5)

(6)

C=WΛWT

(7)

式中：W為n×n的特征向量矩陣；Λ=diag(λ1,λ2, …,λn)是非負(fù)特征值組成的對角矩陣。特征值大小反映了目標(biāo)函數(shù)f在x參數(shù)空間中沿著對應(yīng)特征向量微小擾動的平均變化。將特征值按照降序進(jìn)行排列λ1≥λ2≥…≥λn，特征值越小表示對應(yīng)的特征向量在梯度協(xié)方差矩陣C中的貢獻(xiàn)度越低，降維時可以適當(dāng)舍棄這些向量。

活躍子空間確定基向量維數(shù)有多種方法：通過設(shè)定固定閾值λ0，只保留大于λ0的特征值及對應(yīng)的特征向量；定義κ為前m項特征值占所有特征值的累加比例：

(8)

(9)

(10)

式中：U為活躍子空間的基向量，包含W的前m列；Λ1為對應(yīng)的特征值對角矩陣。確定活躍子空間的基向量后，原輸入的參數(shù)空間可以投影到該活躍子空間：

s=UTx

(11)

此時目標(biāo)函數(shù)f(x)可以近似表示為活躍子空間內(nèi)降維變量s的函數(shù)：

f(x)≈g(UTx)=g(s)

(12)

g(s)的定義域為γ={s=UTx,x∈Rn}∈Rm。在降維變量s的基礎(chǔ)上構(gòu)建低維代理模型：

g(s)≈g*(s)≡R(s;g1,g2,…,gM)

(13)

式中：M為代理模型訓(xùn)練樣本數(shù)；g*(s)≡R為活躍子空間中定義的代理模型，可以降低計算模型的復(fù)雜性和減少計算時間和成本。原函數(shù)和代理模型的誤差可以表示為[34]

(14)

系數(shù)C1取決于x和ρ(x)。由式(14)可以看出，舍棄的特征值之和越小，即累加比例κ的值越大，原函數(shù)和新函數(shù)誤差也越小，即活躍子空間降維法所產(chǎn)生的誤差也就越小。若前m項特征值累加比例κ>90%，則可以選擇前m項特征向量組成基向量用于降低輸入?yún)?shù)空間的維度。

1.3 跨層降維方法

在不確定性縮減法和活躍子空間的基礎(chǔ)上，提出結(jié)合概率盒全局靈敏度和活躍子空間的跨層降維方法。

圖3是跨層降維方法的流程圖，基本步驟如下：

圖3 跨層降維方法流程圖Fig.3 Flow chart of cross-layer method for dimension reduction

步驟1整理系統(tǒng)的不確定性參數(shù)，將其按照隨機(jī)、認(rèn)知和混合不確定性進(jìn)行分類。

步驟2使用概率盒法對初始模型進(jìn)行不確定性分析，得到模型輸出的初始概率盒，計算概率盒的面積。

步驟3縮減待分析的輸入?yún)?shù)的不確定性，對模型進(jìn)行分析得到縮減不確定性后輸出概率盒的面積。

步驟4比較縮減參數(shù)不確定性前后輸出概率盒的面積改變量，計算參數(shù)的靈敏度指標(biāo)S值。

步驟5選取另一輸入?yún)?shù)，重復(fù)步驟3和步驟4，直至完成所有不確定性參數(shù)的靈敏度分析。

步驟6對各輸入?yún)?shù)的S值按照降序排列，將S值較小的輸入?yún)?shù)固定。

步驟7對保留不確定性的輸入?yún)?shù)進(jìn)行抽樣，取樣本數(shù)N，計算模型輸入并得到梯度協(xié)方差矩陣。

步驟8對梯度協(xié)方差矩陣進(jìn)行特征分解，按照特征值大小對特征值進(jìn)行排序，確定保留的活躍子空間的維數(shù)和基向量。

步驟9對輸入?yún)?shù)使用活躍子空間的基向量進(jìn)行降維。

2 NASA多學(xué)科不確定性量化挑戰(zhàn)問題

為了驗證提出的概率盒全局靈敏度和活躍子空間的跨層降維方法，本文以NASA Langley研究中心提出的多學(xué)科不確定性量化挑戰(zhàn)問題為案例[40]。

2.1 挑戰(zhàn)問題描述

挑戰(zhàn)問題的數(shù)學(xué)模型描述了地面基站無線電波操控遙控飛機(jī)的通用傳輸動力學(xué)特性，如圖4所示[40]。要求對不確定性系統(tǒng)模型進(jìn)行綜合評估，解決不確定性量化、全局靈敏度分析、極限狀態(tài)分析和穩(wěn)健優(yōu)化設(shè)計中的一系列關(guān)鍵問題。

圖4 NASA多學(xué)科不確定性量化挑戰(zhàn)問題模型示意圖[40]Fig.4 Model schematic of NASA multidisciplinary uncertainty quantification challenge problem[40]

模型中的不確定性參數(shù)按照類型劃分為隨機(jī)、認(rèn)知以及混合不確定性參數(shù)3類。Ψ表示所研究的系統(tǒng)模型，p表示Ψ中不確定性輸入?yún)?shù)(p∈R21)，d表示設(shè)計參數(shù)(d∈R14)，g表示Ψ性能約束(g∈R8)，g取決于p和d。不確定性輸入?yún)?shù)p、設(shè)計參數(shù)d和約束g之間的關(guān)系可以表示為

x=h(p)

(15)

g=f(x,d)

(16)

式中：x是參數(shù)為p的中間函數(shù)(x∈R5)；g和x假定為連續(xù)可微函數(shù)。式(15)展開為

(17)

x的每個分量都表示系統(tǒng)一個學(xué)科的輸出，x1～x5均為標(biāo)量。若系統(tǒng)Ψ滿足對于g的所有分量均滿足gi≤0,i=1,2,…,8，則可以認(rèn)為該系統(tǒng)是可靠的，符合設(shè)計要求。對于確定設(shè)計參數(shù)d，滿足約束向量g≤0的輸入?yún)?shù)p的集合稱為安全域，反之稱為失效域。g的分量有一個不滿足約束條件，對應(yīng)的參數(shù)p就落在失效域內(nèi)。圖5為挑戰(zhàn)問題系統(tǒng)模型Ψ中各參數(shù)變量的層次關(guān)系示意圖。

圖5 NASA挑戰(zhàn)問題模型參數(shù)層次關(guān)系Fig.5 Hierarchical relationship of model parameters in NASA challenge problem

2.2 全局靈敏度分析

由于系統(tǒng)模型Ψ的輸入?yún)?shù)中包含隨機(jī)、認(rèn)知以及混合不確定性，首先對Ψ進(jìn)行概率盒不確定性分析，研究系統(tǒng)的不確定性傳遞和量化的相關(guān)問題，確定不確定性輸入?yún)?shù)p對學(xué)科輸出x的全局靈敏度，以便后續(xù)活躍子空間降維過程中固定靈敏度低的參數(shù)。

輸入?yún)?shù)p的不確定性通過函數(shù)h傳播，得到x1、x2、x3、x4和x5的概率盒，縮減p中的認(rèn)知和混合不確定性參數(shù)時，對應(yīng)x概率盒的面積也會改變。

由式(17)可知，子學(xué)科Ψ1的輸出x1是p1、p2、p3、p4和p5的函數(shù)，表1給出了參數(shù)p1～p5的具體分布信息，其中p1為單峰Beta分布族組成的混合不確定性參數(shù)，p2是區(qū)間不確定性參數(shù)，p3服從[0,1]均勻分布，p4和p5為相關(guān)的二元正態(tài)分布族組成的混合不確定性參數(shù)。表中:Δ表示參數(shù)的區(qū)間范圍;E(·)、V(·)和ρ分別代表均值、方差和相關(guān)系數(shù)。隨機(jī)不確定性是系統(tǒng)模型的固有屬性，對模型輸出影響較大，為減少計算成本，只對認(rèn)知和混合不確定性參數(shù)按照S值大小排序。

表1 x1對應(yīng)不確定性輸入?yún)?shù)p1～p5

對于子學(xué)科Ψ1，對初始條件下的x1進(jìn)行概率盒不確定性分析。取1 000條x1的樣本累積概率分布曲線，得到初始條件下x1的概率盒，如圖6 所示，此時x1概率盒的面積為0.229。

圖6 初始條件下的x1概率盒Fig.6 Probability box of x1 under initial conditions

使用不確定性縮減法分別縮減參數(shù)p1、p2、p4、p5的不確定性。由于p1為混合不確定性參數(shù)，為計算p1對x1的概率盒全局靈敏度，不確定性縮減法的固定值取為p1的期望值0.7。p2為認(rèn)知不確定性參數(shù)，對應(yīng)的縮減值取為區(qū)間中值0.5。p4～p5服從二維正態(tài)分布，相關(guān)系數(shù)ρ的絕對值在0～1范圍內(nèi)，在計算p4的靈敏度時，取p4為期望值0，p5仍然服從正態(tài)分布，得到p4的靈敏度，同理計算p5的概率盒全局靈敏度。綜合考慮計算精度以及效率成本，結(jié)合x1概率盒邊界情況，每個參數(shù)均取1 000條累積分布函數(shù)曲線，概率盒邊界趨于穩(wěn)定，進(jìn)而得到縮減不確定性后各參數(shù)對應(yīng)的x1概率盒，如圖7所示。運用式(2)計算得到各參數(shù)的靈敏度指標(biāo)S值，計算結(jié)果列于表2之中。由圖7和表2可以看出，子學(xué)科Ψ1的參數(shù)靈敏度排序依次是p1>p5>p4>p2。由靈敏度指標(biāo)S值的計算結(jié)果，p2和p4對x1概率盒的靈敏度均小于0.1，可以將p2和p4視為固定的參數(shù)。

表2 x1概率盒的面積及對應(yīng)輸入?yún)?shù)靈敏度排序

圖7 縮減各參數(shù)不確定性對應(yīng)的x1概率盒Fig.7 Probability boxes of x1 after pinching corresponding uncertain parameters

其余子學(xué)科x2、x3、x4和x5對應(yīng)的16個不確定性輸入?yún)?shù)p6～p21的具體分布信息如表3所示。

表3 x2～x5對應(yīng)不確定性輸入?yún)?shù)p6～p21Table 3 Description of uncertain input parametersp6-p21 corresponding to x2-x5

表4 x2概率盒的面積及對應(yīng)輸入?yún)?shù)靈敏度排序

表5 x3概率盒的面積及對應(yīng)輸入?yún)?shù)靈敏度排序

表6 x4概率盒的面積及對應(yīng)輸入?yún)?shù)靈敏度排序

圖8 縮減參數(shù)p21不確定性對應(yīng)的x5概率盒Fig.8 Probability boxes of x5 after pinching p21

綜合概率盒靈敏度分析的計算結(jié)果，篩選出p1、p5、p6、p7、p12、p16、p17、p18、p21共9個影響程度較大的認(rèn)知或混合不確定性參數(shù)以及p3、p9、p11、p19共4個隨機(jī)不確定性參數(shù)。

2.3 活躍子空間降維

根據(jù)全局靈敏度分析的結(jié)果，共有13個影響較大的不確定性參數(shù)，應(yīng)用概率盒全局靈敏度和活躍子空間的跨層降維方法，采取活躍子空間進(jìn)一步降低不確定性參數(shù)的維度。

為描述系統(tǒng)模型的可靠性，定義如下系統(tǒng)性能度量函數(shù)：

(18)

性能度量函數(shù)ζ描述了系統(tǒng)模型最容易超出約束的g分量，間接反映了系統(tǒng)失效的最大可能性。為構(gòu)建合適的活躍子空間基向量，將ζ作為目標(biāo)函數(shù)。由拉丁超立方抽樣選取300個輸入?yún)?shù)樣本，分別計算對應(yīng)目標(biāo)函數(shù)ζ的值，并根據(jù)式(6) 構(gòu)造ζ的近似梯度協(xié)方差矩陣C。

對梯度協(xié)方差矩陣進(jìn)行特征分解，將特征值按照從大到小排序。表7列出了C的13個特征值，圖9為特征值隨著維度的衰減曲線圖。由表7可知，矩陣C的一維對應(yīng)的特征值為0.063 84，遠(yuǎn)遠(yuǎn)大于二維的特征值1.400×10-17，且一維特征值占所有特征值和的比例κ超過99.99%。由式(8) 和式(14)可知，對于目標(biāo)函數(shù)ζ，活躍子空間維數(shù)為一，選取一維對應(yīng)的特征向量作為活躍子空間的基向量。此時式(11)中基向量，U=[0.256，-0.087 4，-0.010 1，0.015 9，0.025 5，-0.082 6，-0.033 6，0.026 3，-0.019 0，-0.038 2，0.088 8，0.103，0.94 7]。

圖9 梯度協(xié)方差矩陣的特征值衰減曲線Fig.9 Curve of eigenvalues of gradient covariance matrix

表7 梯度協(xié)方差矩陣的特征值Table 7 Eigenvalues of gradient covariance matrix

為進(jìn)一步驗證活躍子空間選取的基向量的準(zhǔn)確性，可以使用統(tǒng)計學(xué)中Bootstrap自助法重復(fù)抽取樣本組成新的梯度協(xié)方差矩陣進(jìn)行活躍子空間的基向量計算。取Bootstrap自助法重復(fù)次數(shù)為1 000，圖10和圖11分別為活躍子空間基向量各分量分布柱狀圖和區(qū)間，由圖可知活躍子空間基向量的各分量分布較為集中，可以認(rèn)為求得的U是準(zhǔn)確的。

圖10 Bootstrap自助法(1 000次)活躍子空間基向量分量分布柱狀圖Fig.10 Histograms of components of active subspace vector using Bootstrap (1 000 times)

圖11 Bootstrap自助法(1 000次)活躍子空間基向量分量分布區(qū)間Fig.11 Ranges of components of active subspace vector using Bootstrap (1 000 times)

確定活躍子空間的基向量后，由式(11)，原13維輸入的參數(shù)空間可以投影到一維的活躍子空間，圖12是活躍子空間降維變量和目標(biāo)函數(shù)ζ在樣本數(shù)據(jù)的對應(yīng)關(guān)系，由圖可以清晰看出目標(biāo)函數(shù)隨著降維變量的變化趨勢，有利于進(jìn)一步不確定性代理模型的構(gòu)建及模型參數(shù)優(yōu)化工作。

圖12 樣本下降維變量和目標(biāo)函數(shù)的關(guān)系Fig.12 Relationship between dimension-reduction variables and objective function based on samples

3 結(jié) 論

針對航空航天系統(tǒng)不確定性建模過程中存在不確定性參數(shù)類型混雜并且維數(shù)眾多問題，本文以概率盒為不確定性分析方法，提出了一種結(jié)合概率盒全局靈敏度和活躍子空間的跨層降維法。主要研究如下：

1) 針對NASA多學(xué)科不確定性量化挑戰(zhàn)問題，使用概率盒不確定性分析方法處理x1～x5混合不確定性量化問題，使得隨機(jī)、認(rèn)知和混合不確定性在模型中有效傳播。

2) 基于概率盒的計算結(jié)果，使用不確定性縮減法對p1～p21認(rèn)知和混合不確定性參數(shù)進(jìn)行全局靈敏度分析，得到對應(yīng)參數(shù)的靈敏度指標(biāo)S值，將輸入?yún)?shù)按照靈敏度進(jìn)行排序。

3) 在全局靈敏度分析的基礎(chǔ)上，固定8個S值小的不確定性參數(shù)，對其余13個輸入變量使用活躍子空間降維法。根據(jù)輸出的梯度協(xié)方差矩陣特征值的大小確定活躍子空間的維數(shù)和降維基向量，將原始的輸入?yún)?shù)空間映射到一維的活躍子空間內(nèi)，降低了系統(tǒng)模型輸入?yún)?shù)的維度。

4) 結(jié)果顯示提出的跨層降維方法能夠?qū)?1維不確定性參數(shù)維數(shù)降低為一維，證明了該方法降低混合不確定性高維復(fù)雜系統(tǒng)輸入?yún)?shù)維度的有效性，顯示出對于處理高維不確定復(fù)雜系統(tǒng)、構(gòu)建精度符合要求且計算效率較高的模型方面的潛力和優(yōu)勢，為進(jìn)一步的模型參數(shù)優(yōu)化奠定基礎(chǔ)。