%--------------------------------------------------------------------------
%?exa130203.m??例13.2.3:?利用小波變換檢測(cè)R波并給出模極大曲線
%--------------------------------------------------------------------------
clear?all;

ecgdata=load（‘ecg.txt‘）;
%plot（ecgdata）;
points=4000;%截取的點(diǎn)數(shù)
range=2:1:32;%尺度范圍
sr=360;%抽樣率
signal=ecgdata（1:points）‘;
figuresubplot（311）
plot（signal）axis?tightgrid?on;

%?對(duì)心電信號(hào)作小波變換，選‘db1‘為分析用小波
wf=‘db1‘;
subplot（312）
ccfs?=?cwt（signalrangewf）;?
imagesc（ccfs）;
%title（‘Continuous?Transform?coefficients.‘）?
colormap（gray（128））;?
%ylabel（‘Scale‘）
swd=ccfs;
ddw=zeros（size（swd））;

%?對(duì)ECG信號(hào)進(jìn)行分段
pddw=ddw;
nddw=ddw;
posw=swd.*（swd>0）;
pdw=（（posw（:1:points-1）-posw（:2:points））<0）;
pddw（:2:points-1）=（（pdw（:1:points-2）-pdw（:2:points-1））>0）;
negw=swd.*（swd<0）;
ndw=（（negw（:1:points-1）-negw（:2:points））>0）;
nddw（:2:points-1）=（（ndw（:1:points-2）-ndw（:2:points-1））>0）;
ddw=pddw|nddw;
ddw（:1）=1;
ddw（:points）=1;
wpeak=ddw.*swd;
wpeak（:1）=wpeak（:1

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????48002??2003-12-03?10:21??心電圖小波\ECGrbojiance\ecg.txt

?????文件???????2756??2004-01-29?17:16??心電圖小波\ECGrbojiance\exa130203.m

?????目錄??????????0??2006-04-21?14:30??心電圖小波\ECGrbojiance

?????目錄??????????0??2009-05-09?22:39??心電圖小波

-----------?---------??----------?-----??----

????????????????50976????????????????????5