【心电信号】基于matlab心电信号PTT+HRV+PRV【含Matlab源码 1551期】

原创

Matlab科研 2021-12-01 11:05:22 ©著作权

文章标签 matlab 图像处理开发语言小波变换参考文献 文章分类 后端开发

©著作权归作者所有：来自51CTO博客作者Matlab科研的原创作品，请联系作者获取转载授权，否则将追究法律责任

一、获取代码方式

获取代码方式1：

完整代码已上传我的资源：【心电信号】基于matlab心电信号PTT+HRV+PRV【含Matlab源码 1551期】

二、心电信号简介

0 引言

心电信号是人类最早研究的生物信号之一, 相比其他生物信号更易于检测, 且具有直观的规律。心电图的准确分析对心脏病的及早治疗有重大的意义。人体是一个复杂精密的系统, 有许多不可抗的外界因素, 得到纯净的心电信号非常困难。可以采用神经网络算法去除心电信号的噪声, 但这种方法存在训练难度大、耗时长的缺点。小波变换在处理非线性、非平稳且奇异点较多的信号时具有一定的优越性, 近年来许多学者使用其对心电信号进行研究。

1 心电信号简介

心电信号由以下几个波段组成, 一个典型的心电图如图1所示。

【心电信号】基于matlab心电信号PTT+HRV+PRV【含Matlab源码 1551期】_图像处理

图1 典型心电图

(1) P波:反映心房肌在除极过程中的电位变化过程;

(2) P-R间期:反映的是激动从窦房结通过房室交界区到心室肌开始除极的时限;

(3) QRS波群:反映心室肌除极过程的电位变化;

(4) T波:代表心室肌复极过程中所引起的电位变化;

(5) S-T段:从QRS波群终点到达T波起点间的一段水平线[2];

(6) Q-T间期:心室从除极到复极的时间[3];

(7) U波:代表动作电位的后电位。

由于心电信号十分微弱, 且低频, 极易受到干扰, 不同的干扰源的噪声虽是随机的, 但来自同一个干扰源的噪声往往具有同一类特征。分析干扰的来源, 针对不同的来源使用合适的处理方法, 是数据采集重点考虑的一个问题。常见干扰有3种: (1) 工频干扰; (2) 基线漂移; (3) 肌电干扰。其中已经证明小波变换在抑制心电信号的工频干扰方面具有较大优势。具体噪声频带如表1所示。

表1 心电信号以及主要噪声频带

【心电信号】基于matlab心电信号PTT+HRV+PRV【含Matlab源码 1551期】_小波变换_02

三、部分源代码

clc;
clear all;
close all;
x=load('day2_0917.txt');
%% ECG signal
y=x(1:95000,1); % ECG signal
figure,plot(y);
title('ECG signal');
xlabel('time');
ylabel('amplitude');
hold on

%% PPG signal
z=x(200:95000,2); % PPG signal
plot(z,'r');
title('PPG signal');
xlabel('time');
ylabel('amplitude');
%% peak detection of ECG
j=1;
n=length(y);
for i=2:n-1
    if y(i)> y(i-1) && y(i)>= y(i+1) && y(i)> 0.45*max(y)
       val(j)= y(i);
       pos(j)=i;
       j=j+1;
     end
end
ecg_peaks=j-1;
ecg_pos=pos./1000;
plot(pos,val,'*r');
title('ECG peak');
%% peak detection of PPG
m=1;
n=length(z);
for i=2:n-1
    if z(i)> z(i-1) && z(i)>= z(i+1) && z(i)> 0.45*max(z)
       val(m)= z(i);
       pos1(m)=i;
       m=m+1;
     end
end
ppg_peaks=m-1;
ppg_pos=pos1./1000;
ppg_val=val;
plot(pos1,val,'*g');
title('ECG & PPG signal');
legend('ECG signal','PPG signal');

%% HRV
j=1;
for i=1:ecg_peaks-1
    e(j)= ecg_pos(i+1)-ecg_pos(i);% gives RR interval
    j=j+1;

end 
hr=60./mean(e); % 60/ mean of RR interval

hrv= (60./e); % 60/ each RR interval
figure,stairs(hrv);
title('HRV');
xlabel('samples');
ylabel('hrv');

%% PRV

end 
pr=60./mean(f); 
prv= 60./f; 
figure,stairs(prv);
title('PRV');
xlabel('samples');
ylabel('prv');
%% PTT
ptt=(ppg_pos-ecg_pos);
figure,stairs(ptt);
title('PTT');
xlabel('ptt');
ylabel('time');

%% notch detection

%%moving average filter
av=smooth(z,150);

%%differentiation
p=100*diff(av,1); % (signal,order of differentiation), 100 to amplify the signal

%%finding peak of the notch on the differentiated signal
np=1;  % notch peak
m=length(p); 
for i=2:m-1
    if p(i)> p(i-1) && p(i)>= p(i+1) 
       val(np)= p(i);
       pos(np)=i;
       np=np+1;
    end
end