數值計算方法朱建新

㈠ 誰有 《數值計算方法 第三版》高等教育出版社 主編朱建新、李有法 課後答案以及 山西師范大學 的歷年考題

主編朱建新、李有法課後答案以及山西師范大學的歷年考題：

有限元法：有限元方法的基礎是變分原理和加權餘量法，其基本求解思想是把計算域劃分為有限個互不重疊的單元，在每個單元內，選擇一些合適的節點作為求解函數的插值點，將微分方程中的變數改寫成由各變數或其導數的節點值與所選用的插值函數組成的線性表達式。

藉助於變分原理或加權餘量法，將微分方程離散求解。採用不同的權函數和插值函數 形式，便構成不同的有限元方法。

在有限元方法中，把計算域離散剖分為有限個互不重疊且相互連接的單元，在每個單元內選擇基函數，用單元基函數的線形組合來逼近單元中的真解，整個計算域上總體的基函數可以看為由每個單元基函數組成的，則整個計算域內的解可以看作是由所有單元 上的近似解構成。

(1)數值計算方法朱建新擴展閱讀：

構造數值積分公式最通常的方法是用積分區間上的n 次插值多項式代替被積函數，由此導出的求積公式稱為插值型求積公式。特別在節點分布等距的情形稱為牛頓-柯茨公式，例如梯形公式與拋物線公式就是最基本的近似公式。但它們的精度較差。

龍貝格演算法是在區間逐次分半過程中，對梯形公式的近似值進行加權平均獲得准確程度較高的積分近似值的一種方法，它具有公式簡練、計算結果准確、使用方便、穩定性好等優點，因此在等距情形宜採用龍貝格求積公式。

㈡ 求數值計算方法 第三版 李有法 朱建新 課後答案

數值計算方法如下：

1、有限元法：有限元方法的基礎是變分原理和加權餘量法，其基本求解思想是把計算域劃分為有限個互不重疊的單元，在每個單元內，選擇一些合適的節點作為求解函數的插值點，將微分方程中的變數改寫成由各變數或其導數的節點值與所選用的插值函數組成的線性表達式。

藉助於變分原理或加權餘量法，將微分方程離散求解。採用不同的權函數和插值函數 形式，便構成不同的有限元方法。

在有限元方法中，把計算域離散剖分為有限個互不重疊且相互連接的單元，在每個單元內選擇基函數，用單元基函數的線形組合來逼近單元中的真解，整個計算域上總體的基函數可以看為由每個單元基函數組成的，則整個計算域內的解可以看作是由所有單元 上的近似解構成。

根據所採用的權函數和插值函數的不同 ，有限元方法也分為多種計算格式。從權函數的選擇來說，有配置法、矩量法、最小二乘法和伽遼金法，從計算單元網格的形狀來劃分，有三角形網格、四邊形網格和多邊形網格，從插值函數的精度來劃分，又分為線性插值函數和高次插值函數等。不同的組合 同樣構成不同的有限元計算格式。

2、多重網格方法：多重網格方法通過在疏密不同的網格層上進行迭代,以平滑不同頻率的誤差分量。具有收斂速度快,精度高等優點。

多重網格法基本原理微分方程的誤差分量可以分為兩大類，一類是頻率變化較緩慢的低頻分量；另一類是頻率高，擺動快的高頻分量。

一般的迭代方法可以迅速地將擺動誤差衰減，但對那些低頻分量，迭代法的效果不是很顯著。高頻分量和低頻分量是相對的，與網格尺度有關，在細網格上被視為低頻的分量，在粗網格上可能為高頻分量。

多重網格方法作為一種快速計算方法，迭代求解由偏微分方程組離散以後組成的代數方程組，其基本原理在於一定的網格最容易消除波長與網格步長相對應的誤差分量。

該方法採用不同尺度的網格，不同疏密的網格消除不同波長的誤差分量，首先在細網格上採用迭代法，當收斂速度變緩慢時暗示誤差已經光滑，則轉移到較粗的網格上消除與該層網格上相對應的較易消除的那些誤差分量，這樣逐層進行下去直到消除各種誤差分量，再逐層返回到細網格上。

3、有限差分方法：有限差分方法(FDM)是計算機數值模擬最早採用的方法，至今仍被廣泛運用。該方法將求解域劃分為差分網格，用有限個網格節點代替連續的求解域。

有限差分法以Taylor級數展開等方法，把控制方程中的導數用網格節點上的函數值的差商代替進行離散，從而建立以網格節點上的值為未知數的代數方程組。該方法是一種直接將微分問題變為代數問題的近似數值解法，數學概念直觀，表達簡單，是發展較早且比較成熟的數值方法。

