方開泰、劉民千、周永道《試驗(yàn)設(shè)計(jì)與建?！氛n件.ppt

《方開泰、劉民千、周永道《試驗(yàn)設(shè)計(jì)與建?！氛n件.ppt》由會(huì)員分享，可在線閱讀，更多相關(guān)《方開泰、劉民千、周永道《試驗(yàn)設(shè)計(jì)與建模》課件.ppt（22頁珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

1、1,第六章計(jì)算機(jī)試驗(yàn),2,6.1引言,計(jì)算機(jī)仿真試驗(yàn)這里屬于某個(gè)函數(shù)族，例如連續(xù)，且平方可積，希望可找到一個(gè)近似模型來近似原模型。,3,6.1引言,近似模型,計(jì)算機(jī)試驗(yàn),,,,4,充滿空間的設(shè)計(jì),完全隨機(jī)抽樣分層隨機(jī)抽樣：超拉丁方抽樣(Latinhypercubesampling)確定性方法：均勻設(shè)計(jì)、minimax或maximin隨機(jī)和確定性混合的方法：隨機(jī)化正交陣列，隨機(jī)化網(wǎng)絡(luò),5,6.2超拉丁方抽樣,設(shè)P=x1,,xn為在X上的試驗(yàn)點(diǎn)，欲通過P來估計(jì)y在X上的均值簡單隨機(jī)抽樣若試驗(yàn)點(diǎn)之間負(fù)相關(guān)，可減少方差,6,A.超拉丁方抽樣,取s個(gè)獨(dú)立的1,,n的隨機(jī)置換j(1),,j(n)，j=1,

2、,s,組成一個(gè)ns矩陣，稱為超拉丁方設(shè)計(jì)(Latinhypercubedesign，LHD)，記為LHD(n,s)，它的第k行j列的元素記為j(k)。取0,1上ns個(gè)均勻分布的獨(dú)立樣本，UijU(0,1),i=1,,n,j=1,,s,記xk=(xk1,,xks)T,其中則設(shè)計(jì)D=x1,,xn即為一個(gè)LHS設(shè)計(jì)，并記為LHS(n,s)。,7,例6.1構(gòu)造超拉丁方抽樣,,8,中點(diǎn)超拉丁方抽樣,,9,LHS的優(yōu)缺點(diǎn),它很容易產(chǎn)生；它可以處理試驗(yàn)次數(shù)n與因素個(gè)數(shù)s較大的問題；與完全隨機(jī)抽樣相比，它估計(jì)y的樣本均值的樣本方差要小有些LHS設(shè)計(jì)會(huì)顯得很不均勻,優(yōu)點(diǎn)：,缺點(diǎn)：,10,B.隨機(jī)化正交列,強(qiáng)度為

3、r的正交列(OA(n,s,q,r))：任意m(mr)列都構(gòu)成完全因子設(shè)計(jì)正交設(shè)計(jì)：r=2U型設(shè)計(jì)：r=1,11,構(gòu)造步驟,選擇合適的正交列OA(n,s,q,r)，記為A，并設(shè)=n/q；對(duì)A的每一列的個(gè)水平為k(k=1,,q)的元素，用(k1)+1,(k1)+2,,(k1)+的一個(gè)隨機(jī)置換代替,產(chǎn)生的超拉丁方設(shè)計(jì)即為隨機(jī)化正交列，記為OH(n,ns)。,12,例6.2.,從正交列OA(8,4,2,3)出發(fā)構(gòu)造OH(8,84)。,13,C.正交超拉丁方設(shè)計(jì),正交超拉丁方設(shè)計(jì)：各因子列之間的相關(guān)系數(shù)為零。二階超拉丁方設(shè)計(jì)：拉丁方設(shè)計(jì)滿足下面條件：(i)設(shè)計(jì)中的任一列與其它列正交，即相關(guān)系數(shù)為零(ii

4、)設(shè)計(jì)中任一列對(duì)應(yīng)的元素平方列及任兩列的元素點(diǎn)乘列與設(shè)計(jì)的所有列正交,14,構(gòu)造LHD(2c+1+1,2c)的方法：,對(duì)c=1,設(shè)對(duì)c1,定義Sc和Tc如下如下得到的Lc即為一個(gè)LHD(2c+1+1,2c)：,15,6.3均勻設(shè)計(jì)在計(jì)算機(jī)試驗(yàn)中的應(yīng)用,考慮下面的偽三階模型其中，x=(x1,x2)，模型系數(shù)R1,S1,R2,S2,A,S為相互獨(dú)立的隨機(jī)選擇，且R1,S1U(0,0.5)，R2,S2U(0.5,1)，AU(0,0.05)，SU(0.04,1)，其中U(a,b)表示在區(qū)間a,b內(nèi)的均勻設(shè)計(jì)。,16,兩次隨機(jī)實(shí)現(xiàn),17,均勻設(shè)計(jì),18,擬合模型,三階回歸模型用最小二乘法估計(jì)系數(shù)可得R2

5、=0.941，s2=0.0002388，,(0.0420,0.2903,0.8484,0.5463,0.5270,1.2069,0.4550,0.2428,0.2610,0.5697).,19,預(yù)測MSE,檢驗(yàn)點(diǎn)：x1=0.02i,x2=0.02i,(i=0,1,,50)構(gòu)成的總共N=2601個(gè)網(wǎng)格點(diǎn)x1,,x2601預(yù)測誤差MSE=0.00022118,20,擬合模型2,逐步回歸法R2=0.919，s2=2.739104，擬合MSE=2.1972104預(yù)測MSE=2.3820104,21,試驗(yàn)點(diǎn)個(gè)數(shù)的影響,22,6.4對(duì)計(jì)算機(jī)試驗(yàn)諸設(shè)計(jì)的注記,完全隨機(jī)抽樣：不穩(wěn)定，效率不高超拉丁方抽樣：仍不夠穩(wěn)定均勻設(shè)計(jì)：穩(wěn)健,沒有一個(gè)設(shè)計(jì)方法永遠(yuǎn)是最好的，也沒有一個(gè)建模方法永遠(yuǎn)高人一籌。計(jì)算機(jī)試驗(yàn)的設(shè)計(jì)和建模是相當(dāng)復(fù)雜的問題，但均勻試驗(yàn)設(shè)計(jì)在大部分情形下均表現(xiàn)優(yōu)秀。,