正交多項式回歸設(shè)計及參數(shù)設(shè)計

1、正交多項式回歸設(shè)計&參數(shù)設(shè)計,正交多項式回歸設(shè)計,定義 將正交試驗法與多項式回歸分析結(jié)合 特點 兼顧了兩種技術(shù)的優(yōu)點，且計算簡單 正交試驗法 可用較少的試驗次數(shù)，獲得能反映全面試驗的情況；通過對試驗結(jié)果的ANOVA，可估計若干因素影響的相對大小及因素間的相互關(guān)系，并利用此種關(guān)系在一定置信度下由各因素的取值去預(yù)測響應(yīng)值的范圍 多項式回歸 可將響應(yīng)值控制在某一區(qū)間內(nèi)，反向確定各個影響因素的取值范圍；尤其當(dāng)影響因素個數(shù)較多，處理數(shù)據(jù)的工作量較大時可大大簡化計算 前提 僅適用于自變量取等間隔數(shù)值的情況,多項式回歸,若x可控，且其水平取值的間距為h（并非均為1），則任一組等距點為x1=a+h，x2=a+

2、2h，xn=a+nh 通過下式轉(zhuǎn)化為一組標準等距點（h=1）1,2,t,n 設(shè)對應(yīng)于xi = t的試驗結(jié)果為yt，則可產(chǎn)生一個k次多項式 設(shè)1（x）、 2（x）、 k（x）分別為x的一、二、及k次多項式，則可見,Cont,k次線性回歸方程的偏回歸系數(shù)由正規(guī)方程組決定 為簡化計算，同時令 （即正交性） 求解偏回歸系數(shù)和截距,Cont,正交多項式 xt的取值為標準等距點，n為因素的水平數(shù) 由于i（x）的值并非均為整數(shù)，則通常引進適當(dāng)?shù)南禂?shù)，使其在幾個整數(shù)點上的值變?yōu)檎麛?shù),Cont,對于給定的n，相應(yīng)的i及i（x）在1，2，n各整數(shù)點的數(shù)值及 均已制表,Cont,Cont,Cont,Cont,Con

3、t,利用正交多項式回歸的計算步驟,根據(jù)因素水平數(shù)n，設(shè)需配一個k次多項式（k5），n個水平至多配（n-1）階的多項式 依次求 回歸方程為,Cont,每次多項式i（x）的系數(shù)bi及相應(yīng)的Bi只與yt及i（x）有關(guān)，而不隨其他各次多項式的增減而變化；在整個回歸中多配一項i（x）將使回歸平方和增加一項biBi，故第i次多項式i（x）的效應(yīng)為Pi=biBi=Bi/si，而回歸平方和則是各次效應(yīng)的和 方差分析表,Cont,為考察甲醛濃度x與縮醛化度y之間的定量關(guān)系，對7種不同甲醛濃度各進行了若干次試驗，測出各種濃度的平均縮醛化度 配一個4次多項式的回歸方程 將x變?yōu)橐唤M標準等距點x（1，2，7） 利用n

4、=7做正交多項式，則回歸方程變?yōu)?Cont,Cont,由于一次項高度顯著，故線性方程為,Cont,由于二次項也顯著，故二次方程為 由于三次項與四次項不顯著，故無需再配,Cont,各個點的回歸計算值與實測值的比較,正交多項式回歸設(shè)計與回歸方程的建立,若用正交表安排試驗，將因素水平的間距取為相等，則可用正交多項式回歸來處理正交試驗的結(jié)果，定量地描述響應(yīng)值與各因素間的關(guān)系 各因素的總效應(yīng)函數(shù)可表示為各因素效應(yīng)之和 式中，P(A)、P(B)、P(C)分別為因素A、B、C等的效應(yīng)函數(shù)，可按正交多項式展開,Cont,各項回歸系數(shù)bk和常數(shù)項b0的計算 式中， ，r為同一因素水平下的重復(fù)試驗次數(shù) 求出各個b

5、和b0，即可建立回歸方程,Cont,對于多因素等水平間距試驗，選用正交表 安排,Cont,標準等距點的轉(zhuǎn)化（式1-6） 按式1-6正交多項式展開，對式1-2（二水平），各因素只需配1項；而對式3-6（四水平），各因素則需配3項（具體展開式見后頁） 由于二水平因素需進行8次試驗，每次試驗進行三次重復(fù)測定，故總的重復(fù)次數(shù)r =8*3次；同理，各個四水平因素的正交試驗次數(shù)則為r =4*3次 計算各自回歸系數(shù),Cont,其余回歸系數(shù)同理計算,Cont,將上述回歸系數(shù)帶入效應(yīng)函數(shù)公式，即可得回歸方程 回歸方程及各回歸系數(shù)的假設(shè)檢驗,Cont,以上為最優(yōu)回歸方程 手工計算過于繁瑣，故需借助軟件實現(xiàn),產(chǎn)品的

