[討論] x,y對t的微分方程

看板MATLAB作者 (新)時間4年前 (2020/03/22 10:13), 4年前編輯推噓2(2015)
留言17則, 2人參與, 4年前最新討論串1/1
問題是這樣的 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

03/22 11:06, 4年前 , 1F
*更 作業用for迴圈無腦加也畫出圖來了 但畢竟是取巧
03/22 11:06, 1F

03/22 11:06, 4年前 , 2F
還是想知道這聯立ode怎麼處理
03/22 11:06, 2F
※ 編輯: Hsinxyzzyx (140.112.250.162 臺灣), 03/22/2020 11:07:44

03/22 12:59, 4年前 , 3F

03/22 12:59, 4年前 , 4F

03/22 13:00, 4年前 , 5F
~部分,裡面應該是t陣列。[t,y]=ode().;
03/22 13:00, 5F

03/22 13:01, 4年前 , 6F
應該能拿t去作圖,y是矩陣,取你要的column 就好了。
03/22 13:01, 6F

03/22 13:02, 4年前 , 7F

03/22 13:02, 4年前 , 8F
e67FTo1.jpg
03/22 13:02, 8F

03/22 13:03, 4年前 , 9F
重傳一次,方程式我是這樣打的。
03/22 13:03, 9F

03/22 13:03, 4年前 , 10F

03/22 13:17, 4年前 , 11F
像這樣,給你參考
03/22 13:17, 11F

03/22 13:17, 4年前 , 12F
dxDt(1)=k1*x1-k2*x1*x2;
03/22 13:17, 12F

03/22 13:17, 4年前 , 13F
dxDt(2)=k3*x1*x2-k4*x2;
03/22 13:17, 13F

03/22 13:17, 4年前 , 14F
where
03/22 13:17, 14F

03/22 13:17, 4年前 , 15F
x1=x
03/22 13:17, 15F

03/22 13:17, 4年前 , 16F
x2=y
03/22 13:17, 16F

03/22 16:14, 4年前 , 17F
大感謝
03/22 16:14, 17F
※ 編輯: Hsinxyzzyx (140.112.250.162 臺灣), 03/22/2020 16:15:13
文章代碼(AID): #1UTid0gm (MATLAB)
文章代碼(AID): #1UTid0gm (MATLAB)