%本程序?qū)崿F(xiàn)基于小波變換的QRS波的檢測(cè)???%
%鄭秀玉??深圳大學(xué)德州儀器DSPs實(shí)驗(yàn)室???%
%2007年7月2日????????????????????????%

clear?all;
close?all;
points=4096;???????level=4;????sr=360;?
%讀入ECG信號(hào)
load?ecgdata.mat;
ecgdata=mydata‘;
g=importdata（‘yy.json‘）;
ecgdata?=g;
plot（ecgdata（1:points））;grid?on;axis?tight;axis（[1points-0.0010.0015]）;
title（‘ECG信號(hào)‘）;

swa=zeros（4points）;
swd=zeros（4points）;
signal=ecgdata（0*4096+1:1*4096）;

%算小波系數(shù)和尺度系數(shù)
for?i=1:points-3
??swa（1i+3）=1/4*signal（i+3-2^0*0）+3/4*signal（i+3-2^0*1）+3/4*signal（i+3-2^0*2）+1/4*signal（i+3-2^0*3）;
???swd（1i+3）=-1/4*signal（i+3-2^0*0）-3/4*signal（i+3-2^0*1）+3/4*signal（i+3-2^0*2）+1/4*signal（i+3-2^0*3）;
end
j=2;
while?j<=level
???for?i=1:points-24
?????swa（ji+24）=1/4*swa（j-1i+24-2^（j-1）*0）+3/4*swa（j-1i+24-2^（j-1）*1）+3/4*swa（j-1i+24-2^（j-1）*2）+1/4*swa（j-1i+24-2^（j-1）*3）;
?????swd（ji+24）=-1/4*swa（j-1i+24-2^（j-1）*0）-3/4*swa（j-1i+24-2^（j-1）*1）+3/4*swa（j-1i+24-2^（j-1）*2）+1/4*swa（j-1i+24-2^（j-1）*3）;
???end
???j=j+1;
end
%畫出原信號(hào)和尺度系數(shù)，小波系數(shù)
figure;
subplot（level11）;?plot（ecgdata（1:points））;?grid?on;axis?tight;
title（‘ECG信號(hào)及其在j=1234尺度下的尺度系數(shù)及小波系數(shù)‘）;
for?i=1:level
????subplot（level+122*（i）+1）;
????plot（swa（i:））;?axis?tight;grid?on;xlabel（‘time‘）;
????ylabel（strcat（‘a(chǎn)???‘num2str（i）））;
????subplot（level+122*（i）+2）;
????plot（swd（i:））;?axis?tight;grid?on;
????ylabel（strcat（‘d???‘num2str（i）））;
end

%畫出原圖及小波系數(shù)
figure;
subplot（level11）;?plot（real（ecgdata（1:points））‘b‘）;?grid?on;axis?tight;
title（‘ECG信號(hào)及其在j=1234尺度下的小波系數(shù)‘）;
for?i=1:level
????subplot（level+11i+1）;
????plot（swd（i:）‘b‘）;?axis?tight;grid?on;
????ylabel（strcat（‘d???‘num2str（i）））;
end

%**************************************求正負(fù)極大值對(duì)*****************************************%
ddw=zeros（size（swd））;
pddw=ddw;
nddw=ddw;
%小波系數(shù)的大于0的點(diǎn)
posw=swd.*（swd>0）;
%斜率大于0
pdw=（（posw（:1:points-1）-posw（:2:points））<0）;
%正極大值點(diǎn)
pddw（:2:points-1）=（（pdw（:1:points-2）-pdw（:2:points-1））>0）;
%小波系數(shù)小于0的點(diǎn)
negw=swd.*（swd<0）;
ndw=（（negw（:1:points-1）-negw（:2:points））>0）;
%負(fù)極大值點(diǎn)
nddw（:2:points-1）=（（ndw（:1:points-2）-ndw（:2:points-1））>0）;
%或運(yùn)算
ddw=pddw|nddw;
ddw（:1）=1;
ddw（:points）=1;
%求出極值點(diǎn)的值其他點(diǎn)置0
wpeak=ddw.*swd;
wpeak（:1）=wpeak（:1）+1e-10;
wpeak（:points）=wpeak（:points）+1e-10;

%畫出各尺度下極值點(diǎn)
figure;
for?i=1:level
????subplot（level1i）;
????plot（wpeak（i:））;?axis?tight;grid?on;
ylabel（strcat（‘j=???‘num2str（i）））;
end
subplot（411）;
title（‘ECG信號(hào)在j=1234尺度下的小波系數(shù)的模極大值點(diǎn)‘）;

interva2=zeros（1points）;
intervaqs=zeros（1points）;
Mj1=wpeak（1:）;
Mj4=wpeak（3:）;

%畫出尺度3極值點(diǎn)
figure;
plot?（Mj4）;
title（‘尺度3下小波系數(shù)的模極大值點(diǎn)‘）;

posi=Mj4.*（Mj4>0）;
%求正極大值的平均
thposi=（max（posi（1:round（points/4）））+max（posi（round（points/4）:2*round（points/4）））+max（posi（2*round（points/4）:3*round（points/4）））+max（posi（3*round（points/4）:4*round（points/4））））/4;
posi=（posi>thposi/3）;
nega=Mj4.*（Mj4<0）;
%求負(fù)極大值的平均
thnega=（min（nega（1:round（points/4）））+min（nega（round（points/4）:2*round（points/4）））+min（nega（2*round（points/4）:

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

?????文件???????4558??2008-03-30?12:34??檢測(cè)QRS超級(jí)好用\100.atr

?????文件????1950000??2008-03-30?12:34??檢測(cè)QRS超級(jí)好用\100.dat

?????文件????????143??2008-03-30?12:34??檢測(cè)QRS超級(jí)好用\100.hea

?????文件???????5680??2018-03-05?17:32??檢測(cè)QRS超級(jí)好用\ecgdata.mat

?????文件???????7178??2018-03-05?16:52??檢測(cè)QRS超級(jí)好用\ecgdetect.m

?????文件???????6646??2008-03-30?12:34??檢測(cè)QRS超級(jí)好用\graduatedemo.asv

?????文件???????5470??2018-02-28?20:28??檢測(cè)QRS超級(jí)好用\rddata.m

?????文件?????165111??2018-03-05?16:25??檢測(cè)QRS超級(jí)好用\yy.json

?????目錄??????????0??2018-03-05?16:47??檢測(cè)QRS超級(jí)好用

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

??????????????2144786????????????????????9