主題
Search

格拉姆-施密特正交化

格拉姆-施密特正交化，也稱為格拉姆-施密特過程，是一個程式，它接受一個非正交的線性無關函式集，並構建一個關於任意區間，相對於任意權重函式 w(x) 的正交基。

將格拉姆-施密特過程應用於函式 1, x, x^2, ... 在區間 [-1,1] 上，使用通常的 L^2 內積 得到 勒讓德多項式（直到常數倍數；Reed and Simon 1972, p. 47）。

給定一個原始的線性無關函式集 {u_n}_(n=0)^infty，令 {psi_n}_(n=0)^infty 表示正交化的（但未歸一化的）函式，{phi_n}_(n=0)^infty 表示正交歸一化的函式，並定義

psi_0(x)=u_0(x)
(1)
phi_0(x)=(psi_0(x))/(sqrt(intpsi_0^2(x)w(x)dx)).
(2)

然後取

psi_1(x)=u_1(x)+a_(10)phi_0(x),
(3)

其中我們要求

intpsi_1phi_0wdx=intu_1phi_0wdx+a_(10)intphi_0^2wdx
(4)
=0.
(5)

根據定義，

intphi_0^2wdx=1,
(6)

所以

a_(10)=-intu_1phi_0wdx.
(7)

因此，第一個正交化函式是

psi_1=u_1(x)-[intu_1phi_0wdx]phi_0,
(8)

相應的歸一化函式是

phi_1=(psi_1(x))/(sqrt(intpsi_1^2wdx)).
(9)

透過數學歸納法，可以得出

phi_i(x)=(psi_i(x))/(sqrt(intpsi_i^2wdx)),
(10)

其中

psi_i(x)=u_i+a_(i0)phi_0+a_(i1)phi_1...+a_(i,i-1)phi_(i-1)
(11)

和

a_(ij)=-intu_iphi_jwdx.
(12)

如果函式被歸一化為 N_j 而不是 1，那麼

int_a^b[phi_j(x)]^2wdx=N_j^2
(13)
phi_i(x)=N_i(psi_i(x))/(sqrt(intpsi_i^2wdx))
(14)
a_(ij)=-(intu_iphi_jwdx)/(N_j^2).
(15)

正交多項式 特別容易使用格拉姆-施密特正交化生成。使用符號

<x_i|x_j>=<x_i|w|x_j>
(16)
=int_a^bx_i(x)x_j(x)w(x)dx,
(17)

其中 w(x) 是一個權重函式，並定義前幾個多項式，

p_0(x)=1
(18)
p_1(x)=[x-(<xp_0|p_0>)/(<p_0|p_0>)]p_0.
(19)

根據定義，p_0 和 p_1 是正交多項式，可以從

<p_0|p_1>=<[x-(<xp_0|p_0>)/(<p_0|p_0>)]p_0>
(20)
=<xp_0>-(<xp_0|p_0>)/(<p_0|p_0>)<p_0>
(21)
=<xp_0>-<xp_0>
(22)
=0.
(23)

現在使用遞推關係

p_(i+1)(x)=[x-(<xp_i|p_i>)/(<p_i|p_i>)]p_i-[(<p_i|p_i>)/(<p_(i-1)|p_(i-1)>)]p_(i-1)
(24)

構造所有更高階的多項式。

為了驗證這個過程確實產生了正交多項式，檢查

<p_(i+1)|p_i>=<[x-(xp_i|p_i)/(p_i|p_i)]p_i|p_i>-<(p_i|p_i)/(p_(i-1)|p_(i-1))p_(i-1)|p_i>
(25)
=<xp_i|p_i>-(<xp_i|p_i>)/(<p_i|p_i>)<p_i|p_i>-(<p_i|p_i>)/(<p_(i-1)|p_(i-1)>)<p_(i-1)|p_i>
(26)
=-(<p_i|p_i>)/(<p_(i-1)|p_(i-1)>)<p_(i-1)|p_i>
(27)
=-(<p_i|p_i>)/(<p_(i-1)|p_(i-1)>)[-(<p_(i-1)|p_(j-1)>)/(<p_(j-2)|p_(j-2)>)<p_(j-2)|p_(j-1)>]
(28)
=...
(29)
=(-1)^j(<p_j|p_j>)/(<p_0|p_0>)<p_0|p_1>
(30)
=0,
(31)

因為 <p_0|p_1>=0。因此，所有的多項式 p_i(x) 都是正交的。

許多常見的數學物理正交多項式可以用這種方式生成。不幸的是，該過程在數值上是不穩定的 (Golub and Van Loan 1996)。

參見

格拉姆行列式, 格拉姆不等式, 格基規約, 正交多項式

用 探索

參考文獻

Arfken, G. "Gram-Schmidt Orthogonalization." §9.3 in Mathematical Methods for Physicists, 3rd ed. Orlando, FL: Academic Press, pp. 516-520, 1985.Cohen, H. A Course in Computational Algebraic Number Theory. New York: Springer-Verlag, 1993.Golub, G. H. and Van Loan, C. F. Matrix Computations, 3rd ed. Baltimore, MD: Johns Hopkins, 1996.Pohst, M. and Zassenhaus, H. "Methods from the Geometry of Numbers." Ch. 3 in Algorithmic Algebraic Number Theory. Cambridge, England: Cambridge University Press, 1989.Reed, M. and Simon, B. Methods of Modern Mathematical Physics, Vol. 1: Functional Analysis. New York: Academic Press, pp. 46-47, 1972.Trefethen, L. N. and Bau, D. III Numerical Linear Algebra. Philadelphia, PA: SIAM, 1997.

在 上被引用

格拉姆-施密特正交化

引用此內容為

Weisstein, Eric W. "格拉姆-施密特正交化。" 來自 Web 資源。 https://mathworld.tw/Gram-SchmidtOrthonormalization.html

主題分類