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

matlab问题requires an initial condition vector of length 4

function du=CQ(t,u)
global delta1 delta2 alpha1 xi1 alpha2 xi2beta2 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=;
=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=;要把你的初始条件写成4行的列向量

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

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

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

回复 4 # lihaitao123 的帖子

>> u0=';
>> tspan=;
>> =ode45(@CQ,tspan,u0,[])这样没问题了,我试过了已经!

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

回复 5 # meiyongyuandeze 的帖子

谢谢了,但是
>> u0=;
>> tspan=;
>> =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的宣告
页: [1]
查看完整版本: matlab问题requires an initial condition vector of length 4