正交試驗

正交試驗

一種實驗設計方法
正交試驗設計(Orthogonalexperimentaldesign)是研究多因素多水平的又一種設計方法,它是根據正交性從全面試驗中挑選出部分有代表性的點進行試驗,這些有代表性的點具備了“均勻分散,齊整可比”的特點,正交試驗設計是分式析因設計的主要方法。是一種高效率、快速、經濟的實驗設計方法。日本著名的統計學家田口玄一将正交試驗選擇的水平組合列成表格,稱為正交表。當析因設計要求的實驗次數太多時,一個非常自然的想法就是從析因設計的水平組合中,選擇一部分有代表性水平組合進行試驗。因此就出現了分式析因設計(fractionalfactorialdesigns),但是對于試驗設計知識較少的實際工作者來說,選擇适當的分式析因設計還是比較困難的。例如作一個三因素三水平的實驗,按全面實驗要求,須進行3^3=27種組合的實驗,且尚未考慮每一組合的重複數。若按L9(3^4)正交表安排實驗,隻需作9次,按L15(3^7)正交表進行15次實驗,顯然大大減少了工作量。因而正交實驗設計在很多領域的研究中已經得到廣泛應用。[1]
  • 中文名:正交試驗
  • 外文名:Orthogonal experimental design
  • 别名:
  • 表達式:
  • 提出者:
  • 适用領域:
  • 類 别:一種實驗設計方法
  • 特 點:“均勻分散,齊整可比”

基本思想

正交試驗設計法,就是使用已經造好了的表格--正交表--來安排試驗并進行數據分析的一種方法。它簡單易行,計算表格化,使用者能夠迅速掌握。下邊通過一個例子來說明正交試驗設計法的基本思想。

[例1]為提高某化工産品的轉化率,選擇了三個有關因素進行條件試驗,反應溫度(

A),反應時間(B),用堿量(C),并确定了它們的試驗範圍:

A:80-90℃

B:90-150分鐘

C:5-7%

試驗目的是搞清楚因子A、B、C對轉化率有什麼影響,哪些是主要的,哪些是次要的,從而确定最适生産條件,即溫度、時間及用堿量各為多少才能使轉化率高。試制定試驗方案。

這裡,對因子A,在試驗範圍内選了三個水平;因子B和C也都取三個水平:

A:A1=80℃,A2=85℃,A3=90℃

B:B1=90分,B2=120分,B3=150分

C:C1=5%,C2=6%,C3=7%

當然,在正交試驗設計中,因子可以是定量的,也可以是定性的。而定量因子各水平間的距離可以相等,也可以不相等。

這個三因子三水平的條件試驗,通常有兩種試驗進行方法:

(Ⅰ)取三因子所有水平之間的組合,即A1B1C1,A1B1C2,A1B2C1,……,A3B3C3,共有

3^3=27次

試驗。用圖表示就是圖1立方體的27個節點。這種試驗法叫做全面試驗法。

全面試驗對各因子與指标間的關系剖析得比較清楚。但試驗次數太多。特别是當因子數目多,每個因子的水平數目也多時。試驗量大得驚人。如選六個因子,每個因子取五個水平時,如欲做全面試驗,則需5^6=15625次試驗,這實際上是不可能實現的。如果應用正交實驗法,隻做25次試驗就行了。而且在某種意義上講,這25次試驗代表了15625次試驗。圖1 全面試驗法取點。

(Ⅱ)簡單對比法,即變化一個因素而固定其他因素,如首先固定B、C于B1、C1,使A變化之:

↗A1

B1C1→A2

↘A3(好結果)

如得出結果A3最好,則固定A于A3,C還是C1,使B變化之:

↗B1

A3C1→B2(好結果)

↘B3

得出結果以B2為最好,則固定B于B2,A于A3,使C變化之:

↗C1

A3B2→C2(好結果)

↘C3

試驗結果以C2最好。于是就認為最好的工藝條件是A3B2C2。

這種方法一般也有一定的效果,但缺點很多。首先這種方法的選點代表性很差,如按上述方法進行試驗,試驗點完全分布在一個角上,而在一個很大的範圍内沒有選點。因此這種試驗方法不全面,所選的工藝條件A3B2C2不一定是27個組合中最好的。其次,用這種方法比較條件好壞時,是把單個的試驗數據拿來,進行數值上的簡單比較,而試驗數據中必然要包含着誤差成分,所以單個數據的簡單比較不能剔除誤差的幹擾,必然造成結論的不穩定。

