zhoutianyue 发表于 2008-3-20 09:36

用ode45求解微分方程时为什么位移及速度值是0?

我用ode45解方程时位移ss和速度vv值结果均为零,不知为什么,附上方程和程序,请指教

y(1,1)=0;y(1,2)=0;beta=0.4488;w=157.1;

N=ceil(20*pi/beta);   
dim(1)=1;
for n=0:1:N                                                      %N=140(相当于t=0.4s)
    q(n+1)=n;
    if mod(n,2)==0;                                                %判断奇数偶数
      tspan=;
      initial=;
       options = odeset('RelTol',1e-6,'AbsTol',1e-14);
             =ode45(@even,tspan,initial,options,q(n+1));
    else
      tspan=;
      initial=;                           %本次初值为上次计算末值
       options = odeset('RelTol',1e-6,'AbsTol',1e-14);
             =ode45(@odd,tspan,initial,options,q(n+1));
    end
   dim=size(y);
   tt(:,n+1)=t;                                                   %提取时间
   ss(:,n+1)=y(:,1);                                             %提取位移
   vv(:,n+1)=y(:,2);                                             %提取速度
end




function dy=even(t,y,q)
    dy=zeros(2,1);
    cca=0.5;ch=81.2;kca=25392;tao6=1e+5;m=3.24;beta=0.4488;w=157.1;
   
dy(1)=y(2);   
dy(2)=(-(ch*((cos(w*t-q*beta)).^2+(cos(w*t-q*beta+2*beta)).^2+(cos(w*t-q*beta+4*beta)).^2+(cos(w*t-q*beta+6*beta)).^2)+cca)*y(2)...      -(tao6*((cos(w*t-q*beta)).^2+(cos(w*t-q*beta+2*beta)).^2+(cos(w*t-q*beta+4*beta)).^2+(cos(w*t-q*beta+6*beta)).^2)+kca)*y(1))/m;



function dy=odd(t,y,q)
    dy=zeros(2,1);
    cca=0.5;ch=81.2;kca=25392;tao6=1e+5;m=3.24;beta=0.4488;w=157.1;
   
dy(1)=y(2);   
dy(2)=(-(ch*((cos(w*t-(q-1)*beta)).^2+(cos(w*t-(q-1)*beta+2*beta)).^2+(cos(w*t-(q-1)*beta+4*beta)).^2)+cca)*y(2)...      -(tao6*((cos(w*t-(q-1)*beta)).^2+(cos(w*t-(q-1)*beta+2*beta)).^2+(cos(w*t-(q-1)*beta+4*beta)).^2)+kca)*y(1))/m;

zhoutianyue 发表于 2008-3-20 14:02

发一张方程的图片
页: [1]
查看完整版本: 用ode45求解微分方程时为什么位移及速度值是0?