正交試驗設計與均勻試驗設計(1)

1、.1 正交試驗設計正交試驗設計 .2 一一.正交試驗設計正交試驗設計 1.1正交試驗設計的基本概念正交試驗設計的基本概念 正交試驗設計正交試驗設計基本概念基本概念 : :就是安排就是安排多因素試驗多因素試驗、尋求、尋求 最優(yōu)水平組合最優(yōu)水平組合的一種高效率試驗設計方法。的一種高效率試驗設計方法。 適用于適用于3 3個或個或3 3個以上的試驗因素。個以上的試驗因素。 正交試驗設計基本特點正交試驗設計基本特點：用部分試驗來代替全面試：用部分試驗來代替全面試 驗，通過對部分試驗結果的分析了解全面試驗。驗，通過對部分試驗結果的分析了解全面試驗。 .3 1.2正交試驗設計的基本原理正交試驗設計的基本原理

2、 eg1: eg1: 考察考察增稠劑用量增稠劑用量、pHpH值值和和殺菌溫度殺菌溫度對豆奶穩(wěn)定對豆奶穩(wěn)定 性的影響。每個因素設置性的影響。每個因素設置3 3個水平進行試驗。個水平進行試驗。 分析分析: : 增稠劑用量定為增稠劑用量定為A A因素因素，設，設5%5%、10%10%、15%( A15%( A1 1 A A2 2、A A3 3) 3) 3水平；水平； pHpH值定為值定為B B因素因素，設，設4.04.0、4.54.5、5.05.0（B B1 1、B B2 2、 B B3 3）3 3個水平；個水平； 殺菌溫度定為殺菌溫度定為C C因素因素，設，設5050、60 60 、70 70 （

3、C C1 1、C C2 2、C C3 3 ）3 3個水平。個水平。 .4 如上例中，如上例中，3個因素的選優(yōu)區(qū)可以用一個立方體個因素的選優(yōu)區(qū)可以用一個立方體 表示（圖表示（圖1），），3個因素各取個因素各取 3個水平，把立方體劃分個水平，把立方體劃分 成成27個格點，反映在個格點，反映在 圖圖1上就是立方體內的上就是立方體內的27個個“.”。 若若27個網(wǎng)格點都試驗，就是全面試驗，其試驗方案如個網(wǎng)格點都試驗，就是全面試驗，其試驗方案如 表表1所示。所示。 .5 圖圖13因素因素3水平全面試驗的立體圖水平全面試驗的立體圖 .6 全面試驗全面試驗各因素的水平之間全部可能組合有各因素的水平之間全部可

4、能組合有27種。種。 .7 正交設計正交設計就是從全面試驗水平組合中挑選出有代就是從全面試驗水平組合中挑選出有代 表性的部分試驗水平組合來進行試驗。圖表性的部分試驗水平組合來進行試驗。圖2 2中標有試驗中標有試驗 號的九個號的九個“( ()”)”，就是利用，就是利用正交表正交表L L9 9(3(34 4) )從從2727個試個試 驗點中挑選出來的驗點中挑選出來的9 9個試驗點。即：個試驗點。即： (1) A(1) A1 1B B1 1C C1 1 (2) A (2) A1 1B B2 2C C2 2 (3) A (3) A1 1B B3 3C C3 3 (4) A (4) A2 2B B1 1

5、C C2 2 (5) A (5) A2 2B B2 2C C3 3 (6) A (6) A2 2B B3 3C C1 1 (7) A (7) A3 3B B1 1C C3 3 (8) A (8) A3 3B B2 2C C1 1 (9) A (9) A3 3B B3 3C C2 2 .8 .9 上述選擇，保證了上述選擇，保證了A A因素的每個水平與因素的每個水平與B B因素、因素、C C因素的各因素的各 個水平在試驗中各搭配一次個水平在試驗中各搭配一次 。對于。對于A A、B B、C 3C 3個因素來個因素來 說，是在說，是在2727個全面試驗點中選擇個全面試驗點中選擇9 9個試驗點，僅是全面

