在通过小波变换检测结构损伤的时候，怎么定义在某个时刻的损伤？？？

猴王

我编译了一个单层框架结构的程序，用的bior6.8进行分解的，加入了噪声，现在要定义一个30s时刻的损伤，但通过分解并没有达到我预期的效果，请高手执教以下是我的执行文件和m文件：
执行文件：
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*0.95;%这里就是用来定义损伤了5%的K
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*0.95;%这里就是用来定义损伤了5%的K
xdot=[x(2);(1/M)*(F-C*x(2)-K*x(1))];