为什么用ode45求解时位移ss和速度vv是零？

zhoutianyue

为什么用ode45求解时位移ss和速度vv是零？是不是计算精度问题？请赐教；
方程
[localimg=640,278]1[/localimg]
程序：

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;



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
    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(@evens,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(@odds,tspan,initial,options,q(n+1));
    end
    tt(:,n+1)=t;                                                                         %取出时间
    dim=size(y);
     ss(:,n+1)=y(:,1);                                                                %取出位移
     vv(:,n+1)=y(:,2);                                                                %取出速度
end

无水1324

简单看了一下，你没有定义odds函数吧？

zhoutianyue

抱歉，odds和evens函数就是程序中的odd和even函数，我发程序时漏掉一个s了，不过把补上后结果仍然是零。