例
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(:))]) |  |  |
周波数応答フィッティング
sys=mrfit(w,mag,order) [num,den]=mrfit(w,mag,order) sys=mrfit(w,mag,order,weight) [num,den]=mrfit(w,mag,order,weight)
周波数ベクトルで正の実数(Hz)
周波数応答の振幅の実数ベクトル (wと同じ大きさ)
整数 (必要な次数, denの次数)
正の実数ベクトル (デフォルト値: ones(w)).
安定な多項式
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(:))]) | | |