matlab问题requires an initial condition vector of length 4

lihaitao123 · 发表于 2011-5-4 21:29

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

function du=CQ(t,u)
global delta1 delta2 alpha1 xi1 alpha2 xi2  beta2 f omega
du=[u(2);
      -delta1*u(1)-alpha1*u(1)^3-xi1*u(2)+alpha2*u(3)+alpha2/2*u(1)*cos(2*omega*t)-beta2*u(4)+1/2*beta2*u(2)-1/2*beta2*u(2)*cos(2*omega*t)+f*sin(omega*t);
      u(4);
      -delta2*u(3)-xi2*u(4)+delta2*u(1)*sin(omega*t)+xi2*u(2)*sin(omega*t)];
u0=[0,0,0,0];
  [t,u]=ode45(@CQ,tspan,u0,[])
  plot(t,u(:,1))
报错是
??? Error using ==> odearguments at 117
Solving CQ requires an initial condition vector of length 4.如何纠正呀。敝人草鸟，望大虾指导！！

meiyongyuandeze · 发表于 2011-5-4 22:10

u0=[0;0;0;0];

复制代码

要把你的初始条件写成4行的列向量

lihaitao123 · 发表于 2011-5-4 22:25

改过u0=[0;0;0;0];之后，报错还是??? Error using ==> odearguments at 117
Solving CQ requires an initial condition vector of length 4.

meiyongyuandeze · 发表于 2011-5-4 22:58

回复 4 # lihaitao123 的帖子

>> u0=[0.00 0 0 0]';
>> tspan=[0 1];
>> [t,u]=ode45(@CQ,tspan,u0,[])

复制代码

这样没问题了，我试过了已经！

lihaitao123 · 发表于 2011-5-5 09:29

回复 5 # meiyongyuandeze 的帖子

谢谢了，但是
>> u0=[0.00 0 0 0];
>> tspan=[0 1];
>> [t,u]=ode45(@CQ,tspan,u0,[])
??? Error using ==> odearguments at 117
Solving CQ requires an initial condition vector of length 4.

meiyongyuandeze · 发表于 2011-5-5 09:35

不知道怎么回事，我自己运行的没问题！

gyq274706322 · 发表于 2011-9-20 08:48

我也有同样的问题，知道怎么回事吗？

ChaChing · 发表于 2011-9-21 00:31

也不知道LZ是否解决了! 好奇试了下!
问题不在u0'或u0! 问题在主程序global的宣告

账号		自动登录	找回密码
密码			我要加入

[编程技巧] matlab问题requires an initial condition vector of length 4

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。