电力线干涉心电图

                                  D =负载(“莎拉·库珀ecgNoise.mat”);
                                 
                                  心电图= D.DD2;%的信号
                                 
                                  L =长度(心电图);%信号长度(样本)
                                 
                                  Fs = 300;%采样频率
                                 
                                  Fn = f / 2;%奈奎斯特频率
                                 
                                  t = (0: l - 1) / Fs;%的时间向量(sec)
                                 
                                  图(1)
                                 
                                  情节(t,心电图)
                                 
                                  网格
                                 
                                  轴(3 ylim [0])
                                 
                                  FTEKG = fft(心电图)/ L;%傅里叶变换
                                 
                                  阵线= linspace(0, 1,修复(L / 2) + 1) * Fn;%频率向量
                                 
                                  4 = 1:长度(艘);%指数向量
                                 
                                  图(2)
                                 
                                  semilogy(阵线、abs (FTEKG (Iv)) * 2)
                                 
                                  网格
                                 
                                  轴(7.5 ylim [0])
                                 
                                  Fs = 300;%采样频率(赫兹)
                                 
                                  Fn = f / 2;%奈奎斯特频率
                                 
                                  Wp = [100] / Fn;%通频带(正常)
                                 
                                  Ws = (110 - 0.5) / Fn;%阻带(正常)
                                 
                                  Rp = 10;%通带波纹(dB)
                                 
                                  Rs = 30;%阻带波动(dB)
                                 
                                  [n, Ws] = cheb2ord (Wp、Ws Rp, Rs);%切比雪夫II型秩序
                                 
                                  [b] = cheby2 (n, Rs, Ws);%传递函数系数
                                 
                                  (sos, g) = tf2sos (b);% Second-Order-Section稳定
                                 
                                  EKG_F = filtfilt (sos, g,心电图);
                                 
                                  图(3)
                                 
                                  情节(t, EKG_F)
                                 
                                  网格
                                 
                                  轴(3 ylim [0])
                                 
                                  标题(星水黾的带通滤波心电图)
                                 
                                  包含(的时间(秒))
                                 
                                  ylabel (的振幅(mV ?))

给了一个很好的为您的心电图结果:

注意,在傅里叶变换没有的情节 50 - 60 赫兹电力线干扰。

我亏本来解释异常ST-T段基线上坡。它出现我是病态的,在临床情况下,我将寄给委员会认证的心脏病专家进行审核。(我是一个认证一般内科医生和生物医学工程师)。 Braunwald的心脏病 10日 Ed。2015年,没有讨论它。它可以是一个人工制品,但不知道它是如何记录,我不知道可以生成它。改变滤波器参数对它没有影响。

登录置评。

在回答这个问题。

类别

信号处理信号处理工具箱数字和模拟滤波器数字滤波器设计

找到更多的在数字滤波器设计在帮助中心和文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!