%本程序?qū)崿F(xiàn)基于小波變換的QRS波的檢測???%
%鄭秀玉??深圳大學(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‘;
plot（ecgdata（1:points））;grid?on;axis?tight;axis（[1points-25]）;
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
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）:3*round（points/4）））+min

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-09-12?19:02??心電r波峰值檢測算法\
?????目錄???????????0??2018-09-12?19:02??心電r波峰值檢測算法\心電信號(hào)R波檢測\
?????文件??????183133??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\A?Real-Time?Microprocessor?QRS?Detector?System?with?a?1-ms?Timing?Accuracy?for?the?Measurement?of?Ambulatory?HRV.pdf
?????目錄???????????0??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\QRSdetect\
?????文件????????4558??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\QRSdetect\100.atr
?????文件?????1950000??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\QRSdetect\100.dat
?????文件?????????143??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\QRSdetect\100.hea
?????文件????????5673??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\QRSdetect\ecgdata.mat
?????文件????????7176??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\QRSdetect\ecgdetect.m
?????文件????????6646??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\QRSdetect\graduatedemo.asv
?????文件????????5452??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\QRSdetect\rddata.m
?????文件??????342760??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\一種高效的QRS波實(shí)時(shí)檢測方法.pdf
?????文件??????214674??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\基于Mexican_hat小波的QRS檢測新方法.pdf
?????文件??????526446??2014-12-08?18:18??心電r波峰值檢測算法\心電信號(hào)R波檢測\基于雙正交樣條小波的QRS波檢測.pdf
?????文件???????????0??2018-09-12?19:02??心電r波峰值檢測算法\心電信號(hào)R波檢測\新建文本文檔.txt