L9(34)是什么意思呢?
字母L表示正交表;數字9表示這張表共有9行,說明用這張表來安排試驗要做9次試驗;數字4表示這張表共有4列,說明用這張表最多可安排4個因素;數字3表示在表中主體部分只出現1,2,3三個數字,它們分別代表因素的3個水平,說明各因素都是3個水平的.
一般的正交表記為Ln(mk),n是表的行數,也就是要安排的試驗次數;k是表中列數,表示因素的個數;m是各因素的水平數。
常見的正交表如下:
(1)2水平正交表——L4(23),L8(27),L12(211),L16(215)等。
這幾張表中的數字2表示各因素都是2水平的;試驗要做的次數分別為4,8,12,16;最多可安排的因素分別為3,7,11,15。
(2)3水平的正交表——L9(34),L27(313)。
這兩張表中的數字3表示各因素都是3水平的,要做的試驗次數分別為9,27;最多可安排的因素分別為4,13。
(3)4水平的正交表——L14(45)。
(4)5水平的正交表——L25(56)。
表2
因素 水平 |
A 焦比 |
|
1:16 1:18 1:14 |
解: 正交表L9(34)
列號 |
|
試驗號 |
|
1 |
|
2 |
|
3 |
|
4 |
|
5 |
|
6 |
|
7 |
|
8 |
|
9 |
|
表3
試驗編號 |
鐵水溫度/℃ |
為了便於分析,將結果與正交表合編,以便利於計算。由於鐵水溫度數值較大,可把每一個鐵水溫度的值減去1350,得到9個較小的數,這樣使計算簡單(表4)。
表4
|
A |
鐵水溫度/℃ |
鐵水溫度值減去1350 |
1 2 3 4 5 6 7 8 9 |
1 1 1 2 2 2 3 3 3 |
1365 1395 1385 1390 1395 1380 1390 1390 1410 |
15 45 35 40 45 30 40 40 60 |
K1 K2 K3 |
|
|
|
K1=K1/3 K2=K2/3 K3=K3/3 |
|
||
極差 |
|
||
優方案 |
|
K1這一行的3個數,分別是因素A,B,C的第1水平所在的試驗中對應的鐵水溫度(減去1350以后)之和。比如對因素A(第1列),它的第l水平安排在第1,2,3號試驗中,對應的鐵水溫度值(減去1350以后)分別為15,45,35,其和為95,記在K1這一行的第1列中。
對於因素B(第2列),它的第1水平安排在第1,4,7號試驗中,對應的鐵水溫度值(減去1350以后)分別為15,40,40,其和為95,記在K1這一行的第2列中。
對於因素C(第3列),它的第l水平安排在第1,6,8號試驗中,對應的鐵水溫度值(減去1350以后)分別為15,30,40,其和為85,記在K1這一行的第3列中。
類似地,K2這一行的3個數,分別是因素A,B,C的第2水平所在的試驗中對應的鐵水溫度(減去1350以后)之和。K3這一行的3個數,分別是因素A,B,C的第3水平所在的試驗中對應的鐵水溫度(減去1350以后)之和.
同一列中,K1,k2,k3這3個數中的最大者減去最小者所得的差叫做極差。
一般地說,各列的極差是不同的,這說明各因素的水平改變時對試驗指標的影響是不同的。極差越大,說明這個因素的水平改變時對試驗指標的影響越大。極差最大的那一列,則那個因素的水平改變時對試驗指標的影響就最大,那個因素就是我們要考慮的主要因素。
這里算出3列的極差分別為15.0,11.6,20.0,顯然第3列即因素C的極差20.0最大。這說明因素C的水平改變時對試驗指標的影響最大,因此因素C是我們要考慮的主要因素。它的3個水平所對應的鐵水溫度(減去1350以后)平均值分別為28.3,48.3,40.O,第2水平所對應的數值48.3最大,所以取它的第2水平最好。
第1列即因素A的極差為15.0,僅次於因素C,它的3個水平所對應的數值分別為31.7,38.3,46.7,第3水平所對應的數值46.7最大,所以取它的第3水平最好。
第2列即因素B的極差為11.6,是3個因素中極差最小的,說明它的水平改變時對試驗指標的影響最小,它的3個水平所對應的數值分別為31.7,43.3,41.7,第2水平所對應的數值43.3最大,所以取它的第2水平最好。
從以上分析可以得出結論:各因素對試驗指標(鐵水溫度)的影響按大小次序來說應當是C(底焦高度)A(焦比)B(風壓);最好的方案應當是C2A3B2,即
可以看出,這里分析出來的最好方案在已經做過的9次試驗中沒有出現,與它比較接近的是第9號試驗.在第9號試驗中只有風壓B不是處在最好水平,而且風壓對鐵水溫度的影響是3個因素中最小的。從實際做出的結果看出,第9號試驗中的鐵水溫度是1410℃,是9次試驗中最高的,這也說明我們找出的最好方案是符合實際的。
為了最終確定上面找出的試驗方案C2A3B2是否為最好方案,可以按這個方案再試驗一次,看是否會得出比第9號試驗更好的結果。若比第9號試驗的效果好,就確定上述方案為最好方案,若不比第9號試驗的效果好,可以取第9號試驗為最好方案。如果出現后一種情況,說明我們的理論分析與實踐有一些差距,最終還是要接受實踐的檢驗。
現將利用正交表安排試驗並分析試驗結果的步驟歸納如下:
(1)明確試驗目的,確定要考核的試驗指標。
(2)根據試驗目的,確定要考察的因素和各因素的水平。要通過對實際問題的具體分析選出主要因素,略去次要因素,這樣可使因素個數少些。如果對問題不太了解,因素個數可適當地多取一些,經過對試驗結果的初步分析,再選出主要因素。因素被確定后,隨之確定各因素的水平數。
(3)選用合適的正交表,安排試驗計划。首先根據各因素的水平選擇相應水平的正交表。同水平的正交表有好幾個,究竟選哪一個要看因素的個數。一般只要正交表中因素的個數比試驗要考察的因素的個數稍大或相等就行了。這樣既能保證達到試驗目的,又使試驗的次數不至於太多,省工省時。
(4)根據安排的計划進行試驗,測定各試驗指標。
(5)對試驗結果進行計算分析,得出合理的結論。
上述方法一般稱為直觀分析法.這種方法比較簡單,計算量不大,是一種很實用的分析方法。
最后再說明一點,這種方法的主要工具是正交表,而在因素及其水平都確定的情況下,正交表並不是惟一的。
2 多指標的分析方法
在上節的問題中,試驗指標只有一個,考察起來比較方便.但在實際問題中,需要考察的指標往往不止一個,可能有兩個、三個,甚至更多,這都是多指標的問題.下面介紹兩種解決多指標試驗的方法:綜合平衡法和綜合評分法.這兩種方法都能找出使每個指標都盡可能好的試驗方案。
2.1
下面通過具體例子來說明這種方法.
例2
解
表2.1
因素 水平 |
水分A(%) |
1 2 3 |
|
表2.2
因素 實驗號 |
1 A |
各指標實驗結果 |
||||
抗壓強度 |
落下強度(0.5m/次) |
裂紋度 |
||||
1 2 3 4 5 6 7 8 9 |
1 1 1 2 2 2 3 3 3 |
11.5 4.5 11.0 7.0 8.0 18.5 9.0 8.0 13.4 |
1.1 3.6 4.6 1.1 1.6 15.1 1.1 4.6 20.2 |
3 4 4 3 2 0 3 2 1 |
||
抗 壓 強 度 |
K1 K2 K3 |
27 33.5 30.4 42.9 |
|
|
||
k1 k2 k3 |
9.0 11.2 10.1 |
|||||
極差 |
2.2 |
列 |
||||
優方案 |
A2 |
1A |
||||
裂 紋 度 |
K1 K2 K3 |
11 5 6 |
落 下 強 度 |
K1 K2 K3 |
9.3 17.8 25.9 |
|
k1 k2 k3 |
3.7 1.7 2.0 |
k1 k2 k3 |
3.1 5.9 8.6 |
|||
極差 |
2.0 |
極差 |
5.5 |
|||
優方案 |
A2 |
優方案 |
A3 |
圖 |
為便於綜合分析,我們將各指標隨因素水平變化的情況用圖形表示出來,畫在圖2.1中(為了看得清楚,將各點用線段連起來,實際上並不是直線)。
從表2.2看出,對抗壓強度和落下強度來講,粒度的極差都是最大的,也就是說粒度是影響最大的因素,從圖2.1看出,顯然取8最好;對裂紋度來講,粒度的極差不是最大,即不是影響最大的因素,但也是取8最好。總之,對3個指標來講,粒度都是取8最好。
(2)鹼度C對各指標的影響
從表2.2看出,對於3個指標,鹼度的極差都不是最大的,也就是說,鹼度不是影響最大的因素,是較次要的因素。從圖2.1看出,對抗壓強度和裂紋度來講,鹼度取1.1最好,對落下強度來講,鹼度取1.3最好,但取1.1也不是太差,對3個指標綜合考慮,鹼度取1.1為好。
從表2.2看出,對裂紋度來講,水分的極差最大,即水分是影響最大的因素.從圖2.1看出,水分取9最好,但對抗壓強度和落下強度來講,水分的極差都是最小的,即是影響最小的因素。從圖2.1看出,對抗壓強度來講,水分取9最好,取7次之;對落下強度來講,水分取7最好,取9次之.對3個指標綜合考慮,應照顧水分對裂紋度的影響,還是取9為好。
下面要介紹的綜合評分法,在一定意義上講,可以克服綜合平衡法的這個缺點。
例3
表2.3
因素 水平 |
A 時間(h) |
1 2 3 |
|
解
。和例1.1一樣,按L9(34)表排出方案(這里有4個因素,正好將表排滿),進行試驗,將得出的試驗結果列入表2.4中。
表2.4
因素 實驗號 |
1 A |
各指標實驗結果 |
綜合評分 |
|
純 |
回收率 |
|||
1 2 3 4 5 6 7 8 9 |
1 1 1 2 2 2 3 3 3 |
17.5 12.0 6.0 8.0 4.5 4.0 8.5 7.0 4.5 |
30.0 41.2 60.0 24.2 51.0 58.4 31.0 20.5 73.5 |
100 89.2 84.0 56.2 69.0 74.4 65.0 48.5 91.5 |
K1 K2 K3 |
273.2 196.6 205.0 |
|
||
k1 k2 k3 |
91.1 65.5 68.3 |
|||
極差 |
25.6 |
|||
優方案 |
A1 |
根據這個算式,算出每個試驗的分數,列在表2.4最右邊。再根據這個分數,用直觀分析法作進一步的分析,整個分析過程都記錄在表2.4中。
3混合水平的正交試驗設計
前兩節介紹的多因素試驗中,各因素的水平數都是相同的,解決這類問題還是比較簡單的。但是在實際問題中,由於具體情況不同,有時各因素的水平數是不相同的,這就是混合水平的多因素試驗問題。
解決混合水平這類問題一般比較復雜。在這里介紹兩個主要的方法:
(1)直接利用混合水平的正交表;
(2)擬水平法——把水平不同的問題化成水平數相同的問題來處理。
混合水平正交表就是各因素的水平數不完全相等的正交表。這種正交表有好多種,比如L8(41×24)就是一個混合水平的正交表,如表3.1。
表3.1
列號 實驗號 |
1 |
1 2 3 4 5 6 7 8 |
1 1 2 2 3 3 4 4 |
表有8行,5列(注意5=1+4),表示用這張表要做8次試驗,最多可安排5個因素,其中一個是4水平的(第1列),4個是2水平的(第2列到第5列)。
例如,第1列中有4個數字1,2,3,4,它們各出現兩次;第2列到第5列中,都只有兩個數字1,2,它們各出現4次。
(2)每兩列各種不同的水平搭配出現的次數是相同的。但要注意一點:每兩列不同水平的搭配的個數是不完全相同的。
比如,第1列是4水平的列,它和其他任何一個2水平的列放在一起,由行組成的不同的數對一共有8個:(1,1),(1,2),(2,1),(2,2),(3,1),(3,2),(4,1),(4,2),它們各出現1次;第2列到第5列都是2水平列,它們之間的任何兩列的不同水平的搭配共有4個:(1,1),(1,2),(2,1),(2,2),它們各出現兩次。
例3.1
表3.2
因素 水平 |
A 品種 |
B 氮肥量(kg) |
C 氮、磷、鉀肥比例 |
D 規格 |
1 2 3 4 |
甲 乙 丙 丁 |
2.5 3.0 |
3:3:1 2:1:2 |
6×6 7×7 |
解
因素 實驗號 |
1 A |
實驗指標(產量) (kg) |
減去200 |
1 2 3 4 5 6 7 8 |
|
195 205 220 225 210 215 185 190 |
-5 5 20 25 10 15 -15 -10 |
K1 K2 K3 K4 |
25 -25 |
|
|
k1 k2 k3 k4 |
22.5 12.5 -12.5 |
||
極差 |
35.0 |
||
優方案 |
A2 |
這里分析計算的方法和例1.1基本上相同。但是要特別注意,由於各因素的水平數不完全相等,各水平出現的次數也不完全相等,因此計算各因素各水平的平均值k1,k2,k3,k4時和例1.1中有些不同。
比如,對於因素A,它有4個水平,每個水平出現兩次,它的各水平的平均值k1,k2,k3,k4是相應的K1,K2,K3,K4分別除以2得到的。
而對於因素B,C,D,它們都只有兩個水平,因此,只有兩個平均值k1,k2,又因為每個水平出現4次,所以它們的平均值k1,k2是相應的K1,K2分別除以4得到的。這樣得出的平均值才是合理的。
從表3.3看出,因素A的極差最大,因此因素A對試驗的影響最大,並且以取2水平為好;因素B的極差僅次於因素A,對試驗的影響比因素A小,也是以取2水平為好;因素C,D的極差都很小,對試驗的影響也就很小,都是以取2水平為好。
總的說來,試驗方案應以A2B2C2D2為好。但這個方案在做過的8個試驗中是沒有的。按理應當照這個方案再試驗一次,從而確定出真正最好的試驗方案。但是,因為農業生產受節氣的制約,只有到第二年再試驗。
事實上,在這里因為因素D的影響很小,這個方案與8個試驗中的第4號試驗A2B2C2D1很接近,從試驗結果看出,第4號試驗是8個試驗中產量最高的,因此完全有理由取第4號試驗作為最好的試驗方案加以推廣。
3.2
表3.4
因素 水平 |
|
1 2 3 |
|
解 這個問題是4個因素的試驗,其中因素C是2水平的,因素A,B,D是3水平的。這種情況沒有合適的混合水平正交表,因此不能用例3.1的方法解決.對這個問題我們可以設想:假若因素C也有3個水平,那么這個問題就變成4因素3水平的問題,因此可以選正交表L9(34)來安排試驗。但是實際上因素C只有兩個水平,不能隨便安排第3個水平。
如何將c變成3水平的因素呢?我們是從第1、第2兩個水平中選一個水平讓它重復一次作為第3水平,這就叫做虛擬水平。
取哪個水平作為第3水平呢?一般來講,都是根據實際經驗,選取一個較好的水平。比如,如果認為第2水平比第1水平好,就選第2水平作為第3水平。這樣因素水平表3.4就變為表3.5的樣子,它比表3.4多了一個虛擬的第3水平(用方框把它圍起來)。
表3.5
|
下面就按L9(34)表安排試驗,測出結果,並進行分析,整個分析過程記錄在表3.6中。
表3.6
因素 實驗號 |
1 |
實驗指標測試結果 |
1 2 3 4 5 6 7 8 9 |
1 2 2 2 3 3 3 |
45 36 12 15 40 15 10 5 47 |
K1 K2 K3 |
|
|
k1 k2 k3 |
|
|
極差 |
|
|
優方案 |
|
這里要注意的是,因素C的“第3水平”實際上就是第2水平,我們把正交表中第3列的C因素的水平安排又重寫一次,兩邊用虛線標出,對應地列在右邊,這一列是真正的水平安排。由於這一列沒有第3水平,因此在求和時並無K3,只出現K1,K2。又因為這里C的第2水平共出現6次,因此平均值k2是K2除以6,即k2=K2/6;C的第l水平出現3次,平均值k1是K1除以3,即k1=K1/3。
因素A,B,D都是3水平的,各水平都出現3次,因此求平均值k1,k2,k3時,都是K1,K2,K3除以3。
從表3.6中的極差看出,因素D對試驗的影響最大,取第3水平最好;其次是因素A,取第3水平為好;再者是因素B,取第l水平為好;因素C的影響最小,取第1水平為好。
總之,這個試驗的最優方案應當是A3B1C1D3。但是這個方案在做過的9個試驗中是沒有的。從試驗結果看,效果最好的是第8號試驗,這個試驗只有因素B不是處在最好情況,而因素B對試驗的影響是最小的。因此我們選出的最優方案是合乎實際的。我們可以按這個方案再試驗一次,看是否會得到比第8號試驗更好的結果,從而確定出真正的最優方案。
從上面的討論可以看出,擬水平法是將水平少的因素歸入水平數多的正交表中的一種處理問題的方法。在沒有合適的混合水平的正交表可用時,擬水平法是一種比較好的處理多因素混合水平試驗的方法.這種方法不僅可以對一個因素虛擬水平,也可以對多個因素虛擬水平,具體做法和上面相同,不再重復。這里要指出的是:虛擬水平以后的表對所有因素來說不具有均衡搭配性質,但是,它具有部分均衡搭配的性質(部分均衡搭配的精確含義這里就不細講了),所以擬水平法仍然保留着正交表的優點。