6、試個試驗點，僅是全面試 驗的三分之一。驗的三分之一。 從圖從圖1 1中可以看到，中可以看到，9 9個試驗點在選優(yōu)區(qū)中分布是均衡的，個試驗點在選優(yōu)區(qū)中分布是均衡的， 在立方體的每個平面上，都恰有在立方體的每個平面上，都恰有3 3個試驗點；在立方體的個試驗點；在立方體的 每條線上也恰有一個試驗點。每條線上也恰有一個試驗點。 9 9個試驗點均衡地分布于整個立方體內，有很強的代表性，個試驗點均衡地分布于整個立方體內，有很強的代表性， 能夠比較全面地反映選優(yōu)區(qū)內的基本情況。能夠比較全面地反映選優(yōu)區(qū)內的基本情況。 .10 1.3 1.3 正交表及其基本性質正交表及其基本性質 1.3.1 1.3.1 正交表

7、正交表 (Latin(Latin方方) ) 表表2 2是一張正交表，記號為是一張正交表，記號為L L8 8(2(27 7) )，其中，其中“L L” 代表正交表；代表正交表；L L右下角的數(shù)字右下角的數(shù)字“8”8”表示有表示有8 8行，用這張行，用這張 正交表安排試驗包含正交表安排試驗包含8 8個處理個處理( (水平組合水平組合) )；括號內的底；括號內的底 數(shù)數(shù)“2” 2” 表示因素的水平數(shù)，括號內表示因素的水平數(shù)，括號內2 2的指數(shù)的指數(shù)“7”7”表表 示有示有7 7列列 ，用這張正交表最多可以安排，用這張正交表最多可以安排7 7個個2 2水平因素。水平因素。 .11 .12 L為正交表符

8、號（為正交表符號（Latin)，n為試驗次數(shù)，即正交表為試驗次數(shù)，即正交表 行數(shù)，行數(shù)，t為因素的水平數(shù)，即為因素的水平數(shù)，即1列中出現(xiàn)的不同數(shù)字列中出現(xiàn)的不同數(shù)字 的個數(shù)，的個數(shù)，q為最多能安排的因素數(shù)，即正交表的列數(shù)。為最多能安排的因素數(shù)，即正交表的列數(shù)。 .13 常用的正交表已由數(shù)學工作者制定出來，供常用的正交表已由數(shù)學工作者制定出來，供 進行正交設計時選用。進行正交設計時選用。2水平正交表除水平正交表除L8(27)外，外， 還有還有L4(23)、L16(215)等；等；3水平正交表有水平正交表有L9(34)、 L27(313) 等等 .14 解答：解答： 3因素，因素，3水平水平 查附

9、表查附表17的（的（4）和（）和（5），有兩個正交表），有兩個正交表 L9（34）和）和L27（313），最終選擇），最終選擇L9（34） .15 1.3.2 正交表的基本性質正交表的基本性質 1.3.2.1 正交性正交性 （1）任一列中，各水平都出現(xiàn)，且出現(xiàn)的次數(shù)相等）任一列中，各水平都出現(xiàn)，且出現(xiàn)的次數(shù)相等 例如例如L8(27)中不同數(shù)字只有中不同數(shù)字只有1和和2，它們各出現(xiàn)，它們各出現(xiàn)4次；次； L9(34)中不同數(shù)字有中不同數(shù)字有1、2和和3，它們各出現(xiàn)，它們各出現(xiàn)3次次。 （2 2）任兩列之間各種不同水平的所有可能組合都出現(xiàn)，且對出）任兩列之間各種不同水平的所有可能組合都出現(xiàn)，且對出

10、 現(xiàn)的次數(shù)相等?，F(xiàn)的次數(shù)相等。 例如例如 L L8 8(2(27 7) )中中(1, 1), (1, 2), (2, 1), (2, 2)(1, 1), (1, 2), (2, 1), (2, 2)各出現(xiàn)兩各出現(xiàn)兩 次；次；L L9 9(3(34 4) ) 中中 (1, 1), (1, 2), (1, 3), (2, 1), (2, 2), (1, 1), (1, 2), (1, 3), (2, 1), (2, 2), (2, 3), (3, 1), (3, 2), (3, 3)(2, 3), (3, 1), (3, 2), (3, 3)各出現(xiàn)各出現(xiàn)1 1次。即每個因素的次。即每個因素的 一個