對於有限差分格式，從格式的精度來劃分，有一階格式、二階格式和高階格式。從差分的空間形式來考慮，可分為中心格式和逆風格式。考慮時間因子的影響，差分格式還可以分為顯格式、隱格式、顯隱交替格式等。

構造差分的方法有多種形式，目前主要採用的是泰勒級數展開方法。其基本的差分表達式主要有三種形式：

一階向前差分、一階向後差分、一階中心差分和二階中心差分等，其中前兩種格式為一階計算精度，後兩種格式為二階計算精度。通過對時間和空間這幾種不同差分格式的組合，可以組合成不同的差分計算格式。

4、有限體積法：有限體積法（Finite Volume Method）又稱為控制體積法。其基本思路是：將計算區域劃分為一系列不重復的控制體積，並使每個網格點周圍有一個控制體積；將待解的微分方程對每一個控制體積積分，便得出一組離散方程。其中的未知數是網格點上的因變數的數值。

為了求出控制體積的積分，必須假定值在網格點之間的變化規律，即假設值的分段的分布的分布剖面。從積分區域的選取方法看來，有限體積法屬於加權剩餘法中的子區域法；從未知解的近似方法看來，有限體積法屬於採用局部近似的離散方法。簡言之，子區域法屬於有限體積發的基本方法。

有限體積法的基本思路易於理解，並能得出直接的物理解釋。離散方程的物理意義，就是因變數在有限大小的控制體積中的守恆原理，如同微分方程表示因變數在無限小的控 制體積中的守恆原理一樣。

限體積法得出的離散方程，要求因變數的積分守恆對任意一組控制體積都得到滿足，對整個計算區域，自然也得到滿足。這是有限體積法吸引人的優點。有一些離散方法，例如有限差分法，僅當網格極其細密時，離散方程才滿足積分守恆。

而有限體積法即使在粗網格情況下，也顯示出准確的積分守恆。就離散方法而言，有限體積法可視作有限單元法和有限差分法的中間物。有限單元法必須假定值在網格點之間的變化規律（既插值函數），並將其作為近似解。

有限差分法只考慮網格點上的數值而不考慮值在網格點之間如何變化。有限體積法只尋求的結點值 ，這與有限差分法相類似；但有限體積法在尋求控制體積的積分時，必須假定值在網格點之間的分布，這又與有限單元法相類似。

在有限體積法中，插值函數只用於計算控制體積的積分，得出離散方程之後，便可忘掉插值函數；如果需要的話，可以對微分方程 中不同的項採取不同的插值函數。

5、近似求解的誤差估計方法：近似求解的誤差估計方法共有三大類:單元餘量法,通量投射法及外推法。

單元餘量法廣泛地用於以FEM離散的誤差估計之中,它主要是估計精確運算元的餘量,而不是整套控制方程的全局誤差。

這樣就必須假定周圍的單元誤差並不相互耦合,誤差計算採用逐節點演算法進行。單元餘量法的各種不同做法主要來自對單元誤差方程的邊界條件的不同處理辦法。基於此,該方法能夠有效處理局部的殘餘量,並能成功地用於網格優化程序。

通量投射法的基本原理來自一個很簡單的事實:精確求解偏微分方程不可能有不連續的微分,而近似求解卻可以存在微分的不連續,這樣產生的誤差即來自微分本身,即誤差為系統的光滑求解與不光滑求解之差。該方法與單元餘量法一樣,對節點誤差採用能量范數,故也能成功地用於網格優化程序。

單元餘量法及通量投射法都局限於局部的誤差計算(採用能量范數),誤差方程的全局特性沒有考慮。另外計算的可行性(指誤差估計方程的計算時間應小於近似求解計算時間)不能在這兩種方法中體現,因為獲得的誤差方程數量,階數與流場控制方程相同。

外推是指採用後向數值誤差估計思想由精確解推出近似解的誤差值。各類文獻中較多地採用Richardson外推方法來估計截斷誤差。無論是低階還是高階格式,隨著網格的加密數值計算結果都會趨近於准確解。但由於計算機內存與計算時間的限制,實際上不能採用這種網格無限加密的辦法。

6、多尺度計算方法：近年來發展的多尺度計算方法包括均勻化方法、非均勻化多尺度方法、以及小波數值均勻化方法、多尺度有限體積法、多尺度有限元法等。

