簡易正交試驗設計方法

1、2.5 正交試驗設計方法正交試驗設計方法 一、試驗設計方法概述一、試驗設計方法概述 試驗設計是數(shù)理統(tǒng)計學的一個重要的分支。多數(shù)數(shù)理統(tǒng)計方法主要用于分析已經(jīng)得到的數(shù)據(jù)，而試驗設計卻是用于決定數(shù)據(jù)收集的方法。試驗設計方法主要討論如何合理地安排試驗以及試驗所得的數(shù)據(jù)如何分析等。 例例1 某化工廠想提高某化工產(chǎn)品的質(zhì)量和產(chǎn)量，對工藝中三個主要因素各按三個水平進行試驗（見表1）。試驗的目的是為提高合格產(chǎn)品的產(chǎn)量，尋求最適宜的操作條件。 對此實例該如何進行試驗方案的設計呢？水平 因素 溫度 壓力Pa 加堿量kg 符號 T p m 123 T1 (80 )T2(100)T3(120) p1(5.0) p2(

2、6.0) p3(7.0) m1(2.0)m2(2.5)m3(3.0) 表表1 因素水平因素水平 很容易想到的是全面搭配法方案（如圖1所示）：圖1 全面搭配法方案 此方案數(shù)據(jù)點分布的均勻性極好，因素和水平的搭配十分全面，唯一的缺點是實驗次數(shù)多達3327次（指數(shù)3代表3個因素，底數(shù)3代表每因素有3個水平）。因素、水平數(shù)愈多，則實驗次數(shù)就愈多，例如，做一個6因素3水平的試驗，就需36729次實驗，顯然難以做到。那么采用簡單比較法方案又如何呢？ 先固定T1和p1，只改變m，觀察因素m不同水平的影響，做了如圖2（1）所示的三次實驗，發(fā)現(xiàn) mm2時的實驗效果最好（好的用 表示），合格產(chǎn)品的產(chǎn)量最高，因此認

3、為在后面的實驗中因素m應取m2水平。圖2 簡單比較法方案 固定T1和m2，改變p的三次實驗如圖2（2）所示，發(fā)現(xiàn)pp時的實驗效果最好，因此認為因素p應取p水平。 固定p和m2，改變T 的三次實驗如圖2（3）所示，發(fā)現(xiàn)因素T 宜取T2水平。l 因此可以引出結(jié)論：為提高合格產(chǎn)品的產(chǎn)量，最適宜的操作條件為T2pm2。與全面搭配法方案相比，簡單比較法方案的優(yōu)點是實驗的次數(shù)少，只需做9次實驗。但必須指出，簡單比較法方案的試驗結(jié)果是不可靠的。需要尋找一種合適的試驗設計方法。l 在改變m值（或p值，或T值）的三次實驗中，說m2（或p或T2 ）水平最好是有條件的。在T T，p p時，m2 水平不是最好的可能性

4、是有的。l 在改變m的三次實驗中，固定T T，p p 應該說也是可以的，是隨意的，故在此方案中數(shù)據(jù)點的分布的均勻性是毫無保障的。l 用這種方法比較條件好壞時，只是對單個的試驗數(shù)據(jù)進行數(shù)值上的簡單比較，不能排除必然存在的試驗數(shù)據(jù)誤差的干擾。l試驗設計方法常用的術語定義如下：l 試驗指標：試驗指標：指作為試驗研究過程的因變量，常為試驗結(jié)果特征的量（如得率、純度等）。例1的試驗指標為合格產(chǎn)品的產(chǎn)量。l 因素：因素：指作試驗研究過程的自變量，常常是造成試驗指標按某種規(guī)律發(fā)生變化的那些原因。如例1的溫度、壓力、堿的用量。l 水平：水平：指試驗中因素所處的具體狀態(tài)或情況，又稱為等級。如例1的溫度有3個水平

5、。溫度用T表示，下標1、2、3表示因素的不同水平，分別記為T1、T2、T3。 二、正交試驗設計方法的優(yōu)點和特點二、正交試驗設計方法的優(yōu)點和特點l用正交表安排多因素試驗的方法，稱為正交試驗設計法。其優(yōu)點為： 完成試驗要求所需的實驗次數(shù)少。 數(shù)據(jù)點的分布很均勻。 可用相應的極差分析方法、方差分析方法、回歸分析方法等對試驗結(jié)果進行分析，引出許多有價值的結(jié)論。l運用正交試驗設計方法，不僅兼有上述兩個方案的優(yōu)點，而且實驗次數(shù)少，數(shù)據(jù)點分布均勻，結(jié)論的可靠性較好。l正交試驗設計方法是用正交表來安排試驗的。對于例1適用的正交表是L9（34），其試驗安排見表2。 表表2 試驗安排表試驗安排表試驗號 列號 12

