H は SISO であるため、S はスカラーです。H が MIMO であった場合、S の次元は H の I/O 次元と一致します。
S
S =
1000000
元のシステムを、元のシステムと因数分解したシステム間の差分と比較し、G と S が H = G'*S*G を満たすことを確認します。差分が非常に小さいため、sigmaplot は警告をスローします。
Hf = G'*S*G;
sigmaplot(H,H-Hf)
Warning: The frequency response has poor relative accuracy. This may be because the response is nearly zero or infinite at all frequencies, or because the state-space realization is ill conditioned. Use the "prescale" command to investigate further.
[G,S] = spectralfact(H) 構文を使用する場合、H = G'*S*G。S の次元は H および G の I/O の次元と一致します。
[G,S] = spectralfact(F,R) 構文を使用する場合、G'*S*G = F'*R*F。各次元に沿った S のサイズは F の出力の数と一致します。
ヒント
spectralfact は H が自己共役であると仮定します。H が自己共役でない場合、spectralfact は H = G'*S*G を満たさない G と S を返します。したがって、spectralfact を使用する前に入力モデルが実際に自己共役であることを確認します。H を確認する 1 つの方法は、特異値プロットで H を H - H' と比較することです。
sigmaplot(H,H-H')
H が自己共役の場合、プロット上の H - H' ラインが H ラインよりもはるかに下に配置されます。
You can also select a web site from the following list:
How to Get Best Site Performance
Select the China site (in Chinese or English) for best site performance. Other MathWorks country sites are not optimized for visits from your location.