frfit
周波数応答フィッティング
呼び出し手順
sys=frfit(w,fresp,order) [num,den]=frfit(w,fresp,order) sys=frfit(w,fresp,order,weight) [num,den]=frfit(w,fresp,order,weight)
引数
- 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])
参照
| Report an issue | ||
| << Identification | Identification | lattn >> |