移动平均滤波器的频率响应定量分析

主旨

在生理信号处理中，经常有滤除低频信号的需求，例如在分析心率相关问题时排除掉呼吸等因素造成的基线漂移等低频干扰。在这个场景的工程实践中，移动平均滤波器由于设计实现简单，经常被用到。但是移动平均滤波器的频谱分析缺乏类似于FIR或IIR相关定量资料，在工作中时常凭经验确定阶数，严谨性不足。本文试图采用理论分析+数值计算方法给出移动平均滤波器频谱的一般规律，用于指导日常工作

代码位置

本文中的数值计算使用Matlab实现，代码和输出表格放在如下位置
https://github.com/wangyaobupt/MovingAverageFilter

移动平均滤波器介绍

考虑时间上无限长度的输入信号 $x(n),n=-\infty, .. , 0,1,...\infty$
移动平均滤波器的输出为

$$y(n) = \frac{1}{2m+1}\sum_{k=n-m}^{n+m}x(n-k)$$
其中m为正整数，2m+1称为滤波器阶数

工程上不可能有无限长度信号，对于输入信号$x(n),n=0,1,..N$，当$k < 0$时取