數值逼近
數值逼近
數值逼近,泛指數學計算問題的近似解法。
目錄
泛指數學計算問題的近似解法。狹義的理解則專指對函數的逼近,即對於給定的較廣泛的函數類F中的函數,從較小的子類H中尋求在某種意義下ƒ的一個近似函數h(x),以便於計算和處理。∏.Л.切比雪夫和K.(T.W.)外爾斯特拉斯曾於19世紀中後期做了奠基性工作。函數逼近的主要內容有,對於某些特定的被逼近函數類F與逼近函數類H,討論逼近的可能性,最佳逼近的存在性、特徵、惟一性、誤差估計以及演演算法等。它是現代數值分析的基本組成部分,除自身具有獨立學科分支的意義外,還可用於構造數值積分、求函數零點、解微分方程和積分方程的近似方法。
設被逼近函數,逼近函數類記作H嶅,定義兩個函數ƒ與g之間的距離為 式中為取定的權函數。當時,通常取ω(x)呏1,此時⑴簡化為
這種度量下的逼近稱為一致逼近;另一種重要情形是p=2的度量,稱為均方逼近或平方逼近。
最佳逼近
若滿足
則稱h為距離度量(1)的意義下ƒ 在H 中的最佳逼近。對於和∞,相應的h分別稱為ƒ 在H 中的最佳平方逼近和最佳一致逼近,后一種情形又稱切比雪夫逼近或極小極大逼近,它是由切比雪夫在1854年首先開始研究的。
多項式逼近
指H 取作多項式類的情形。關於用多項式一致逼近連續函數到任意精度的可能性問題,外爾斯特拉斯於 1885 年以定理的形式給出肯定的答案:若,則對於任何,都存在代數多項式p(x),使。關於用三角多項式一致逼近周期連續函數到任意精度的可能性問題,他也給出平行的結果。該定理本身及其各種不同的證明和推廣對逼近論的研究和發展有重要的影響。
最佳一致多項式逼近
取為次數不大於n的多項式集合。若滿足
(2)
則稱p 為ƒ 的次數不大於 n的最佳一致多項式逼近,稱En(ƒ)為極小極大偏差。這樣的多項式p是存在和惟一的。它的特徵可表述為如下的交錯定理:設,則p(x)滿足(2)當且僅當上存在一組分點(稱為偏差點組或交錯點組)
,
使得
成立。這個定理除了理論上的意義外,還是構造最佳一致多項式逼近演演算法的依據。
列梅茲演演算法 以交錯定理為基礎的尋求最佳一致多項式逼近的一種典型方法。主要步驟如下:
①選取初始偏差點組
通常取作上n 1次切比雪夫多項式的極值點,即
② 解關於和en的線性代數方程組
求得逼近多項式和在偏差點組上均衡了的偏差量|en|。
③在上求一點x使
④ 若,則p(x)已是最佳;否則以x替換某個與之鄰近的xj使同號,然後回到②重新開始。當且異號時,則保留x1而去掉xn 2並將諸偏差點按序重新編號。對於出現在右端的這種情形亦做類似處理。經過若干次循環即可得到足夠精確的結果。逼近階 指極小極大偏差En(ƒ)當n增長時的下降速度,它與被逼近函數的光滑性質有著內在的聯繫。D.傑克森於1911年做了開創性研究。以
表示函數g(x)的連續模,LipM(α)表示所有滿足條件
的函數g(x)的集合,E奱(ƒ)表示周期為2π的連續函數ƒ用次數不大於n 的最佳一致三角多項式逼近的極小極大偏差。傑克森的基本結果可表述如下:
若,則
,
當時,有
;
若,則
,
當時,有
,
式中A為絕對常數。
從相反方向的研究,即從序列{E奱(ƒ)}或{En(ƒ)}的遞減速度來推斷ƒ的光滑性質,是С.Η.伯恩斯坦1912年的工作。他的結果幾乎就是傑克森定理的逆定理,只是的情形稍有差別。
平方逼近 採用時的距離度量(1),被逼近函數ƒ可以屬於比連續函數類更廣的函數類,即所有使存在的ƒ 之集合,記作L嵣。定義L嵣中兩個函數ƒ 與g 的內積為
。
當時,則稱ƒ 與g 正交。設φ0,φ1,…,φn為L嵣中一組線性無關的函數,它們的所有線性組合所構成的函數集合記作φ,則每個ƒ∈L嵣在φ 中的最佳平方逼近
(3)
存在且惟一,其特徵為與每個都正交,諸係數сk由 (4)確定,其中
為格拉姆行列式。當諸φj兩兩正交時,сk簡化為
,
並稱為ƒ 的廣義傅里葉係數,相應的(3)稱為ƒ 按正交函數系{φ0,φ1,…,φn}的級數展開式。一種重要的特殊情形是切比雪夫級數展開式。
切比雪夫級數展開式 即,的情形。係數具體表示為
H.L.勒貝格曾證明:若,則
此即表明,對於連續函數類而言,切比雪夫級數展開式的部分和是最佳一致多項式逼近的很好近似,並且由於它簡便易行,在數值逼近的實踐中廣為採用。
有理逼近 指逼近函數類取作 的情形,其中p(x),Q(x)為多項式,嬠p,嬠Q表示它們的次數。當時,ƒ在中的最佳一致逼近R*存在且惟一,其特徵(充要條件)為在上有一組點 使得誤差R*-ƒ在這些點上達到其最大絕對值且符號正負交替變化,即
實踐中有各種各樣尋求最佳一致有理逼近的數值方法,其中效果令人滿意者有基於上述特徵的列梅茲演演算法,還有加權極小極大演演算法,又稱勞勃演演算法。前者類同於尋求最佳一致多項式逼近的相應演演算法,只是將那裡的多項式p(x)代之以有理函數R(x),初始偏差點組中包含n m 2個點,以及在偏差點組上等化偏差時要解一個非線性方程組。後者的主要步驟可表述為:
① 選取初始值。
② 確定和,使
為極小。為避免得到零解,可固定pk或Qk中的係數之一為非零常數,例如取b0呏1。經過若干次迭代,求得的即為足夠精確的結果。
帕德逼近 一種特殊類型的有理逼近,被逼近的函數由形式冪級數定義。設
若存在多項式
滿足條件以及pn/Qm不可約且規範條件為,則稱外)的(n,m)階帕德逼近,並簡記為存在,則必惟一,其係數pj,qj滿足方程組
這個方程組稱為帕德方程組,其中若當(5)非奇異時,雅可比給出的顯式解
式中若出現求和號的下指標超過上指標,則規定該和數為零。這個顯式解儘管(當m 較大時)在計算上並不便利,但它對於研究帕德逼近的代數性質有重要的作用。帕德將諸[n/m]依自然順序排列成形如
的表格,被稱為帕德表,並研究了它的結構性質。表中相鄰近的一些元素之間還存在著若干內在聯繫,可用於構造帕德逼近的各種遞推演演算法。實踐表明,帕德表中主對角線上及其兩側的元素,即在n m相等的諸元素中通常有更好的逼近效果。帕德逼近是函數的泰勒級數的自然引伸,已證明它是最佳局部有理切比雪夫逼近,它在數值分析的一些領域中以及物理學和化學的某些計算問題中有著各種應用。