11、水平與另一因素的各個水平所有可能組合次數(shù)相等，表明一個水平與另一因素的各個水平所有可能組合次數(shù)相等，表明 任意兩列各個數(shù)字之間的搭配是均勻的。任意兩列各個數(shù)字之間的搭配是均勻的。 .16 .17 .18 根據(jù)以上特性，我們用正交表安排的試根據(jù)以上特性，我們用正交表安排的試 驗，具有驗，具有均衡分散均衡分散和和整齊可比整齊可比的特點。的特點。 .19 1.4 1.4 正交表的類別正交表的類別 1、等水平正交表等水平正交表 各列水平數(shù)相同的正交表各列水平數(shù)相同的正交表 稱為等水平正交表。如稱為等水平正交表。如L4(23)、L8(27)、L12(211)等等 各列中的水平為各列中的水平為2，稱為，稱

12、為2水平正交表；水平正交表；L9(34)、 L27(313)等各列水平為等各列水平為3，稱為，稱為3水平正交表。水平正交表。 2、混合水平正交表混合水平正交表 各列水平數(shù)不完全相同各列水平數(shù)不完全相同 的正交表稱為混合水平正交表。如的正交表稱為混合水平正交表。如L8(424)表中有表中有 一列的水平數(shù)為一列的水平數(shù)為4，有，有4列水平數(shù)為列水平數(shù)為2。也就是說該表。也就是說該表 可以安排一個可以安排一個4水平因素和水平因素和4個個2水平因素。再如水平因素。再如 L16(4423)，L16(4212)等都混合水平正交表。等都混合水平正交表。 .20 .21 .22 試驗結果方差分析試驗結果方差分

13、析 列方差分析表，列方差分析表， 進行進行F F 檢驗檢驗 計算各列偏差平方和、計算各列偏差平方和、 自由度自由度 分析檢驗結果，分析檢驗結果， 寫出結論寫出結論 .23 et2：為提高山楂原料的利用率，研究酶法液化工藝：為提高山楂原料的利用率，研究酶法液化工藝 制造山楂原汁，擬通過正交試驗來尋找酶法液制造山楂原汁，擬通過正交試驗來尋找酶法液 化的最佳工藝條件?；淖罴压に嚄l件。 .24 選因素選因素:根據(jù)專業(yè)知識、以往的研究結論和經(jīng)驗，從影響試驗指根據(jù)專業(yè)知識、以往的研究結論和經(jīng)驗，從影響試驗指 標的諸多因素中，通過因果分析篩選出需要考察的試驗因素。一標的諸多因素中，通過因果分析篩選出需要考

14、察的試驗因素。一 般確定試驗因素時，應以對試驗指標影響大的因素、尚未考察過般確定試驗因素時，應以對試驗指標影響大的因素、尚未考察過 的因素、尚未完全掌握其規(guī)律的因素為先。的因素、尚未完全掌握其規(guī)律的因素為先。 定水平定水平:試驗因素選定后，根據(jù)所掌握的信息資料和相關知識，試驗因素選定后，根據(jù)所掌握的信息資料和相關知識， 確定每個因素的水平，一般以確定每個因素的水平，一般以2-42-4個水平為宜。對主要考察的試個水平為宜。對主要考察的試 驗因素，可以多取水平，但不宜過多（驗因素，可以多取水平，但不宜過多（6 6），否則試驗次數(shù)驟），否則試驗次數(shù)驟 增。因素的水平間距，應根據(jù)專業(yè)知識和已有的資料，

