Question

matlab7.0解常微分方程组并绘图

麻烦问下~我想画dn_1=(a*p^2+b*(1-p)^2).*n_1+(1-n_1).*n_1+c_2*n_2
dn_2=(c*q^2+d*(1-q)^2).*n_2+(a_2-n_2).*n_2+c_1*n_1
dp=p.*(1-p).*(a*p-b*(1-p))+c_2*(q-p).*n_2./n_1
dq=q.*(1-q).*(c*q-d*(1-q))+c_1*(p-q).*n_1./n_2
其中n1,n2的图。a=-1;b=-1;c=-5;d=-5;c_1=0.2;c_2=1;编程如下：
function dx=dot(t,x)
dx=[
(a*x(3)^2+b*(1-x(3))^2).*x(1)+(1-x(1)).*x(1)+c_2*x(2);
(c*x(4)^2+d*(1-x(4))^2).*x(2)+(a_2-x(2)).*x(2)+c_1*x(1);
x(3).*(1-x(3)).*(a*x(3)-b*(1-x(3)))+c_2*(x(4)-x(3)).*x(2)./x(1);
x(4).*(1-x(4)).*(c*x(4)-d*(1-x(4)))+c_1*(x(3)-x(4)).*x(1)./x(2)
];

>> a=-1;b=-1;c=-5;d=-5;c_1=0.2;c_2=1;
x0=[100;50;0.5;0.5];
[t,x]=ode45(@dot,[0,100],x0);
plot(x(:,1),x(:,2));
但是运行不出来，哪里出错了？麻烦了

Answer 1

问题

1.(c*x(4)^2+d*(1-x(4))^2).*x(2)+(a_2-x(2)).*x(2)+c_1*x(1); a_2没定义

2 a=-1;b=-1;c=-5;d=-5;c_1=0.2;c_2=1; 并没有传递给ode函数dx

可以这样写

 a=-1;b=-1;c=-5;d=-5;c_1=0.2;c_2=1; 
odefun=@(t,x)[ 
(a*x(3)^2+b*(1-x(3))^2).*x(1)+(1-x(1)).*x(1)+c_2*x(2); 
(c*x(4)^2+d*(1-x(4))^2).*x(2)+(c_2-x(2)).*x(2)+c_1*x(1); 
x(3).*(1-x(3)).*(a*x(3)-b*(1-x(3)))+c_2*(x(4)-x(3)).*x(2)./x(1); 
x(4).*(1-x(4)).*(c*x(4)-d*(1-x(4)))+c_1*(x(3)-x(4)).*x(1)./x(2) 
]; 
x0=[100;50;0.5;0.5]; 
[t,x]=ode45(odefun,[0,100],x0); 
plot(x(:,1),x(:,2));