sindy209 发表于 2008-4-11 16:03

求助ode45求解问题！

clear all
lb=inline('','t','x');

p=-1.5*xa(:,1)-xa(:,4).*xa(:,1).*(xa(:,1).^4+2*xa(:,5)+1);
a=-1.5-xa(:,4).*(5*xa(:,1)^4+2*xa(:,5)+1);
b=-xa(:,1)^5-2*xa(:,1).*xa(:,1);
c=-2*xa(:,1).*xa(:,4);
r1=xa(:,1).^2;
r2=xa(:,1).^2+xa(:,2).^2;
q1=xa(:,1).^4+xa(:,1).^2.*2xa(:,5)+xa(:,1).^2-0.001*xa(:,4);
q2=q1+(xa(:,2)-p).^2*(1+a^2)*(2*xa(:,5)+1)+p^2*(r2^2+(a*r1)^2)
gama=-2*p+(-1.5-xa(:,4).*(5*xa(:,1)^4+2*xa(:,5)+1))*xa(:,2)+0.5*(-xa(:,1)^5-2*xa(:,1).*xa(:,1)).*q2+(-2*xa(:,1).*xa(:,4)).*xa(:,5)-xa(:,4).*p*(1+a^2)*(2*xa(:,5)+1)-xa(:,4).*q2*(r2^2+(a*r1)^2);
options=odeset('RelTol',1e-4,'AbsTol',1e-4);
=ode45(lb,,,options);
figure(1)
plot(t2,xa(:,1),'k:')

高手帮忙看下错在那！谢谢了！

[ 本帖最后由 无水1324 于 2008-4-14 18:25 编辑 ]

无水1324 发表于 2008-4-14 18:24

回复 楼主 的帖子

提示的错误是什么？
页: [1]
查看完整版本: 求助ode45求解问题！