intfilt

插FIR滤波器设计

描述

b= intfilt (pα设计的线性相位FIR滤波器,其执行理想使用最近带限内插2*p非零样本,在一个序列中使用时,交错-1连续零每的初始带宽α乘以奈奎斯特频率。返回的过滤器b是否与所使用的相同插补

b= intfilt (ñ“拉格朗日”设计一个FIR滤波器执行ñ交错序列上的三阶拉格朗日多项式插值-1连续零每样本。

例子

全部收缩

设计一个数字内插滤波器由七到上采样的信号,利用该带限方法。指定的0.5和使用“bandlimitedness”因子 2 × 2 样品中的插值。

upfac = 7;阿尔法= 0.5;H1 = intfilt(upfac,2,阿尔法);

当原始信号带宽限制的过滤效果最好α乘以奈奎斯特频率。通过生成200个高斯随机数,并使用40阶FIR低通滤波器对序列进行滤波,产生带限噪声信号。重置随机数生成器以获得可重复的结果。

lowp = FIR1(40,阿尔法);RNG(“默认”)X =滤波器(lowp,1,randn(200,1));

由每对的样品的之间插入零增加的信号的采样率X

XR =上采样(X,upfac);

使用过滤函数来产生内插信号。

y =过滤器(h1 1 xr);

补偿滤波器引入的延迟。绘制原始信号和插值信号。

延迟=平均值(grpdelay(H1));Y(1:延迟)= [];茎(1:upfac:upfac *长度(x)的,x)的保持情节(Y)XLIM([400 700])

intfilt还执行拉格朗日多项式插值。

  • 一阶多项式插值就是线性插值,它是由一个三角形滤波器完成的。

  • 零阶插值是通过一个移动平均滤波器完成的,它类似于采样保持显示的输出。

内插的原始信号和叠加的结果。

H2 = intfilt(upfac,1,“拉格朗日”);Y2 =滤波器(h2,1,XR);Y2(1:地板(平均值(grpdelay(H2))))= [];图(Y2)保持

输入参数

全部收缩

样品数,指定为正整数标量。intfilt设计使用的序列的线性相位FIR滤波器穿插-1连续零每样本。

非零样本数,指定为正整数标量。intfilt设计了一个线性相位FIR滤波器,利用最接近的2*进行带限插值p非零的样本。

过渡带宽的逆量度,指定为一个标量。α是成反比滤波器的过渡带宽,同时也影响到阻带不关心区域的带宽。指定α允许你指定你的输入信号多少奈奎斯特区间占据。这是有益的要插补信号,因为它可以让你增加过渡带宽,不影响更好的阻带衰减插值和结果对于给定的p。如果你设置α1,假设你的信号占据整个奈奎斯特区间。设置α到小于1允许不关心的区域在阻带。例如,如果您的输入占用了一半的奈奎斯特间隔,您可以将alpha设置为0.5。

拉格朗日多项式指定为正整数标量。FIR滤波器执行ñ交错序列上的三阶拉格朗日多项式插值-1连续零每样本。如果两个ñ均为偶数,滤波器的设计不是线性相位。

多项式内插法,指定为“拉格朗日。

输出参数

全部收缩

滤波器系数,返回作为载体。要点b为FIR滤波器的系数。如果α指定,它假设的原始带宽为α乘以奈奎斯特频率。b是长度为2 **p-1。

为了ñ阶拉格朗日多项式插值,b长度(ñ+1)*对于ñ甚至,和长度(ñ+1)*L-1对于ñ奇。

算法

带限方法使用firls设计了一个插值FIR滤波器。多项式方法是利用拉格朗日多项式插值公式对等间距样本进行插值,从而构造出合适的滤波器。这两种滤波器基本上都是低通的,增益为在通带。

扩展功能

之前介绍过的R2006a