小波分析如何确定时间序列的奇异性？如何通过matlab实现？

(1)小波模极大值重构MATLAB代码_天天_新浪博客/s/blog _ 6 c 00 b 0e 30100 u 1 S5 . html

功能

);

signal = ECG data(1:points)'+offset；

%信号的小波变换，给出了总体情况和详细的分步波形。

[swa，swd] =

swt(信号、电平、Lo_D、Hi _ D)；

图；

支线剧情(等级，1，1)；

plot(真实(信号))；网格打开；轴紧；

对于i=1:level

支线剧情(等级+1，2，2 *(I)+1)；

plot(swa(i，)；轴

紧绷；网格打开；xlabel(' time ')；

ylabel(strcat('a

，num 2 str(I)))；

支线剧情(等级+1，2，2*(一)+2)；

plot(swd(i，))；轴

紧绷；网格打开；

ylabel(strcat('d '，num 2 str(I)))；

结束

找出小波变换的模极大值及其位置的%。

ddw =零(大小(swd))；

pddw = ddw

nddw = ddw

posw=swd。*(社署& gt0);

pdw=((posw(:，1:点-1)-posw(:，2:点))& lt0);

pddw(:，2:points-1)=((pdw(:，1:points-2)-pdw(:，2:points-1))& gt；0);

negw=swd。*(社会福利署& lt0);

ndw=((negw(:，1:点-1)-negw(:，2:点))& gt0);

nddw(:，2:点-1)=((ndw(:，1:点-2)-ndw(:，2:点-1))& gt；0);

ddw = pddw | nddw

ddw(:，1)= 1；

ddw(:，点数)= 1；

wpeak=ddw。*社署；

wpeak(:，1)=wpeak(:，1)+1e-10；

wpeak(:，点数)=wpeak(:，点数)+1e-10；

%给出了各阶段小波变换模极大值的波形。

图；

对于i=1:level

支线剧情(等级，1，I)；

plot(wpeak(i，)；轴紧；格子

开；

ylabel(strcat('j= '，num 2 str(I)))；

结束

注意:运行此程序时，必须添加要处理的信号，即程序的红色部分。