作业帮matlab 微分方程求解作图这个方程怎么解,还有图象怎么出来~dN1/dt=[r1*(1-N1/K1)-m*N2]*N1;dN2/dt=[r2*(1-N2/K2)-n*N1]*N2;其中r1=0.01,r2=0.05,K1=1000,K2=1000,m=9e-6,n=4e-5;N1(0)=0,N2(0)=0.
来源：学生作业帮助网 编辑：作业帮 时间：2024/04/28 13:39:08
![matlab 微分方程求解作图这个方程怎么解,还有图象怎么出来~dN1/dt=[r1*(1-N1/K1)-m*N2]*N1;dN2/dt=[r2*(1-N2/K2)-n*N1]*N2;其中r1=0.01,r2=0.05,K1=1000,K2=1000,m=9e-6,n=4e-5;N1(0)=0,N2(0)=0.](/uploads/image/z/1157112-0-2.jpg?t=matlab+%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E6%B1%82%E8%A7%A3%E4%BD%9C%E5%9B%BE%E8%BF%99%E4%B8%AA%E6%96%B9%E7%A8%8B%E6%80%8E%E4%B9%88%E8%A7%A3%2C%E8%BF%98%E6%9C%89%E5%9B%BE%E8%B1%A1%E6%80%8E%E4%B9%88%E5%87%BA%E6%9D%A5%7EdN1%2Fdt%3D%5Br1%2A%281-N1%2FK1%29-m%2AN2%5D%2AN1%3BdN2%2Fdt%3D%5Br2%2A%281-N2%2FK2%29-n%2AN1%5D%2AN2%3B%E5%85%B6%E4%B8%ADr1%3D0.01%2Cr2%3D0.05%2CK1%3D1000%2CK2%3D1000%2Cm%3D9e-6%2Cn%3D4e-5%3BN1%280%29%3D0%2CN2%280%29%3D0.)
matlab 微分方程求解作图这个方程怎么解,还有图象怎么出来~dN1/dt=[r1*(1-N1/K1)-m*N2]*N1;dN2/dt=[r2*(1-N2/K2)-n*N1]*N2;其中r1=0.01,r2=0.05,K1=1000,K2=1000,m=9e-6,n=4e-5;N1(0)=0,N2(0)=0.
matlab 微分方程求解作图
这个方程怎么解,还有图象怎么出来~
dN1/dt=[r1*(1-N1/K1)-m*N2]*N1;
dN2/dt=[r2*(1-N2/K2)-n*N1]*N2;
其中r1=0.01,r2=0.05,K1=1000,K2=1000,m=9e-6,n=4e-5;
N1(0)=0,N2(0)=0.
matlab 微分方程求解作图这个方程怎么解,还有图象怎么出来~dN1/dt=[r1*(1-N1/K1)-m*N2]*N1;dN2/dt=[r2*(1-N2/K2)-n*N1]*N2;其中r1=0.01,r2=0.05,K1=1000,K2=1000,m=9e-6,n=4e-5;N1(0)=0,N2(0)=0.
【1】用解析法
r1=0.01,r2=0.05,K1=1000,K2=1000,m=9e-6,n=4e-5
[N1,N2]=dsolve('DN1=(r1*(1-N1/K1)-m*N2)*N1','DN2=(r2*(1-N2/K2)-n*N1)*N2','N1,N2','N1(0)=0,N2(0)=0')
运行结果：
N1 =0
N2 =0
【2】数值法
function haha
clear all
clc
N0 = [0 0];%初值不为0时,图形有变化
[t,N] = ode45(@f,[0:0.2:1],N0)
plot(t,N)
function dNdt = f(t,N)
r1=0.01;r2=0.05;K1=1000;K2=1000;m=9e-6;n=4e-5;
dN1dt=(r1*(1-N(1)/K1)-m*N(2))*N(1);
dN2dt=(r2*(1-N(2)/K2)-n*N(1))*N(2);
dNdt = [dN1dt;dN2dt];
运行结果：
t =
0
0.2000
0.4000
0.6000
0.8000
1.0000
N =
0 0
0 0
0 0
0 0
0 0
0 0
再检查一下方程或初始条件吧,两种方法结果一样.
如果初值不为0,N的结果不是0.