freq

周波数応答

呼び出し手順

[x]=freq(A,B,C [,D],f)
[x]=freq(NUM,DEN,f)

パラメータ

A, B, C, D: 適当な次元nxn, nxp, mxn, mxpの実数行列 .
NUM,DEN: 次元 mxpの多項式行列
x: 実数または複素数の行列

説明

x=freq(A,B,C [,D],f) は次の関係を満たす実数または複素数の mxp*t 行列を返します.:

x(:,k*p:(k+1)*p)= C*inv(f(k)*eye()-A)*B + D.

ここで, f が虚軸上または単位円上に値を有する場合, xは連続系または離散系における (A,B,C,D)の周波数応答となります.

x=freq(NUM,DEN,f) は, x の k*(p-1)+1列から k*p列が行列NUM(f(k))./DEN(f(k))を含むような実数または複素数の行列xを返します.

例

s=poly(0,'s');
sys=(s+1)/(s^3-5*s+4)
rep=freq(sys("num"),sys("den"),[0,0.9,1.1,2,3,10,20])
[horner(sys,0),horner(sys,20)]
//
Sys=tf2ss(sys);
[A,B,C,D]=abcd(Sys);
freq(A,B,C,[0,0.9,1.1,2,3,10,20])

説明

Report an issue
<< dbphi	Frequency Domain	freson >>