6、34因素 溫度 壓力Pa 加堿量kg 符號 T p m 123456789 1（T1）1（T1）1（T1）2（T2）2（T2）2（T2）3（T3）3（T3）3（T3） 1（p1）2（p2）3（p3）1（p1）2（p2）3（p3）1（p1）2（p2）3（p3） 1（m1）2（m2）3（m3）2（m2）3（m3）1（m1）3（m3）1（m1）2（m2） 123312231 所有的正交表與L9（34）正交表一樣，都具有以下兩個特點：（1）在每一列中，各個不同的數(shù)字出現(xiàn)的次數(shù)相同。在表L9（34）中，每一列有三個水平，水平1、2、3都是各出現(xiàn)3次。（2）表中任意兩列并列在一起形成若干個數(shù)字對，不同數(shù)字

7、對出現(xiàn)的次數(shù)也都相同。在表L9（34）中，任意兩列并列在一起形成的數(shù)字對共有9個：（1,1），（1,2），（1,3），（2,1），（2,2），（2,3），（3,1），（3,2），（3,3），每一個數(shù)字對各出現(xiàn)一次。l這兩個特點稱為正交性。正是由于正交表具有上述特點，就保證了用正交表安排的試驗方案中因素水平是均衡搭配的，數(shù)據(jù)點的分布是均勻的。因素、水平數(shù)愈多，運用正交試驗設計方法，愈發(fā)能顯示出它的優(yōu)越性，如上述提到的6因素3水平試驗，用全面搭配方案需729次，若用正交表L27（313）來安排，則只需做27次試驗。 l在化工生產(chǎn)中， 因素之間常有交互作用。 如果上述的因素T的數(shù)值和水平發(fā)生變化時，

8、試驗指標隨因素p變化的規(guī)律也發(fā)生變化，或反過來，因素p的數(shù)值和水平發(fā)生變化時，試驗指標隨因素T變化的規(guī)律也發(fā)生變化。這種情況稱為因素T、p間有交互作用，記為Tp 。三、正交表三、正交表l使用正交設計方法進行試驗方案的設計，就必須用到正交表。正交表可查閱有關參考書。（如數(shù)學手冊） 1、各列水平數(shù)均相同的正交表各列水平數(shù)均相同的正交表l 各列水平數(shù)均相同的正交表，也稱單一水平正交表。這類正交表名稱的寫法舉例如下： 各列水平數(shù)均為2的常用正交表有：L4（23），L8（27），L12（211），L16（215），L20（219），L32（231）。 各列水平數(shù)均為3的常用正交表有：L9（34），L27

9、（313） 各列水平數(shù)均為4的常用正交表有：L16（45） 各列水平數(shù)均為5的常用正交表有：L25（56）L 9（3 4）正交表的列數(shù)每一列的水平數(shù)實驗的次數(shù)正交表的代號 2、混合水平正交表、混合水平正交表l各列水平數(shù)不相同的正交表，叫混合水平正交表，下面就是一個混合水平正交表名稱的寫法：l L 8（4124）常簡寫為L 8（424）。此混合水平正交表含有1 個4水平列，4個2水平列，共有145列。 L 8（4124）2水平列的列數(shù)為44水平列的列數(shù)為1實驗的次數(shù)正交表的代號3、選擇正交表的基本原則、選擇正交表的基本原則l一般都是先確定試驗的因素、水平和交互作用，然后選擇適用的L表。在確定因素

10、的水平數(shù)時，主要因素宜多安排幾個水平，次要因素可少安排幾個水平。 （1）先看水平數(shù)。若各因素全是2水平，就選用L(2)表；若各因素全是3水平，就選L(3)表。若各因素的水平數(shù)不相同，就選擇適用的混合水平表。 （2）每一個交互作用在正交表中應占一列或兩列。要看所選的正交表是否足夠大，能否容納得下所考慮的因素和交互作用。為了對試驗結(jié)果進行方差分析或回歸分析，還必須至少留一個空白列，作為“誤差”列，在極差分析中要作為“其他因素”列處理。 （3）要看試驗精度的要求。若要求高，則宜取實驗次數(shù)多的L表。 （4）若試驗費用很昂貴，或試驗的經(jīng)費很有限，或人力和時間都比較緊張，則不宜選實驗次數(shù)太多的L表。 （5

