breakpoint2.m

用MATLAB小波工具箱实现的信号奇异点的检测

load nearbrk;
x=nearbrk;
%使用DB4小波进行2层分解
[c,l]=wavedec(x,2,'db4');
subplot(4,1,1);
plot(x);
ylabel('x');
%对分解的第2层低频系数进行重构
a=wrcoef('a',c,l,'db4',2);
subplot(4,1,2);
plot(a);
ylabel('a2');
for i=1:2
    %对分解的第2层到第1层的高频系数分别进行重构
    d=wrcoef('d',c,l,'db6',3-i);
    subplot(4,1,i+2);
    plot(d);
    ylabel(['d',num2str(3-i)]);
end