frfit

引数

w

周波数 (Hz)のベクトル(正の実数)

fresp

周波数応答の複素ベクトル (wと同じ大きさ)

order

整数 (指定する次数, denの次数)

weight

正の実数ベクトル (デフォルト値: ones(w)).

num,den

安定な多項式

説明

sys=frfit(w,fresp,order,weight) は, 周波数応答G(w(i))がfresp(i)に一致する, すなわち,freq(num,den,%i*w)がfrespに 近くなるような,指定した次数orderの 双安定の伝達関数 G(s)=sys=num/denを返します. weight(i)はw(i)に関する重みです.

例

w=0.01:0.01:2;s=poly(0,'s');
G=syslin('c',2*(s^2+0.1*s+2), (s^2+s+1)*(s^2+0.3*s+1));
fresp=repfreq(G,w);
Gid=frfit(w,fresp,4);
frespfit=repfreq(Gid,w);
bode(w,[fresp;frespfit])

参照

• frep2tf — 周波数応答から伝達関数実現を得る
• factors — 実数因数分解
• cepstrum — ケプストラム計算
• mrfit — 周波数応答フィッティング
• freq — 周波数応答
• calfrq — 離散的な周波数応答