discontinuites_detection1.m

MATLAB 小波分析高级技术——书中的一些应用程序

load freqbrk;
x= freqbrk;
% 对信号进行傅里叶变换
f = fft(x,1024);
f = abs(f);
figure(1);
subplot(211);
plot(x);
subplot(212);
plot(f);
% 使用DB6小波进行6层分解
[c,l]=wavedec(x,6,'db6');
figure(2);
subplot(8,1,1);
plot(x);
ylabel('x');
% 对分解的第6层低频系数进行重构
a=wrcoef('a',c,l,'db6',6);
subplot(8,1,2);
plot(a);
ylabel('a6');
for i=1:6
    % 对分解的第6层到第1层的高频系数分别进行重构
    d=wrcoef('d',c,l,'db6',7-i);
    subplot(8,1,i+2);
    plot(d);
    ylabel(['d',num2str(7-i)]);
end