《同濟大學-工程數學-第1章-數值分析與科學計算引論》由會員分享，可在線閱讀，更多相關《同濟大學-工程數學-第1章-數值分析與科學計算引論（40頁珍藏版）》請在裝配圖網上搜索。

1、課程簡介,科學和工程計算是工程類碩士研究生的一門應用性很強的 重要基礎課程，是計算機科學的重要內容。 科學計算是工程實踐的重要工具，本課程主要研究用計算機求解各種數學問題的 數值計算方法 及其理論，簡稱數值計算方法或數值分析。,計算方法,科學計算突破了實驗與 理論方法的局限性 科學計算已和理論、實驗 并列為三大科學方法 （There are three great branches of science：theory，experiment， and computation.）,科學計算的重要性：,1、已經測得在某處海洋不同深度處的水溫如下： 深度（M） 466 741 950 1422 163

2、4 水溫（oC）7.04 4.28 3.40 2.54 2.13 根據這些數據，希望合理地估計出其它深度（如500米，600米，1000米）處的水溫,插值法!,應用問題舉例：,1.梁志偉，馬旭東等.一種基于雙層插值的路徑規(guī)劃及跟蹤算法.機器人2010.11 2.任茂棟，梁晉等.數字圖像相關法中的優(yōu)化插值濾波器.西安交通大學學報2014.7 3.周峰，趙春宇等.基于時域線性插值的信號周期 計算方法及誤差分析儀器儀表學報20101.8 4.史再峰等.基于邊緣方向插值的視頻縮放算法及電路設計.吉林大學學報2009.4,相關領域的應用文獻,2、鋁制波紋瓦的長度問題,建筑上用的一種鋁制波紋瓦是用一種機器

3、將一塊平整的鋁板壓制而成的.,假若要求波紋瓦長48英寸,每個波紋的高度(從中心線)為1英寸,且每個波紋以近似2英寸為一個周期. 求制做一塊波紋瓦所需鋁板的長度L.,這個問題就是要求由函數f(x)=sin x給定的曲線從x=0到x=48英寸間的弧長L. 由微積分學我們知道,所求的弧長可表示為:,上述積分稱為第二類橢圓積分,它不能用普通方法來計算.,數值積分!,主要應用于物理、航天等領域; 神經網絡 1.許少華，王穎等一種基于數值積分的過程神經元 網絡訓練算法.計算機科學.2010.11 2.李盼池，施光堯.基于數值積分的離散過程神經網 絡算法及應用;系統(tǒng)工程理論與實踐2013.12 使用智能優(yōu)化

4、算法求解數值積分問題 1.梁莉莉，韋修喜.云自適應粒子群優(yōu)化算法在 數值積分中的應用.計算機工程與應用 2012.6 2.賴志柱，張云艷.基于遺傳算法求任意函數的數值積分2014.2,全球定位系統(tǒng)：在地球的任何一個位置，至少可以同時收到4顆以上衛(wèi)星發(fā)射的信號,3、全球定位系統(tǒng)(Global Positioning System, GPS),表示地球上一個接收點R的當前位置，衛(wèi)星Si的位置為 ，則得到下列非線性方程組,非線性方程組的數值方法!,10,2020年6月20日,國與國之間和地區(qū)之間的種族歧視、民族矛盾、利益沖突、歷史遺留問題等原因造成了局部戰(zhàn)爭和地區(qū)性武裝沖突時有發(fā)生，有的長期處于敵對

5、狀態(tài)，必然會導致敵對雙方的軍備競賽，軍事裝備現已成為決定戰(zhàn)爭勝負的重要因素 軍事裝備: 軍事實力的總和，主要包括武器裝備、電子信息裝備、軍事兵力、軍事費用等,現代戰(zhàn)爭的特點是多兵種的協(xié)同作戰(zhàn)，根據不同兵種的特點,在不同的區(qū)域參加戰(zhàn)斗，都對戰(zhàn)爭的結果產生一定的影響,4.戰(zhàn)爭的預測與評估問題,11,2020年6月20日,4.戰(zhàn)爭的預測與評估問題,現在要求建立數學模型討論的問題： (1) 分析研究引起軍備競賽的因素，并就諸多因素之間的相互關系進行討論； (2) 在多兵種的作戰(zhàn)條件下，對作戰(zhàn)雙方的戰(zhàn)勢進行評估分析. （3）分析研究作戰(zhàn)雙方的兵力消耗，并預測初始總兵力和戰(zhàn)斗力變化對作戰(zhàn)結果的影響。,12

6、,2020年6月20日,4.戰(zhàn)爭的預測與評估問題,2. 模型的假設,13,2020年6月20日,3. 模型的建立與求解,4.戰(zhàn)爭的預測與評估問題,14,2020年6月20日,4.戰(zhàn)爭的預測與評估問題,微分方程模型！,第一章 緒論,1.1 計算方法的意義 1.2 誤差及有關概念 1.3 數值計算中必須注意的幾個原則,主要內容：,本課程主要內容包括線性方程組的數值求解、非線性方程求根、插值與逼近、數值積分、常微分方程數值解等.,1.1 計算方法的意義,特點：,現 實 世 界,研究 對象,測量 數據,數學模型的建立,數值分析,程序設計,測量 誤差,模型誤差,截斷誤差(方法誤差),舍入誤差,上機計算

