Question

画一个微分方程组的x，y，z分别随a的取值变化的图像

function dy=myfun(t,y,a,A,b,B,f,F,c,C,e,E,g,G)
dy=zeros(3,1);
dy(1)=y(1)*(1-y(1))*(y(2)*y(3)*(a-A+b*B)+y(3)*(C-c)-C-y(2)*A);
dy(2)=y(2)*(1-y(2))*(y(1)*y(3)*(1-b)*B-(1-y(3))*e);
dy(3)=y(3)*(1-y(3))*((f-F)+(1-y(1)*y(2))*(g-G)-E);

参数：A=2,b=0.5,B=2.5,c=2.2,C=3,e=3.6,E=2.5,f=3.8,F=3,g=2.5,G=1.8，其中a的变化范围在[1,8]初值y(1)=0.3，y(2)=0.5，y(3)=0.4,。

我想请问一下如何在MATLAB上画图，使得X轴是参数a，Y轴分别是y(1)，y(2)，y(3)。简单来说就是描述参数a在区间[1,8]上与y(1)，y(2)，y(3)的关系，而不是时间t与y(1)，y(2)，y(3)的关系.

Answer 1

你的微分方程组应该这样写才比较合理。即

function dy=myfun(a,y)

A=2,b=0.5,B=2.5,c=2.2,C=3,e=3.6,E=2.5,f=3.8,F=3,g=2.5,G=1.8

dy=zeros(3,1);

dy(1)=y(1)*(1-y(1))*(y(2)*y(3)*(a-A+b*B)+y(3)*(C-c)-C-y(2)*A);

dy(2)=y(2)*(1-y(2))*(y(1)*y(3)*(1-b)*B-(1-y(3))*e);

dy(3)=y(3)*(1-y(3))*((f-F)+(1-y(1)*y(2))*(g-G)-E);

end

然后用ode（）函数来求解。求得结果的图形如下。

追问

上述描绘的图形是这样的，是不是固定t的某个值，然后再看y(1)、y(2)、y(3)与参数a的关系？如果是的话要怎么做，不是的话又该怎么做？麻烦了。

追答

计算后，取某个t值，然后绘出a一y1（a），。。。图像。

Answer 2

这样的话，你把t取一个固定的值，就可以话和a的图像了

追问

请问t需要如何固定，试了很多遍都不行