一階常系數(shù)線性齊次微分方程組的求解_第1頁
已閱讀1頁,還剩1頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、<p>  一階常系數(shù)線性齊次微分方程組的求解</p><p>  【模型準(zhǔn)備】一只蟲子在平面直角坐標(biāo)系內(nèi)爬行. 開始時位于點P0(1, 0)處. 如果知道蟲子在點P(x, y)處沿x軸正向的速率為4x 5y, 沿y軸正向的速率為2x 3y. 如何確定蟲子爬行的軌跡的參數(shù)方程? </p><p>  圖31 蟲子爬行的軌跡</p><p>  【模型假

2、設(shè)】設(shè)t時刻蟲子所處位置的坐標(biāo)為(x(t), y(t)). </p><p>  【模型構(gòu)成】由已知條件和上述假設(shè)可知</p><p>  而且(x(0), y(0)) = (1, 0). </p><p>  現(xiàn)要由此得出蟲子爬行的軌跡的參數(shù)方程. </p><p>  【模型求解】令A(yù) =, 則|EA| == (+1)(2). 可見A的特

3、征值為1 = 1, 2 = 2. </p><p>  (EA)x = 0的一個基礎(chǔ)解系為: 1 = (1, 1)T; </p><p>  (2EA)x = 0的一個基礎(chǔ)解系為: 2 = (5, 2)T. </p><p>  令P = (1, 2), 則P1AP =. </p><p>  記X =, Y =, 并且作線性變換X = PY

4、, 則Y = P1X, </p><p>  = P1= P1AX = P1APY =Y, </p><p><b>  即</b></p><p><b>  =, </b></p><p>  故u = c1et, v = c2e2t, 即Y =. 因而</p><p>

5、  = Y|t =0 = P1X|t =0 ==. </p><p><b>  于是</b></p><p>  Y =, X = PY ==.</p><p>  這就是說, 蟲子爬行的軌跡的參數(shù)方程為</p><p>  如果在Matlab命令窗口輸入以下命令</p><p>  >&

6、gt;ezplot('-2/3*exp(-t)+5/3*exp(2*t)','-2/3*exp(-t)+2/3*exp(2*t)',[0,1]) </p><p>  >> grid on;</p><p>  >> axis([0, 12, 0, 5])</p><p>  Matlab執(zhí)行后得</p

7、><p>  圖32 Matlab繪制的蟲子爬行軌跡</p><p>  【模型分析】從圖32可以看出蟲子爬行的軌跡接近一條直線.</p><p><b>  Matlab實驗題</b></p><p>  一只蟲子在平面直角坐標(biāo)系內(nèi)爬行. 開始時位于點P0(0, 1)處. 如果知道蟲子在點P(x, y)處沿x軸正向的速率為

溫馨提示

  • 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. 眾賞文庫僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論