該方法通過對單胞問題的求解，把細觀尺度的信息映射到宏觀尺度上，從而推導出宏觀尺度上的均勻化等式，即可在宏觀尺度上求解原問題。均勻化方法在很多科學和工程應用中取得了巨大成功，但這種方法建立在系數細觀結構周期性假設的基礎上，因此應用范圍受到了很大限制。

鄂維南等提出的非均勻化多尺度方法，是構造多尺度計算方法的一般框架。該方法有兩個重要的組成部分：基於宏觀變數的整體宏觀格式和由微觀模型來估計缺少的宏觀數據，多尺度問題的解通過這兩部分共同得到。

該方法基於多分辨分析，在細尺度上建立原方程的離散運算元，然後對離散運算元進行小波變換，得到了大尺度上的數值均勻化運算元。此方法在大尺度上解方程，大大地減小了計算時間。

該法在宏觀尺度上進行網格剖分，然後通過在每個單元里求解細觀尺度的方程(構造線性或者振盪的邊界條件)來獲得基函數。從而把細觀尺度的信息反應到有限元法的基函數里，使宏觀尺度的解包含了細觀尺度的信息。但多尺度有限元方法在構造基函數時需要較大的計算量。

藉助於變分原理或加權餘量法，將微分方程離散求解。採用不同的權函數和插值函數 形式，便構成不同的有限元方法。

在有限元方法中，把計算域離散剖分為有限個互不重疊且相互連接的單元，在每個單元內選擇基函數，用單元基函數的線形組合來逼近單元中的真解，整個計算域上總體的基函數可以看為由每個單元基函數組成的，則整個計算域內的解可以看作是由所有單元 上的近似解構成。

根據所採用的權函數和插值函數的不同 ，有限元方法也分為多種計算格式。從權函數的選擇來說，有配置法、矩量法、最小二乘法和伽遼金法，從計算單元網格的形狀來劃分，有三角形網格、四邊形網格和多邊形網格，從插值函數的精度來劃分，又分為線性插值函數和高次插值函數等。不同的組合 同樣構成不同的有限元計算格式。

2、多重網格方法：多重網格方法通過在疏密不同的網格層上進行迭代,以平滑不同頻率的誤差分量。具有收斂速度快,精度高等優點。

多重網格法基本原理微分方程的誤差分量可以分為兩大類，一類是頻率變化較緩慢的低頻分量；另一類是頻率高，擺動快的高頻分量。

一般的迭代方法可以迅速地將擺動誤差衰減，但對那些低頻分量，迭代法的效果不是很顯著。高頻分量和低頻分量是相對的，與網格尺度有關，在細網格上被視為低頻的分量，在粗網格上可能為高頻分量。

多重網格方法作為一種快速計算方法，迭代求解由偏微分方程組離散以後組成的代數方程組，其基本原理在於一定的網格最容易消除波長與網格步長相對應的誤差分量。

該方法採用不同尺度的網格，不同疏密的網格消除不同波長的誤差分量，首先在細網格上採用迭代法，當收斂速度變緩慢時暗示誤差已經光滑，則轉移到較粗的網格上消除與該層網格上相對應的較易消除的那些誤差分量，這樣逐層進行下去直到消除各種誤差分量，再逐層返回到細網格上。

3、有限差分方法：有限差分方法(FDM)是計算機數值模擬最早採用的方法，至今仍被廣泛運用。該方法將求解域劃分為差分網格，用有限個網格節點代替連續的求解域。

有限差分法以Taylor級數展開等方法，把控制方程中的導數用網格節點上的函數值的差商代替進行離散，從而建立以網格節點上的值為未知數的代數方程組。該方法是一種直接將微分問題變為代數問題的近似數值解法，數學概念直觀，表達簡單，是發展較早且比較成熟的數值方法。

對於有限差分格式，從格式的精度來劃分，有一階格式、二階格式和高階格式。從差分的空間形式來考慮，可分為中心格式和逆風格式。考慮時間因子的影響，差分格式還可以分為顯格式、隱格式、顯隱交替格式等。

構造差分的方法有多種形式，目前主要採用的是泰勒級數展開方法。其基本的差分表達式主要有三種形式：

一階向前差分、一階向後差分、一階中心差分和二階中心差分等，其中前兩種格式為一階計算精度，後兩種格式為二階計算精度。通過對時間和空間這幾種不同差分格式的組合，可以組合成不同的差分計算格式。

4、有限體積法：有限體積法（Finite Volume Method）又稱為控制體積法。其基本思路是：將計算區域劃分為一系列不重復的控制體積，並使每個網格點周圍有一個控制體積；將待解的微分方程對每一個控制體積積分，便得出一組離散方程。其中的未知數是網格點上的因變數的數值。

