[討論] 聯立一階微分方程處理

[Hsinxyzzyx (新)]

問題是這樣的
http://i.imgur.com/N2e46Am.jpg
其中k都是常數
然後有給起始值x=500   y=200

我知道一般只有x對t的微方
用ode45之類的指令就可以處理

像這種因為有x和y對t又相互關聯的
要怎麼對t作圖啊

抱歉小弟很菜只會處理沒有連立的
以前沒處理過類似問題
請各位大大給點指點

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.112.250.162 (臺灣)
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1584843200.A.AB0.html
※ 編輯: Hsinxyzzyx (140.112.250.162 臺灣), 03/22/2020 10:19:23
※ 編輯: Hsinxyzzyx (140.112.250.162 臺灣), 03/22/2020 10:20:24
→ Hsinxyzzyx: *更 作業用for迴圈無腦加也畫出圖來了 但畢竟是取巧
→ Hsinxyzzyx: 還是想知道這聯立ode怎麼處理
※ 編輯: Hsinxyzzyx (140.112.250.162 臺灣), 03/22/2020 11:07:44
推 Howard1648: http://i.imgur.com/ml9t8xU.jpg
→ Howard1648: http://i.imgur.com/T9zyX06.jpg
→ Howard1648: ~部分，裡面應該是t陣列。[t,y]=ode().;
→ Howard1648: 應該能拿t去作圖，y是矩陣，取你要的column 就好了。
→ Howard1648: http://i.imgur.com/w7LDWQT.jpg http://i.imgur.com/
→ Howard1648: e67FTo1.jpg
→ Howard1648: 重傳一次，方程式我是這樣打的。
→ Howard1648: http://i.imgur.com/vIsCMR5.jpg
推 Howard1648: 像這樣，給你參考
→ Howard1648: dxDt(1)=k1*x1-k2*x1*x2；
→ Howard1648: dxDt(2)=k3*x1*x2-k4*x2；
→ Howard1648: where
→ Howard1648: x1=x
→ Howard1648: x2=y
→ Hsinxyzzyx: 大感謝
※ 編輯: Hsinxyzzyx (140.112.250.162 臺灣), 03/22/2020 16:15:13