作者zupo (朝著目標前進)
標題Re: [討論] 聯立一階微分方程處理
時間Sun Apr 5 09:31:08 2020
※ 引述《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://webptt.com/m.aspx?n=bbs/MATLAB/M.1586050275.A.530.html