為了求出控制體積的積分，必須假定值在網格點之間的變化規律，即假設值的分段的分布的分布剖面。從積分區域的選取方法看來，有限體積法屬於加權剩餘法中的子區域法；從未知解的近似方法看來，有限體積法屬於採用局部近似的離散方法。簡言之，子區域法屬於有限體積發的基本方法。

有限體積法的基本思路易於理解，並能得出直接的物理解釋。離散方程的物理意義，就是因變數在有限大小的控制體積中的守恆原理，如同微分方程表示因變數在無限小的控 制體積中的守恆原理一樣。

限體積法得出的離散方程，要求因變數的積分守恆對任意一組控制體積都得到滿足，對整個計算區域，自然也得到滿足。這是有限體積法吸引人的優點。有一些離散方法，例如有限差分法，僅當網格極其細密時，離散方程才滿足積分守恆。

而有限體積法即使在粗網格情況下，也顯示出准確的積分守恆。就離散方法而言，有限體積法可視作有限單元法和有限差分法的中間物。有限單元法必須假定值在網格點之間的變化規律（既插值函數），並將其作為近似解。

有限差分法只考慮網格點上的數值而不考慮值在網格點之間如何變化。有限體積法只尋求的結點值 ，這與有限差分法相類似；但有限體積法在尋求控制體積的積分時，必須假定值在網格點之間的分布，這又與有限單元法相類似。

在有限體積法中，插值函數只用於計算控制體積的積分，得出離散方程之後，便可忘掉插值函數；如果需要的話，可以對微分方程 中不同的項採取不同的插值函數。

5、近似求解的誤差估計方法：近似求解的誤差估計方法共有三大類:單元餘量法,通量投射法及外推法。

單元餘量法廣泛地用於以FEM離散的誤差估計之中,它主要是估計精確運算元的餘量,而不是整套控制方程的全局誤差。

這樣就必須假定周圍的單元誤差並不相互耦合,誤差計算採用逐節點演算法進行。單元餘量法的各種不同做法主要來自對單元誤差方程的邊界條件的不同處理辦法。基於此,該方法能夠有效處理局部的殘餘量,並能成功地用於網格優化程序。

通量投射法的基本原理來自一個很簡單的事實:精確求解偏微分方程不可能有不連續的微分,而近似求解卻可以存在微分的不連續,這樣產生的誤差即來自微分本身,即誤差為系統的光滑求解與不光滑求解之差。該方法與單元餘量法一樣,對節點誤差採用能量范數,故也能成功地用於網格優化程序。

單元餘量法及通量投射法都局限於局部的誤差計算(採用能量范數),誤差方程的全局特性沒有考慮。另外計算的可行性(指誤差估計方程的計算時間應小於近似求解計算時間)不能在這兩種方法中體現,因為獲得的誤差方程數量,階數與流場控制方程相同。

外推是指採用後向數值誤差估計思想由精確解推出近似解的誤差值。各類文獻中較多地採用Richardson外推方法來估計截斷誤差。無論是低階還是高階格式,隨著網格的加密數值計算結果都會趨近於准確解。但由於計算機內存與計算時間的限制,實際上不能採用這種網格無限加密的辦法。

6、多尺度計算方法：近年來發展的多尺度計算方法包括均勻化方法、非均勻化多尺度方法、以及小波數值均勻化方法、多尺度有限體積法、多尺度有限元法等。

該方法通過對單胞問題的求解，把細觀尺度的信息映射到宏觀尺度上，從而推導出宏觀尺度上的均勻化等式，即可在宏觀尺度上求解原問題。均勻化方法在很多科學和工程應用中取得了巨大成功，但這種方法建立在系數細觀結構周期性假設的基礎上，因此應用范圍受到了很大限制。

鄂維南等提出的非均勻化多尺度方法，是構造多尺度計算方法的一般框架。該方法有兩個重要的組成部分：基於宏觀變數的整體宏觀格式和由微觀模型來估計缺少的宏觀數據，多尺度問題的解通過這兩部分共同得到。

該方法基於多分辨分析，在細尺度上建立原方程的離散運算元，然後對離散運算元進行小波變換，得到了大尺度上的數值均勻化運算元。此方法在大尺度上解方程，大大地減小了計算時間。

該法在宏觀尺度上進行網格剖分，然後通過在每個單元里求解細觀尺度的方程(構造線性或者振盪的邊界條件)來獲得基函數。從而把細觀尺度的信息反應到有限元法的基函數里，使宏觀尺度的解包含了細觀尺度的信息。但多尺度有限元方法在構造基函數時需要較大的計算量。