一维提升小波程序

llbpp

小弟初学，想做一维振动信号处理的提升小波程序?可做论坛里提升程序主要都是做图像处理的，不太明白。那位大哥能发个一维的？

yazhoulong

clear all;
close all;
load d:\我的文档\MATLAB\驱动端轴承故障\IR0071.mat   %读数据序列

x=X106_DE_time(1:1000) ;              
fs=12000;                              %采样率
dt=1/fs ;                           %采样间隔时间1/fs
n=length(x);                       %得到序列的长度
t=[0:n-1]*dt;

figure(1)
subplot(311)
plot(x)
xlabel('时间 t')
ylabel('幅值 A')
% 提升小波
lshaar=liftwave('haar')%得到HAAR小波的提升方案
els={'p',[-0.125 0.125],0}%添加ELS到提升方案中
lsnew=addlift(lshaar,els)
%[ca,cd]=lwt(x,lsnew)%进行单层提升小波分解
xdec=lwt(x,lsnew,2)
ca1=lwtcoef('ca',xdec,lsnew,2,1)
ca2=lwtcoef('ca',xdec,lsnew,2,2)
cd1=lwtcoef('cd',xdec,lsnew,2,1)
cd2=lwtcoef('cd',xdec,lsnew,2,2)
subplot(323)
plot(ca1)
title('第一层近似信号')
subplot(324)
plot(ca2)