first_break_point_detection.m

db1小波对任意信号进行分解

%信号的奇异性检测
%两类：一种是信号在某一时刻其幅值发生突变，引起信号的非连续，第一种类型的间断点
%另一种是信号图形光滑，幅值没有突变，但是信号的一阶导数不连续且有突变产生，第二种类型的间断点

%第一种类型间断点检测实例

%装载原始信号，并用wavedec函数对原始信号进行六级分解
clear all
load freqbrk;
s=freqbrk;
lx=length(s);
[c,l]=wavedec(s,6,'db5');
subplot(811);plot(s);
title('用db5小波对信号进行六级分解');Ylabel('s');
%=====================================
%用wrcoef函数从[c,l]中提取信号的第六层低频系数并显示
a6=wrcoef('a',c,l,'db5',6);
subplot(812);plot(a6);Ylabel('a6');
%=====================================
%用wrcoef函数从[c,l]中提取信号的各层高频系数并显示
for i=1:6
    decmp=wrcoef('d',c,l,'db5',7-i);
    subplot(8,1,i+2);
    plot(decmp);
    Ylabel(['d',num2str(7-i)]);
end