正交多項式
![正交多項式回歸](/img/7/f80/wZwpmL1AjNygDO0QzN2YjN1UTM1QDN5MjM5ADMwAjMwUzL0czL1AzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
對於定義在區間上[a,b]的一個函式系 ,如果其中任何兩個函式在此區間上的積分為零,而他們之中每個函式自乘的積分不等於零,即
![正交多項式回歸](/img/c/914/wZwpmL3UDNwIjN1gjN2YjN1UTM1QDN5MjM5ADMwAjMwUzL4YzLygzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/3/9be/wZwpmL0ADOxUTOzITN2YjN1UTM1QDN5MjM5ADMwAjMwUzLyUzL1gzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
則稱此函式係為在此區間上關於權函式 的正交函式系。
![正交多項式回歸](/img/6/f7f/wZwpmL2ETO4EzNxkDN2YjN1UTM1QDN5MjM5ADMwAjMwUzL5QzL4YzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
當 時稱之為規範的正交函式系;當此函式系中每一個函式均為多項式時稱之為正交多項式(系)。
數學模型
設變數y和x的n組觀測數據服從以下k次多項式
![正交多項式回歸](/img/a/209/wZwpmLzITNzMTN4ETN2YjN1UTM1QDN5MjM5ADMwAjMwUzLxUzLzQzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
令
![正交多項式回歸](/img/5/63c/wZwpmL4gzMxUzN5cjN2YjN1UTM1QDN5MjM5ADMwAjMwUzL3YzL1MzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/2/d67/wZwpmL3UTNwMDN1QzN2YjN1UTM1QDN5MjM5ADMwAjMwUzL0czLxMzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/f/3f4/wZwpmL4QDM2IDN3MjMzEzM1UTM1QDN5MjM5ADMwAjMwUzLzIzL3QzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
分別是x的一次,二次,...,k次多項式, 是一些適當選擇的常數,如何選擇將在下面討論(i=1,2,…,n)。將(2)式帶入(1)式,則有
![正交多項式回歸](/img/6/76f/wZwpmL3UzM5EDMzQzN2YjN1UTM1QDN5MjM5ADMwAjMwUzL0czL4EzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/3/118/wZwpmLwIDN2ADN3EzN2YjN1UTM1QDN5MjM5ADMwAjMwUzLxczL1gzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/8/8aa/wZwpmL4ADNykTM2cTN2YjN1UTM1QDN5MjM5ADMwAjMwUzL3UzLxAzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
比較(3)和(1)式可知,二者係數間存在簡單的函式關係,只要求出 ,就可以求出 。
![正交多項式回歸](/img/2/d67/wZwpmL3UTNwMDN1QzN2YjN1UTM1QDN5MjM5ADMwAjMwUzL0czLxMzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
若把 看作新的自變數,則(3)式就成為一個k元線性模型,其結構矩陣為
![正交多項式回歸](/img/c/96c/wZwpmLzcTM3EjN2AzN2YjN1UTM1QDN5MjM5ADMwAjMwUzLwczL2MzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
正規方程為
![正交多項式回歸](/img/b/0cf/wZwpmLygzN5EDOzUzN2YjN1UTM1QDN5MjM5ADMwAjMwUzL1czL3QzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/3/618/wZwpmL0ITN1UDOwgTN2YjN1UTM1QDN5MjM5ADMwAjMwUzL4UzLyQzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
其中
![正交多項式回歸](/img/f/96d/wZwpmL0gDM0MDMzEzN2YjN1UTM1QDN5MjM5ADMwAjMwUzLxczL4UzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/3/9a0/wZwpmLwQjNzUDM4kDN2YjN1UTM1QDN5MjM5ADMwAjMwUzL5QzL1YzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/1/289/wZwpmLwczN1czNzgzN2YjN1UTM1QDN5MjM5ADMwAjMwUzL4czL0YzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/3/31f/wZwpmL3QTMxYDO5gTN2YjN1UTM1QDN5MjM5ADMwAjMwUzL4UzL4MzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
在多項式回歸中遇到的困難是解正規方程係數矩陣的工作量太大,可以使其對角線上的元素不為零,而其餘元素均為零,從而簡化計算,而且同時消去了係數間的相關性。
![正交多項式回歸](/img/2/d67/wZwpmL3UTNwMDN1QzN2YjN1UTM1QDN5MjM5ADMwAjMwUzL0czLxMzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/5/f17/wZwpmL3YDNyQzN1cDN2YjN1UTM1QDN5MjM5ADMwAjMwUzL3QzL1czLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
對於 我們可以通過選擇係數 使得
![正交多項式回歸](/img/0/9d8/wZwpmL1gDNxQTN2QzN2YjN1UTM1QDN5MjM5ADMwAjMwUzL0czL3EzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/6/8be/wZwpmLyYjN5kzMxcTN2YjN1UTM1QDN5MjM5ADMwAjMwUzL3UzL2UzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
從而使
![正交多項式回歸](/img/2/b00/wZwpmLyQTM4YTO4cjN2YjN1UTM1QDN5MjM5ADMwAjMwUzL3YzL0IzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/f/fde/wZwpmL4UzNwYTO0MjN2YjN1UTM1QDN5MjM5ADMwAjMwUzLzYzL0MzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/0/920/wZwpmLxgjN1kTNzUzN2YjN1UTM1QDN5MjM5ADMwAjMwUzL1czLyUzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/c/6ac/wZwpmLxcjN1YDOyAjN2YjN1UTM1QDN5MjM5ADMwAjMwUzLwYzLyYzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
則正規方程組為
![正交多項式回歸](/img/d/4db/wZwpmLzMTOzkjMxcDN2YjN1UTM1QDN5MjM5ADMwAjMwUzL3QzLxAzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
回歸係數為
![正交多項式回歸](/img/1/1ef/wZwpmL1MTOyYzN2AzN2YjN1UTM1QDN5MjM5ADMwAjMwUzLwczLwMzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/d/e16/wZwpmL4MDMwkjN3ETN2YjN1UTM1QDN5MjM5ADMwAjMwUzLxUzLyczLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
滿足(7)和(8)式的多項式組…我們稱之為 正交多項式。
在正交多項式回歸中自變數的選擇是等間隔的,設間隔為h,x=a, 則
![正交多項式回歸](/img/3/532/wZwpmLzEDM3gTM5kTN2YjN1UTM1QDN5MjM5ADMwAjMwUzL5UzL4gzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
若令
![正交多項式回歸](/img/2/5f8/wZwpmL3czM4YTO5MjN2YjN1UTM1QDN5MjM5ADMwAjMwUzLzYzLwEzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
則
![正交多項式回歸](/img/1/ae0/wZwpmLxgzN5MjN5QTN2YjN1UTM1QDN5MjM5ADMwAjMwUzL0UzLyIzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/7/929/wZwpmL2AjN2gjM4UjN2YjN1UTM1QDN5MjM5ADMwAjMwUzL1YzLwQzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
是1至n的正整數。
![正交多項式回歸](/img/4/b5b/wZwpmL0cjM2QzMzEzN2YjN1UTM1QDN5MjM5ADMwAjMwUzLxczL0gzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
![正交多項式回歸](/img/6/236/wZwpmL4cDOzIzMyYzN2YjN1UTM1QDN5MjM5ADMwAjMwUzL2czLxUzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
為簡化問題,用 代替x作為自變數。在條件許可時,取自變數x=1,x=2,…,x=n。當x=1,x=2,…,x=n時有 ,這時驗證以下多項式是正交的,即
![正交多項式回歸](/img/3/682/wZwpmLzUTNwEjN0cjN2YjN1UTM1QDN5MjM5ADMwAjMwUzL3YzL0YzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/d/2cb/wZwpmL2ETO4czNwYjMzEzM1UTM1QDN5MjM5ADMwAjMwUzL2IzLxIzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
顯然,當x取正整數時, 不一定是整數,為了克服這給計算上帶來的困難,取
![正交多項式回歸](/img/5/d3a/wZwpmL0UzNyUDN3QjN2YjN1UTM1QDN5MjM5ADMwAjMwUzL0YzL2EzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/8/754/wZwpmL0ADO3gTNwAzNwMzM1UTM1QDN5MjM5ADMwAjMwUzLwczLzYzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/0/10f/wZwpmL4MDNyMTO5ATN2YjN1UTM1QDN5MjM5ADMwAjMwUzLwUzL2MzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/f/060/wZwpmL3EDO5gTN0kjNwYjN1UTM1QDN5MjM5ADMwAjMwUzL5YzL1UzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
![正交多項式回歸](/img/0/10f/wZwpmL4MDNyMTO5ATN2YjN1UTM1QDN5MjM5ADMwAjMwUzLwUzL2MzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/0/10f/wZwpmL4MDNyMTO5ATN2YjN1UTM1QDN5MjM5ADMwAjMwUzLwUzL2MzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
使x取正整數時 是整數。可以驗證用正交多項式 代替 所求得的回歸方程與用正交多項式 所求得的回歸方程是完全一樣的。
![正交多項式回歸](/img/f/060/wZwpmL3EDO5gTN0kjNwYjN1UTM1QDN5MjM5ADMwAjMwUzL5YzL1UzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
對於正交多項式 有
![正交多項式回歸](/img/c/18b/wZwpmLxQzNxQTMzcDN2YjN1UTM1QDN5MjM5ADMwAjMwUzL3QzL0YzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
![正交多項式回歸](/img/1/8fb/wZwpmLycTNzETMyAjN2YjN1UTM1QDN5MjM5ADMwAjMwUzLwYzL3AzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
![正交多項式回歸](/img/0/109/wZwpmLxczMzczN1cDM0YTN1UTM1QDN5MjM5ADMwAjMwUzL3AzL4YzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
不同的n相對應的 在 時的值以及S值都已製成正交多項式表,根據正交多項式表,可以計算出回歸方程的係數。令
![正交多項式回歸](/img/d/48f/wZwpmL1UzNzYDNxcDN2YjN1UTM1QDN5MjM5ADMwAjMwUzL3QzLzUzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
則
![正交多項式回歸](/img/0/fa3/wZwpmL0gjN4gDM0cDN2YjN1UTM1QDN5MjM5ADMwAjMwUzL3QzL2UzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
回歸方程為
![正交多項式回歸](/img/6/270/wZwpmLyADNxMTM4MjN2YjN1UTM1QDN5MjM5ADMwAjMwUzLzYzLwUzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
由於正交多項式回歸係數之間不存在相關性,因此某一項如果不顯著,只要將它剔除即可,而不必對整個回歸方程重新計算。
回歸係數的顯著性檢驗
正交多項式回歸方程與回歸係數的顯著性檢驗可利用正交多項式的性質按表1進行。經檢驗不顯著的高次項可以剔除,將其效應併入殘差平方和,自由度也同時併入,如果對回歸方程精度不滿意,可以增加高次項,而已經計算出的結果不必重算。
表1 正交多項式回歸方差分析表
方差來源 | 平方和 | 自由度 | 均方 | F | 顯著性 |
回歸 | ![]() | k | ![]() | ![]() | |
一次項 ![]() | ![]() | 1 | ![]() | ![]() | |
二次項 ![]() | ![]() | 1 | ![]() | ![]() | |
... | ... | ... | ... | ... | ... |
k次項 ![]() | ![]() | 1 | ![]() | ![]() | |
殘差 | ![]() | n-k-1 | ![]() | ||
總計 | ![]() | n-1 |
程式框圖
1.數學模型:
![正交多項式回歸](/img/f/a60/wZwpmL4UjMwIjNzITN2YjN1UTM1QDN5MjM5ADMwAjMwUzLyUzL3IzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
2.變數及數組說明:
J-正確讀入數據的控制變數;
N-試驗組數;
M-所取正交多項式項數;
X(I)-存自變數數值;
Y(I)-存因變數數值;
Z(I)-存Y(I)的平方項;
![正交多項式回歸](/img/e/cbe/wZwpmLzIDO2MDN5EzNwMzM1UTM1QDN5MjM5ADMwAjMwUzLxczLwAzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
E(I,1)-存在正交多項式一次項 ;
![正交多項式回歸](/img/9/d47/wZwpmL1QzN2kzM3MzMxMzM1UTM1QDN5MjM5ADMwAjMwUzLzMzL2IzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
E(I,2)-存在正交多項式二次項 ;
![正交多項式回歸](/img/1/4eb/wZwpmLyUDO1QDMzcTN2YjN1UTM1QDN5MjM5ADMwAjMwUzL3UzL3czLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
E(I,3)-存在正交多項式三次項 (其中I=1,…N);
S(J)-結構矩陣逆矩陣元素 J=1,2,3;
B(J)-常數項矩陣B J=1,2,3;
D(J)-回歸係數 J=0,1,2,3;
Q(J)-偏回歸平方和 J=0,1,2,3;
S-剩餘平方和;
S-標準離差;
S-總平方和;
F(J)-F檢驗值。
3.程式框圖:
見如下圖1
![圖1](/img/9/184/wZwpmL1UjMzkTMxgDO2YjN1UTM1QDN5MjM5ADMwAjMwUzL4gzL1IzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)