主要内容

普龙尼

滤波器设计的Prony方法

描述

例子

b一个= prony(h矿房aord返回具有脉冲响应的因果有理传递函数的分子和分母系数h,分子阶矿房,分母顺序aord

例子

全部折叠

用四阶IIR模型拟合低通滤波器的脉冲响应。绘制原始和prony设计的脉冲响应。

D = designfilt(“lowpassiir”“NumeratorOrder”,4,“DenominatorOrder”,4,...“HalfPowerFrequency”, 0.2,“DesignMethod”“黄油”);H = filter(d,[1 0 (1,31)]);Bord = 4;Aord = 4;[b,a] = prony(h, board,aord);Subplot (2,1,1) stem(impz(b,a,length(h)))标题“Prony设计的冲动反应”副图(2,1,2)干(h)标题“输入脉冲响应”

图中包含2个轴对象。标题为Impulse Response with Prony Design的Axes对象1包含一个stem类型的对象。标题为Input Impulse Response的Axes对象2包含一个stem类型的对象。

将10阶FIR模型拟合到高通滤波器的脉冲响应中。绘制原始和prony设计的频率响应。响应匹配精度高。

D = designfilt(“highpassfir”“FilterOrder”10“CutoffFrequency”, 0.8);H = filter(d,[1 0 (1,31)]);Bord = 10;Aord = 0;[b,a] = prony(h, board,aord);FVT = fvtool(b,a,d);传奇(fvt“普龙尼”“原始”

图形过滤可视化工具-幅度响应(dB)包含一个轴对象和其他类型的uitoolbar, uimenu对象。标题为Magnitude Response (dB)的axis对象包含2个类型为line的对象。这些物件分别代表Prony、Original。

输入参数

全部折叠

脉冲响应,用矢量表示。

例子:0.5) impz (fir1(20日)给出了20阶FIR滤波器在归一化截止频率下的脉冲响应π/ 2rad /样品。

数据类型:|
复数支持:金宝app是的

分子和分母的顺序,指定为正整数标量。如果的长度h小于max (矿房aord时,函数将脉冲响应填充为零。

  • 如果你想要一个全极传递函数,请指定矿房作为0

  • 如果你想要一个全零传递函数,请指定aord作为0

数据类型:|

输出参数

全部折叠

传递函数系数,作为向量返回。b长度矿房+ 1而且一个长度aord+ 1

更多关于

全部折叠

传递函数

传递函数是脉冲响应的z变换hn

H z n h n z n

一个有理传递函数多项式的比值在吗z1.这个方程描述了一个分子阶的因果有理传递函数分母的顺序p

H z B z 一个 z k 0 b k z k 1 + l 1 p 一个 l z l

在哪里一个[0] = 1

参考文献

[1]托马斯·W·帕克斯和c·西德尼·伯勒斯。数字滤波器设计.美国纽约:Wiley-Interscience, 1987。

扩展功能

C/ c++代码生成
使用MATLAB®Coder™生成C和c++代码。

R2006a之前介绍