正交試驗設計
正交試驗設計是利用正交表來安排與分析多因素試驗的一種設計方法。它是由試驗因素的 全部水平組合中,挑選部分有代表性的水平組合進行試驗的,通過對這部分試驗結果的分 析瞭解全面試驗的情況,找出最優的水平組合。可以通過代表性很強的少數次試驗,摸清 各個因素對試驗指標的影響情況,確定出因素的主次順序,找出較好的生產條件或最佳參 數組合。
<!–more–>
1. 正交表的性質
1.1. 正交性
- 均勻分散性:任一列中,不同數字出現的次數相等;
- 整齊可比性:任兩列中,同一橫行所組成的數字對出現的次數相等;
- 可置換性:
- 正交表各列的地位是平等的,表中各列之間可以相互置換,稱爲列間置換;
- 正交表各行之間也可相互置換,稱列間置換;
- 正交表中同一列的水平數字也可以相互置換,稱水平置換。
1.2. 分散均勻性
- 任一列的各水平都出現,使得部分試驗中包含所有因素的所有水平;
- 任意兩列間的所有組合全部出現,使任意兩因素間都是全面試驗。
1.3. 綜合可比性
- 任一列各水平出現的次數都相等;
- 任兩列間所有可能的組合出現的次數都相等,因此使任一因素各水平的試驗條件相同。
2. 拉丁方陣
拉丁方陣(Latin square)是一種 \(n\times n\) 的方陣,在這種 $n× n$的方陣裏, 恰有 $n$種不同的元素,每一種不同的元素在同一行或同一列裏只出現一次。如下面的三階 拉丁方。
\begin{bmatrix} 1 & 2 & 3\\ 2 & 3 & 1\\ 3 & 1 & 2\\ \end{bmatrix}
2.1. 標準型
當一個拉丁方陣的第一行與第一列的元素按順序排列時,此爲這個拉丁方陣的標準型,英語 稱爲"reduced Latin square, normalized Latin square, or Latin square in standard form"。
2.2. 構造拉丁方
同階的拉丁方有多個,構造其中的一個有一種通用的方法。
若 \(n\) 爲偶數,則構造第一行爲\(1, 2, n, 3, n-1, 4, n-2, \dots\),之後的每一行都在 上一行對應列元素上加一併對\(n\) 取模,即可表示爲如下矩陣。
\begin{bmatrix} 1 & 2 & n & 3 & n-1 & 4 & n-2 & \cdots\\ 2 & 3 & n+1 & 4 & n & 5 & n-1 & \cdots\\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots & \vdots\\ i & i+1 & i+n-1 & i+2 & i+n-2 & i+3 & i+n-3 & \cdots\\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots & \vdots\\ n & 1 & n-1 & 2 & n-2 & 3 & n-3 & \cdots\\ \end{bmatrix} \mathrm{mod}\ n
若 \(n\) 爲奇數,則先按偶數法則構造一個拉丁方,然後將上述模式左右對調。
2.3. 正交拉丁方
設有兩個階數相同且爲 $n$的拉丁方陣\(A_1=(a^{(1)}_{i,j})_{n\times n}, A_2=(a^{(2)}_{i,j})_{n\times n}\),其中將所有放置位置相同的元素組合成一個元組,組 合成一個新的矩陣\(((a^{(1)}_{i,j}, a^{(2)}_{i,j}))_{n\times n}\)。當這個新的矩陣 $((a(1)i,j, a(2)i,j))n× n$中每一個元素互不相同時,拉丁方陣 \(A_1\) 和 $A2$是互相正交的。此時,\(A_1\) 和 \(A_2\) 即爲一對正交拉丁方。而在階數固 定的情況下,所有兩兩正交的拉丁方所成的集合稱爲正交拉丁方族。如下面的兩個3 階拉丁 方。
\begin{matrix} A = \begin{bmatrix} 1 & 2 & 3\\ 2 & 3 & 1\\ 3 & 1 & 2\\ \end{bmatrix} & B = \begin{bmatrix} 1 & 2 & 3\\ 3 & 1 & 2\\ 2 & 3 & 1\\ \end{bmatrix} \end{matrix}
2.4. 素數階拉丁方
若 \(n\) 階拉丁方存在 \(r\) 個兩兩正交的拉丁方,那麼 \(r\leq n-1\)。素數階拉丁方(以 \(3, 5, 7, \dots\) 爲階)是完全的,即 $n$階正交拉丁方族中有 \(n-1\) 個拉丁方。
3. 構造正交表
因爲素數階拉丁方都可構造完全的正交拉丁方族,則可繼續構造正交表。若把素數 \(n\) 作爲 水平數,則以 \(n\) 階拉丁方構造的正交拉丁方有 \(n-1\) 個。以 $n-1$作爲因素數構造正交 表。
以 3 階拉丁方爲例,3 階正交拉丁方族包含 2 個正交拉丁方,並將對應位置元素合併成二 元有序數對。
\begin{matrix} \begin{bmatrix} 1 & 2 & 3\\ 2 & 3 & 1\\ 3 & 1 & 2\\ \end{bmatrix} & \begin{bmatrix} 1 & 2 & 3\\ 3 & 1 & 2\\ 2 & 3 & 1\\ \end{bmatrix} & = & \begin{bmatrix} 11 & 22 & 33\\ 23 & 31 & 12\\ 32 & 13 & 21\\ \end{bmatrix} & \end{matrix}
以第一個元素作爲因素 A,第二個元素作爲因素 B,然後與因素 C、D 均勻搭配,構造 4 因素 3 水平正交表。
C, A, B, D | 1 | 2 | 3 |
---|---|---|---|
1 | 11 | 22 | 33 |
2 | 23 | 31 | 12 |
3 | 32 | 13 | 21 |
展開後爲
試驗號 | A | B | C | D |
---|---|---|---|---|
1 | 1 | 1 | 1 | 1 |
2 | 1 | 2 | 2 | 2 |
3 | 1 | 3 | 3 | 3 |
4 | 2 | 2 | 3 | 1 |
5 | 2 | 3 | 1 | 2 |
6 | 2 | 1 | 2 | 3 |
7 | 3 | 3 | 2 | 1 |
8 | 3 | 1 | 3 | 2 |
9 | 3 | 2 | 1 | 3 |
若只與一個因素 C 搭配,則可構造 3 因素 3 水平正交表。
C, A, B | |||
---|---|---|---|
1 | 11 | 22 | 33 |
2 | 23 | 31 | 12 |
3 | 32 | 13 | 21 |
展開後爲
試驗號 | A | B | C |
---|---|---|---|
1 | 1 | 1 | 1 |
2 | 1 | 2 | 2 |
3 | 1 | 3 | 3 |
4 | 2 | 2 | 3 |
5 | 2 | 3 | 1 |
6 | 2 | 1 | 2 |
7 | 3 | 3 | 2 |
8 | 3 | 1 | 3 |
9 | 3 | 2 | 1 |
由此可見,以素數 \(n\) 爲水平數設計正交試驗,可滿足 $[n-1, n+1]$個因素的設計,並且 總共需要 \(n^2\) 次試驗。即對於 3 因素 3 水平與 4 因素 3水平正交試驗具有相同的試驗 組數,同爲 9 組。可類似地構造出\(L_9(3^3)\),\(L_9(3^4)\),\(L_{25}(5^5)\), \(L_{25}(5^6)\),\(L_{49}(7^7)\),$L49(78)$等。