[工学]9-10第7章FIR设计.ppt

实验二、有限冲击响应滤波器（FIR）算法实验。 滤波器性能要求： 通带截止频率10kHz，阻带截止频率22kHz，阻带 衰减75dB，采样频率50kHz。 设计步骤： 1）首先计算过渡带宽度 2）选择窗函数 3）确定窗函数参数 4）确定滤波器系数 * 令m=N-1/2-n * 了解线性相位FIR滤波器的各种幅频特性可以 * * * 第七章 FIR数字滤波器的设计方法 优点：含有零极点，可以利用模拟滤波器的设计， 且模拟滤波器设计有大量图表，设计简单。 缺点：相位的非线性。 IIR滤波器： FIR滤波器： 优点：严格线性相位；滤波器总是稳定的；可以利用FFT 快速处理； 缺点：取得很好的衰减特性，FIR滤波器阶次要比IIR高。 7.1 FIR滤波器满足线性相位的条件 其线性相位满足： 单位冲激响应满足偶对称： 单位冲激响应满足奇对称： 其线性相位满足： 或 7.2 线性相位FIR滤波器的频率特性 系统函数： 频率特性： 所以 则幅频特性和相频特性按照奇对称和偶对称不同有两种情况： 偶对称时： 奇对称时： 频率特性为： 情况1： 由h(n)的奇偶两种对称性以及N等于奇偶不同情况，可以有四种频率特性与之对应。 N为奇数时，有一个中间项，其余两两对称，因此： 其中 特 点： 对 ω＝0， π，2π 呈 偶 对 称。 clear; clc; omega = linspace(0, 2*pi, 1000); h = [4 -1 5 3 9 3 5 -1 4]; or h = [2 4 -2 -1 -2 -1 -2 4 2]; N = length(h); a = zeros(1, (N-1)/2+1); a(1) = h((N-1)/2+1); 仿 真 程 序 ： for i=2:(N-1)/2+1 a(i) = 2*h((N-1)/2-(i-2)); end Homega = zeros(1, length(omega)); for n =1 : (N-1)/2+1 Homega(:) = Homega(:) + a(n) * cos(n*omega(:)) ; end plot(omega, Homega); 情况2： 其中 频率特性： N为奇数时，中间项为0，其余两两对称，因此： clear; clc; omega = linspace(0, 2*pi, 1000); h = [2 4 -2 -1 0 1 2 -4 -2]; or h = [2 -2 3 -3 0 3 -3 2 -2]; N = length(h); c = zeros(1, (N-1)/2); 仿 真 程 序 ： 特 点：当ω＝0，π，2π时， ； 当ω＝0，π， 2π 时， 成 奇 对 称； 固定的 相移 for i=1:(N-1)/2 c(i) = 2*h((N-1)/2-i+1); end Homega = zeros(1, length(omega)); for n =1 : (N-1)/2 Homega(:) = Homega(:) + c(n) * sin(n*omega(:)) ; end plot(omega, Homega); 情况3： 其中 频率特性： N为偶数时，不存在中间项，其余两两对称，因此： 特 点：当ω=π时， ， 是奇对称； 当ω=0，2π 时， 是偶 对 称 clear; clc; omega = linspace(0, 2*pi, 1000); h = [-3 -4 -1 -2 4 4 -2 -1 -4 -3]; or h = [3 -1 2 4 5 5 4 2 -1 3]; N = length(h); b = zeros(1, N/2); 仿 真 程 序 ： for i=1:N/2 b(i) = 2*h(N/2-i+1); end Homega = zeros(1, length(omega)); for n =1 : N/2 Homega(:) = Homega(:) + b(n) * cos((n-0.5)*omega(:)) ; end plot(omega, Homega); 情况4： 其中 频