


下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、專業(yè)序號(hào)姓名日期實(shí)驗(yàn) 3 常微分方程數(shù)值解【實(shí)驗(yàn)?zāi)康摹? .掌握用MATLAB求微分方程初值問題數(shù)值解的方法;2 .通過實(shí)例學(xué)習(xí)微分方程模型解決簡(jiǎn)化的實(shí)際問題;3 .了解歐拉方法和龍格庫(kù)塔方法的基本思想。【實(shí)驗(yàn)內(nèi)容】用歐拉方法和龍格庫(kù)塔方法求下列微分方程初值問題的數(shù)值解,畫出解的圖形,對(duì)結(jié)果進(jìn)行分析比較yy2x,x(1)W(0 x1),精確解y=3e-2x-2;y(0)a-22y=x-y,(2)(0_x_10).y(0)=0或丫(0)=1【解】:手工分析怎樣求解【計(jì)算機(jī)求解】:怎樣設(shè)計(jì)程序?流程圖?變量說明?能否將某算法設(shè)計(jì)成具有形式參數(shù)的函數(shù)形式?【程序如下】:functionf=f(x,y
2、)f=y+2*x;clc;clear;a=0;b=1;%求解區(qū)間x1,y_r=ode45(f,ab,1);%調(diào)用龍格庫(kù)塔求解函數(shù)求解數(shù)值解;%以下利用Euler方法求解y(1)=1;N=100;h=(b-a)/N;x=a:h:b;fori=1:Ny(i+1)=y(i)+h*f(x(i),y(i);endfigure(1)plot(x1,y_r,r*,x,y,b+,x,3*exp(x)-2*x-2,k-);%數(shù)值解與真解圖title(數(shù)值解與真解圖);legend(RK4,Euler,真解);xlabel(x);ylabel(y);figure(2)plot(x1,abs(y_r-(3*exp(x1)-2*x1-2),k-);%龍格庫(kù)塔方法的誤差title(龍格庫(kù)塔方法的誤差)xlabel(x);ylabel(Error);figure(3)plot(x,abs(y-(3*exp(x)-2*x-2),r-)%Euler方法的誤差title(Euler方法的誤差)xlabel(x);ylabel(Error);【運(yùn)行結(jié)果如下】4.50040.0350.030.0250020.0150010.0050【結(jié)果分析】:龍格庫(kù)塔方法和Euler方法求解常微分方程都能獲得比較好的數(shù)值解,相比較而言龍格庫(kù)塔方法的數(shù)值解的精度遠(yuǎn)遠(yuǎn)要比Euler方法的數(shù)值解的精度
溫馨提示
- 1. 本站所有資源如無(wú)特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 影響跨境電商增長(zhǎng)的全球貿(mào)易壁壘分析
- 遼寧省遼西重點(diǎn)高中2025屆高三下學(xué)期模擬預(yù)測(cè)試題 政治 含答案
- 統(tǒng)籌推進(jìn)教師教育能力提升的背景意義及必要性
- 白酒行業(yè)創(chuàng)新驅(qū)動(dòng)與轉(zhuǎn)型路徑
- 多元化學(xué)習(xí)模式在語(yǔ)文教學(xué)中的應(yīng)用
- 國(guó)際兒童節(jié)課件4
- 智能健美操設(shè)備的設(shè)計(jì)與應(yīng)用前景
- 新能源與抽水蓄能的綜合利用方案
- 智游新紀(jì)元模板
- 電商節(jié)購(gòu)物金融攻略
- FZT 73013-2017 針織泳裝行業(yè)標(biāo)準(zhǔn)
- 藥品生產(chǎn)監(jiān)督管理辦法培訓(xùn)課件
- 口腔門診納入定點(diǎn)后醫(yī)療保障基金的預(yù)測(cè)性分析報(bào)告
- 血液科疾病常見癥狀護(hù)理培訓(xùn)課件
- 遣送押解方案
- 2024-2025年全國(guó)初中化學(xué)競(jìng)賽試卷及答案
- 廢水處理作業(yè)風(fēng)險(xiǎn)點(diǎn)和安全防范應(yīng)急措施
- 零星工程維修投標(biāo)方案技術(shù)標(biāo)
- 鄉(xiāng)村醫(yī)生從業(yè)管理?xiàng)l例全面解讀
- 《軟件工程》課程設(shè)計(jì)學(xué)校訂購(gòu)系統(tǒng)
- 財(cái)經(jīng)素養(yǎng)知識(shí)考試題庫(kù)(濃縮500題)
評(píng)論
0/150
提交評(píng)論