下載本文檔
版權(quán)說(shuō)明:本文檔由用戶(hù)提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、實(shí)驗(yàn)五用matlab求解常微分方程1. 微分方程的概念未知的函數(shù)以及它的某些階的導(dǎo)數(shù)連同自變量都由一已知方程聯(lián)系在一起的方程稱(chēng)為 微分方程。如果未知函數(shù)是一元函數(shù),稱(chēng)為常微分方程0常微分方程的一般形式為F(心八八嚴(yán))=0如果未知函數(shù)是多元函數(shù),成為偏微分方程。聯(lián)系一些未知函數(shù)的一組微分方程組稱(chēng)為 微分方程組微分方程中出現(xiàn)的未知函數(shù)的導(dǎo)數(shù)的最高階解數(shù)稱(chēng)為微分方程的階。若方程中 未知函數(shù)及貝各階導(dǎo)數(shù)都是一次的,稱(chēng)為線(xiàn)性常微分方程,一般表示為yg + (f+ +(r)y+a 八(t)y = b(t)若上式中的系數(shù)( = 1,2,均與f無(wú)關(guān),稱(chēng)之為常系數(shù)。2. 常微分方程的解析解dy=y +1有些微分
2、方程可宜接通過(guò)積分求解例如,一解常系數(shù)常微分方程力可化為dy=(itJy+1,兩邊積分可得通解為一1苴中C為任意常數(shù)有些常微分方程可用一些技巧,如分離變量法,積分因子法,常數(shù)變異法,降階法等可化為可積分的方程而求得解析解 線(xiàn)性常微分方程的解滿(mǎn)足疊加原理,從而他們的求解可歸結(jié)為求一個(gè)特解和柑應(yīng)齊次微 分方程的通解一階變系數(shù)線(xiàn)性微分方程總可用這一思路求得顯式解。高階線(xiàn)性常系數(shù)微分 方程可用特征根法求得柑應(yīng)齊次微分方程的基本解,再用常數(shù)變異法求特解。一階常微分方程與髙階微分方程可以互化,已給一個(gè)階方程宀Twy,嚴(yán))設(shè)兒=”匕=,兒=ySi,可將上式化為一階方程組*1 -*2兒T =兒.兒=/(人片*
3、2,兒)反過(guò)來(lái),在許多情況下 一階微分方程組也可化為髙階方程。所以一階微分方程組與高 階常微分方程的理論與方法在許多方而是相通的,一階常系數(shù)線(xiàn)性微分方程組也可用特征根 法求解。3.微分方程的數(shù)值解法除常系數(shù)線(xiàn)性微分方程可用特征根法求解,少數(shù)特殊方程可用初等積分法求解外,大 部分微分方程無(wú)限世界,應(yīng)用中主要依靠數(shù)值解法.考虎一階常微分方程初值問(wèn)題卩(/) = /(人曲)0 v/s=dsolve(,Df=f+g,jDg=g-f,Jf (0)=1,Jg(0)=l,) y (s f jl%s 是一個(gè)結(jié)構(gòu)3implif y (s g)結(jié)果為ans =exp(t)*cos(i)+exp(l)*sin(l)
4、ans =-exp(t)*sin(t)+exp(i)*cost)例2求解微分方程yJ-y+r + h y(0) = 1,先求解析解,再求數(shù)值解,并進(jìn)行比較。由clear ;s=dsolve(,Dy=-y+t+l,Jy(0)=1,Jt,)3implif y (s)可得解析解為=+廣1下而再求其數(shù)值解先編寫(xiě)M文件funS.m%M 函數(shù) funS.mfunction f=fun8(t,y)f=-y+t+l;再用命令clear ; close ; t=0 : 0.1:1;y=t+exp (-t); plot (t, y); %化解析解的圖形hold on; %保留已經(jīng)畫(huà)好的圖形,如果下面再畫(huà)圖,兩個(gè)圖
5、形和并在-起 t,y=ode45*fun8 0,1,1);plot(t,y, *ro*); %iSj數(shù)值解圖形,用紅色小圈畫(huà)xlabel ( P , ylabel ( 了)結(jié)果見(jiàn)圖71圖1661解析解與數(shù)值解 由圖1661可見(jiàn),解析解和數(shù)值解吻合得很好。例3求方程mg sin 3. &(0) =O(0) = 0的數(shù)值解不妨取/ = hg=9&(0) = 15 則上而方程可化為0= 98sm0 6*(0) = 15, 6(0) = 0先看看有沒(méi)有解析解運(yùn)行MATLAB代碼clear ;3=d3olve (* D2y=9.8 *3in (y) *y(0)=15 *Dy(0)=0)3implif y
6、 (s)知原方程沒(méi)有解析解下而求數(shù)值解令門(mén)可將原方程化為如下方程組yi-vi y2=9.8sin(y|)X (0) = 15,2 (0) = 0建立M文件fun9m如下%M 文件 fun9 mfunction f=fun9(t,y)f=y(2), 9.8*3inyl)%f 向量必須為列向量運(yùn)行MATLAB代碼clear; close;t,y=ode45(* funS *, 0,10, 15,0);plot (t, y (: z 1); %iSJ &隨時(shí)間變化圖,y (:2)則衣示&的值xlabel ( P , ylabelyl J 結(jié)果見(jiàn)圖7216.51615.515910圖7.2數(shù)值解由圖72可見(jiàn),&隨時(shí)間f周期變化。L求下列微分方程的解析解y + 2j/ 3卩=界” y-3y = 20 sinx(4)y + ay = sinx (a 0) 刃-宀1 = Q(5)才必+ 2(5?_期2妙=0, y(6)y + /-i-y=GQSX, 7PS-OCT)+ y = / + os卷 7h + 2尹+才=0 7y=2$-0_32=12.求方程(l + x-)y= 2xy y(0) = L y(0) = 3的解析
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 2023三年級(jí)語(yǔ)文下冊(cè) 第一單元 2 燕子配套說(shuō)課稿 新人教版
- 2024-2025學(xué)年高中語(yǔ)文 名著導(dǎo)讀 莎士比亞戲劇說(shuō)課稿 新人教版必修4
- 9古詩(shī)三首清明說(shuō)課稿2023-2024學(xué)年統(tǒng)編版語(yǔ)文三年級(jí)下冊(cè)
- Unit 4 Natural Disasters Reading for Writing 說(shuō)課稿-2024-2025學(xué)年高中英語(yǔ)人教版(2019)必修第一冊(cè)
- Unit 2 lconic Attractions Learning About Language (1)說(shuō)課稿 2023-2024學(xué)年高中英語(yǔ)人教版選擇性第四冊(cè)
- 2025主體信用評(píng)級(jí)合同
- 2025吊頂勞務(wù)承包合同
- 19《夜宿山寺》(說(shuō)課稿)2024-2025學(xué)年部編版語(yǔ)文二年級(jí)上冊(cè)
- 2024-2025學(xué)年高中生物 第一章 人體的內(nèi)環(huán)境與穩(wěn)態(tài) 專(zhuān)題1.2 內(nèi)環(huán)境穩(wěn)態(tài)的重要性說(shuō)課稿(基礎(chǔ)版)新人教版必修3001
- 7《壓歲錢(qián)的使用與思考》(說(shuō)課稿)-2023-2024學(xué)年四年級(jí)下冊(cè)綜合實(shí)踐活動(dòng)長(zhǎng)春版
- 放射科護(hù)理常規(guī)
- 2024年大宗貿(mào)易合作共贏協(xié)議書(shū)模板
- 初中數(shù)學(xué)教學(xué)經(jīng)驗(yàn)分享
- 新聞?dòng)浾咦C600道考試題-附標(biāo)準(zhǔn)答案
- 2024年公開(kāi)招聘人員報(bào)名資格審查表
- TSG ZF001-2006《安全閥安全技術(shù)監(jiān)察規(guī)程》
- 長(zhǎng)螺旋鉆孔壓灌樁工程勞務(wù)清包合同(范本)
- 中考語(yǔ)文二輪復(fù)習(xí):記敘文閱讀物象的作用(含練習(xí)題及答案)
- 老年外科患者圍手術(shù)期營(yíng)養(yǎng)支持中國(guó)專(zhuān)家共識(shí)(2024版)
- 子宮畸形的超聲診斷
- 2024年1月高考適應(yīng)性測(cè)試“九省聯(lián)考”數(shù)學(xué) 試題(學(xué)生版+解析版)
評(píng)論
0/150
提交評(píng)論