周波数応答フィッティング
sys=frfit(w,fresp,order) [num,den]=frfit(w,fresp,order) sys=frfit(w,fresp,order,weight) [num,den]=frfit(w,fresp,order,weight)
周波数 (Hz)のベクトル(正の実数)
周波数応答の複素ベクトル (w
と同じ大きさ)
整数 (指定する次数, den
の次数)
正の実数ベクトル (デフォルト値: ones(w)
).
安定な多項式
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]) | ![]() | ![]() |