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(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=';xp0=';tspan=;options=odeset('RelTol',le-5);=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老说出错呢 没人理啊 中神通 发表于 2013-6-9 21:41 static/image/common/back.gif
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 ...
；后不必加.. 出错的内容贴上来，便于分析吧 缺少报错内容!
我想高手都懒得复製LZ的问题吧! 中神通 发表于 2013-6-9 21:41 static/image/common/back.gif
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上去 看一下{:{31}:}
页:
[1]