15、盡可能把增。因素的水平間距，應根據(jù)專業(yè)知識和已有的資料，盡可能把 水平值取在理想?yún)^(qū)域。水平值取在理想?yún)^(qū)域。 （2） 選因素、定水平，列因素水平表選因素、定水平，列因素水平表 .25 水平 試驗因素 加水量 （mL/100g） A 加酶量 （mL/100g） B 酶解溫度 （） C 酶解時間 （h） D 1101201.5 2504352.5 3907503.5 .26 一般情況下，試驗因素的水平數(shù)應等于正交一般情況下，試驗因素的水平數(shù)應等于正交 表中的水平數(shù)；因素個數(shù)（包括交互作用）應不表中的水平數(shù)；因素個數(shù)（包括交互作用）應不 大于正交表的列數(shù)。大于正交表的列數(shù)。 此例有此例有4 4個個3

16、3水平因素，可以選用水平因素，可以選用L L9 9(3(3 4 4) )或 或 L L27 27(3 (313 13) ) ；因本試驗僅考察四個因素對液化率的；因本試驗僅考察四個因素對液化率的 影響效果，不考察因素間的交互作用，故宜選用影響效果，不考察因素間的交互作用，故宜選用 L L9 9（3 34 4）正交表。）正交表。 .27 把正交表中安排各因素的列中的每個水把正交表中安排各因素的列中的每個水 平數(shù)字換成該因素的實際水平值，便形成了正平數(shù)字換成該因素的實際水平值，便形成了正 交試驗方案（表交試驗方案（表5 5）。）。 .28 試驗號 因 素 ABCD 11111 21222 31333

17、 42123 52231 62312 73132 83213 93321 試驗結果 （液化率 %） 0 17 24 12 47 28 1 18 42 .29 2.2 2.2 試驗結果分析試驗結果分析- - 分清各因素及其交互作用的主次順序，分清哪個是主要因素，分清各因素及其交互作用的主次順序，分清哪個是主要因素， 哪個是次要因素；哪個是次要因素； 判斷因素對試驗指標影響的顯著程度；判斷因素對試驗指標影響的顯著程度； 找出試驗因素的優(yōu)水平和試驗范圍內的最優(yōu)組合，即試驗因找出試驗因素的優(yōu)水平和試驗范圍內的最優(yōu)組合，即試驗因 素各取什么水平時，試驗指標最好；素各取什么水平時，試驗指標最好； 分析因素

18、與試驗指標之間的關系，即當因素變化時，試驗指分析因素與試驗指標之間的關系，即當因素變化時，試驗指 標是如何變化的。找出指標隨因素變化的規(guī)律和趨勢，為進標是如何變化的。找出指標隨因素變化的規(guī)律和趨勢，為進 一步試驗指明方向；一步試驗指明方向； 了解各因素之間的交互作用情況；了解各因素之間的交互作用情況； 估計試驗誤差的大小。估計試驗誤差的大小。 .30 .31 Kjm為第為第j列因素列因素m水平所對應的試驗指標和，水平所對應的試驗指標和，kjm為為 Kjm平均值。由平均值。由kjm大小可以判斷第大小可以判斷第j列因素優(yōu)水平列因素優(yōu)水平 和優(yōu)組合。和優(yōu)組合。 Rj為第為第j列因素的列因素的kjm的

19、極差，反映了第的極差，反映了第j列因素水平列因素水平 波動時，試驗指標的變動幅度。波動時，試驗指標的變動幅度。Rj越大，說明該越大，說明該 因素對試驗指標的影響越大。根據(jù)因素對試驗指標的影響越大。根據(jù)Rj大小，可以大小，可以 判斷因素的主次順序。判斷因素的主次順序。 .32 試驗號 因 素 ABCD 11111 21222 31333 42123 52231 62312 73132 83213 93321 試驗結果 （液化率 %） 0 17 24 12 47 28 1 18 42 .33 .34 .35 .36 .37 極差分析法簡單明了，計算工作量少便于推廣普及，極差分析法簡單明了，計算工作

