作者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/cn.aspx?n=bbs/MATLAB/M.1586050275.A.530.html