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

看板MATLAB作者 (朝著目標前進)時間4年前 (2020/04/05 09:31), 編輯推噓0(000)
留言0則, 0人參與, 最新討論串1/1
※ 引述《Hsinxyzzyx (新)》之銘言: : 問題是這樣的 : http://i.imgur.com/N2e46Am.jpg
: 其中k都是常數 : 然後有給起始值x=500 y=200 : 我知道一般只有x對t的微方 : 用ode45之類的指令就可以處理 : 像這種因為有x和y對t又相互關聯的 : 要怎麼對t作圖啊 : 抱歉小弟很菜只會處理沒有連立的 : 以前沒處理過類似問題 : 請各位大大給點指點 好久沒解,順便複習 本題目因為係數差異頗大,所以多用了一個解題器解  看看是否會因係數差異值大而有所影響 ----------------------------------- function pttex148 [t,x1] = ode45(@xy,[0,20],[500;200]); figure(1) plot(t,x1) legend('x','y') xlabel('t') ylabel('x,y') title('一階聯立ODE') [t,x2] = ode15s(@xy,[0,20],[500;200]); figure(2) plot(t,x2) legend('x','y') xlabel('t') ylabel('x,y') title('一階聯立ODE') function f = xy(t,x) % x(1) = x;f(1) = x'; % x(1) = y;f(2) = y'; k1 = 0.02; k2 = 0.00004; k3 = 0.0004; k4 = 0.04; f(1) = k1.*x(1) - k2.*x(1).*x(2); f(2) = k3.*x(1).*x(2) - k4.*x(2); f = f'; ------------------------------------- 看起來是沒有影響的 -- 1.MATLAB programming 2.ASPEN process simulation package 3.FORTRN programming 4.Advance Engineering Mathematics 5.Process Control Theory 6.Chemical Engineering Basic Theory(Kinetic.thermodynamics.transport) 7.Numerical Method and Analysis 8.MATLAB Toolbox.Simulink system basic design -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 123.18.162.138 (越南) ※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1586050275.A.530.html
文章代碼(AID): #1UYJJZKm (MATLAB)
文章代碼(AID): #1UYJJZKm (MATLAB)