7-2第七章FIR设计-2概述.ppt

例1：用矩形窗设计设计FIR低通滤波器： N=21，;例2：用Hanning窗设计线性相位FIR高通滤波器， N=51，;给定指标窗函数设计法;目的：设计出h(n);高通：;窗函数法设计步骤;例3：设计一个数字FIR低通滤波器，其技术指标： ;FIR滤波器设计完毕的标识是：h(n) 求出。;注意： 1、理想滤波器的单位冲激响应需移位(N-1)/2，不管N是奇数还是偶数，以保证对称，保证线性相位。其线性相位为： 用的是第一类线性相位条件，即偶对称：Φ(ω)=-(N-1)/2 2、窗函数本身均定义在0~N-1之间，所以，窗函数直接乘以理想滤波器函数上，不用移位! 3、第一类线性相位数字高通滤波器、数字带阻滤波器，N必须为奇数。若计算出N=M（偶数），取N=M+1。 4、理想滤波器的频域表示, 必须写上相位(不是1): 这里 ;;当z=-1时，H(z)=0;例：因果线性相位FIR系统，单位脉冲响应h[n]是实序列，N=7， h[n]= -h[6-n]，已知系统函数的一个零点是0.5ejπ/4，写出系统函数。 ;四种类型广义线性相位FIR滤波器的应用限制;作业： 7-4 7-14（4）并交换ωp和ωs设计高通滤波器 讲解： 7-6;设计步骤 按阻带衰减要求选定窗函数 计算过渡带，选较小的一个 计算N： A与所选窗有关 计算两个理想低通滤波器的截止频率 写出两个理想截止频率的低通滤波器的N点单位抽样响应 理想带通滤波器的单位抽样响应为上述两个低通滤波器之差 生成所选窗的窗函数 用窗函数与理想带通滤波器相乘，得到实际带通滤波器的单位抽样响应h(n) 计算h(n)的幅频响应，用函数freqz(hd,[1],w) 验证Rp、As是否满足指标，并做调整，记录实际的通带衰减和阻带衰减 绘出带通滤波器的单位抽样响应图和幅频特性图 See lab3.m ;DTMF;频率的选择是为了避免谐波。;Max and max instead of comprison.;检波电路;半波整流;全波整流 直流稳压电源……;另一种方法用Goertzel算法;;要求： 按一定采样率仿真0-9双音多频信号，播放拨号音，显示时域波形、频域波形。 任意给定某个双音多频信号，通过Goertzel算法（或其他算法）识别出拨号数字，并显示出来。 考虑到线路质量，在双音多频信号上加入噪声，再??识别。 这里涉及到滤波、DFT、系统响应、比较大小的编程。 能做到哪步就做到哪步，仔细在实验报告上写明你的实验过程。;　　DTMF系统可以用模拟手段或数字手段实现。当用数字手段实现时，可以通过软件来模拟合成DTMF信号。在检测端，可以用8个数字滤波器检测DTMF信号， 也可以用其软件实现信号的检测。数字信号的产生可以用查表法实现。电话音频信号在进行数字处理时，其采样周期为125 μs（采样频率为8 kHz）， 在这段时间内，也可以用计算两个频率的正弦波信号并将它们相加的方法来实现。考虑到DSP芯片的存储容量有限，显然用增加计算复杂度来换取减少存储容量的方法更好些。检测端收到数据后， 对数据进行数字处理以检测出接收的信号。这里的检测可用滤波器组实现，也可用FFT算法实现, 当频率成分较少时，前者更有效。 ;2. Goertzel算法 N点数据的DFT公式为; 我们知道，相位因子WNk具有周期性，且WN－kＮ 。利用这一点， 在式（8-1）的两边同时乘以它， 得到如下表达式： ;　　这个滤波器在n=N处的输出就是DFT在频率点ωk=2πk/N的值，就此可以得到该滤波器的系统函数为 ;　　对于每一个k，W-kN为一个固定值，它可以预先计算出来并作为常数项使用, 但该值一般为复数，给计算带来不便。通过对滤波器系统函数的变换： ;　　这就是Goertzel算法。每输入一个样值，分别用迭代法计算出各个vk(n)，每次计算只需进行两次实数乘和两次实数加。 当n=N时，用式（8-3）计算出各yk(n)的值，这些值就是所要求的各点的DFT值。 　　另外，考虑到最后对信号进行判断时需要幅度值|X(k)|或幅度的平方值|X(k)|2，因此对式（8-3）进行变换： ;　　3. DTMF信号产生与检测算法的实现 　　在检测信号时，为避免信号前后的0值带来的对DFT计算的影响，同时考虑到 DTMF 信号的高频和低频分量的最小间隔为73 Hz， 故选取每200点为一段， 分别用Goertzel算法计算每一段的8个频率点的DFT的平方值，然后设定一个阈值对得到的值进行判断（若某一频率点的幅度平方值大于阈