心电信号的特征提取、分析与处理

心电信号的特征提取、分析与处理*

数据来源：MIT-BIH数据库（可从以下数据中任选两组进行实验）
给出4组不同病例的心电信号数据，分别命名为“100-2-3”，“105-2-3”，“109-2-3”，“111-2-3”，每组数据以“.mat”形式存储。（在文章最后面附带四组数据库的压缩包）
每组数据长度N=2048，采样率fs=360Hz(硬件采集心电信号时，每秒采集360个点。注意设计FIR,IIR时可能用到该采样率。). 即2048点对应时间约为5.69s（）

##内容

（1）谱分析： 取两段心电信号数据（不同病例），对该数据进行频谱分析（幅度谱、相位谱、功率谱）；
（2）相关分析：分别计算两段心电信号的均值、方差、自相关函数与互相关函数；分析两段信号的相干函数曲线
对于相关函数进行循环相关函数与线性相关函数的对比；
（3）数字滤波器设计：
取一段心电信号，添加白噪声。分别作出原始信号、加噪信号的图；作出原始信号、加噪信号的自相关曲线，分析估计心电信号的准周期；取一段心电信号，添加高频噪声（1k-2khz），如高频正弦信号，频率自己选择。结合（1）中得出的结论，即ECG的主要能量分布结果，设计数字滤波器（IIR或FIR），去除高频噪声。（也可直接设计数字滤波器去除基线漂移）要求：对滤波器的截止频率的设置要给出说明；
（4）维纳滤波器去除工频干扰：
取一段心电信号，添加频率为50Hz的高斯白噪声（工频干扰）。设计维纳滤波器，分析滤波结果。
##matlab代码实现

clear
Fs=360;N=2048;
load('109-2-3.mat');
x1=y;
load('111-2-3.mat');
x2=y;
figu