matlab中Filter Design和Simulink仿真的实操与训练2

上期讲了模拟滤波器的设计与Simulink仿真实验，今天讲解下数字滤波器的设计与Simulink仿真实验。
今天是采用三个滤波器将包含随机噪声的信号中分离出两个有用信号，这两个有用信号分别为Sine Wave1（0.7mHz毫赫兹）和Sine Wave2（5uHz微赫兹），为了避免在单位转换中的你晕头转向，这里给出进率1Hz=1000mHz=1000000uHz。这两个信号频率之间相差约140倍，相差很大，所以滤波起来还是有点难度的，下面开始我们的工作。
本实验使用了三个滤波器分别是IIR1、IIR2和IIR3，其中IIR1和IIR3都是单极点低通滤波器（通带截至频率4.09mHz，这是我项目里的参数我直接拿来用了，如果读者很难理解的话直接将它视为4mHz吧，影响不大），IIR2为三级点低通椭圆滤波器（通带截至频率为0.03mHz），这里我们放入两个有用信号分别是Sine Wave1（0.7mHz）和Sine Wave2（5uHz），并且加入一个随机噪声信号，通过滤波器的设置将两个有用信号分离出来。
通过上面的叙述显而易见，IIR1滤波器通带带宽为4.09mHz（<4.09mHz的信号都可以通过），所以两个有用信号Sine Wave1（0.7mHz）和Sine Wave2（5uHz） 都可以通过，但是高频随机噪声（IIR1：>4.09mHz的信号不能通过）就会被滤除。
这样经过第一个滤波器滤除了随机噪声，留下两个有用信号的叠加。如何分离两个有用信号？先滤除高频段信号得到低频段信号，再用原始信号减去已经滤出的低频段有用信号就得到高频段的有用信号了。
具体来说的话首先确定IIR2滤波器通带截止频率，它要求通过两个信号中频率较低的那一个，滤除频率较高的那一个。我们将IIR2滤波器的通带设为0.03mHz（<0.03mHz的信号可通过），这样频率较低的Sine Wave2（5uHz） 信号可以通过滤波器IIR2，而频率较高的信号Sine Wave1（0.7mHz） 就被滤除了。
被滤除了不就没了？先别急。在通过IIR1滤波器的时候我们滤除了随机噪声，得到了两个信号的叠加信号。在通过IIR2滤波器的时候我们得到了两个有用信号中频率较低的那一个