普通高等教育“十一五”***規劃教材/中國科學技術大學數學教學叢書:數值計算方法與算法(第三版) pdf epub mobi txt 電子書 下載
産品特色
編輯推薦
適讀人群 :計算機是數值計算方法最常用的計算工具,隨著計算機技術的迅速發展和普及,計算方法課程已成為所有理工科學生的必修課程。 中國科學技術大學考研指定圖書,內容經典,高屋建瓴。
內容簡介
本書覆蓋瞭計算方法最基本的內容,內容的順序為“插值”、“數值微分和數值積分”、“麯綫擬閤”、“非綫性方程求根”、“解綫性方程組的直接法”、“解綫性方程組的迭代法”、 “計算矩陣特徵值和特徵嚮量”和“常微分方程數值解”。最後一章給齣用符號計算語言Mathematica做各章計算方法的例題。
本書參考瞭國內外多本計算方法教材,例如,由教育部高等教育司推薦的國外優秀信息科學與技術係列教學用書,Richard L. Burden數值分析(Numerical Analysis),並吸取瞭他們的優點,例如,給齣大部分方法對應的算法。通過算法縮短數學方法和計算機實現的距離。
本書例題豐富,通過典型例題幫助學生進一步理解計算對象、計算公式、限定條件和計算步驟。學習計算方法中的逼近和迭代等數學思想,掌握常用的數值方法,獲取近似計算的能力,激發學生的學習興趣,擴大學生數值計算的知識麵,並能觸類旁通地應用到各自的科研和技術領域中,培養學生的數學綜閤分析能力和計算能力。
作者簡介
張韻華,中國,漢,1972.4至1975.12就讀於中國科學技術大學數學係計算數學專業,畢業後留校任教至今,1991.10獲中國科學技術大學計算機係軟件碩士學位。1991.10獲中國科學技術大學計算機係軟件碩士學位。從事數值計算方法、符號計算係統的和計算機輔助教學的應用。已發錶相關論文10多篇和著作7部。
目錄
第三版目錄
第1章 插 值
1.1拉格朗日(Lagrange)插值多項式
1.1.1綫性插值
1.1.2 二次插值
1.1.3 n次拉格朗日插值多項式
1.2牛頓(Newton) 插值多項式
1.2.1 差商及其計算
1.2.2 Newton插值
*1.3 Hermite插值
1.4 三次樣條函數
1.4.1分段插值
1.4.2 三次樣條插值的M關係式
1.4.3 三次樣條插值的m關係式
習題1
第2章?最小二乘擬閤
2.1擬閤函數
2.2多項式擬閤
2.3矛盾方程組
習題2
第3章 非綫性方程求解
3.1 迭代法
3.1.1實根的對分法
3.1.2 不動點迭代
3.1.2 不動點迭代
3.2牛頓迭代法
3.3弦截法
3.4求解非綫性方程組的牛頓方法
習題3
第4章?求解綫性方程組的直接法
4.1 Gauss消元法
4.1.1 Gauss順序消元法
4.1.2 Gauss列主元消元法
4.2 直接分解法
4.2.1 Doolittle分解
4.2.2 Crout分解
4.2.3 特殊綫性方程組
習題4
第5章 求解綫性方程組的迭代方法
5.1 簡單(Jacobi)迭代
5.1.1 Jacobi迭代計算公式
5.1.2 Jacobi 迭代收斂條件
5.2高斯-賽德爾(Gauss-Seidel)迭代
5.2.1高斯-賽德爾迭代計算
5.2.2 高斯-賽德爾迭代矩陣
5.3 鬆弛迭代
5.3.1鬆弛迭代計算公式
5.3.2 鬆弛迭代矩陣
習題5
第6章數值積分和數值微分
6.1 牛頓-柯特斯數值積分
6.1.1 插值型數值積分
6.1.2 牛頓-柯特斯(Newton-Cote’s)積分
6.2 復化數值積分
6.2.1 復化梯形積分
6.2.2復化Simpson積分
6.2.3自動控製誤差的復化積分
6.2.4龍貝格 (Romberg )積分
*6.3 重積分計算
*6.4 高斯(Gauss)型積分
6.5 數值微分
6.5.1 差商與數值微分
6.5.2插值型數值微分
習 題 6
第7章 常微分方程數值解
7.1 歐拉(Euler)公式
7.1.1 基於數值微商的歐拉公式
*7.1.2 歐拉公式的收斂性
7.1.3 基於數值積分的近似公式
7.2 龍格-庫塔方法
7.2.1二階龍格-庫塔方法
7.2.2 四階龍格-庫塔格式
7.3 綫性多步法
7.4 常微分方程組的數值解法
7.4.1 一階常微分方程組的數值解法
7.4.2 高階常微分方程數值方法
*7.5常微分方程的穩定性
習題7
第8章 計算矩陣的特徵值和特徵嚮量
8.1 冪法
8.1.1 冪法計算
8.1.2 冪法的規範運算
8.2反冪法
*8.3 實對稱矩陣的Jacobi方法
*8.4 QR方法簡介
8.4.1 QR方法初步
8.4.2矩陣的QR分解
習題8
附錄1 上機作業題
附錄2 C語言程序示例
程序1 構造Newton插值多項式
程序2 弦截法求根
程序3 超鬆弛迭代求解方程組
程序4 龍貝格( Romberg )積分算法
程序5 四階龍格—庫塔法求解常微分方程初值問題
附錄3 在符號語言Mathematica中做題
1插值
2麯綫擬閤
3求解非綫性方程
4求解綫性方程組
5數值積分
6常第1章 插 值
第2章?最小二乘擬閤
第3章 非綫性方程求解
第4章?求解綫性方程組的直接法
第5章 求解綫性方程組的迭代方法
第6章數值積分和數值微分
第7章 常微分方程數值解
第8章 計算矩陣的特徵值和特徵嚮量
附錄1 上機作業題
附錄2 C語言程序示例
附錄3 在符號語言Mathematica中做題
精彩書摘
緒論
0.1數值計算方法與算法
數值計算方法,是一種研究數學問題的數值近似解方法,是在計算機上使用的解數學問題的方法,簡稱計算方法。它的計算對象是那些在理論上有解而又無法用手工計算的數學問題,以及沒有解析解的數學問題。例如,解一個有300個未知量的綫性方程組;計算6階矩陣的全部特徵值。
在科學研究和工程技術中都要用到各種計算方法。例如,在航天航空、地質勘探、汽車製造、橋梁設計、天氣預報和漢字字樣設計中都有計算方法的蹤影。在70年代,大多數學校僅在數學係的計算數學專業和計算機係開設計算方法這門課程。隨著計算機技術的迅速發展和普及,現在計算方法課程幾乎已成為所有理工科學生的必修課程。
計算方法是一門理論性和實踐性都很強的學科,計算方法既有數學類課程中理論上的抽象性和嚴謹性,又有實用性和實驗性的技術特徵。計算方法的前提課程是微積分,綫性代數,常微分方程和一門計算機語言。
大多數人學習計算方法的目的是為瞭使用方法,在學習計算方法中,在套用計算公式、修改計算公式和創建計算公式中,都需要不同程度的專業知識和數學基礎。要注重學習計算方法中的逼近和迭代等數學思想和常用手法,獲取近似計算的能力,並能觸類旁通地應用到各個領域中。一些有創造力的工程師不僅擅長使用某些計算方法,並能創建齣簡便有效的計算方法。例如,樣條函數、快速富裏葉變換和有限元方法都是有創造力的工程師們創建的,再由數學傢們完善這些方法的理論基礎,並從理論上進行提高和推廣。
從方法的計算公式到在計算機上實際運行,兩者之間還有距離,這是數學能力與計算機應用技術能力之間的距離,還與計算機的運行環境和編程工具有關,為瞭縮小兩者之間的距離,本文將給齣部分計算公式的算法描述。用算法容易準確而簡便地描述計算公式,在算法中能簡潔地錶達計算公式中的“循環”和“迭代”等操作。有瞭方法的算法,將它轉化成C或PASCAL等語言的程序上機運行也就容易瞭。
在學習計算方法過程中,如果能用某種語言編製該方法的程序並運行通過,那麼有利於準確而深刻地掌握該方法的計算步驟和過程。
本教材中提供瞭部分上機作業題,在平時作業中布置一些上機編程題目,其目的是通過編程上機,加深對方法實施的理解和體會,訓練和提高數學與計算機應用能力和水平。
0.2 誤差與有效數字
絕對誤差與絕對誤差界
近似計算必然産生誤差,誤差錶示精確值與近似值的距離。
定義0.1 設 為精確值(或準確值), 是 的一個近似值,稱 為近似值 的絕對誤差或誤差。
絕對誤差 = 精確值 ? 近似值
誤差 的值可正可負,如果得不到精確值 ,也就算不齣絕對誤差 的值。常用限製誤差絕對值的範圍 描述和控製誤差的範圍。
定義0.2 如果精確值 與近似值 的誤差的絕對值不超過某正數 ,即
稱 為絕對誤差限或誤差限。
精確值 也可錶示為: 。通常,在誤差允許的範圍內的近似值 即認為是精確值,這也是計算中控製循環中止的常用手段。
例0.1 若經四捨五入得到 ,對於數 , 的近似值都是 ,
即第四位小數大於5時,必然進位到第三位小數;第四位小數小於5時,必然捨去。
它的誤差限是:
若 ,則它的誤差限是:
相對誤差與相對誤差限
在很多情況下,絕對誤差並不能全麵的反映近似程度。例如,某電器公司兩次進貨的某型號電風扇分彆為1000颱和2000颱,其中開箱不閤格電風扇分彆為8和12(絕對誤差的值)。不閤格率分彆為8/1000=0.8%和12/2000=0.6 %(相對誤差的值),這說明該電風扇的質量有所提高。我們把近似值與準確值的比值定義為相對誤差。
定義0.3 設 為精確值, 是 的一個近似值
稱 為近似值 的相對誤差。
在實際計算中,有時得不到精確值 ,當 較小時 可用近似值 代替,即
相對誤差 = 或 相對誤差 =
相對誤差 的值也可正可負,與絕對誤差一樣不易計算,常用相對誤差限控製相對誤差的範圍。
定義0.4 如果有正數 使得 ,則稱 為 的相對誤差限。
産生誤差的因素很多,産生誤差的原因主要有
原始誤差
由客觀存在的模型抽象到物理模型産生的誤差。包括模型誤差和原始數據誤差。
截斷誤差
用有限項近似無限項時,由截取函數的部分項而産生的誤差,稱為截斷誤差。
例如: ,在計算中用
的截斷誤差
捨入誤差
在數值計算中,通常都按有限位進行運算。例如,按照4捨5入的原則,2 / 3=0.666667
或2 / 3=0.667,由捨入産生的誤差,稱為捨入誤差。
在實際計算中的數據通常是近似值,它們由觀察、估計或一些計算而得到,這些數在計算機錶示後也會帶來進一步誤差,即誤差的積纍和傳播。關於誤差的傳播似乎沒有多少統一的理論,通常積纍誤差的界是以通例分析為基礎而建立的。
有效位數
定義0.5 當 的誤差限為某一位的半個單位,則這一位到第一個非零位的位數稱為 的有效位數。
例如, = 12.34 , = 0.004067均有4位有效數字,而3.00與3.0000分彆有3位和5位有效位數。
有效位的多少直接影響到近似值的絕對誤差和相對誤差,因此,在計算中也應注意保持一定的有效位數。
數值計算的近似計算免不瞭有誤差相隨,隻能盡量約束和控製誤差。
選擇收斂的穩定的方法
對同一問題選擇不同的數值計算方法,可能得到不同的計算結果。在計算方法中,除瞭給齣方法的數值計算公式,還要討論計算公式的收斂性、穩定性和截斷誤差的特性。選擇收斂性要求低、穩定性好的方法是約束誤差擴張最重要的措施。例如,樣條插值函數比高次多項式的效果好的多,是構造插值函數的首選方法。
提高數值計算精度
數值在計算機中存放的位數稱為字長。有限位的字長是帶來捨入誤差和抑製數值計算精度的根源。對同一種方法,在字長大的計算機上的計算效果要比在字長小的計算機上優越。
同一計算問題,簡化計算步驟、減少運算次數、控製除法中分母的值等措施都會約束和減少捨入誤差.
例如:將多項式錶達式 ,改寫為
在計算機上,用同一種數值計算方法對數據選用不同的數值類型,有時會直接影響到計算效果。例如,對病態的綫性方程組,采用單精度數據的Gauss消元方法,其數據解大大失真,而用雙精度數據Gauss列主元消元方法卻可得到滿意的數值解。
0.3 矩陣和嚮量範數
0.3.1嚮量範數
1. 嚮量範數的定義
在一維空間中,實軸上任意兩點 的距離用兩點坐標差的絕對值 錶示. 絕對值是單變量的一種度量距離的定義.
範數是在廣義長度意義下,對函數、嚮量和矩陣的一種度量定義. 任何對象的範數值都是一個非負實數. 使用範數可以測量兩個函數、嚮量或矩陣之間的距離. 嚮量範數是度量嚮量長度的一種定義形式. 範數有多種定義形式,隻要滿足定義1.1中的三個條件即可定義一個範數.
對任一嚮量 ,按照一個規則確定一個非負實數與它對應,記該實數為 ,若 滿足下麵三個性質:
(1) 任取 ,有 ,當且僅當 時, (非負性)
(2) 任取 , 有 (齊次性)
(3) 任取 ,有 (三角不等式)
那麼稱實數 為嚮量 的範數.
定義 0.6 嚮量 的 範數( 範數)定義為
, (0.1)
其中,經常使用三種 嚮量範數是 .
1-範數(曼哈頓範數)
2-範數(歐幾裏得範數)
或寫成
-範數
注:
例0.2 計算嚮量 的嚮量範數.
,
例0.3 設A是一個正定矩陣,對任何嚮量 ,定義函數 ,
是一種嚮量範數.
例0.4 當 , 不是嚮量範數.
證明 取 則
,
不是嚮量範數.
2. 不同嚮量範數的關係
同一嚮量,在不同的範數定義下,得到不同的範數值. 定理0.1給齣有限維綫性空間 中任意嚮量範數都是等價的.
定理0.1 若 是 上兩種不同的範數定義,則必存在 ,使 ,均有
(0.2)
或 (證明略)
可以驗證,對於嚮量的1、2和 範數有下列等價關係
,
例0.5 圖示 中嚮量1範數、2範數、4範數和 範數的單位“圓”.
圖0-1範數的單位“圓”
嚮量的極限
嚮量範數的定義提供瞭度量兩個嚮量的距離標準,即可定義嚮量的極限和收斂概念瞭.
定義 0.7 設 為 上嚮量序列,若存在嚮量??
有 ,則稱嚮量列 是收斂的, 稱為該嚮量序列的極限.
由嚮量範數的等價性,嚮量序列是否收斂與選取哪種範數無關. 嚮量的極限是通過它的所有分量的極限定義的。不論選取那種範數,嚮量序列 收斂的充分必要條件為其序列的每個分量收斂,即 存在. 若 ,則 就是嚮量序列 的極限. 在數值計算中,當迭代的嚮量序列中相鄰兩個嚮量的誤差 給定精度時,視 為極限嚮量 .
0.3.2矩陣範數
1. 矩陣範數定義
設 ,記方陣 的範數為 ,矩陣範數滿足下列性質:
(1) 當且僅當 時, (非負性)
(2) (齊次性)
(3) 對於任意兩個同階矩陣 有
(三角不等式)
(4)設 為同階矩陣,則 (相容性)
(5)設 為 階陣,對 ,
恒有
隻要滿足(1) (2) (3)就可以定義一個矩陣範數。矩陣範數可用嚮量範數定義.
定義 0.8 設 ,定義矩陣範數 (0.3)
下麵簡化矩陣範數的(0.3)式的定義.
,設 且 ,則
即:
這樣,在 上 的選取範圍由一張平麵壓縮到單位圓周上;在 上選取範圍由三維空間壓縮到單位球麵上.
定義0.9 設 是 上的一個嚮量範數,則由
(0.4)
定義的實值函數 是一個矩陣範數。
這類範數稱為算子範數,誘導範數或從屬範數。幾何直觀上,矩陣範數是矩陣對嚮量的最大拉伸。
2. 常用矩陣範數
對應於嚮量的三種範數,相應的三種矩陣範數形式為:
(列和範數)
(行和範數)
,其中 , 是 的特徵值
* 證明: 矩陣的範數是 上滿足 嚮量範數 的上確界,那麼,找到這個上確界也就找到瞭矩陣的範數.
(1) 任取 設 ,則
即
設極大值在 列達到,有 ,取 , 除第 個分量為1外,其餘分量均為0,
於是有 .
由定義和 ,故 . 因此有
(2) 任取 設 ,則
即
另一方麵設極大值在k行達到,取
這裏
於是 . 故
(3) 為半正定對稱矩陣,具有非負特徵值,並具有n個相互正交的單位特徵嚮量.
設 的特徵值為 ,相應的特徵嚮量為 ,其中 為相互正交的
普通高等教育“十一五”***規劃教材/中國科學技術大學數學教學叢書:數值計算方法與算法(第三版) 下載 mobi epub pdf txt 電子書
普通高等教育“十一五”***規劃教材/中國科學技術大學數學教學叢書:數值計算方法與算法(第三版) pdf epub mobi txt 電子書 下載