创建或加载信号

在Live编辑器中，将一个有噪声的心电图(ECG)信号加载到MATLAB®工作区中。图数据。

负载noisyecg情节(noisyECG_withTrend)

心电信号有噪声。有几种噪声源会影响信号，包括运动伪影、高频噪声和电源干扰。交互式设计一个滤波器来去除信号中的噪声。在住编辑器选项卡中,扩大任务列表并选择设计滤波器打开任务。

利用Kaiser窗设计低通FIR滤波器

要去除高频噪声，首先选择一种低通滤波器冷杉筛选并指定订单作为10．频率、幅度和算法参数的可用选项取决于所选滤波器响应类型和滤波器顺序。

低通滤波器从输入信号中去除超出指定阈值的不需要的频率内容。在指定频率参数部分中,选择截止频率(6 db)从频率约束列表。当知道抽样率时，可以进行选择赫兹从频率的单位列表。一个采样率选项出现，您可以从工作空间中的变量中选择一个抽样速率。在本例中，采样速率未知，因此指定的归一化截止频率0.3rad /样品。

对于FIR低通滤波器，在指定级参数部分，您可以指定约束来控制通带波纹和阻带衰减的量。选择通带波纹和阻带衰减从大小限制列表。幅值约束和滤波器顺序也会影响滤波器的过渡宽度。

任务选择一个冷杉contrained最小二乘默认根据指定的频率和震级参数设计算法。保留设计选项的默认设置。

在显示过滤器响应部分中,选择级&阶段而且群时延使设计的滤波器响应可视化。在幅值图中，可以看到阻带的衰减水平为60 dB。群延迟图显示延迟为5个样本，且滤波器为线性相位。

的下面的箭头显示过滤器响应部分，以显示为设计的过滤器生成的代码。您可以在命令行上复制并粘贴代码，以手动编辑筛选器设计规范。

将所设计的滤波器应用于有噪声的心电信号。考虑滤波器引入的延迟，并绘制结果。

负载designedFilterfilteredECG =过滤器(designedFilter noisyECG_withTrend);延迟= grpdelay (designedFilter);mdelay =意味着(延迟);filteredECG (1: mdelay) = [];情节(noisyECG_withTrend (1: end-mdelay))在情节(filteredECG)传说([“原始”，“过滤”)举行从

设计等纹波带阻FIR滤波器

像心电图监护仪这样的医疗设备会受到电磁干扰的影响。电源的工作频率一般为50hz或60hz。在这个例子中，一个60赫兹的正弦信号被添加到一个从MIT-BIH心律失常数据库中提取的心电信号中作为噪声[1]．采样速率为360赫兹。要去除噪音，打开设计滤波器设计了一种最小阶带阻FIR滤波器。将默认筛选器名称更改为bandstop60Hz．

指定频率的单位作为赫兹．要指定抽样速率，请输入一个值或从列表中选择一个抽样速率变量。要出现在列表中，必须在工作空间中保存一个样本速率变量。创建一个变量,fs，令它等于360赫兹,然后选择fs从采样率列表。指定通带和阻带频率值，以衰减55 - 65hz之间的频率为一个10hz的陷波滤波器的中心为60hz。

fs = 360;

设置通带波纹2 (dB)来0．5,增加阻带衰减(dB)来80．

任务默认采用等纹波设计方法。显示滤波器的幅度和相位响应。

您也可以选择过滤信息从显示过滤器响应节以查看有关所设计过滤器的其他详细信息。

负载ecg60Hz进入工作区。mat文件包含原始心电信号加上附加噪声(ecg60)和滤波后的信号(ecgFilt)．将这两个信号画出来以可视化滤波结果。

负载ecg60Hzt = 0:1 / fs:(长度(ecg60) 1) / fs;情节(t) [ecg60 ecgFilt])传说([“原始”；“过滤”])