版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認(rèn)領(lǐng)
文檔簡介
1、Lecture: numerical methods using Matlab(3)The school of information and electronic engineering, NBUTLecturer: B. TanNumerical Differentiation數(shù)值微分 為計算f(x),首先對連續(xù)函數(shù)進行離散化,x1x2x3x4x5f(x1)f(x2)f(x3)f(x4)f(x5)hNumerical Differentiation數(shù)值微分 定義計算機只能近似處理無窮小hh足夠小,近似計算微分滿足誤差需求Numerical Differentiation數(shù)值微分 h足夠小M
2、atlab的最小數(shù)epseps=2.2204e-016取h=eps(sin(pi/3+eps)-sin(pi/3)/epsans = 0.5000計算sin(pi/3)的導(dǎo)數(shù)。Numerical Differentiation數(shù)值微分 微分方程Eulers methodNumerical Differentiation數(shù)值微分 運用遞推迭代算法可以計算上式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),遞推迭代計算一次優(yōu)點:一個變量y,節(jié)約空間缺點:只能得到最后的y,中間值沒有存儲。Numerical Differentiation數(shù)值微分 分析解為:將分析解和數(shù)值解畫在一個圖中,對比觀察。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米寬的河道如下圖所示,為了測量其流量需要知道河道的截面積。為此從一端開始每個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單位:米估計河道的截面積,假定流速1m/s計算流量。若沿河床鋪設(shè)一條光纜,估計鋪設(shè)的長度。Numerical Integration 數(shù)值積分 問題分析:由于沒有確定的解析表達式,無法計算精度可控的遞推算法。采用的方法確定為:5米的間距過大,對測量數(shù)據(jù)進行數(shù)據(jù)插值,再用trapz語句直接計算。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ù)據(jù)Numerical Integration 數(shù)值積分 S=abs(trapz(x1,-y1)S = 339.4320為了測量精度的提高,方法是更密集的測量河道深度。對數(shù)據(jù)進行插值,能夠提高計算的精度嗎?Numerical Integration 數(shù)值積分 dxdydl計算光纜長度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)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 課件:《躺在波浪上看書》
- 2022高考英語唐山市路南區(qū)閱讀短文復(fù)習(xí)單詞暑假練習(xí)(12)及答案
- 【原創(chuàng)】江蘇省2020-2021學(xué)年高二第一學(xué)期第五次周練語文試題
- 山東省德州市2025屆高三上學(xué)期期中考試數(shù)學(xué)試卷(含解析)
- 【名師一號】2020-2021學(xué)年高中地理湘教版必修二-雙基限時練3
- 2024-2025學(xué)年部編版歷史七年級下冊期中綜合評估卷(第1-11課)(含答案)
- 【名師一號】2020-2021學(xué)年高中生物必修三:第二章-動物和人體生命活動的調(diào)節(jié)-單元檢測
- 《金版學(xué)案》2022屆高考數(shù)學(xué)理科一輪復(fù)習(xí)課時作業(yè)-3-8解三角形的應(yīng)用-
- 一年級數(shù)學(xué)計算題專項練習(xí)1000題匯編
- 【2020秋備課】高中生物學(xué)案新人教版必修3-2.1-通過神經(jīng)系統(tǒng)的調(diào)節(jié)
- 人體解剖上肢血管
- 精品解析浙教版科學(xué) 九年級上冊 3.43 簡單機械之機械效率 同步練習(xí)
- 六棱塊護坡施工方案
- 機械制圖課件(完整版)
- 夸美紐斯-大教學(xué)論-文本細(xì)讀
- 《行政組織學(xué)小抄》word版
- 日立多聯(lián)機系統(tǒng)調(diào)試培訓(xùn)教材
- (完整版)環(huán)境科學(xué)與工程-專業(yè)英語詞匯必備(免費)
- 交通管理與控制課件(全)全書教學(xué)教程完整版電子教案最全幻燈片
- 小學(xué)鋼筆字寫字課教案(20課時完整版)
- 紅金簡約風(fēng)教師退休歡送會PPT通用模板
評論
0/150
提交評論