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

zhoutianyue

我用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=[n*beta/w,(n+1)*beta/w];
        initial=[y(dim(1),1);y(dim(1),2)];
       options = odeset('RelTol',1e-6,'AbsTol',1e-14);
             [t,y]=ode45(@even,tspan,initial,options,q(n+1));
    else
        tspan=[n*beta/w,(n+1)*beta/w];
        initial=[y(dim(1),1);y(dim(1),2)];                           %本次初值为上次计算末值
       options = odeset('RelTol',1e-6,'AbsTol',1e-14);
             [t,y]=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

发一张方程的图片