约束最小二乘线性相位FIR低通和高通滤波器设计
dp, b = fircls1 (n, ds)
dp, b = fircls1 (n, ds,“高”)
dp, b = fircls1 (n, ds, wt)
dp, b = fircls1 (n, ds, wt,“高”)
dp, b = fircls1 (n, ds, wp, ws, k)
dp, b = fircls1 (n, ds, wp, ws, k,“高”)
dp, b = fircls1 (n, ds,…,”design_flag
”)
dp, b = fircls1 (n, ds)
生成低通FIR滤波器b
,在那里n + 1
为滤波器长度,我们
为0 - 1范围内的归一化截止频率(其中1为奈奎斯特频率),dp
最大通带偏离1(通带纹波),和ds
为从0(阻带纹波)的最大阻带偏差。
dp, b = fircls1 (n, ds,“高”)
生成高通FIR滤波器b
.fircls1
对于高通配置,总是使用一个均匀的过滤器顺序。这是因为对于奇阶,奈奎斯特频率处的频率响应必然为0。如果你指定一个奇值n
,fircls1
把它加1。
dp, b = fircls1 (n, ds, wt)
和
dp, b = fircls1 (n, ds, wt,“高”)
指定一个频率wt
上面(wt
>我们
)或更低(为wt
<我们
)保证滤波器满足给定的频带准则。这将帮助您设计满足通带或阻带边缘要求的滤波器。有四种情况:
低通滤波器:
0
<wt
<我们
<1
:滤波器的振幅在dp
1除以频率范围0
<ω<wt
.
0
<我们
<wt
<1
:滤波器的振幅在ds
在频率范围内为0wt
<ω<1
.
高通滤波:
0
<wt
<我们
<1
:滤波器的振幅在ds
在频率范围内为00
<ω<wt
.
0
<我们
<wt
<1
:滤波器的振幅在dp
1除以频率范围wt
<ω<1
.
dp, b = fircls1 (n, ds, wp, ws, k)
生成低通FIR滤波器b
带加权函数,其中n + 1
为滤波器长度,我们
为归一化截止频率,dp
最大通带偏离1(通带纹波),和ds
为从0(阻带纹波)的最大阻带偏差。wp
是L2权函数的通带边和吗ws
是L2权函数的止带边,在哪里wp
<我们
<ws
.k
为通频带L2误差/阻频带L2误差的比值
dp, b = fircls1 (n, ds, wp, ws, k,“高”)
生成高通FIR滤波器b
带加权函数,其中ws
<我们
<wp
.
dp, b = fircls1 (n, ds,…,”
使您可以监视过滤器的设计,其中design_flag
”)'
design_flag
'
可以
“跟踪”
,用于在设计中使用的设计表的文本显示
“阴谋”
,表示滤波器的幅值、群延迟、零点和极点。在每个迭代步骤更新所有的图。图上的O是新迭代的估计极值,X是前一次迭代的估计极值,其中的极值是滤波器波纹的峰值(最大值和最小值)。只有具有相应的O和X的波纹才相等。
“两个”
,用于文本显示和情节
在设计中非常窄带滤波器用小dp
和ds
,则可能不存在满足规格要求的给定长度的过滤器。
fircls1
采用迭代最小二乘算法得到等纹波响应。该算法是一种多重交换算法,每次迭代使用拉格朗日乘子和库恩-塔克条件。
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页。