prony

Prony 法によるフィルター設計

ページ内をすべて折りたたむ

構文

[b,a] = prony(h,bord,aord)

説明

例

[b,a] = prony(h,bord,aord) からは、インパルス応答が h、分子次数が bord、分母次数が aord である因果性有理伝達関数の分子係数と分母係数が返されます。

例

すべて折りたたむ

Prony 法を使用したフィルター応答

ライブスクリプトを開く

ローパスフィルターのインパルス応答に、4 次の IIR モデルを近似させます。元のインパルス応答と Prony 法で設計したインパルス応答をプロットします。

d = designfilt('lowpassiir','NumeratorOrder',4,'DenominatorOrder',4, ...
    'HalfPowerFrequency',0.2,'DesignMethod','butter');

h = filter(d,[1 zeros(1,31)]);
bord = 4;
aord = 4;
[b,a] = prony(h,bord,aord);

subplot(2,1,1) 
stem(impz(b,a,length(h)))
title 'Impulse Response with Prony Design'

subplot(2,1,2)
stem(h)
title 'Input Impulse Response'

Figure contains 2 axes objects. Axes object 1 with title Impulse Response with Prony Design contains an object of type stem. Axes object 2 with title Input Impulse Response contains an object of type stem.

ハイパスフィルターのインパルス応答に、10 次の FIR モデルを近似させます。元の周波数応答と Prony 法で設計した周波数応答をプロットします。応答は高精度で一致します。

d = designfilt('highpassfir','FilterOrder',10,'CutoffFrequency',0.8);

h = filter(d,[1 zeros(1,31)]);
bord = 10;
aord = 0;
[b,a] = prony(h,bord,aord);

fvt = fvtool(b,a,d);
legend(fvt,'Prony','Original')

{"String":"Figure Figure 1: Magnitude Response (dB) contains an axes object. The axes object with title Magnitude Response (dB) contains 2 objects of type line. These objects represent Prony, Original.","Tex":"Magnitude Response (dB)","LaTex":[]}