MATLAB在实时信号处理中的作用:性能仿真,设计合理可能的算法;产生仿真系数;与实际DSP结果进行比较、验证。MATLAB不适用于实时信号处理,因为MATLAB是解释型语言,程序运行起来太慢,并且MATLAB要在计算机上才能运行,因此计算机是不适合于实时信号处理的。
MATLAB中的数字信号处理工具:MATLA中有很多用于信号处理的函数,还有一些用于信号处理的工具,例如FDAtool ,Blockset Fixed point , Embedded Target for TI C6000 DSP等等。
滤波器设计
一、题目
滤波器设计,要求通带分别为50%和20%两种,波纹1%,阻带65db,对过渡带没有要求。各自分别用IIR和FIR设计,且比较系数分别取浮点、16位、12位定点的差别。
二、题目分析
采用FDAtool设计为低通滤波器,采样频率设置为4000Hz。
50%通带时,通带截止频率为Fpass=1000Hz,阻带截止频率为Fstop=1100Hz,过渡带设计为100Hz。 通带波纹1%,即通带最大衰减为Apass=-20lg(1-1%)=0.087,留一定余量取0.06。阻带衰减为Astop=65dB,取68dB。
通带为20%时,通带截止频率为400Hz, 阻带截止频率为Fstop=500Hz,过渡带设计为100Hz。 通带波纹1%,即通带最大衰减为Apass=-20lg(1-1%)=0.087,留一定余量取0.06。阻带衰减为Astop=65dB,取68dB。计算精度默认为双精度。
IIR选择切比雪夫II型,FIR选择等波纹法。
三、设计过程
1.IIR滤波器设计参数
通带50%
通带截止频率:Fpass=1000Hz
阻带截止频率:Fstop=1100Hz
通带最大衰减:Apass=0.06
阻带最小衰减:Astop=68dB
- 双精度
![4f62c7b3c1db7f3b10936521f89e1678.png](https://i-blog.csdnimg.cn/blog_migrate/e73d6ab94ce0eae960bbbb7f2fdb3748.jpeg)
- 定点12位
![e8ebe61fcc5539f4d1677e7b0b1a5bda.png](https://i-blog.csdnimg.cn/blog_migrate/6cb8275726cebb41f42ce813d937009e.jpeg)
- 定点16位
![48240810148f1cde1286296dba44f66e.png](https://i-blog.csdnimg.cn/blog_migrate/d041198a9111104a36ecfb592bd42a18.jpeg)
分析:50%的通带,定点和浮点型对滤波器的性能并没有多大的影响,原因是他们的极点都保持在单位圆内部,如下图所示。
![8d59eac9384265902dc3cb83ee78b36e.png](https://i-blog.csdnimg.cn/blog_migrate/8de7f7bb970bf3fcb57f3cf8ff8c75b5.png)
双精度
![82d315fa27dde056899db04f3cd7e757.png](https://i-blog.csdnimg.cn/blog_migrate/b976c40652a9fd9e9f4eda67f47235f3.png)
16位定点
通带20%
通带截止频率:Fpass=400Hz
阻带截止频率:Fstop=500Hz
通带最大衰减:Apass=0.06
阻带最小衰减:Astop=68dB
(1)双精度
![05b5cf352db4db9fe57ad549c10d7180.png](https://i-blog.csdnimg.cn/blog_migrate/f73a03940f8d43eac5ef11eb6acd29ff.jpeg)
(2)定点12位
![c54a3d7e9478f39ee2ab620642edb226.png](https://i-blog.csdnimg.cn/blog_migrate/09c6099f57806433c204a725f430a128.jpeg)
(3)定点16位
![81e234152bdf3d164ee788e51c1ce65f.png](https://i-blog.csdnimg.cn/blog_migrate/2fafa865e6ae4b3a24fecc7974d24f27.jpeg)
分析:20%的通带与50%的通带相比,20%所需要的滤波器阶数要高一些。定点和浮点型对滤波器的性能并没有多大的影响,原因是他们的极点都保持在单位圆内部,如下图所示。
![47d7099eba55790f03e2f19c6daae53a.png](https://i-blog.csdnimg.cn/blog_migrate/2f3c1cc85dbc883fc778c9cd1a23f564.png)
双精度极点分布
![dc21a9cfedbbd1f255524c816a0d573e.png](https://i-blog.csdnimg.cn/blog_migrate/060134741d08ca782913641cc60876ee.png)
定点16位极点分布
2.FIR滤波器设计参数
通带50%
通带截止频率:Fpass=1000Hz
阻带截止频率:Fstop=1100Hz
通带最大衰减:Apass=0.06
阻带最小衰减:Astop=68dB
(1)双精度
![bd0cc7471fa4ff99e6ba39b0dcd0d4d1.png](https://i-blog.csdnimg.cn/blog_migrate/e9a27856cbfe33646152c8c21869a428.jpeg)
(2)定点12位
![50f5cfd7a6c996a246d813eee9d7f3f6.png](https://i-blog.csdnimg.cn/blog_migrate/acfd25e44f3c63dd50d79e8abb5d0d01.jpeg)
(3)定点16位
![951caccb13849d5bb6c86feedb82384e.png](https://i-blog.csdnimg.cn/blog_migrate/cd8078d98e1f07f35f4cf54ed7de88d9.jpeg)
分析:滤波器达到相同的性能,FIR的阶数要比IIR的阶数要高很多。12位量化效果非常差,在阻带部分已经不满足我们需要的最小衰减了。16位的效果稍微好一点,但是还是有一些误差。
通带20%
通带截止频率:Fpass=400Hz
阻带截止频率:Fstop=500Hz
通带最大衰减:Apass=0.06
阻带最小衰减:Astop=68dB
- 双精度
![91c7feb9a28237f2f919de609335c255.png](https://i-blog.csdnimg.cn/blog_migrate/883ca0a4ceecbf88fced199768df0b86.jpeg)
- 定点12位
![daa1f78ef66c1c621538f56d8862d5c9.png](https://i-blog.csdnimg.cn/blog_migrate/2e2528410cbbd80da754770bf40ad9c4.jpeg)
- 定点16位
![f19d0f882b56eea49fb49cba87b241c6.png](https://i-blog.csdnimg.cn/blog_migrate/13640883664bee0f2af8e6c3569144ef.jpeg)
分析:20%的定点12位效果也是非常差,阻带最小衰减已经不满足我们的要求了,定点16位勉强能满足我们的要求,但是还是有略微的误差,原因是量化具有误差。