哈工程數值計算方法

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

數值計算方法如下：

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、多尺度計算方法：近年來發展的多尺度計算方法包括均勻化方法、非均勻化多尺度方法、以及小波數值均勻化方法、多尺度有限體積法、多尺度有限元法等。

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

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

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

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

2. 計算機專業本科的《數值計算方法》都講了哪些內容

《數值計算方法》是數學類專業（如信息與計算專業、數學與應用數學專業）的專業基礎課，主要包括數值逼近、數值代數和微分方程數值解三個部分。隨著學分制改革的推進，該課程也可作為學校部分工科專業學生的選修課。以前我校面向部分工科專業學生開設的《計算方法》課程的大部分內容都包含在《數值計算方法》課程中。

隨著計算機技術的發展和科學技術的進步，科學計算的應用范圍已擴大到許多的學科領域，已經形成了一些邊緣學科。例如，計算物理、計算力學、計算化學等。目前，實驗、理論和計算已經成為了人們進行科學活動的三大方法。對從事工程與科學技術工作的人員，學習和掌握《數值計算方法》是非常必要的。

數值計算方法是數學的一個分支，但它又不象純數學那樣只研究數學本身的理論，而是把數學理論與計算方法緊密結合，既有純數學高度抽象性的特點，又有應用的廣泛性與實際試驗的高度技術性的特點，是一門與計算機使用密切結合的實用性很強的數學課程，著重研究數學問題的數值方法及其理論。

http://www1.snut.e.cn/math/2007/reseach/math_web/

3. 工程數值方法的MATLAB程序編寫

數值積分法是求定積分的近似值的數值方法。即用被積函數的有限個抽樣值的離散或加權平均近似值代替定積分的值。

求某函數的定積分時，在多數情況下，被積函數的原函數很難用初等函數表達出來，另外，許多實際問題中的被積函數往往是列表函數或其他形式的非連續函數，對這類函數的定積分，也不能用不定積分方法求解。對微積分學作出傑出貢獻的數學大師，如I.牛頓、L.歐拉、C.F.高斯、拉格朗日等人都在數值積分這個領域作出了各自的貢獻，並奠定了這個分支的理論基礎。

數值積分法也是計算機模擬中常用的一種方法。在已知函數的微分方程時，求解函數下一時刻的值，我們主要有歐拉法、梯形法和龍格庫塔法。

歐拉法，這些方法中精度最低的，程序相對簡單。歐拉法的表達式可以寫成下面的形式：

圖中振盪最大的是歐拉法，梯形法和四階龍格庫塔法精度差不多。四階龍格庫塔法精度最高。

4. 傳統的數值計算方法包括哪些內容現在的數值計算方法包括哪些內容

隨著計算機和計算方法的飛速發展，幾乎所有學科都走向定量化和精確化，從而產生了一系列計算性的學科分支，如計算物理、計算化學、計算生物學、計算地質學、計算氣象學和計算材料學等，計算數學中的數值計算方法則是解決「計算」問題的橋梁和工具。我們知道，計算能力是計算工具和計算方法的效率的乘積，提高計算方法的效率與提高計算機硬體的效率同樣重要。科學計算已用到科學技術和社會生活的各個領域中。
數值計算方法，是一種研究並解決數學問題的數值近似解方法， 是在計算機上使用的解數學問題的方法，簡稱計算方法。
在科學研究和工程技術中都要用到各種計算方法。 例如，在航天航空、地質勘探、汽車製造、橋梁設計、 天氣預報和漢字字樣設計中都有計算方法的蹤影。
計算方法既有數學類課程中理論上的抽象性和嚴謹性，又有實用性和實驗性的技術特徵， 計算方法是一門理論性和實踐性都很強的學科。 在70年代，大多數學校僅在數學系的計算數學專業和計算機系開設計算方法這門課程。 隨著計算機技術的迅速發展和普及， 現在計算方法課程幾乎已成為所有理工科學生的必修課程。
計算方法的計算對象是微積分，線性代數，常微分方程中的數學問題。 內容包括：插值和擬合、數值微分和數值積分、求解線性方程組的直接法和迭代法、 計算矩陣特徵值和特徵向量和常微分方程數值解等問題。

5. 標准偏差計算公式是什麼

1、樣本標准偏差：

(5)哈工程數值計算方法擴展閱讀

標准差能反映一個數據集的離散程度，標准偏差越小，這些值偏離平均值就越少，反之亦然。標准偏差的大小可通過標准偏差與平均值的倍率關系來衡量。平均數相同的兩個數據集，標准差未必相同。

例如，A、B兩組各有6位學生參加同一次語文測驗，A組的分數為95、85、75、65、55、45，B組的分數為73、72、71、69、68、67。這兩組的平均數都是70，但A組的標准差應該是17.078分,B組的標准差應該是2.160分，說明A組學生之間的差距要比B組學生之間的差距大得多。