mrfit

引数

w

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

mag

周波数応答の振幅の実数ベクトル (wと同じ大きさ)

order

整数 (必要な次数, denの次数)

weight

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

num,den

安定な多項式

説明

sys=mrfit(w,mag,order,weight) は, 周波数応答の振幅abs(G(w(i)))が mag(i)に一致する,すなわち, abs(freq(num,den,%i*w))がmagに近い, 指定された次数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)); // syslin('c',Num,Den);
fresp=repfreq(G,w);
mag=abs(fresp);
Gid=mrfit(w,mag,4);
frespfit=repfreq(Gid,w);
plot2d([w',w'],[mag(:),abs(frespfit(:))])

参照

• cepstrum — ケプストラム計算
• frfit — 周波数応答フィッティング
• freq — 周波数応答
• calfrq — 離散的な周波数応答