怎么在信号中插入损伤信号？

猴王

我仿真了一个加速度信号，我想在信号中插入一个损伤信号，怎么插入，在什么地方插入，损伤发生在30S，这是我的信号模拟程序，请高手指教：
  function xdot=danziyoudukuangjia(t,x)
F=20*sin(25*t);
  M=31.2;
  C=0.01;
  K=67.51;
xdot=[x(2);(1/M)*(F-C*x(2)-K*x(1))];
执行文件
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=20*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]);%单层框架加速度信号输出


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


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

%提取结构的低频和高频信号
ca3=appcoef(c,l,'bior6.8',5);%提取第三尺度系数的低频
[cd1,cd2,cd3,cd4,cd5]=detcoef(c,l,[1,2,3,4,5]);%提取第一、二、三尺度系数的高频

%重构信号的低频和高频部分
a3=wrcoef('a',c,l,'bior6.8',3);
d1=wrcoef('d',c,l,'bior6.8',1);
d2=wrcoef('d',c,l,'bior6.8',2);
d3=wrcoef('d',c,l,'bior6.8',3);
d4=wrcoef('d',c,l,'bior6.8',4);
d5=wrcoef('d',c,l,'bior6.8',5);

%显示多尺度一维信号的分解结果


subplot(611);plot(t,a3);title('信号的低频细节部分');%AXIS([0 60 -5 5]);

subplot(612);plot(t,d1);title('Detailsd1');%AXIS([0 60 -5 5]);

subplot(613);plot(t,d2);title('Detailsd2');%AXIS([0 60 -5 5]);

subplot(614);plot(t,d3);title('Detailsd3');%AXIS([0 60 -5 5]);

subplot(615);plot(t,d4);title('Detailsd4');%AXIS([0 60 -5 5]);

subplot(616);plot(t,d5);title('Detailsd5');%AXIS([0 60 -5 5]);

[ 本帖最后由 xinyuxf 于 2007-6-23 13:24 编辑 ]

VibrationMaster

首先要定义什么叫损伤信号？ 我从来没有听过这个名词;@P

huangyong87

楼主想必欲生成一个单自由度体系刚度突降型损伤的响应信号吧？在30S处将刚度变小，并保证位移和速度的连续性，解方程是否能够得到呢？直接向完好结构响应加入表示损伤的尖峰信号，恐怕不行。因为，损伤后响应频率也会降低。

猴王

是不是可以这么理解，我如果对速度信号进行小波多尺度分析，然后对其导数进行连续性的分析，即分析其速度函数的导数，就是加速度曲线的连续性，也就是小波分析中的分析第二类间断点的问题，但现在的问题是这样的，我定义损伤，损伤应该怎么来加速到我原来的无损的信号当中去，还请楼上的指导，其实以后的分析不一定分析加速度，我也可以通过小波来分析速度，然后再用小波来分析加速度，这都是可以的，但现在的关键是 损伤我不知道该怎么加到无损信号当中去，或者说我不知道损伤信号是不是单纯的定义一个冲击信号就可以了呢