matlab 常微分方程的画图问题dA/dt=0.00875+7.5A^2/(0.2^2+A^2)-4.1^(-8)AR-10^(-4)AdR/dt=0.025+2.5A^5/(0.12^5+A^5)-10^(-4)R需要绘出A,R与时间的二维图,及A,R间的关系曲线图不胜感激!

来源:学生作业帮助网 编辑:作业帮 时间:2024/04/28 10:51:19

matlab 常微分方程的画图问题dA/dt=0.00875+7.5A^2/(0.2^2+A^2)-4.1^(-8)AR-10^(-4)AdR/dt=0.025+2.5A^5/(0.12^5+A^5)-10^(-4)R需要绘出A,R与时间的二维图,及A,R间的关系曲线图不胜感激!
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 常微分方程的画图问题dA/dt=0.00875+7.5A^2/(0.2^2+A^2)-4.1^(-8)AR-10^(-4)AdR/dt=0.025+2.5A^5/(0.12^5+A^5)-10^(-4)R需要绘出A,R与时间的二维图,及A,R间的关系曲线图不胜感激!
这个你没给初值条件啊 我假设是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))