加载并绘制时间域内的心电图信号，并实施Q因子为1的陷波滤波器以去除50 Hz频率研究（Matlab代码实现）

 欢迎来到本博客❤️❤️

博主优势：博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

本文目录如下：

目录

1 概述

2 运行结果

3 参考文献

4 Matlab代码实现

---

1 概述

陷波滤波器是一种带阻滤波器，用于抑制特定频率的信号。为了去除50 Hz频率，您可以设计一个带通滤波器，将该频率陷波。研究方法分为五个步骤，分别如下：

1. 设计一个一阶带通滤波器：可以使用巴特沃斯或切比雪夫等滤波器类型。选择适当的截止频率来包含50 Hz频率。

2. 确定滤波器的Q因子：Q因子表示滤波器的带宽与中心频率之比。Q因子为1意味着带宽等于中心频率。对于陷波滤波器，Q因子越大，滤波器的带宽越窄。

3. 根据所选的滤波器类型和Q因子，计算滤波器的参数：根据所选的滤波器设计方法和Q因子，计算滤波器的系数或传递函数。

4. 实施滤波器：将心电图信号传递通过滤波器以去除50 Hz频率。可以使用信号处理工具或编程语言（MATLAB）来实施滤波器。

5. 检查滤波效果：绘制滤波后的心电图信号，并检查50 Hz频率的抑制效果。

2 运行结果

 

 

 

 

部分代码：

%% Task 2-a
figure, subplot 211  % study useage of subplot under help section 
plot(t,ecg); title('Raw ECG Data plotting ')             
xlabel('time')
ylabel('amplitude')
legend(' ORIGINAL ECG SIGNAL')
subplot 212
plot(t1,ecg_notch,'r'); title('Filtered ECG signal ')             
xlabel('time')
ylabel('amplitude')
legend(' Flitered ECG SIGNAL')

%% Task 2-b
figure, subplot 121  % study useage of subplot under help section 
plot(t,ecg); title('Raw ECG Data plotting ')             
xlabel('time')
ylabel('amplitude')
legend(' ORIGINAL ECG SIGNAL')
subplot 122
plot(t1,ecg_notch,'r'); title('Filtered ECG signal ')             
xlabel('time')
ylabel('amplitude')
legend(' Flitered ECG SIGNAL')
%% Task 2-c
figure, plot(t(1:201),ecg(1:201),'r'); title('Data plotting for 0 to 0.804 time frame')             
xlabel('time')
ylabel('amplitude')
hold on
plot(t1(1:201),ecg_notch(1:201),'g');             
legend('ORIGINAL ECG SIGNAL',' Flitered ECG SIGNAL')
hold off

3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]刘东辉,孙晓云,何品纯,等.数字滤波及陷波器在心电图信号检测中的应用[J].工业仪表与自动化装置, 2002(2):3.DOI:10.3969/j.issn.1000-0682.2002.02.016.

[2]江峰.心电信号滤波与QRS波形定位[J].武汉工业学院, 2012.

[3]刘东辉孙晓云何品纯高慧泉何钺.数字滤波及陷波器在心电图信号检测中的应用[J].工业仪表与自动化装置, 2002, 000(002):46-48.

4 Matlab代码实现