ode15i 解方程，求大神

中神通

1.7726*x''-5.01*x'-2000*x=

6.52*                              sqrt(x'^2+y'^2)+40*cos100t-40sin100t

    1.7726*y''-5.01*y'  -  2000*y=40sin100t+40cos100t-14.11148

function dy=li9(t,x,xp);

dy=[xp(1)-x(2);...

  xp(2)*1.7726-x(2)*5.01-x(1)*2000-sqrt(x(2)^2+x(4)^2)*6.52+40*cos(t*100)-40*sin(t*100);...

    xp(3)-x(4);...

   x(4)*1.7726-x(4)*5.01-x(3)*2000-sin(100*t)*40-cos(100*t)*40-14.11148] ；

保存li9文件，然后

clear all;

x0=[0.000275,0.000275 ]';

xp0=[0,0]';

tspan=[0,100];

options=odeset('RelTol',le-5);

[t,x]=ode15i('li9',tspan,x0,xp0,options);

figure;

plot(t,x(:,1),'k-');

hold on;

plot(t,x(:,3),'k-');

legend('x 曲线','y 曲线');

xlabel('时间t');

ylabel('转子位移');

中神通

1.7726*x''-5.01*x'-2000*x=6.52* sqrt(x'^2+y'^2)+40*cos100t-40sin100t

1.7726*y''-5.01*y'  -  2000*y=40sin100t+40cos100t-14.11148

这是我的方程为什么matlab老说出错呢

中神通

没人理啊

fanxinga

中神通 发表于 2013-6-9 21:41

                               
登录/注册后可看大图

1.7726*x''-5.01*x'-2000*x=6.52* sqrt(x'^2+y'^2)+40*cos100t-40sin100t

1.7726*y''-5.01*y'  -  2000 ...

；后不必加.. 出错的内容贴上来，便于分析吧

ChaChing

缺少报错内容!
我想高手都懒得复製LZ的问题吧!

牛小贱

中神通 发表于 2013-6-9 21:41

                               
登录/注册后可看大图

1.7726*x''-5.01*x'-2000*x=6.52* sqrt(x'^2+y'^2)+40*cos100t-40sin100t

1.7726*y''-5.01*y'  -  2000 ...

出错在什么地方？能不能把出错的内容copy上去 看一下