breakpoint1.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