20、量少便于推廣普及， 但不能區(qū)分因素各水平間對應的試驗結果的差異究竟但不能區(qū)分因素各水平間對應的試驗結果的差異究竟 是由于因素水平不同引起的，還是由于試驗誤差引起是由于因素水平不同引起的，還是由于試驗誤差引起 的，無法估計試驗誤差的大小。此外，各因素對試驗的，無法估計試驗誤差的大小。此外，各因素對試驗 結果的影響大小無法給以精確的數(shù)量估計，不能提出結果的影響大小無法給以精確的數(shù)量估計，不能提出 一個標準來判斷所考察因素作用是否顯著。為了彌補一個標準來判斷所考察因素作用是否顯著。為了彌補 極差分析的缺陷，可采用極差分析的缺陷，可采用方差分析方差分析。 2.3 正交試驗結果的方差分析正交試驗結果的方

21、差分析 .38 方差分析：將數(shù)據(jù)的總變異分解成方差分析：將數(shù)據(jù)的總變異分解成因素引起的因素引起的 變異變異和和誤差引起的變異誤差引起的變異兩部分，構造兩部分，構造F統(tǒng)計量統(tǒng)計量，作，作F 檢驗，即可判斷因素作用是否檢驗，即可判斷因素作用是否顯著顯著。 .39 空列（誤差）因素 SSSSSST 空列（誤列（因素T dfdfdf 誤差 誤差 誤差 因素 因素 因素 df SS MS df SS MS ， .40 誤差 因素 因素 MS MS F .41 誤差 因素 因素 MS MS F .42 處理處理 號 號 第第1 1列（列（A A） 第第2 2列 列 第第3 3列列 第第4 4列列 試驗結果

22、試驗結果yiyi 1 11 11 11 11 1y1y1 2 21 12 22 22 2y2y2 3 31 13 33 33 3y3y3 4 42 21 12 23 3y4y4 5 52 22 23 31 1y5y5 6 62 23 31 12 2y6y6 7 73 31 13 32 2y7y7 8 83 32 21 13 3y8y8 9 93 33 32 21 1y9y9 因素因素重復重復1 1重復重復2 2重復重復3 3 A A1 1y1y1y2y2y3y3 A A2 2y4y4y5y5y6y6 A A3 3y7y7y8y8y9y9 和 y1+y2+y3K K1 1 y4+y5+y6K K

23、2 2 y7+y8+y9K K3 3 .43 表頭設計表頭設計A AB B試驗數(shù)據(jù)試驗數(shù)據(jù) 列號列號1 12 2k kx xi ix xi i2 2 試驗號試驗號 1 11 1x x1 1x x1 12 2 2 21 1x x2 2x x2 22 2 n nm mx xn nx xn n2 2 K K1j 1jK K1111K K1212K K1k1k K K2j 2jK K2121K K2222K K2k2k K Kmj mjK Km1m1K Km2m2K Kmkmk K K1j 1j2 2K K11112 2K K12122 2K K1k1k2 2 K K2j 2j2 2K K21212

24、2K K22222 2K K2k2k2 2 K Kmj mj2 2K Km1m12 2K Km2m22 2K Kmkmk2 2 SSSSj jSSSS1 1SSSS2 2SSSSk k CTQSS r 1 Q CTQSS xQ n T CTxT jj i j TT n i i T 2 n i i m 1 2 ij 1 2 1 K .44 n x xSS n 1i 2 i n 1i 2 iT ）（ ），（ ）（ k.21j n x K r 1 SS n 1i 2 i m 1i 2 ijj ），（）（k.21j K-K n 1 SS 2 2j1jj .45 n-1dfT 為因素水平個數(shù)，m j 1

25、mdf .46 均勻試驗設計均勻試驗設計 .47 正交設計正交設計: 可使試驗點可使試驗點“均勻分散、整齊可比均勻分散、整齊可比”。 為保證為保證“整齊可比性整齊可比性”，使試驗設計的均勻性受到，使試驗設計的均勻性受到 一一 定限制，使試驗點的代表性還不夠強，試驗次數(shù)不能充定限制，使試驗點的代表性還不夠強，試驗次數(shù)不能充 分地少。分地少。 均勻設計均勻設計:可以用較少的試驗次數(shù)，安排多因素、多水平的析因可以用較少的試驗次數(shù)，安排多因素、多水平的析因 試試驗，是在驗，是在均均勻性勻性的度量下的度量下最好最好的析因試驗設計方法。的析因試驗設計方法。 可以使試驗點在試驗范圍內充分地可以使試驗點在試驗