6、三次設(shè)計,一種產(chǎn)品的質(zhì)量，是設(shè)計質(zhì)量、制造質(zhì)量和使用質(zhì)量的總和；而前者是形成產(chǎn)品質(zhì)量關(guān)鍵的第一步，以防“先天不足，后患無窮” 美國質(zhì)量管理專家JurAn博士認為，設(shè)計質(zhì)量占整個產(chǎn)品質(zhì)量的比率為60% 日本質(zhì)量工程專家田口玄一博士認為，設(shè)計質(zhì)量（包括產(chǎn)品設(shè)計與工藝設(shè)計）占整個產(chǎn)品質(zhì)量的比率為70%左右 田口博士在20世紀70年代首次提出“三次設(shè)計”的概念。認為產(chǎn)品的質(zhì)量，首先是設(shè)計出來的，其次才是創(chuàng)造和檢驗出來的,Cont,概念 設(shè)計某種性能的產(chǎn)品以滿足市場的需求，一般要經(jīng)歷三個階段，系統(tǒng)設(shè)計、參數(shù)設(shè)計與容差設(shè)計 本質(zhì) 在專業(yè)設(shè)計的基礎(chǔ)上，用正交設(shè)計的方法選擇最佳組合和最合理的容差范圍，盡量用價

7、格低廉的、低等級的零部件進行組裝的優(yōu)化設(shè)計方法,Cont,系統(tǒng)設(shè)計（System Design） 又稱基礎(chǔ)設(shè)計，為三次設(shè)計的第一步 專業(yè)技術(shù)人員運用專業(yè)知識與技術(shù)，利用系統(tǒng)工程的思想與方法，對產(chǎn)品的結(jié)構(gòu)、性能、壽命、材料等進行綜合設(shè)計，以探討產(chǎn)品如何最經(jīng)濟、最合理地滿足市場用戶需求 該階段質(zhì)量的好壞完全取決于專業(yè)技術(shù)人員的能力和經(jīng)驗；而試驗設(shè)計方法在此階段不起任何作用 故需要數(shù)學(xué)領(lǐng)域?qū)＜遗浜蠈I(yè)技術(shù)人員，盡力尋找產(chǎn)品性能指標與系統(tǒng)中各參數(shù)（即因素）之間的函數(shù)關(guān)系，才會為后續(xù)的參數(shù)設(shè)計帶來極大方便,Cont,參數(shù)設(shè)計（Parameter Design） 在系統(tǒng)設(shè)計的基礎(chǔ)上，對影響產(chǎn)品特征的各項參

8、數(shù)及水平，運用試驗設(shè)計的技術(shù)方法，目的在于尋找試驗因素的最優(yōu)參數(shù)組合 常見分類 “直接擇優(yōu)” 若性能指標具有某個確定的目標值，或具有越大（或越?。┰胶弥匦?，試驗?zāi)康膬H在于使指標值達到或接近優(yōu)良狀態(tài) “穩(wěn)定性擇優(yōu)” 通過良好的參數(shù)組合，來提高產(chǎn)品指標的穩(wěn)定性。一般地，按照一組優(yōu)良參數(shù)設(shè)計制造的產(chǎn)品，將表現(xiàn)在其目標特性的波動性小，穩(wěn)定性好。,Cont,本質(zhì) 乃產(chǎn)品質(zhì)量的優(yōu)化設(shè)計，屬于三次設(shè)計的核心階段；主要利用非線性設(shè)計的性質(zhì)減少輸出響應(yīng)的波動，減少質(zhì)量損失 步驟 分析、明確問題的提出，選擇因素及其水平 選擇正交設(shè)計表，確定試驗方案 進行試驗，測出所需的響應(yīng)特征值 進行數(shù)據(jù)的統(tǒng)計分析 確定最佳方

9、案,Cont,容差設(shè)計（Tolerance Design） 對產(chǎn)品的質(zhì)量、成本和市場等問題綜合考慮，做出規(guī)劃 基本思想 參數(shù)設(shè)計只決定各因素的參數(shù)中心值，但其實各參數(shù)均由一個公差（公差的一半即容差）。明顯地，當(dāng)各因素的參數(shù)在公差范圍內(nèi)變動時，產(chǎn)品的性能也隨之改變 要使產(chǎn)品的性能接近目標值，各因素的公差應(yīng)該多少為宜？這就是第三次設(shè)計需要考慮的重點問題 一般地，除了使產(chǎn)品滿足容差的要求，還需要考慮主要影響因素，是否可以用波動范圍較小的一級品、二級品替代三級品，從而提高產(chǎn)品質(zhì)量，降低質(zhì)量損失，增加經(jīng)濟收益，但同時成本會上升，因此需要衡量得失利弊,Cont,三次設(shè)計的特點 基礎(chǔ)設(shè)計 不用試驗設(shè)計 參數(shù)