簡單對比法的最大優點就是試驗次數少,例如六因子五水平試驗,在不重複時,隻用5+(6-1)×(5-1)=5+5×4=25次試驗就可以了。

考慮兼顧這兩種試驗方法的優點,從全面試驗的點中選擇具有典型性、代表性的點,使試驗點在試驗範圍内分布得很均勻,能反映全面情況。但我們又希望試驗點盡量地少,為此還要具體考慮一些問題。

如上例,對應于A有A1、A2、A3三個平面,對應于B、C也各有三個平面,共九個平面。則這九個平面上的試驗點都應當一樣多,即對每個因子的每個水平都要同等看待。具體來說,每個平面上都有三行、三列,要求在每行、每列上的點一樣多。這樣,作出如圖2所示的設計,試驗點用⊙表示。我們看到,在9個平面中每個平面上都恰好有三個點而每個平面的每行每列都有一個點,而且隻有一個點,總共九個點。這樣的試驗方案,試驗點的分布很均勻,試驗次數也不多。

當因子數和水平數都不太大時,尚可通過作圖的辦法來選擇分布很均勻的試驗點。但是因子數和水平數多了,作圖的方法就不行了。

試驗工作者在長期的工作中總結出一套辦法,創造出所謂的正交表。按照正交表來安排試驗,既能使試驗點分布得很均勻,又能減少試驗次數,圖2正交試驗設計圖例而且計算分析簡單,能夠清晰地闡明試驗條件與指标之間的關系。用正交表來安排試驗及分析試驗結果,這種方法叫正交試驗設計法。

正交表的性質

(1)每列中不同數字出現的次數是相等的,如L9(),每列中不同的數字是1,2,3,它們各出現3次;

(2)在任意兩列中,将同一行的兩個數字看成有序數對時,每種數對出現的次數是相等的,如L9(),有序數對共有9個:(1,1),(1,2),(1,3),(2,1),(2,2),(2,3),(3,1),(3,2),(3,3),它們各出現一次。

由于正交表有這兩條性質,用它來安排試驗時,各因素的各種水平的搭配是均衡的。

正交表

為了叙述方便,用L代表正交表,常用的有L8(),L9(),L16(),L8(4×),L12(),等等。此符号各數字的意義如下  :

L8()

7為此表列的數目(最多可安排的因子數)

2為因子的水平數

8為此表行的數目(試驗次數)

L16(2×)

有7列是3水平的

有1列是2水平的

L16(2×)的數字告訴我們,用它來安排試驗,做16個試驗最多可以考察一個2水平因子和7個3水平因子。

在行數為mn型的正交表中(m,n是正整數),試驗次數(行數)=Σ(每列水平數-1)+1

如L8(),

8=7×(2-1)+l

利用上述關系式可以從所要考察的因子水平數來決定最低的試驗次數,進而選擇合适的正交表。比如要考察五個3水平因子及一個2水平因子,則起碼的試驗次數為5×(3-1)+1×(2-1)+1=12(次)

這就是說,要在行數不小于13,既有2水平列又有3水平列的正交表中選擇,L16(2×)适合。

正交表具有兩條性質:(1)每一列中各數字出現的次數都一樣多。(2)任何兩列所構成的各有序數對出現的次數都一樣多。所以稱之謂正交表。

例如在L9()中(見表1),各列中的1、2、3都各自出現3次;任何兩列,例如第3、4列,所構成的有序數對從上向下共有九種,既沒有重複也沒有遺漏。其他任何兩列所構成的有序數對也是這九種各出現一次。這反映了試驗點分布的均勻性。

方案設計

安排試驗時,隻要把所考察的每一個因子任意地對應于正交表的一列(一個因子對應一列,不能讓兩個因子對應同一列),然後把每列的數字"翻譯"成所對應因子的水平。這樣,每一行的各水平組合就構成了一個試驗條件(不考慮沒安排因子的列)。

對于[例1],因子A、B、C都是三水平的,試驗次數要不少于

3×(3-1)+1=7(次)