11、）按原來考慮的因素、水平和交互作用去選擇正交表，若無正好適用的正交表可選，簡便且可行的辦法是適當修改原定的水平數(shù)。 （6）對某因素或某交互作用的影響是否確實存在沒有把握的情況下，選擇L表時常為該選大表還是選小表而猶豫。若條件許可，應盡量選用大表，讓影響存在的可能性較大的因素和交互作用各占適當?shù)牧?。某因素或某交互作用的影響是否真的存在，留到方差分析進行顯著性檢驗時再做結(jié)論。這樣既可以減少試驗的工作量，又不致于漏掉重要的信息。 4、正交表的表頭設計、正交表的表頭設計l所謂表頭設計，就是確定試驗所考慮的因素和交互作用，在正交表中該放在哪一列的問題。 （1）有交互作用時，表頭設計必須嚴格地按規(guī)定辦事。

12、此處不討論，可以查閱有關書籍。 （2）若試驗不考慮交互作用，則表頭設計可以是任意的。如在例1中，對L 9（3 4）表頭設計，表3所列的各種方案都是可用的。但是正交表的構(gòu)造是組合數(shù)學問題，必須滿足（二、）中所述的特點。對試驗之初不考慮交互作用而選用較大的正交表，空列較多時，最好仍與有交互作用時一樣，按規(guī)定進行表頭設計。只不過將有交互作用的列先視為空列，待試驗結(jié)束后再加以判定。列號1234方案1234T空mPpT空mmpT空空mpT表表3 L 9（3 4）表頭設計方案）表頭設計方案 四、正交試驗的操作方法四、正交試驗的操作方法 （1）分區(qū)組。對于一批試驗，如果要使用幾臺不同的機器，或要使用幾種原料

13、來進行，為了防止機器或原料的不同而帶來誤差，從而干擾試驗的分析，可在開始做實驗之前，用L表中未排因素和交互作用的一個空白列來安排機器或原料。 與此類似，若試驗指標的檢驗需要幾個人（或幾臺機器）來做，為了消除不同人（或儀器）檢驗的水平不同給試驗分析帶來干擾，也可采用在L表中用一空白列來安排的辦法。這樣一種作法叫做分區(qū)組法。 （2）因素水平表排列順序的隨機化。如在例1中，每個因素的水平序號從小到大時，因素的數(shù)值總是按由小到大或由大到小的順序排列。按正交表做試驗時，所有的1水平要碰在一起，而這種極端的情況有時是不希望出現(xiàn)的，有時也沒有實際意義。因此在排列因素水平表時，最好不要簡單地按因素數(shù)值由小到大

14、或由大到小的順序排列。從理論上講，最好能使用一種叫做隨機化的方法。所謂隨機化就是采用抽簽或者查隨機數(shù)值表的辦法，來決定排列的順序。 （3）試驗進行的次序沒必要完全按照正交表上試驗號碼的順序。為減少試驗中由于先后實驗操作熟練的程度不勻帶來的誤差干擾，理論上推薦用抽簽的辦法來決定試驗的次序。 （4）在確定每一個實驗的實驗條件時，只需考慮所確定的幾個因素和分區(qū)組該如何取值，而不要（其實也無法）考慮交互作用列和誤差列怎么辦的問題。交互作用列和誤差列的取值問題由實驗本身的客觀規(guī)律來確定，它們對指標影響的大小在方差分析時給出。 （5）做實驗時，要力求嚴格控制實驗條件。這個問題在因素各水平下的數(shù)值差別不大時

15、更為重要。例如，例1中的因素（加堿量）m的三個水平：m1=2.0，m2=2.5，m3=3.0，在以m=m2=2.5為條件的某一個實驗中，就必須嚴格認真地讓m2=2.5。若因為粗心和不負責任，造成m2=2.2或造成m2=3.0，那就將使整個試驗失去正交試驗設計方法的特點，使極差和方差分析方法的應用喪失了必要的前提條件，因而得不到正確的試驗結(jié)果。 五、正交試驗結(jié)果分析方法五、正交試驗結(jié)果分析方法 正交試驗方法之所以能得到科技工作者的重視并在實踐中得到廣泛的應用，其原因不僅在于能使試驗的次數(shù)減少，而且能夠用相應的方法對試驗結(jié)果進行分析并引出許多有價值的結(jié)論。因此，用正交試驗法進行實驗，如果不對試驗結(jié)

