基于小波变换的方法识别ECG信号的中的PT波程序

该算法，是通过小波变化的方法，识别出了心电信号的中RR波，同时通过求取极值来确定PT波的位置。希望对那些在学习人工智能医疗的这块的有所帮助！

clc;
load ‘ECG_100.MAT‘ ;
E=sig（:2）;
E=E‘;
n=size（E）;
s=E（1:2000）;%取部分数据
%wavelet
[cl]=wavedec（E3‘db5‘）;
ca3=appcoef（cl‘db5‘3）;
cd1=detcoef（cl1）;
cd2=detcoef（cl2）;
cd3=detcoef（cl3）;

%使用stein的无偏似然估计原理进行选择各层的阈值
%cD1cD2cD3为各层小波系数， 
%‘rigrsure’为无偏似然估计阈值类型
thr1=thselect（cd1‘rigrsure‘）;
thr2=thselect（cd2‘rigrsure‘）;
thr3=thselect（cd3‘rigrsure‘）;
%各层的阈值
TR=[thr1thr2thr3];
%‘s‘为软阈值;‘h‘硬阈值。
SORH=‘s‘;
%---------去噪----------------
%XC为去噪后信号
%[CXCLXC]为的小波分解结构
%PERF0和PERF2是恢复和压缩的范数百分比。
%‘lvd‘为允许设置各层的阈值
%‘gbl‘为固定阈值。
%3为阈值的长度
[XCCXCLXCPERF0PERF2]=wdencmp（‘lvd‘E‘db5‘3TRSORH）;
%---------去噪效果衡量（SNR越大效果越好MSE越小越好）
%选取信号的长度
N=n（2）;
x=E;
y=XC;
F=0;
M=0;
for ii=1:N
   m（ii）=（x（ii）-y（ii））^2;
   t（ii）=y（ii）^2;
   f（ii）=t（ii）/m（ii）;
    F=F+

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     240184  2001-05-14 16:43  ECG denoising based on wavelet transform\ECG_100.MAT

     文件       1187  2019-03-14 18:21  ECG denoising based on wavelet transform\ECG_wavelet.m

     文件       9119  2019-03-12 19:59  ECG denoising based on wavelet transform\locationP_QRS_T.m

     目录          0  2019-03-12 19:55  ECG denoising based on wavelet transform

----------- ---------  ---------- -----  ----

               250490                    4