关于MATLAB的解方程绘图问题

关于MATLAB的解方程绘图问题
如何用MATLAB绘出X3与r3的关系图?
K1=11/5;K2=5;K3=15.3;r1=1.1;r2=-0.5;
a12=-0.5;a13=-0.1;a21=0.5;a23=-0.1;a31=-1.43;a32=0;
[X1,X2,X3]=dsolve('DX1=X1*(r1*(1-X1/K1)+a12*X2+a13*X3','DX2=X2*(r2*(1-X2/K2)+a21*X1+a23*X3)','DX3=X3*(r3*(1-X3/K3)+a31*X1+a32*X2','0.33,0.6,1.8','t');
plot(r3,X3,'.');xlabel('r3');ylabel('X3');
这样做哪里错了呢?

没看懂题目

 

你是想变r3然后求出坦哪X3再作图吧

 

但你清启没有给r3的取值范围啊,还是时间t的取值范围是多少?X的初值条件等

 

我随便设一个:r3=0.1,时间范让正码围为10

 

代码:

clc

K1=11/5;K2=5;K3=15.3;r1=1.1;r2=-0.5;
a12=-0.5;a13=-0.1;a21=0.5;a23=-0.1;a31=-1.43;a32=0;
r3=0.1;

f=@(t,X)([X(1)*(r1*(1-X(1)/K1)+a12*X(2)+a13*X(3));X(2)*(r2*(1-X(2)/K2)+a21*X(1)+a23*X(3));X(3)*(r3*(1-X(3)/K3)+a31*X(1)+a32*X(2))]);
ode45(f,[0 10],[0.33 0.6 1.8])

 

 

结果为:

 


← 上一篇: 下一篇: →