单支撑的wrcoef函数具体实施的意义是什么,还有单支撑函数到底有什么样的作用呢,如果是尺度问题的话,我师兄同样应用的是bior6.8小波函数,在三个尺度上进行分解,就可以得到小波的奇异点,我为什么得不到呢,如果您需要我可以把师兄的论文给您看看,因为现在时间很紧张,所以还要多请您帮忙,关于信号的奇异性我是这么来看的,首先是信号本身就具有奇异性,第二就是一阶或者二阶信号具有奇异性,所以还要请你帮忙啊,以下是我的执行文件
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;
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,3,'bior6.8');%对第一信号进行3尺度一维离散小波分解,采用墨西哥小帽函数
%提取结构的低频和高频信号
ca3=appcoef(c,l,'bior6.8',3);%提取第三尺度系数的低频
[cd1,cd2,cd3]=detcoef(c,l,[1,2,3]);%提取第一、二、三尺度系数的高频
%重构信号的低频和高频部分
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);
%显示多尺度一维信号的分解结果
subplot(411);plot(t,a3);title('信号的低频细节部分');AXIS([0 60 -5 5]);
subplot(412);plot(t,d1);title('Detailsd1');AXIS([0 60 -5 5]);
subplot(413);plot(t,d2);title('Detailsd2');AXIS([0 60 -5 5]);
subplot(414);plot(t,d3);title('Detailsd3');AXIS([0 60 -5 5]); |