Scilab Home page | Wiki | Bug tracker | Forge | Mailing list archives | ATOMS | File exchange
Please login or create an account
Change language to: English - Français - Português -

Please note that the recommended version of Scilab is 6.0.1. This page might be outdated.
See the recommended documentation of this function

Scilab help >> Signal Processing > yulewalk

yulewalk

最小二乗フィルタを設計

呼び出し手順

Hz = yulewalk(N,frq,mag)

パラメータ

N

整数 (指定するフィルタの次数)

frq

実数行ベクトル (非降順), 周波数.

mag

非負の実数行ベクトル (frqと同じ大きさ), 指定する振幅.

Hz

フィルタ B(z)/A(z)

説明

Hz = yulewalk(N,frq,mag) は N次IIRフィルタを見つけます.

n-1         n-2            
B(z)   b(1)z     + b(2)z    + .... + b(n)
H(z)= ---- = ---------------------------------
n-1       n-2
A(z)    z   + a(2)z    + .... + a(n)

は,ベクトルfrqおよびmagで指定された周波数応答に振幅が 一致します. 周波数frqは 0.0および1.0の間とする必要があり, 1.0はサンプルレートの半分に対応します. これらは昇順で, 0.0 から始まり, 1.0で終わる筆王が あります.

f=[0,0.4,0.4,0.6,0.6,1];H=[0,0,1,1,0,0];Hz=yulewalk(8,f,H);
fs=1000;fhz = f*fs/2;  
clf(0);xset('window',0);plot2d(fhz',H');
xtitle('Desired Frequency Response (Magnitude)')
[frq,repf]=repfreq(Hz,0:0.001:0.5);
clf(1);xset('window',1);plot2d(fs*frq',abs(repf'));
xtitle('Obtained Frequency Response (Magnitude)')
Scilab Enterprises
Copyright (c) 2011-2017 (Scilab Enterprises)
Copyright (c) 1989-2012 (INRIA)
Copyright (c) 1989-2007 (ENPC)
with contributors
Last updated:
Thu Mar 03 11:00:56 CET 2011