Scilab Website | Contribute with GitLab | Mailing list archives | ATOMS toolboxes
Scilab Online Help
2024.1.0 - 日本語


mrfit

周波数応答フィッティング

呼び出し手順

sys=mrfit(w,mag,order)

[num,den]=mrfit(w,mag,order)

sys=mrfit(w,mag,order,weight)

[num,den]=mrfit(w,mag,order,weight)

引数

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 — 離散的な周波数応答
Report an issue
<< lattp Identification phc >>

Copyright (c) 2022-2024 (Dassault Systèmes)
Copyright (c) 2017-2022 (ESI Group)
Copyright (c) 2011-2017 (Scilab Enterprises)
Copyright (c) 1989-2012 (INRIA)
Copyright (c) 1989-2007 (ENPC)
with contributors
Last updated:
Mon Jun 17 17:54:18 CEST 2024