|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
%改写自xjzou的示范编程
clear all
%6維降階得12維
n=12
%Omega为先前定义的6*6系数矩阵
O=[Omega,zeros(6);zeros(6),zeros(6)]
%D1为先前定义的6*6系数矩阵
D=[D1,zeros(6);(-1)*eye(6),zeros(6)]
%q为强制外力
q=inline('[sin(4.14*t);sin(4.14*t);sin(4.14*t);sin(4.14*t);sin(4.14*t);sin(4.14*t);0;0;0;0;0;0]','t')
motion=inline('[z(n+1:end,1); q-O*[z(1);z(2);z(3);z(4);z(5);z(6);z(1);z(2);z(3);z(4);z(5);z(6)]-D*[z(7);z(8);z(9);z(10);z(11);z(12);0;0;0;0;0;0]]',...
't','z','flag','n','O','D','q')
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5])
[t,z]=ode45(motion,[0,20],[0;0;0;0;0;0;0;0;0;0;0;0],options,n,q,D,O)
??? Error using ==> inline.feval
Not enough inputs to inline function.
Error in ==> funfun\private\odearguments at 110
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.
Error in ==> ode45 at 171
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...
我应该已经将要用到的变数系数都定义完成
为何依旧显示如上错误讯息呢
先感谢各位看完 |
|