最佳平方逼近多項(xiàng)式課件

最佳平方逼近多項(xiàng)式最佳平方逼近多項(xiàng)式本節(jié)內(nèi)容1.內(nèi)積空間2.兩類(lèi)特殊的函數(shù)族3.函數(shù)的最佳平方逼近4.舉例5.MATLAB程序?qū)崿F(xiàn)§5.2最佳平方逼近多項(xiàng)式本節(jié)內(nèi)容§5.2最佳平方逼近多項(xiàng)式1.內(nèi)積空間權(quán)函數(shù)：考慮到在區(qū)間[a,b]上各點(diǎn)的函數(shù)值比重不同,常引進(jìn)加權(quán)形式的定義，設(shè)在區(qū)間[a,b]上的非負(fù)函數(shù)滿(mǎn)足條件：1）存在；2）對(duì)非負(fù)的連續(xù)函數(shù)，若則在[a,b]上，，即不恒為0。就稱(chēng)為[a,b]上的權(quán)函數(shù)。它的物理意義可以解釋為密度函數(shù)。1.內(nèi)積空間權(quán)函數(shù)：考慮到在區(qū)間[a,b]上各點(diǎn)的函數(shù)1.內(nèi)積空間內(nèi)積：設(shè)是[a,b]上的權(quán)函數(shù)，則稱(chēng)積分

為函數(shù)與在[a,b]上的內(nèi)積,有下列性質(zhì)：1）2）為常數(shù);3）4）當(dāng)且僅當(dāng)時(shí)，1.內(nèi)積空間內(nèi)積：設(shè)是[1.內(nèi)積空間內(nèi)積空間：滿(mǎn)足內(nèi)積定義的函數(shù)空間稱(chēng)為內(nèi)積空間。如在連續(xù)函數(shù)空間上定義了內(nèi)積就形成了一個(gè)內(nèi)積空間。向量的模：在n維歐氏空間中，內(nèi)積就是兩向量的數(shù)量積，即向量的模（范數(shù))的定義為：1.內(nèi)積空間內(nèi)積空間：滿(mǎn)足內(nèi)積定義的函數(shù)空間稱(chēng)為內(nèi)積空間。如1.內(nèi)積空間歐式范數(shù)：若,則量稱(chēng)為的歐式范數(shù)。對(duì)任何,有以下結(jié)論：（1），又稱(chēng)柯西-施瓦茨不等式；（2），又稱(chēng)三角不等式；（3），又稱(chēng)平行四邊形定律。1.內(nèi)積空間歐式范數(shù)：若,則量稱(chēng)為的2.兩類(lèi)特殊的函數(shù)族正交：若為[a,b]上的權(quán)函數(shù)且滿(mǎn)足則稱(chēng)與在[a,b]上帶權(quán)正交。正交函數(shù)族：若函數(shù)族滿(mǎn)足關(guān)系則稱(chēng)是[a,b]上帶權(quán)的正交函數(shù)族；若，則稱(chēng)為標(biāo)準(zhǔn)正交函數(shù)族。2.兩類(lèi)特殊的函數(shù)族正交：若2.兩類(lèi)特殊的函數(shù)族可以證明，三角函數(shù)族滿(mǎn)足上述條件，是在上的正交函數(shù)族。線(xiàn)性無(wú)關(guān)：若函數(shù)在區(qū)間[a,b]上連續(xù)，如果當(dāng)且僅當(dāng)時(shí)成立，則稱(chēng)在[a,b]上是線(xiàn)性無(wú)關(guān)的。2.兩類(lèi)特殊的函數(shù)族可以證明，三角函數(shù)族線(xiàn)性無(wú)關(guān)：若函2.兩類(lèi)特殊的函數(shù)族線(xiàn)性無(wú)關(guān)函數(shù)族：若函數(shù)族中的任何有限個(gè)線(xiàn)性無(wú)關(guān)，則稱(chēng)為線(xiàn)性無(wú)關(guān)函數(shù)族。充要條件：

