◆ 俞鐘行 / 文

響應(yīng)曲面與混料（也稱“配方”）設(shè)計，每一種都是在DOE（試驗設(shè)計）領(lǐng)域特色鮮明且重要的品類，在六西格瑪管理（第二版）[1]中基于“五步法”作了詳細介紹，所用軟件為Minitab。本文基于excel和因素趨勢圖的方法[2]來分析同例的數(shù)據(jù)，得到結(jié)果相同但在綜合成本上有明顯下降。

1 響應(yīng)曲面法

1.1原例7-8“粘合劑生產(chǎn)條件優(yōu)化問題”簡解：在粘合劑生產(chǎn)中，經(jīng)過因子的篩選，得關(guān)鍵因子A(溫度、temp)，低水平200度、高水平300度；關(guān)鍵因子B（時間、time），低水平40秒、高水平70秒。為求最大（粘性、stick ）Y，做全因子試驗22，又在中心點250度、55秒也做了3次試驗，試驗結(jié)果見表1。

原例后來又補做4次星號試驗，構(gòu)成11個點的CCC（中心復合序貫設(shè)計），最后得到回歸方程：y=-377.43+0.056A+0.109B-0.000002A2-0.000051B2，相關(guān)系數(shù)R=0.996，標準誤差s=1.269。最優(yōu)解為：A(溫度)=2 6 0.5 6度、B（時間）=67.98秒時，粘性Y=50.68。

1.2用“因素趨勢法”求解：先做常規(guī)的極差分析，并畫因素趨勢圖，見表2、圖1。

表1 全因子試驗22+3次中心點試驗

表2 22+3試驗的極差分析

如果直接對表1用excel“數(shù)據(jù)分析”中“回歸”模塊分析，擬合效果很差，如復相關(guān)系數(shù)Multiple R= 0.563599、標準誤差s= 9.917526、方程p 值=0.465611等。但是圖1提示我們：兩個因子強度相當，且都呈現(xiàn)還算比較規(guī)范的“角”狀，根據(jù)經(jīng)驗它們都是“二次項”。經(jīng)操作，如果把A、B、A2和B2都用來做回歸擬合的話，由于總的試驗次數(shù)才7次，得不到正確結(jié)果。最后插入A2和B2，刪去B，按表3情況回歸，得回歸結(jié)果見圖2。

圖1 粘度Y的因素趨勢圖

表3 22+3試驗（回歸前變項）

圖2 經(jīng)變項后所得的回歸結(jié)果

總的來說，這個回歸結(jié)果相當令人滿意，與原例（多做4次試驗）相比回歸系數(shù)基本相等，只標準差大了7%。圖2中最下表中各項的P-value也全都小于0.01。得到的回歸方程為：

現(xiàn)在用excel的“規(guī)劃求解”對方程選優(yōu)（望大）。結(jié)果是，無論初始值設(shè)為0，還是設(shè)為表1中實測結(jié)果最好的A=250、B=55，最后得到的最優(yōu)方案都是A=258.5276度、B=70秒、Y=53.36508?？偟膩砜?，這個結(jié)果至少不比原方案差，但在試驗次數(shù)及數(shù)據(jù)分析成本上都顯著降低了。

表4 四因子二階單純形格子點混料設(shè)計

2 混料設(shè)計（配方設(shè)計）

2.1原例7-11“新型高強度合金鋼配方設(shè)計問題”：在新型高強度合金冶煉技術(shù)中，關(guān)鍵是添加劑的配比問題。添加劑由A、B、C及D共4種成分構(gòu)成。采用四因子二階單純形格子點設(shè)計，共進行10次試驗，其試驗方案和結(jié)果見表4。

原例用Minitab軟件，最后求得當A=0.6636、D=0.3364時，可得最大強度145.77。

2.2用“因素趨勢法”求解：先做常規(guī)的極差分析，并畫因素趨勢圖，見表5、圖3。

在因素趨勢圖的提示下挑選交互作用項時，宜優(yōu)先考慮至少其中之一是最強因子。如這里A是最強因子，而且它的走勢也與B和D明顯交叉。當然這也是經(jīng)excel“回歸”模塊的實際操作驗證的。最終插入3個交互作用項AB、AD和BC，如表6所示。

對表6用excel的“回歸”模塊分析，結(jié)果為圖4，各個指標都令人滿意，僅BC的P-value略大于0.05。

最后就是用excel的“規(guī)劃求解”對方案選優(yōu)，但要注意“混料設(shè)計”的特殊性。假設(shè)在excel電子表格的A62：C67處設(shè)置好規(guī)劃求解區(qū)域的標識（參見表7，但目前應(yīng)認為表7中的數(shù)字都是沒有的）。注意在C67處鍵入約束條件：“=SUM（C63:C66）”，然后在A63處輸入回歸方程：y=142.152A+136.3186 B+133.2402C+132.4069D+11.05 882AB+30.88235AD-7.11765BC（系數(shù)從圖4最下表的Coefficients列獲得，方程中的A～D就是電子表格的C63：C66）。在規(guī)劃求解的對話框里的設(shè)置見圖5，注意到約束條件里有個$C$67=1就是針對混料設(shè)計約束條件的。

表5 四因子二階單純形格子點混料設(shè)計極差分析

圖3 混料設(shè)計因素趨勢圖

無論是初始條件設(shè)為0，還是設(shè)為表6中實測最佳的第4次條件A=0.5、B=0、C=0和D=0.5，規(guī)劃求解后的結(jié)果都如表7所示。應(yīng)當認為這與原例用Minitab所得結(jié)果是一樣的。

表6 插項后的excel電子表格界面

圖4 經(jīng)插項后的回歸結(jié)果

圖5 規(guī)劃求解的設(shè)置

表7 規(guī)劃求解的界面