計(jì)算方法第一章誤差.ppt

數(shù)值計(jì)算方法,授課老師：涂泳秋 Email：tyq_,數(shù)值計(jì)算方法,上課時(shí)間： 1-18周周三上午三、四節(jié) 上課地點(diǎn)： 第4，8，12，16周：1班4#機(jī)房，2班5#機(jī)房 其余周次：A1-403 考試方式： 閉卷 成績(jī)計(jì)算方法： 筆試60%，平時(shí)20%，上機(jī)20% 實(shí)驗(yàn)工具：MatLab,數(shù)值分析,數(shù)值分析是數(shù)學(xué)與計(jì)算機(jī)技術(shù)結(jié)合的 一門學(xué)科，是利用計(jì)算機(jī)解決數(shù)學(xué)問(wèn)題的理論和方法，是計(jì)算數(shù)學(xué)的一個(gè)重要分支。,現(xiàn)代復(fù)雜工程技術(shù)問(wèn)題的解決步驟,工程問(wèn)題,上機(jī)計(jì)算,數(shù)學(xué)模型,結(jié)果分析,問(wèn)題解答,設(shè)計(jì)算法,數(shù)值分析涉及的主要內(nèi)容,計(jì)算機(jī)只能進(jìn)行加減乘除四則運(yùn)算和一些簡(jiǎn)單的函數(shù)計(jì)算（即使函數(shù)也是通過(guò)數(shù)值分析方法處理，轉(zhuǎn)化為四則運(yùn)算而形成的小型軟件包） 數(shù)值代數(shù)：求解線性方程組和非線性方程組的解法，分直接方法和間接方法 插值和數(shù)值逼近。離散的點(diǎn)上的函數(shù)值，想辦法得到點(diǎn)之間的值 數(shù)值微分和數(shù)值積分。很多函數(shù)無(wú)法求出積分，利用數(shù)值方法求解 常微分方程和偏微分方程的數(shù)值解法,數(shù)值分析需要考慮哪些問(wèn)題,計(jì)算速度 例如：求解一個(gè)20階線性方程組，20個(gè)未知量，用加減消原法需3000次乘法運(yùn)算，用行列式求解需進(jìn)行9.7*1020次運(yùn)算，如果用每秒1億次乘法運(yùn)算的計(jì)算機(jī)要30萬(wàn)年。 說(shuō)明了算法方法的重要性 2. 存儲(chǔ)量。大型問(wèn)題有必要考慮 例如算法所需要保留的中間結(jié)果比較少，則可以省下為保留中間結(jié)果所需要的額外的存儲(chǔ)空間。,數(shù)值分析需要考慮哪些問(wèn)題,數(shù)值穩(wěn)定性 在大量計(jì)算中，舍入誤差是積累還是能控制，這與算法有關(guān)。 例：一元二次方程,其精確解為,如用求根公式：,和字長(zhǎng)為8位的計(jì)算器求解，有,又,數(shù)值分析需要考慮哪些問(wèn)題,誤差的來(lái)源和基本概念,誤差的來(lái)源 模型誤差：在建立數(shù)學(xué)模型過(guò)程中，不可能將所有因素都考慮，必然要進(jìn)行必要的簡(jiǎn)化，帶來(lái)了與實(shí)際問(wèn)題的誤差。不是數(shù)值方法考慮的問(wèn)題 測(cè)量誤差：測(cè)量已知參數(shù)時(shí)，數(shù)據(jù)帶來(lái)的誤差。也不是數(shù)值方法考慮的問(wèn)題 截?cái)嗾`差：在設(shè)計(jì)算法時(shí)，必然要近似處理，尋求簡(jiǎn)化。這是計(jì)算數(shù)學(xué)考慮問(wèn)題 舍入誤差：計(jì)算機(jī)的字長(zhǎng)是有限的，每一步運(yùn)算均需四舍五入，由此產(chǎn)生的誤差稱舍入誤差。 數(shù)值分析主要討論截?cái)嗾`差。測(cè)量誤差看成初始的舍入誤差，數(shù)值分析也要從整體上討論舍入誤差的影響,誤差的基本概念,絕對(duì)誤差和絕對(duì)誤差限 X*是精確值，x是它的一個(gè)近似值，稱e=x-x*是近似值x的絕對(duì)誤差，簡(jiǎn)稱誤差。絕對(duì)誤差可正可負(fù)，是有量綱的。誤差是無(wú)法計(jì)算的，但可以估計(jì)出它的一個(gè)上界。即 ：,相對(duì)誤差和相對(duì)誤差限,相對(duì)誤差和相對(duì)誤差限,記作,相對(duì)誤差是個(gè)相對(duì)數(shù)，是無(wú)量綱的，也可正,可負(fù)。相對(duì)誤差的估計(jì),，即：,實(shí)際計(jì)算中，x*未知，用x代替，兩者的差為：,例：重力加速度常數(shù)g。,兩者均有三位有效數(shù)字,后者的絕對(duì)誤差大，而相對(duì)誤差分別為,兩者相等，與量綱的選取無(wú)關(guān),有效數(shù)字,若近似值x的絕對(duì)誤差限是某一位數(shù)的半個(gè)單位，則說(shuō) x 精確到該位，若從該位到 x的左面第一位非零數(shù)字一共有n位，則稱近似值x有n位有效數(shù)字。,例：,3.14有三位有效數(shù)字,；3.1416有五位有效數(shù)字，,誤差限為0.00005,又例：0.003529是四位有效數(shù)字，0.00352900,是六位有效數(shù)字。前者的誤差限為,，,后者為,，寫成標(biāo)準(zhǔn)的浮點(diǎn)數(shù)為：,當(dāng),其中，有效數(shù)字的個(gè)數(shù)是l,即,有效數(shù)字：由絕對(duì)誤差決定,有效數(shù)字與相對(duì)誤差,此定理說(shuō)明，相對(duì)誤差限是由有效數(shù)字決定。,數(shù)值計(jì)算中誤差的傳播,1. 對(duì)函數(shù)的計(jì)算：,數(shù)值計(jì)算中誤差的傳播,2. 對(duì)多元函數(shù)的計(jì)算：,數(shù)值計(jì)算中誤差的傳播,3. 四則運(yùn)算中誤差的傳播：,證明第二條如下：,數(shù)值計(jì)算中的若干準(zhǔn)則,關(guān)于數(shù)值穩(wěn)定性的算法 一個(gè)程序往往需要進(jìn)行大量的四則運(yùn)算才能得出結(jié)果，每一步的運(yùn)算均會(huì)產(chǎn)生舍入誤差。在運(yùn)算過(guò)程中，舍入誤差能控制在某個(gè)范圍內(nèi)的算法稱之為數(shù)值穩(wěn)定的算法，否則就稱之為不穩(wěn)定的算法。,數(shù)值計(jì)算中的若干準(zhǔn)則,結(jié)果只有一位有效數(shù)字，有效數(shù)字大量損失，造成,相對(duì)誤差擴(kuò)大。這是由兩個(gè)比較接近的數(shù)相減造成的,結(jié)果仍然有四位有效數(shù)字。這說(shuō)明算法設(shè)計(jì)的重要性,2. 注意避免兩個(gè)相近數(shù)的減法,數(shù)值計(jì)算中的若干準(zhǔn)則,3. 簡(jiǎn)化計(jì)算步驟，減少運(yùn)算次數(shù),例：計(jì)算 的值。,只需14次乘法。,又如：,采用“秦九韶算法”,例：計(jì)算多項(xiàng)式,只需n次乘法和n次加法。,一般要注意，能在循環(huán)外