小波分析 检测 跳变点

%wavlet break detection
%清除以前的数据
clear;
close all;
%载入原始信号的波形数据
load freqbrk;%幅度明显跳变，第一类跳变点
%显示数据的基本信号
whos;
figure;
plot(freqbrk);
%自定义坐标轴
xlabel('时间'); ylabel('幅值');
%自定义标题
title('频率突变信号');

%对信号进行傅里叶变换
f=fft(freqbrk);
figure;
%显示处理后的信号图像
plot(abs(f));
%自定义标题
title('傅里叶变换后的信号示意图');
%对原始信号进行三层小波分解
[d,a]=wavedec(freqbrk,3,'db5');
%重构第三层近似系数
a3=wrcoef('a',d,a,'db5',3);
%重构第1层到第3层细节系数
d3=wrcoef('d',d,a,'db5',3);
d2=wrcoef('d',d,a,'db5',2);
d1=wrcoef('d',d,a,'db5',1);

figure;
subplot(411); plot(a3); ylabel('近似信号 a3'); %显示各层小波系数
subplot(412); plot(d3); ylabel('近似信号 d3'); %显示各层小波系数
subplot(413); plot(d2); ylabel('