|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
复制动态方程:5 I$ z$ b4 ^5 D7 q6 k$ I4 x+ r& J: w
function dy=myfun(t,y,a,A,b,B,f,F,c,C,e,E,g,G)
5 x* z! w6 C9 B9 i* \) A1 Wdy=zeros(3,1);
, _/ R6 X3 E$ X Udy(1)=y(1)*(1-y(1))*(y(2)*y(3)*(a-A+b*B)+y(3)*(C-c)-C-y(2)*A);6 t3 |& n4 d _/ X! w) X
dy(2)=y(2)*(1-y(2))*(y(1)*y(3)*(1-b)*B-(1-y(3))*e);
3 R$ X' `' `/ \5 O# y+ @" ndy(3)=y(3)*(1-y(3))*((f-F)+(1-y(1)*y(2))*(g-G)-E);
0 M( O0 K+ X' s( q* B( I
. q) v P: [$ ~; s A* g% M1 @参数: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,。
2 j$ C: D4 _9 Z; v' J+ Q+ Y X1 i
5 B( Y# r. L+ M* |* S% ]我想请问一下如何在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)的关系.
& J" g: z8 a1 P$ ?- a( b
1 G" g. }- m! Y5 s, s2 _8 w |
|