程序问题

猴王

执行文件，这个程序主要是对一个单层框架的加速度信号进行的墨西哥小帽的分解，总感觉不对，请高手指教：
clc;clear all;
close all;
t0=0;
   tf=60;
   x0=[0,0]';

   [t,x]=ode23('danziyoudukuangjia',[0:0.005:60],x0);
   d=x(:,1);
   v=x(:,2);
   
t=0:0.005:60;

F=0.02*sin(25*t);
M=31.2;
C=0.01;
K=67.51;
s1=(1/M)*(F-C*v.'-K*d.');
s=s1+randn(1,12001);

subplot(311);plot(t,s);
AXIS([0 60 -10 10]);%单层框架加速度信号输出


%下面进行离散的单尺度小波变换并生成，各尺度上的信号


[ca1,cd1]=dwt(s,'bior6.8');%对第一信号进行单尺度一维离散小波分解，采用墨西哥小帽函数




%重构信号的低频和高频部分
a1=upcoef('a','ca1','bior6.8',1,12001);
d1=upcoef('d','cd1','bior6.8',1,12001);
%a3=wrcoef('a',c,l,'db6',3);
%d1=wrcoef('d',c,l,'db6',1);
%d2=wrcoef('d',c,l,'db6',2);
%d3=wrcoef('d',c,l,'db6',3);
%显示多尺度一维信号的分解结果

subplot(312);plot(a1);
title('信号的低频细节部分');

subplot(313);plot(d1);
title('d1细节部分');


m文件：
  function xdot=danziyoudukuangjia(t,x)
F=0.02*sin(25*t);
  M=31.2;
  C=0.01;
  K=67.51;
xdot=[x(2);(1/M)*(F-C*x(2)-K*x(1))];