基于ARM微处理器和ML2035的低频信号发生器设计
0 引言
数字滤波器作为语音与图象处理、模式识别、雷达信号处理、频谱分析等应用中最基本的处理部件,现已成为最常用的工具之一。它既能满足滤波器对幅度和相位特性的严格要求,又能避免模拟滤波器所无法克服的电压漂移、温度漂移和噪声等问题。而对于具有线性相位特性的滤波问题,
设计时一般都选择FIR滤波器。
相对于窗函数法和频率
设计法,在将理想频率响应和实际频率响应之间的加权逼近误差均匀地分散到滤波器的整个通带和阻带最小化和最大误差这个意义上来说,Chebyshev逼近法可以被视为最佳的
设计准则。
1设计原理
1.1 FIR数字滤波器
对于长度为N、输入为x(n)、输出为y(n)的FIR滤波器,其输出函数可用差分方程表示为:
事实上,具有线性相位的滤波器都具有对称性或反对称性,即单位样本响应可满足条件:
1.2 Chebyshev逼近法
事实上,具有线性相位的滤波器都具有对称性或反对称性,即单位样本响应可满足条件:
1.2 Chebyshev逼近法
(1)线性相位FIR滤波器的四种情况
根据单位样本响应的对称性或反对称性,以及滤波器长度的奇偶性,其线性相位FIR滤波器有以下四种情形:
情形1:单位样本响应具有对称性,即h(n)=h(N-1-n),且N为奇数;
情形2:单位样本响应具有对称性,即h(n)=h(N-1-n),且N为偶数;
情形3:单位样本响应具有反对称性,即h(n)=-h(N-1-n),且N为奇数;
情形4:单位样本响应具有反对称性,即h(n)=-h(N-1-h),且N为偶数。
情形2:单位样本响应具有对称性,即h(n)=h(N-1-n),且N为偶数;
情形3:单位样本响应具有反对称性,即h(n)=-h(N-1-n),且N为奇数;
情形4:单位样本响应具有反对称性,即h(n)=-h(N-1-h),且N为偶数。
(2)误差函数E(ω)
若定义实值理想频率响应Hd(ω)在通带内为1,在阻带内为0;同时定义加权函数W(ω)在通带内为δ2/δ1(δ1为通带波纹,δ1为阻带波纹),阻带内为1。则可将加权逼近误差E(ω)定义为:
如误差函数已知,则Chebyshev逼近只需确定滤波器参数{α(k)},然后使其逼近频带E(ω)上的最大绝对值最小化。即要找到下式的解:
该问题的解法已由Parks和McClellan解决,称之为Remez交换算法。该算法是建立在交错定理的基础上的。图1所示是Remez算法的流程图。
若定义实值理想频率响应Hd(ω)在通带内为1,在阻带内为0;同时定义加权函数W(ω)在通带内为δ2/δ1(δ1为通带波纹,δ1为阻带波纹),阻带内为1。则可将加权逼近误差E(ω)定义为:
如误差函数已知,则Chebyshev逼近只需确定滤波器参数{α(k)},然后使其逼近频带E(ω)上的最大绝对值最小化。即要找到下式的解:
该问题的解法已由Parks和McClellan解决,称之为Remez交换算法。该算法是建立在交错定理的基础上的。图1所示是Remez算法的流程图。
评论