二、生理信号处理 ——1.心电信号(含Matlab代码及数据)

本文适合快速了解心电信号，并能够进行数据的滤波处理。

一. 心电数据预处理(消除工频干扰、基线漂移)

* 心电数据及rdmat函数见文章底部

1. 导入心电数据

## 心电图导入及读取
clc;
[TIME,M,Fs,siginfo]=rdmat('100m');# 通过读取函数ramat对心电图进行处理
Fs=1500;# 采样频率
plot(TIME,M);

2. 选取前1000个数据点进行后续处理

ecg = M(1:1000);
TIME = TIME(1:1000);
plot(TIME,ecg)

3. 去除心电信号中的50Hz频率干扰(国内通常为50Hz)——陷波滤波器

## Butterworth 陷波滤波器-去除工频干扰
# 设计滤波器
d = designfilt('bandstopiir','FilterOrder',2, ...'HalfPowerFrequency1',49,'HalfPowerFrequency2',51, ...'DesignMethod','butter','SampleRate',Fs);
#应用滤波器去除50Hz干扰
ecg_50 = filtfilt(d,ecg);

4. 去除心电信号中的低频信号(5Hz以下)——消除基线漂移

##bandpass滤波器-解决基线漂移
fmaxd_1=5;# 截止频率为5Hz
fmaxn_1=fmaxd_1/(Fs/2);
[B,A]=butter(1,fmaxn_1,'low');
ecg_low=filtfilt(B,A,ecg_50);# 通过5Hz低通滤波器的信号
ecg1=ecg_50-ecg_low; # 减去5Hz低频信号
plot(TIME,ecg1)
xlabel('t(s)');
ylabel('mv');
title('期望信号');

5. 添加随机噪声

* wgn函数：随机添加白噪声

##添加随机噪声信号
Noise_White = (0.1*wgn(1,length(TIME),2))';
plot(Noise_White)
title('噪声信号');

6. 噪声干扰后信号

##噪声干扰后信号
Mix_Signal = ecg1 + Noise_White;
plot(TIME,Mix_Signal)
xlabel('t(s)');
ylabel('mv');
title('噪声干扰后信号');

二. 滤波处理(维纳滤波)

* h为维纳滤波运行后得到的数值，运行下文中维纳滤波器设计部分可得：

一、信号处理 ——3.维纳滤波（含matlab代码)_蛮有趣的_的博客-CSDN博客

##维纳滤波
Signal_Filter = filter(h,1,Mix_Signal);# 将输入信号通过维纳滤波器
figure(1)
plot(TIME,Signal_Filter)
xlabel('t(s)');
ylabel('mv');
title('维纳滤波后信号');

三. 均方误差

* 引入噪声后信号相对原信号的统计均方误差:

mse1 = 0.0153

* 滤波后的信号相对原信号的统计均方误差:

mse2 = 0.0069

应用维纳滤波后，信号更接近于原始信号。

* 心电数据及rdmat函数

链接：https://download.csdn.net/download/qq_41958946/86841684

心电数据原始来源：MIT-BIH Arrhythmia Database（mitdb）

二、生理信号处理 ——1.心电信号(含Matlab代码及数据)相关推荐

【数字信号处理】基于DFT的滤波系列3之插值滤波(含MATLAB代码)
四.基于DFT的(理想)滤波例2:一个"警告" "理想DFT滤波器"虽然简单.有效,但可能会导致意想不到的问题.在博客 [数字信号处理]基于DFT的滤波系列2 ...
ResNet50模型识别二维化的心电信号——以MIT-BIH心律失常数据库为例
ResNet50模型识别二维化的心电信号--以MIT-BIH心律失常数据库为例一维信号二维化格拉米角场实现步骤马尔可夫变迁场实现步骤代码实现 ResNet50模型介绍代码实现总结问题 ...
COI实验室技能：常见的图像增强算法（含MATLAB代码）
COI实验室技能:常见的图像增强算法(含MATLAB代码) 成像系统采集到的图片有时不一定具有很好的图像质量,往往需要进行一些图像增强的操作.本文主要整理了几种常见的图像增强算法,包括:同态滤波. ...
浣熊优化算法(COA)(含MATLAB代码)
先做一个声明:文章是由我的个人公众号中的推送直接复制粘贴而来,因此对智能优化算法感兴趣的朋友,可关注我的个人公众号:启发式算法讨论.我会不定期在公众号里分享不同的智能优化算法,经典的,或者是近几年提出 ...
matlab去除大跳变,一种快速灰度校正算法（处理亮度不均等情况）（含MATLAB代码）...
一种快速灰度校正算法(处理亮度不均等情况)(含MATLAB代码) 一种快速灰度校正算法(处理亮度不均等情况)(含MATLAB代码) 文章目录前言一.MATLAB代码二.结果示例总结前言方法 ...
回声消除AEC算法（含Matlab代码）
基于自适应滤波器的回声消除AEC算法(含Matlab代码) 摘要自适应滤波器声学回波抵消 AEC算法解析 LMS算法 NLMS算法 VSNLMS算法 APLMS算法 LMS-Newton算法 PF ...
四维空间的二维线框投影可视化（附matlab代码）
四维空间的二维线框投影可视化(附matlab代码) 1 三维空间在2维屏幕上的投影 1.1平行投影 1.2透视投影 2 四维空间在2维屏幕上的投影 2.1 四维空间与三维空间的一些区别 2.2 四维空 ...
灰狼优化(GWO)算法(含MATLAB代码)
先做一个声明:文章是由我的个人公众号中的推送直接复制粘贴而来,因此对智能优化算法感兴趣的朋友,可关注我的个人公众号:启发式算法讨论. 文献[1]中的灰狼优化(Grey Wolf Optimizer, ...
灰狼优化(GWO)算法的8种中文变体(含MATLAB代码)
先做一个声明:文章是由我的个人公众号中的推送直接复制粘贴而来,因此对智能优化算法感兴趣的朋友,可关注我的个人公众号:启发式算法讨论.我会不定期在公众号里分享不同的智能优化算法,经典的,或者是近几年提出 ...