7、求得結果,1.2 誤差及有關概念(error),1.2.1 誤差來源,19,例如，用泰勒(Taylor)多項式,近似代替函數 ，,則數值方法的截斷誤差(truncation error)是,在0與x之間。,20,產生的誤差,用3.14159近似代替 ，,就是舍入誤差.,例如，,有了計算公式后，在用計算機做數值計算時，還要受計算機字長的限制，計算過程又可能產生新的誤差，這種誤差為舍入誤差(roundoff error).,若能根據測量工具或計算情況估計出誤差絕對值的一個上界，即,1.2.2 絕對誤差與相對誤差,設 為準確值，,為 的一個近似值，,通常準確值x 是未知的，,因此誤差 也是未知的.,

8、為近似值的絕對誤差(absolute error) ，,定義1,稱,簡稱誤差.,則 叫做近似值的誤差限，,它總是正數.,例如，用毫米刻度的米尺測量一長度 ，讀出和該長度接近的刻度 ，,是 的近似值，,它的誤差限是 ，,于是,如讀出的長度為 ，,則有 .,雖然從這個不等式不能知道準確的 是多少，但可知,結果說明 在區(qū)間 內.,對于一般情形 ，,即,也可以表示為,需要注意的是誤差限的大小并不能完全表示近似值的好壞.,例如，有兩個量 ，,則,把近似值的誤差 與準確值 的比值,稱為近似值 的相對誤差(relative error) ，,記作 .,上例中 與 的相對誤差限分別為,可見 近似 的程度比 近

9、似 的程度好.,根據定義，,當準確值x 位數比較多時，常常按四舍五入的原則得 到x的前幾位近似值 ，,1.2.3 準確位數與有效數字,如取 作為 的近似值，,取 ，,按這個定義，,就有3位有效數字，,就有5位有效數字.,按定義，,187.93, 0.037856, 8.0000, 2.7183.,的5位有效數字近似數是8.0000，而不是8，,例1 按四舍五入原則寫出下列各數具有5位有效數字的 近似數：187.9325, 0.03785551, 8.000033, 2.7182818.,上述各數具有5位有效數字的近似數分別是,因為8只有1位有效數字.,注意：,（1-1）,其中 是0到9中的一個

10、數字， 為整數，則此時 有n位有效數字,（1-2）,且絕對誤差限,若近似值 表示為,在 相同的情況下， 越大則 越小， 故有效位數越多，絕對誤差限越小.,其相對誤差限為,有效位數越多，相對誤差限越小.,1.3 數值計算中必須注意的幾個原則,前面的誤差分析只適用于簡單情形，一個工程或 科學計算問題往往要運算千萬次，由于每步運算都有誤差， 如果每步都做誤差分析是不可能的，也不科學.,因為誤差積累有正有負，絕對值有大有小，都按最壞 情況估計誤差限得到的結果比實際誤差大得多，這種保守 的誤差估計不反映實際誤差積累.到目前為止，誤差分析只 是定性分析，定量分析尚沒有好的方法。,要避免絕對值很小的數做除數

11、,要避免兩相近數相減,要防止大數“吃掉”小數,注意簡化計算步驟，減少運算次數,數值計算中通常不采用數值不穩(wěn)定算法，在設計算法時 還應盡量避免誤差危害，防止有效數字損失，注意這樣幾個 原則：,例如：,對分母做微小變化：,分母只有0.0001的變化，計算結果變化很大。 避免很小的數做除數。,解,只有一位有效數字.,則具有3位有效數字.,若改用,此例說明，可通過改變計算公式避免或減少有效數字 的損失.,類似地，如果 和 很接近時，由,用右邊算式有效數字就不損失.,也應該用右端算式代替左端.,當 很大時，,一般情況，當 時，可用泰勒展開,取右端的有限項近似左端.,如果無法改變算式，則采用增加有效位數進

12、行運算；,在計算機上則采用雙倍字長運算，但這要增加機器計算時間和多占內存單元.,例3 已知 a0, a1, a2 , an, x, 計算多項式： 直接計算：運算量（乘） 秦九韶算法（1247年）（horner方法1819）：,運算量n,例4 利用公式,的前 項和，可計算 的近似值（令 ）.,若要精確到 ，需要對 N=100000項求和，此時不但計 算量大，舍入誤差的積累也很嚴重.,若改用,取 ，只要計算前10項之和，其截斷誤差便小于 .,1.了解數值分析研究的對象及特點；了解誤差的來源及分類； 2.掌握誤差與有效數字的概念，掌握數值運算的誤差估計方法； 3.了解數值算法的穩(wěn)定性；掌握避免誤差危害的若干原則。,本章總結,