26、范圍內充分地均勻分散均勻分散，不僅可大大，不僅可大大 減少試驗點，而且仍能得到反映試驗體系主要特征的試驗結減少試驗點，而且仍能得到反映試驗體系主要特征的試驗結 果。果。 .48 .49 下面通過下面通過制藥工業(yè)中的一個實例制藥工業(yè)中的一個實例來說來說 明均勻試驗設計方法。明均勻試驗設計方法。 例例1.11.1 ：阿魏酸的制備 根據(jù)試驗目的，確定以根據(jù)試驗目的，確定以阿魏酸產(chǎn)量阿魏酸產(chǎn)量作為試驗作為試驗 指標指標Y Y。 阿魏酸是某些藥品的主要成分，在制備阿魏酸是某些藥品的主要成分，在制備 過程中，我們想提高阿魏酸產(chǎn)量。過程中，我們想提高阿魏酸產(chǎn)量。 .50 全面交叉試驗要全面交叉試驗要N=73

27、=343次，太多了。次，太多了。 建議使用均勻設計。查閱均勻設計表。建議使用均勻設計。查閱均勻設計表。 如何安排試驗如何安排試驗? 經(jīng)過資料查閱，分析研究后經(jīng)過資料查閱，分析研究后 1.選出影響阿魏酸產(chǎn)量的試驗因素選出影響阿魏酸產(chǎn)量的試驗因素 2.確定試驗因素為：確定試驗因素為： 原料配比：原料配比：1.0-3.4 吡啶總量：吡啶總量：10-28 反應時間：反應時間：0.5-3.5 3.確定每個因素相應的水平數(shù)為確定每個因素相應的水平數(shù)為7。 .51 第第1步：步： 列出試驗因素水平表列出試驗因素水平表 表表 1 試驗因素水平表試驗因素水平表 .52 第第2步步: : 選擇相應的均勻設計表選擇

28、相應的均勻設計表 均勻設計表格式見下，其含義為：均勻設計表格式見下，其含義為：uniformity(U) Un(qs) 均勻設計 試驗總次數(shù) 因素水平數(shù) 因素個數(shù) .53 例如: )7( 4 7 U )9( 4 9 U 表表 2: 表表 3: .54 每個均勻設計表都有一個每個均勻設計表都有一個使用表使用表，它將建議我們如何選，它將建議我們如何選 擇適當?shù)牧邪才旁囼炓蛩?，進行試驗設計，這樣可以減少擇適當?shù)牧邪才旁囼炓蛩?，進行試驗設計，這樣可以減少“ 試驗偏差試驗偏差”。其中。其中偏差偏差為均勻性的度量值，數(shù)值小的設為均勻性的度量值，數(shù)值小的設 計表示均勻性好。例如計表示均勻性好。例如 U7 (

29、74)的使用表為：的使用表為： 因素數(shù)列號偏差 21, 30.2398 31, 2, 30.3721 41, 2, 3, 40.4760 No.1234 11236 22465 33624 44153 55312 66541 77777 No.123 1123 2246 3362 4415 5531 6654 7777 ) 4 7( 7 U表表 1.1.4:表表1.1.2: .55 第第3步步: 應用選擇的 UD-表安排試驗，設計試驗方案 No.123 1123 2246 3362 4415 5531 6654 7777 1. 將 x1, x2和 x3放入均勻設 計表的1,和3列； x1 x2 x3 2用x1的個水平（值）替 代第一列的1到 7； 1.0 1.4 1.8 2.2 2.6 3.0 3.4 3. 對第二列，第三列做同樣 的替代； 13 1.5 19 3.0 25 1.0 10 2.5 16 0.5 22 2.0 28 3.5 4. 按設計的方案進行試驗，得 到個結果，將其放入最后一 列。 y 0.330 0.366 0.294 0.476 0.209 0.451 0.482 表表 1.1.5: 試驗方案試驗方案