约束最小二乘FIR多频带滤波器设计
b = fircls (n、f、amp, lo)
fircls (n、f, amp,,,”design_flag
”)
b = fircls (n、f、amp, lo)
生成一个长度n + 1
线性相位FIR滤波器b
.该滤波器的频率-幅度特性与向量给出的特性相匹配f
和amp
:
f
为转换频率的向量,取值范围为0 ~ 1,其中1为奈奎斯特频率。第一点f
必须0
最后一点1
.频率点必须按递增顺序排列。
amp
是描述频率响应的分段常数期望振幅的向量。的长度amp
等于响应中的频带数应该等于长度(f) 1
.
向上
和罗
向量的长度和amp
.它们定义了每个频带的频率响应的上下限。
fircls
对于通频带在奈奎斯特频率(即高通滤波器和带阻滤波器)的配置,总是使用一个均匀的滤波器顺序。这是因为对于奇阶,奈奎斯特频率处的频率响应必然为0。如果你指定一个奇值n
,fircls
把它加1。
fircls (n、f, amp,,,”
使您可以监视过滤器的设计,其中design_flag
”)'
design_flag
'
可以
“跟踪”
,以文本形式显示每个迭代步骤的设计错误。
“阴谋”
,用于显示滤波器的全频带幅度响应和每个子频带幅度响应的放大视图的图集合。在每个迭代步骤更新所有的图。图上的O是新迭代的估计极值,X是前一次迭代的估计极值,其中的极值是滤波器波纹的峰值(最大值和最小值)。只有具有相应的O和X的波纹才相等。
“两个”
,用于文本显示和情节。
请注意
通常,阻带中较低的值将被指定为负值。通过设置罗
等于0
在阻带中,可以得到一个非负的频率响应幅值。这样的滤波器可以被频谱分解得到最小相位滤波器。
fircls
采用迭代最小二乘算法得到等纹波响应。该算法是一种多重交换算法,每次迭代使用拉格朗日乘子和库恩-塔克条件。
Selesnick, i.w., M. Lang和c.s. Burrus。无指定过渡带的FIR滤波器的约束最小二乘设计。1995年声学、语音和信号处理国际会议论文集。1995年第2卷,第1260-1263页。
Selesnick, i.w., M. Lang和c.s. Burrus。无指定过渡带的FIR滤波器的约束最小二乘设计。IEEE®信号处理汇刊.1996年第44卷第8期,1879-1892页。