sindy209 发表于 2008-4-9 16:33

程序运行结果有错求助

clear all

%the non-zero point is (sqrt(35*3),sqrt(35*3),-35),(-sqrt(35*3),-sqrt(35*3),-35)

ch=inline...
('','t','y');
options=odeset('RelTol',1e-4,'AbsTol',1e-4);
=ode45(ch,,,options);
figure(1)
%plot the system states x,y,z
plot(t,ya(:,1)+sqrt(35*3),'k-',t,ya(:,2)+sqrt(35*3),'k-.',t,ya(:,3)-35,'k:')
xlabel('t/s');
ylabel('x,y,z');


结果有错，望高人指点！

[ 本帖最后由 eight 于 2008-4-9 16:35 编辑 ]

eight 发表于 2008-4-9 16:34

原帖由 sindy209 于 2008-4-9 16:33 发表 http://www.chinavib.com/forum/images/common/back.gif
clear all

%the non-zero point is (sqrt(35*3),sqrt(35*3),-35),(-sqrt(35*3),-sqrt(35*3),-35)

ch=inline...
('新人发帖前先看看置顶的会员守则，你这样提问对解答毫无帮助

donkeyxu 发表于 2008-4-9 17:36

改下区间就可以了：

clear all

%the non-zero point is (sqrt(35*3),sqrt(35*3),-35),(-sqrt(35*3),-sqrt(35*3),-35)

ch=inline...
('','t','y');
options=odeset('RelTol',1e-4,'AbsTol',1e-4);
=ode45(ch,[-0.05 0.05],,options);
figure(1)
%plot the system states x,y,z
plot(t,ya(:,1)+sqrt(35*3),'k-',t,ya(:,2)+sqrt(35*3),'k-.',t,ya(:,3)-35,'k:')
xlabel('t/s');
ylabel('x,y,z');
页: [1]
查看完整版本: 程序运行结果有错求助