在通过小波变换检测结构损伤的时候，怎么定义在某个时刻的损伤？？？
我编译了一个单层框架结构的程序，用的bior6.8进行分解的，加入了噪声，现在要定义一个30s时刻的损伤，但通过分解并没有达到我预期的效果，请高手执教以下是我的执行文件和m文件：执行文件：
clc;clear all;
close all;
t0=0;
tf=60;
x0=';
=ode23('danziyoudukuangjia',,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();%单层框架加速度信号输出
%下面进行离散的单尺度小波变换并生成，各尺度上的信号
=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=;
页:
[1]