16、果進行認真的分析，并引出應該引出的結(jié)論，那就失去用正交試驗法的意義和價值。1、極差分析方法、極差分析方法 下面以表4為例討論L4（23）正交試驗結(jié)果的極差分析方法。極差指的是各列中各水平對應的試驗指標平均值的最大值與最小值之差。從表4的計算結(jié)果可知，用極差法分析正交試驗結(jié)果可引出以下幾個結(jié)論： （1）在試驗范圍內(nèi)，各列對試驗指標的影響從大到小的排隊。某列的極差最大，表示該列的數(shù)值在試驗范圍內(nèi)變化時，使試驗指標數(shù)值的變化最大。所以各列對試驗指標的影響從大到小的排隊，就是各列極差D的數(shù)值從大到小的排隊。 （2）試驗指標隨各因素的變化趨勢。為了能更直觀地看到變化趨勢，常將計算結(jié)果繪制成圖。 （3）使

17、試驗指標最好的適宜的操作條件（適宜的因素水平搭配）。 （4）可對所得結(jié)論和進一步的研究方向進行討論 。 表表4 L4（23）正交試驗計算）正交試驗計算 注：注：j 第第 j 列列“1”水平所對應的試驗指標的數(shù)值之和；水平所對應的試驗指標的數(shù)值之和； j 第第 j 列列“2”水平所對應的試驗指標的數(shù)值之和；水平所對應的試驗指標的數(shù)值之和； kj 第第 j 列同一水平出現(xiàn)的次數(shù)。等于試驗的次數(shù)（列同一水平出現(xiàn)的次數(shù)。等于試驗的次數(shù)（n）除以第）除以第 j 列的水平數(shù)。列的水平數(shù)。 j / kj 第第 j 列列“1”水平所對應的試驗指標的平均值；水平所對應的試驗指標的平均值； j / kj 第第 j

18、 列列“2”水平所對應的試驗指標的平均值；水平所對應的試驗指標的平均值； Dj 第第 j 列的極差。等于第列的極差。等于第 j 列各水平對應的試驗指標平均值中的最大值減最列各水平對應的試驗指標平均值中的最大值減最 小值，小值， 即即 Djmaxj / kj ，j / kj ， -min j / kj ，j / kj ， 列 號 123試驗指標yi 試驗號 123n4 1122 1212 1221 y1y2y3y4 jjkjj / kjj / kj極差（Dj） 1y1y21y3y4k121/ k11/ k1max -min 2y1y32y2y4k222/ k22/ k2max -min 3y1y

19、43y2y3k323/ k33/ k3max -min 2、方差分析方法、方差分析方法（1）計算公式和項目計算公式和項目 l 試驗指標的加和值 ，試驗指標的平均 值 ，以第j列為例： j “1”水平所對應的試驗指標的數(shù)值之和 j “2”水平所對應的試驗指標的數(shù)值之和 kj 同一水平出現(xiàn)的次數(shù)。等于試驗的次數(shù)除以第j列的水平數(shù) j / kj “1”水平所對應的試驗指標的平均值 j / kj “2”水平所對應的試驗指標的平均值 以上7項的計算方法同極差法（見表4）。 niiy1niiyny11 偏差平方和偏差平方和 fj 自由度。自由度。fj 第第j 列的水平數(shù)列的水平數(shù)1。 Vj 方差。方差。V

20、j Sj fj 。 Ve 誤差列的方差。誤差列的方差。Ve Se fe 。式中，。式中，e為正交表的誤差列。為正交表的誤差列。 Fj 方差之比方差之比 Fj Vj Ve 。 查查F分布數(shù)值表（分布數(shù)值表（F分布數(shù)值表請查閱有關參考書）做顯著性檢驗。分布數(shù)值表請查閱有關參考書）做顯著性檢驗。 總的偏差平方和總的偏差平方和 總的偏差平方和等于各列的偏差平方和之和。即總的偏差平方和等于各列的偏差平方和之和。即 式中，式中，m為正交表的列數(shù)。為正交表的列數(shù)。 若誤差列由若誤差列由5個單列組成，則誤差列的偏差平方和個單列組成，則誤差列的偏差平方和Se等于等于5個單列的個單列的偏差平方和之和，即：偏差平方

21、和之和，即：Se Se1 Se2 Se3 Se4 Se5 ；也可用；也可用Se S總總 S 來計算，其中來計算，其中S 為安排有因素或交互作用的各列的偏差平為安排有因素或交互作用的各列的偏差平方和之和。方和之和。 222IIIIIIykkykkykkSjjjjjjjjjjmjjSS1總niiyyS12總（2）可引出的結(jié)論）可引出的結(jié)論 與極差法相比，方差分析方法可以多引出一個結(jié)論：各列對試驗指標的影響是否顯著，在什么水平上顯著。在數(shù)理統(tǒng)計上，這是一個很重要的問題。顯著性檢驗強調(diào)試驗在分析每列對指標影響中所起的作用。如果某列對指標影響不顯著，那么，討論試驗指標隨它的變化趨勢是毫無意義的。因為在某