10、設(shè)計 不需要增加單位制造費就能提高產(chǎn)品質(zhì)量，是 最便宜的設(shè)計 容差設(shè)計 只能在參數(shù)設(shè)計之后進行，且產(chǎn)品的單位制造 費用很高 穩(wěn)定性擇優(yōu)設(shè)計就集中在如何有效地進行參數(shù)設(shè)計階段,參數(shù)設(shè)計的基本思想（1）,產(chǎn)品的質(zhì)量用產(chǎn)品帶來的損失大小來衡量 “理想質(zhì)量” 用戶總期望購買的產(chǎn)品在給定的使用條件下和使用壽命期內(nèi)，既能達到既定的目標性能，又無有害的副作用，這也是用戶評價質(zhì)量的參考點 實際上，由于產(chǎn)品間的差異、使用環(huán)境的不同或產(chǎn)品提前老化等諸多原因，“理想質(zhì)量”會有不同程度的偏離，自然也就給供需雙方，乃至社會帶來不同程度的損失 一般認為，理想產(chǎn)品的損失為0，損失越小，則質(zhì)量越高，需方人數(shù)將越多，供方利潤也

11、越豐厚,參數(shù)設(shè)計的基本思想（2）,產(chǎn)品的損失用質(zhì)量損失函數(shù)來度量 產(chǎn)品的特性值y不僅與目標值m之間存在差異（由表示容差），同時還可能因為來自生產(chǎn)條件、使用環(huán)境即時間因素等多方面的干擾而發(fā)生變化，產(chǎn)生波動，造成損失 當(dāng)產(chǎn)品的y落在m- ，m+的公差范圍內(nèi)，即為合格品 反之，則為不合格（需報廢、降級或返修）或次品（需方） 設(shè)質(zhì)量波動超過作為不合格處理時，給供方造成的損失為A元；而作為次品出廠，給需方造成的損失為D元；D A 若質(zhì)量波動未超過，仍會因為少許的波動，導(dǎo)致y偏離目標值m，也會造成一定的損失，但其必然小于A,Cont,設(shè)產(chǎn)品特性值為y的產(chǎn)品，出廠給需方帶來的損失為L（y） 當(dāng)y = m時，

12、L（m）= 0取極小值，其導(dǎo)數(shù)必然也為0 當(dāng)y偏離m達到時， L（m-）= L（m+）= D 圍繞目標值m，對損失函數(shù)L（y）做Taylor展開，并只取到二次項，可得以下近似表達式 可見，損失函數(shù)的第一項即為第三項，即質(zhì)量損失函數(shù)在y = m附近可近似地表達為一條二次曲線,Cont,質(zhì)量損失函數(shù)為,Cont,質(zhì)量損失函數(shù)的確定，關(guān)鍵在于m與k k可由功能界限0及其相應(yīng)損失D0確定 功能界限是指產(chǎn)品失去功能時偏離中心的界限值，即當(dāng) y - m 0 時，產(chǎn)品可正常使用，發(fā)揮功能；反之，當(dāng) y - m 0 時，產(chǎn)品則完全喪失功能，不能使用 即質(zhì)量損失函數(shù)為,Cont,某一電視機的電源電壓，輸入交流電

13、壓220V，輸出直流電壓y的目標值m0為110V。用戶在使用過程中，其穩(wěn)壓電源的直流輸出電壓y與目標值m0之間產(chǎn)生10V的偏差，即喪失功能。此時，用戶因維修、廢棄、重構(gòu)等而蒙受損失，假設(shè)平均損失40元。請估計損失函數(shù)。,Cont,若生產(chǎn)出某一電視，其輸出電壓105V，偏離5V。試問該產(chǎn)品能否作為合格品出廠呢？并估計該產(chǎn)品將給用戶帶來的損失。 若廠家能支付2元對產(chǎn)品進行返修，使輸出電壓達到110V，此時輸出電壓的波動損失則為0 反之，若廠家只考慮自身利益，為節(jié)省2元返修費，則講給用戶和社會增加（10 - 2）= 8元的損失,參數(shù)設(shè)計的基本思想（3）,產(chǎn)品質(zhì)量特性y是隨機變量，由于事先無法得知確切