可考慮選用L9()。因子A、B、C可任意地對應于L9()的某三列,例如A、B、C分别放在l、2、3列,然後試驗按行進行,順序不限,每一行中各因素的水平組合就是每一次的試驗條件,從上到下就是這個正交試驗的方案,見表2。這個試驗方案的幾何解釋正好是圖2。

三個3水平的因子,做全面試驗需要3*3*3=27次試驗,現用L9()來設計試驗方案,隻要做9次,工作量減少了2/3,而在一定意義上代表了27次試驗.。

再看一個用L9()安排四個3水平因子的例子。

[例2]某礦物氣體還原試驗中,要考慮還原時間(A)、還原溫度(B)、氣體流速(C)、還原氣體比例(D)這四個因子對全鐵含量X〔越高越好)、金屬化率Y(越高越好)、二氧化钛含量Z(越低越好)這三項指标的影響。希望通過試驗找出主要影響因素,确定最适工藝條件。

首先根據專業知識以确定各因子的水平:

時間:A1=3(小時),A2=4(小時),A3=5(小時)

溫度:B1=1000(℃),B2=1100(℃),B3=1200(℃)

流速:Cl=600(毫升/分),C2=400(毫升/分),

C3=800(毫升/分)

CO:H2:D1=1:2,D2=2:1,D3=1:1

這是四因子3水平的多指标(X、Y、Z)問題,如果做全面試驗需3^4=81次試驗,而用L9( )來做隻要9次。具體安排如表3。同全面試驗比較,工作量少了8/9。由于縮短了試驗周期,可以提高試驗精度,時間越長誤差幹擾越大。并且對于多指标問題,采用簡單對比法,往往顧此失彼,最适工藝條件很難找;而應用正交表來設計試驗時可對各指标通盤考慮,結論明确可靠。

數據分析

正交表的另一個好處是簡化了試驗數據的計算分析。還是以[例1]為例來說明。按照表2的試驗方案進行試驗,測得9個轉化率數據。

通過9次試驗,我們可以得兩類收獲。

第一類收獲是拿到手的結果。第9号試驗的轉化率為64,在所做過的試驗中最好,可取用之。因為通過L9()已經把試驗條件均衡地打散到不同的部位,代表性是好的。假如沒有漏掉另外的重要因素,選用的水平變化範圍也合适的話,那麼,這9次試驗中最好的結果在全體可能的結果中也應該是相當好的了,所以不要輕易放過。

第二類收獲是認識和展望。9次試驗在全體可能的條件中(遠不止3^3=27個組合,在試驗範圍内還可以取更多的水平組合)隻是一小部分,所以還可能擴大。精益求精。尋求更好的條件。利用正交表的計算分折,分辨出主次因素,預測更好的水平組合,為進一步的試驗提供有份量的依據。

其中I、Ⅱ、Ⅲ分别為各對應列(因子)上1、2、3水平效應的估計值,其計算式是:

Ⅰi(Ⅱi,Ⅲi)=第i列上對應水平1(2,3)的數據和

K1為1水平數據的綜合平均=Ⅰ/水平1的重複次數

Si為變動平方和=

[例1]的轉化率試驗數據與計算分析見表4。

先考慮溫度對轉比率的影響。但單個拿出不同溫度的數據是不能比較的,因為造成數據差異的原因除溫度外還有其他因素。但從整體上看,80℃時三種反應時間和三種用堿量全遇到了,85℃時、90℃時也是如此。這樣,對于每種溫度下的三個數據的綜合數來說,反應時間與加堿量處于完全平等狀态,這時溫度就具有可比性。所以算得三個溫度下三次試驗的轉化率之和:

80℃:ⅠA=x1+x2+x3=31+54+38=123;

85℃:ⅡA=x4+x5+x6=53+49+42=144;

90℃:ⅢA=x7+x8+x9=57+62+64=183。

分别填在A列下的Ⅰ、Ⅱ、Ⅲ三行。再分别除以3,表示80℃、85℃、90℃時綜合平均意義下的轉化率,填入下三行K1、K2、K3。R行稱為極差,表明因子對結果的影響幅度。

同樣地,為了比較反應時間;用堿量對轉化率的影響,也先算出同一水平下的數據和IB、ⅡB、ⅢB,IC、ⅡC、ⅢC,再計算其平均值和極差。都填入表4中;