在[a,b]上線(xiàn)性無(wú)關(guān)的充要條件是它的Gramer行列式，其中2.兩類(lèi)特殊的函數(shù)族線(xiàn)性無(wú)關(guān)函數(shù)族：若函數(shù)族充要條件：3.函數(shù)的最佳平方逼近最佳平方逼近函數(shù)：對(duì)于及中的一個(gè)子集若存在使下式成立:則稱(chēng)是在子集中的最佳平方逼近函數(shù)，其中是一組線(xiàn)性無(wú)關(guān)函數(shù)族，函數(shù)3.函數(shù)的最佳平方逼近最佳平方逼近函數(shù)：對(duì)于3.函數(shù)的最佳平方逼近對(duì)函數(shù)s*(x)的求解：等價(jià)于求以下多元函數(shù)

的最小值。令則引入內(nèi)積定義，可得即3.函數(shù)的最佳平方逼近對(duì)函數(shù)s*(x)的求解：等價(jià)于求以下多3.函數(shù)的最佳平方逼近上式是關(guān)于的線(xiàn)性方程組，稱(chēng)為法方程。用矩陣形式可表示為簡(jiǎn)記為。其中3.函數(shù)的最佳平方逼近上式是關(guān)于的線(xiàn)3.函數(shù)的最佳平方逼近由于線(xiàn)性無(wú)關(guān)，故其系數(shù)矩陣H的行列式非奇異，即，該法方程有唯一解為

則最佳平方逼近函數(shù)為令，則平方誤差3.函數(shù)的最佳平方逼近由于線(xiàn)性無(wú)關(guān)，3.函數(shù)的最佳平方逼近

特別地，取，，求其最佳平方逼近多項(xiàng)式。此時(shí)，3.函數(shù)的最佳平方逼近特別地，取3.函數(shù)的最佳平方逼近又稱(chēng)為希爾伯特矩陣。則方程

的唯一解即為所求多項(xiàng)式s*(x)的系數(shù)。3.函數(shù)的最佳平方逼近又稱(chēng)為希爾伯特矩陣。則方程4.舉例1.求在區(qū)間[0,1]上的一次最佳平方逼近多項(xiàng)式。解：取由得則由4.舉例1.求在區(qū)間[0,1]上的一次最佳4.舉例得解得：故所求一次最佳平方逼近多項(xiàng)式為：所求最佳一次逼近多項(xiàng)式為：4.舉例得解得：故所求一次最佳平方逼近多項(xiàng)式為：所求最佳一次4.舉例4.舉例4.舉例Matlab求定積分（int函數(shù)）d0=(2*2^(1/2))/5-(6*ellipticF(asin(1/(3/2+(3^(1/2)*i)/2)^(1/2)),-(3/2+(3^(1/2)*i)/2)/(-3/2+(3^(1/2)*i)/2))*(-1/(-3/2+(3^(1/2)*i)/2))^(1/2))/5+(6*(3/2+(3^(1/2)*i)/2)*(2/(3/2+(3^(1/2)*i)/2))^(1/2)*((-1/2+(3^(1/2)*i)/2)/(3/2+(3^(1/2)*i)/2))^(1/2)*((1/2+(3^(1/2)*i)/2)/(-3/2+(3^(1/2)*i)/2))^(1/2)*ellipticF(asin((2/(3/2+(3^(1/2)*i)/2))^(1/2)),-(3/2+(3^(1/2)*i)/2)/(-3/2+(3^(1/2)*i)/2))*(-1/(2*(-1/2+(3^(1/2)*i)/2)*(1/2+(3^(1/2)*i)/2)))^(1/2))/54.舉例Matlab求定積分（int函數(shù)）d0=4.舉例4.舉例4.舉例4.舉例4.舉例二次4.舉例二次4.舉例三次4.舉例三次4.舉例四次4.舉例四次4.舉例4.舉例4.舉例4.舉例4.舉例4.舉例4.舉例4.舉例5.MATLAB編程實(shí)現(xiàn)functionA=ZJPFBJ(f,n,a,b)C=zeros(n+1,n+1);var=findsym(f);f=f/varfori=1:n+1C(1,i)=(power(b,i)-power(a,i))/i;f=f*var;d(i,1)=int(sym(f),var,a,b);endfori=2:n+1C(i,1:n)=C(i-1,2:n+1);f1=power(b,n+i);f2=power(a,n+i);C(i,n+1)=(f1-f2)/(n+i);endA=C\d;A=real(double(A));end5.MATLAB編