作业帮 > 综合 > 作业

matlab 常微分方程的画图问题

来源:学生作业帮 编辑:搜狗做题网作业帮 分类:综合作业 时间:2024/04/29 02:47:15
matlab 常微分方程的画图问题
dA/dt=0.00875+7.5A^2/(0.2^2+A^2)-4.1^(-8)AR-10^(-4)A
dR/dt=0.025+2.5A^5/(0.12^5+A^5)-10^(-4)R
需要绘出A,R与时间的二维图,及A,R间的关系曲线图
不胜感激!
matlab 常微分方程的画图问题
这个你没给初值条件啊 我假设是x0=[0,0]; x(1)是A,x(2)是R
程序如下
首先定义这个方程
function dx = funct(t,x)
dx = zeros(2,1);
dx(1)=0.00875+7.5*x(1)^2/(0.2^2+x(1)^2)-4.1^(-8)*x(1)*x(2)-10^(-4)*x(1);
dx(2)=0.025+2.5*x(1)^5/(0.12^5+x(1)^5)-10^(-4)*x(2);
然后再写M文件:
clear all;
close all;
clc
%x0=[0.0001;0 0.0000001];
x0=[0;0];
[t,x]=ode45('funct',[0,10],x0);
plot(x(:,1),x(:,2))
再问: 谢谢,但如果是要看长时间的曲线变化呢?比如时间段为0-30hour,参考里的峰值大概A在13*10^3, R在2*10^4左右。。没用过matlab,不知道这样说清不清楚,非常感谢~~
再答: 那横轴就是时间t了是吗?
再问: 不好意思,之前没看以为没人回了。。横轴是时间t。。