由此分别得出結論:溫度越高轉化率越好,以90℃為最好,但可以進一步探索溫度更好的情況。反應時間以120分轉化率最高。用堿量以6%轉化率最高。所以最适水平是A3B2C2。

正交試驗的方差分析

(一)假設檢驗

在數理統計中假設檢驗的思想方法是:提出一個假設,把它與數據進行對照,判斷是否舍棄它 [5]。其判斷步驟如下:

(1)設假設H0正确,得到一個理論結論,設此結論為R0;

(2)再根據試驗得出一個試驗結論,與理論結論相對應,設為R1;

(3)比較R0與R1:若R0與R1沒有大的差異,則沒有理由懷疑H0,從而判定為:"不舍棄H。"(采用H。);若R0與R1有較大差異,則可以懷疑H0,此時判定為:"舍棄H0"。

但是,R1/R0比值為多少才能舍棄H0呢?為确定這個量的界限,需要利用數理統計中F分布的理論。

若yl服從自由度為φ1的χ2分布,y2服從自由度為φ2的χ2分布,并且yl、y2相互獨立,則(y1/φ1)/(y2/φ2)服從自由度為(φ1,φ2)的F分布。F分布是連續分布,分布模數是兩個自由度(φ1,φ2)。稱φ1為分子自由度,稱φ2為分母自由度。在自由度為(φ1,φ2)的F分布中,某點右側面積為p,也就是F比此值大的概率為p,把這個值寫為(p)。若檢驗的顯著性水平(或危險率)給定為α時,則可以把 (α)作為臨界值來檢驗假設。

這裡,Se/σ2服從自由度為φe,的χ2分布;當H。成立,σ2=0時,SA/σ2也服從自由度為φA的χ2分布;又SA與Se相互成立,所以(SA/(φAσ2)/Se/(φeσ2))=VA/Ve服從自由度為(φA,φe)的F分布。這就是假定H。正确時的理論結論R。。而試驗結論Rl要與理論結論R。相比較。由給定的顯著性水平,通常是α=0.05;分子自由度φ1=φA=a-1,分母自由度φ2=φe=a(n-1);查F分布表得出 (α)。所以H。:α1=α2=……=αa=0(σA2=0)的檢驗是:(顯著性水平α)

FA=VA/Ve>(α)→ 舍棄H0;

FA=VA/Ve≤(α)→ 不舍棄H0;

通常,(α)一般性地表示成Fα(φA,φB)。

假設因子A對試驗結果的影響不顯著,那麼A的兩個水平的效應該表現為相等或相近,即假設H0:α1=α2=0。如果因子A顯著,則舍棄假設。

為了判斷因子A是否顯著,首先要計算比值顯然,這個比值越大,因子A對指标的影響越顯著;反之,因子A就不顯著。在給定置信度α後,如α=0.05,查F分布表,自由度φA是因子A的,自由度φe是誤差的,其臨界值Fα(φA,φe),如果FA>Fα(φA,φe)就舍棄假設,可以認為因子A是顯著的;如果FA≤Fα(φA,φe)就沒有理由否定假設,而隻能認為因子A是不顯著的。因為按照F分布表的物理念義,F值小于Fα(φA,φe)的概率是95%,即有95%的機會出現小于Fα(φA,φe)的F值,既然出現了這種情況,就有了95%的把握,所以就沒有理由否定假設,隻能接受假設,認為因子A不顯著。另一方面,F值大于Fα(φA,φe)的概率是5%,也就是隻有5%的機會出現大于Fα(φA,φe)的F值,這是小概率事件,如果小概率事件居然發生了,則可認為情況異常,假設不可信,必須否定假設,因子A是顯著的。對其他因子的顯著性檢驗完全類似。

(二)方差分析表

由總平方和與各因素平方和即可求得誤差平方和,亦稱剩餘平方和。是總平方和減各因素平方和所得。如正交表有一空列,則該列的平方和就是誤差平方和。但在正交表飽和試驗的情況下,即所有各列全部排滿時,誤差平方和一般用各因素平方和中幾個最小的平方和之和來代替,同時,這幾個因素不再作進一步的分析。

自由度:φT=試驗次數一1

φA,B…=水平數一1

φA×B=φA×φB

φe=φT-φA-φB-……-φD

相關詞條

相關搜索

其它詞條