22、列對指標的影響不顯著時，即使從表中的數(shù)據(jù)可以看出該列水平變化時，對應的試驗指標的數(shù)值與在以某種“規(guī)律”發(fā)生變化，但那很可能是由于實驗誤差所致，將它作為客觀規(guī)律是不可靠的。有了各列的顯著性檢驗之后，最后應將影響不顯著的交互作用列與原來的“誤差列”合并起來。組成新的“誤差列”，重新檢驗各列的顯著性。六、正交試驗方法在化工原理實驗中的六、正交試驗方法在化工原理實驗中的 應用舉例應用舉例 例例2 為提高真空吸濾裝置的生產(chǎn)能力，請用正交試驗方法確定恒壓過濾的最佳操作條件。 影響實驗的主要因素和水平見表5（a）。表中p為過濾壓強差；T為漿液溫度；w為漿液質(zhì)量分數(shù)；M為過濾介質(zhì)（材質(zhì)屬多孔陶瓷）。 解：解：

23、（1）試驗指標的確定：恒壓過濾常數(shù)K（m2/s） （2）選正交表：根據(jù)表5（a）的因素和水平，可選用L 8（424）表。 因素 壓強差kPa 溫度 質(zhì)量分數(shù) 過濾介質(zhì) 符號 p T w M 水平 1234 2.943.924.905.88 （室溫）18（室溫15）33 稀（約5%）濃（約10%） G2*G3* * G2 、G3為過濾漏斗的型號。過濾介質(zhì)孔徑：G2 為3050m、G3為1630m。 表表5（a） 過濾實驗因素和水平過濾實驗因素和水平 （3）制定實驗方案：按選定的正交表，應完成8次實驗。實驗方案見表5（b）。 （4）實驗結(jié)果：將所計算出的恒壓過濾常數(shù)K（m2/s）列于表5（b）。

24、列號 j1 23456因素 p T w M e K（m2/s） 試驗號 水 平 12345678 11223344 12121212 12122121 12211221 12212112 4.011042.931045.211045.551044.831041.021035.111041.10103 表表5（b）正交試驗的試驗方案和實驗結(jié)果）正交試驗的試驗方案和實驗結(jié)果 （5）指標K的極差分析和方差分析： 分析結(jié)果見表5（c）。以第2列為例說明計算過程： 2 4.011045.211044.831045.11104 1.92103 2 2.931045.551041.021031.101032

25、.97103 k24 2/ k21.92103 44.79104 2/ k22.9710347.42104 D27.42104 4.791042.63104 K4.88103 6.11104 S2k2（2/ k2 ）2k2（2/ k2 ）2 4（4.791046.11104 ）2 4（7.421046.11104 ）2 1.38107 KKK f2第二列的水平數(shù)第二列的水平數(shù)1211 V2S2f21.3810711.38107 SeS5k5（5/ k5 ）2k5（5/ k5 ）2 4（6.221046.11104 ）2 4（5.991046.11104 ）2 1.06109 fef51 VeS

26、efe1.0610911.06109 F2 =V2Ve1.381071.06109130.2 查查F 分布數(shù)值表分布數(shù)值表可知：可知： F（0.01，f11，f21）4052 F2 F（0.05，f11，f21）161.4 F2 F（0.10，f11，f21）39.9 F2 F（0.25，f11，f21）5.83 F2 （其中：（其中：f1 為分子的自由度，為分子的自由度，f2 分母的自由度）分母的自由度） 所以第二列對試驗指標的影響在所以第二列對試驗指標的影響在0.10水平上顯著。水平上顯著。 練習：對剩余列作級差和方差分析練習：對剩余列作級差和方差分析KK 表表5（c） K的極差分析和方差

27、分析的極差分析和方差分析列號 j1 23456因素 p T w M e K（m2/s） 項目 jjjjkjj / kjj / kjj / kjj / kjDjSjfjVjFjF0.01F0.05F0.10F0.25顯著性 6.941041.081031.501031.6110323.471045.381047.521048.061034.591042.6510738.8410883.45403215.753.68.202*(0.10) 1.921032.9710344.791047.421042.631041.3810711.38107130.24052161.439.95.832*(0.10) 3.041031.8410347.611044.611043.001041.8010711.80107169.84052161.4