主要内容

集成滤波器

插值FIR滤波器设计

描述

实例

B= intfilt (L,P,阿尔法)设计了一个线性相位FIR滤波器,该滤波器使用最近邻滤波器进行理想的带限插值2.*P非零采样,用于与L-每隔1个连续零L样本,假设原始带宽限制为阿尔法乘以奈奎斯特频率。返回的滤波器B内页

B= intfilt (L,N,“拉格朗日”)设计一个FIR滤波器执行N三阶拉格朗日多项式的插值L-每隔1个连续零L样品。

例子

全部崩溃

设计一个数字内插滤波器,使用带限方法将信号采样增加7倍。指定“带限”因子为0.5,并使用 2. × 2. 插值中的采样。

upfac=7;α=0.5;h1=intfilt(upfac,2,alpha);

当原始信号的带宽限制为时,滤波器工作最佳阿尔法乘以奈奎斯特频率。通过生成200个高斯随机数和用40阶FIR低通滤波器过滤序列来创建一个限带噪声信号。重置随机数生成器以获得可重复的结果。

lowp=fir1(40,α);rng(“默认”)x=滤波器(低电平,1,随机数(200,1));

通过在信号的每对采样之间插入零来增加信号的采样率x

xr=上采样(x,upfac);

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

y =过滤器(h1 1 xr);

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

延迟=平均值(GRP延迟(h1));y(1:延迟)=[];阀杆(1:upfac:upfac*长度(x),x)保持地块(y)xlim([400 700])

图中包含一个轴对象。轴对象包含两个类型为stem、line的对象。

集成滤波器还执行拉格朗日多项式插值。

  • 一阶多项式内插是线性内插,用三角滤波器实现。

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

对原始信号进行插值并叠加结果。

h2=intfilt(upfac,1,“拉格朗日”); y2=滤波器(h2,1,xr);y2(1:楼层(平均值(grpdelay(h2)))=[];绘图(y2)暂停

图中包含一个轴对象。轴对象包含3个类型为stem、line的对象。

输入参数

全部崩溃

样本数,指定为正整数标量。集成滤波器设计了一个线性相位FIR滤波器,该滤波器使用了一个插入L-每隔1个连续零L样品。

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

转换带宽的反向度量,指定为标量。阿尔法与滤波器的过渡带宽成反比,并影响阻带中无关区域的带宽。指定阿尔法允许您指定输入信号占用的奈奎斯特间隔的多少。这对要插值的信号是有益的,因为它允许您在不影响插值的情况下增加过渡带宽,并在给定的时间内产生更好的阻带衰减LP.如果你设置阿尔法如果设置为1,则假定您的信号占用整个奈奎斯特间隔。设置阿尔法小于1允许停止带中的“不关心”区域。例如,如果您的输入占用了一半的Nyquist间隔,您可以将alpha设置为0.5。

指定为正整数标量的拉格朗日多项式的阶数。FIR滤波器执行N三阶拉格朗日多项式的插值L-每隔1个连续零L样品。如果两者都有NL均为偶数时,设计的滤波器相位不是线性的。

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

输出参数

全部崩溃

过滤器系数,作为向量返回。元素B为FIR滤波器的系数。如果阿尔法时,假设原始带宽限制为阿尔法乘以奈奎斯特频率。B长度是2 *L*P-1.

对于N次拉格朗日多项式插值,B长度(N+1)*L对于N均匀,长度(N+1)*l-1对于N古怪的

算法

带宽限制方法使用冷杉设计一个插值FIR滤波器。多项式法使用拉格朗日多项式插值公式对等间距样本构造适当的滤波器。这两种滤波器基本上都是低通的,增益为L在通带中。

扩展能力

之前介绍过的R2006a