版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
1、Lecture: numerical methods using Matlab(3)The school of information and electronic engineering, NBUTLecturer: B. TanNumerical Differentiation數(shù)值微分 為計(jì)算f(x),首先對連續(xù)函數(shù)進(jìn)行離散化,x1x2x3x4x5f(x1)f(x2)f(x3)f(x4)f(x5)hNumerical Differentiation數(shù)值微分 定義計(jì)算機(jī)只能近似處理無窮小hh足夠小,近似計(jì)算微分滿足誤差需求Numerical Differentiation數(shù)值微分 h足夠小M
2、atlab的最小數(shù)epseps=2.2204e-016取h=eps(sin(pi/3+eps)-sin(pi/3)/epsans = 0.5000計(jì)算sin(pi/3)的導(dǎo)數(shù)。Numerical Differentiation數(shù)值微分 微分方程Eulers methodNumerical Differentiation數(shù)值微分 運(yùn)用遞推迭代算法可以計(jì)算上式Numerical Differentiation數(shù)值微分 Example:初始條件 y(0)=1求解y when t from 0 to 5初始條件:to=0,y(0)=1Numerical Differentiation數(shù)值微分 賦初值:
3、dt = 0.1;t = 0;y = 1;如何編程:遞推迭代可在循環(huán)體內(nèi)完成Numerical Differentiation數(shù)值微分 while t=5y = y + dt * (1-t)*y;t=t+dt;end每次循環(huán),遞推迭代計(jì)算一次優(yōu)點(diǎn):一個(gè)變量y,節(jié)約空間缺點(diǎn):只能得到最后的y,中間值沒有存儲(chǔ)。Numerical Differentiation數(shù)值微分 分析解為:將分析解和數(shù)值解畫在一個(gè)圖中,對比觀察。Numerical Differentiation數(shù)值微分 dt = 0.1;t = 0;y = 1;while t=error h=h/2; x=a:h:b; I2=trapz(x
4、,fun(x); di=I2-I1; I1=I2;endfprintf(%.6fn,I2)quadrature4.mNumerical Integration 數(shù)值積分 例二:一條100米寬的河道如下圖所示,為了測量其流量需要知道河道的截面積。為此從一端開始每個(gè)5米測量出河床的深度如表所示:0 x1x2x3100Numerical Integration 數(shù)值積分 坐標(biāo)x1x2x3x4x5x6x7x8x9x10深度2.412.962.152.653.124.235.126.215.684.22坐標(biāo)X11x12x13x14x15x16x17x18x19x20深度3.913.262.852.353
5、.023.634.123.462.080單位:米估計(jì)河道的截面積,假定流速1m/s計(jì)算流量。若沿河床鋪設(shè)一條光纜,估計(jì)鋪設(shè)的長度。Numerical Integration 數(shù)值積分 問題分析:由于沒有確定的解析表達(dá)式,無法計(jì)算精度可控的遞推算法。采用的方法確定為:5米的間距過大,對測量數(shù)據(jù)進(jìn)行數(shù)據(jù)插值,再用trapz語句直接計(jì)算。Numerical Integration 數(shù)值積分 %quadrature3.my=0,2.41,2.96,2.15,2.65,3.12,4.23,5.12,6.21,5.68,4.22,. 3.91,3.26,2.85,2.35,3.02,3.63,4.12,3.46,2.08,0;x=0:5:100;x1=0:100;y1=interp1(x,y,x1,spline);畫圖語句略實(shí)線為插值后的數(shù)據(jù)Numerical Integration 數(shù)值積分 S=abs(trapz(x1,-y1)S = 339.4320為了測量精度的提高,方法是更密集的測量河道深度。對數(shù)據(jù)進(jìn)行插值,能夠提高計(jì)算的精度嗎?Numerical Integration 數(shù)值積分 dxdydl計(jì)算光纜長度dl=sqrt(diff(x1).2+diff(y1).2);l=sum(dl)l = 102.2249Numerical Integration 數(shù)值積分 d
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 《龍湖葵花寶典培訓(xùn)》課件
- 經(jīng)濟(jì)效益的年度跟蹤計(jì)劃
- 《數(shù)據(jù)圖示餅》課件
- 沿海工程防護(hù)設(shè)備采購合同三篇
- 內(nèi)部激勵(lì)措施的年度優(yōu)化計(jì)劃
- 《通信技術(shù)原理》課件
- 冷拔鋼相關(guān)行業(yè)投資方案
- 合結(jié)鋼行業(yè)相關(guān)投資計(jì)劃提議
- 食品加工合同三篇
- 《液壓與氣動(dòng)》課件 1氣動(dòng)系統(tǒng)概述
- 《高等數(shù)學(xué)》課程思政教學(xué)大綱
- 國際疾病分類手術(shù)碼(ICD-9-CM-3)使用手冊
- 物資部對標(biāo)管理實(shí)施方案
- 上海工程技術(shù)大學(xué)《管理信息系統(tǒng)》 ~學(xué)年 第 一 學(xué)期 實(shí)驗(yàn)報(bào)告
- 職工醫(yī)?;鶖?shù)調(diào)整對比明細(xì)表Excel模板
- 送你一個(gè)字評語2022
- 放射科優(yōu)質(zhì)護(hù)理服務(wù)PPT學(xué)習(xí)教案
- GB_T 22627-2022水處理劑 聚氯化鋁_(高清-最新版)
- 教學(xué)團(tuán)隊(duì)建設(shè)總結(jié)報(bào)告
- 看守所釋放證明書
- 魚骨圖-PPT模板
評論
0/150
提交評論