14、值，故損失函數(shù)L（y）也是隨機變量 對隨機變量最好的評定方式為其數(shù)學(xué)期望（平均值），即EL = E L（y）為平均損失，該值越小，則產(chǎn)品質(zhì)量越高 可見，要減少平均損失，需盡量減少方差和絕對偏差；減少后者相較前者容易一些,Cont,參數(shù)設(shè)計的兩步法 穩(wěn)健設(shè)計減少y的波動，使產(chǎn)品質(zhì)量較為穩(wěn)定 靈敏度設(shè)計尋找“調(diào)節(jié)因子”，減少偏差，使y的均值盡量向目標值靠近,參數(shù)設(shè)計的基本思想（4-1）,進行穩(wěn)健設(shè)計時，需要明確兩個指標，一是產(chǎn)品的質(zhì)量特性y，二是表征y波動的指標（信噪比） 望目特性 具備望目特性的y應(yīng)該是連續(xù)測定的且不為負的任意值，其目標值m是非零的有限值，且y越接近m越好 信噪比SN 本質(zhì)為變異

15、系數(shù)平方的倒數(shù)，SN越大，則質(zhì)量波動越小 記y1，y2，yn是y的n個觀察，則SN估計值為（dB分貝）,參數(shù)設(shè)計的基本思想（4-2）,望小特性 具備望小特性的y應(yīng)該是連續(xù)測定的且不為負的任意值，同時y應(yīng)越接近理想值0越好 由于望小特性沒有目標值，故其參數(shù)設(shè)計只有穩(wěn)健設(shè)計 此時，使用m = 0的平方損失函數(shù) 其平均損失為 其信噪比為 可見，信噪比越大，則等價于平均損失越小,參數(shù)設(shè)計的基本思想（4-3）,望大特性 具備望大特性的y應(yīng)該是連續(xù)測定的且不為負的任意值，同時y應(yīng)取值越大越好 同樣望大特性沒有目標值，故其參數(shù)設(shè)計只有穩(wěn)健設(shè)計 若y為望大特性，則（1/y）為望小特性，利用前頁望小特性的公式，

16、也可估計出望大特性的信噪比,穩(wěn)健設(shè)計,目的旨在減少y的波動，尋找使得質(zhì)量波動達到最小的可控因素的水平組合 舉例說明穩(wěn)健設(shè)計與分析的全過程 欲設(shè)計一個電源電路，能把交流110V變?yōu)橹绷?15V。工程中已有成熟的電路可以使用，但其中有兩個電子元器件的參數(shù)需要確定 因素A（電阻） 因素B（電流放大倍數(shù)） 兩因素各選5個水平，并安排11次試驗,Cont,三組參數(shù)搭配可使輸出為直流115V ：A=200 B=800 ：A=250 B=500 ：A=350 B=260 需要選擇一組參數(shù)搭配，使得偏差和波動均較??！,Cont,尋找波動源 由專業(yè)經(jīng)驗與實際情況判斷 A 電阻的標稱值有10%的波動 B 電流放大

17、倍數(shù)的標稱值有50%的波動 直觀分析 當(dāng)A固定在250時，y是B的線性函數(shù) 當(dāng)B固定在500時，y是A的非線性函數(shù),Cont,從上圖可見，要使y波動小，電阻A應(yīng)選大一些為好 綜上，選擇參數(shù)搭配（A=350 B=260） 注意事項 一般地，線性效應(yīng)有助于將y的均值平移至目標值m；而非線性效應(yīng)則有助于降低噪聲的干擾 面臨包含多個因素的復(fù)雜試驗中，噪聲因素也多種多樣時，直觀分析是非常困難的,Cont,一般方案借助正交設(shè)計 以參數(shù)搭配（A=200 B=800）為例，在噪聲因素A和B的干擾下，取其最大干擾，形成噪聲因素水平表,Cont,由y與B之間的線性效應(yīng)，得出y = yA+（B-500）/20 以2

18、號試驗為例 將A=180，yA = 89，B=800帶入上式 ，得y = 104 估計y的均值，方差的無偏估計和信噪比SN,Cont,類似地，對參數(shù)搭配和同樣可以獲得相應(yīng)的SN,穩(wěn)健設(shè)計的基本步驟,明確參數(shù)設(shè)計的問題 確定影響質(zhì)量特性y的因素 可控因素 3-8個為宜，各因素的設(shè)定為2-3個水平為宜 噪聲因素 外部噪聲、內(nèi)部噪聲、產(chǎn)品間噪聲 內(nèi)外表設(shè)計 可控因素在內(nèi)，噪聲因素在外 內(nèi)表中的每個試驗點各對應(yīng)一張外表 進行試驗，獲得每個試驗結(jié)果yij 內(nèi)n外m，共nm個結(jié)果 統(tǒng)計分析 一個外表數(shù)據(jù)計算一個對應(yīng)的信噪比 將各個信噪比放入內(nèi)表，按正交設(shè)計的分析處理 驗證試驗,Cont,靈敏度設(shè)計,目的 在不增加y波動的前提下，設(shè)法將均值調(diào)節(jié)到目標值 適用 望目特性的y 設(shè)隨機變量y的均值期望為，則其平方則為y的靈敏度 其無偏估計為 參照信噪比的做法，則y的靈敏度為,Cont,靈敏