声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1391|回复: 0

[综合] 在通过小波变换检测结构损伤的时候,怎么定义在某个时刻的损伤???

[复制链接]
发表于 2007-6-9 10:51 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我编译了一个单层框架结构的程序,用的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))];
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-30 14:30 , Processed in 0.088324 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表