一维提升小波程序
小弟初学,想做一维振动信号处理的提升小波程序?可做论坛里提升程序主要都是做图像处理的,不太明白。那位大哥能发个一维的? 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=*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)
%=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)
页:
[1]