马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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))]; |