線形時不変なシステムのH∞ノルムを計算する

MATLAB/simulink

線形時不変なシステムのH∞ノルムを計算する。適当な伝達関数

\begin{align}
G(s)=\frac{1}{(s+1)(s+2)}
\end{align}

を定義する。H∞ノルムの定義は

\begin{align}
\| G(j \omega) \|_{\infty} = \max_{\omega \in(-\infty,\infty)}|G(j \omega)|
\end{align}

である。一応、MATLABの関数を使えば

hinfnorm(G)

で同様の結果を計算できる。

コード

s=tf('s');
j=sqrt(-1);
N=1;

omega=linspace(0,10^(N),10000);
G=1/((s+1)*(s+2));
G0=zeros(size(omega));

for i=1:1:length(omega)
    G0(1,i)=norm(1./((j*omega(1,i)+1).*(j*omega(1,i)+2)));
end

figure
plot(omega,G0)
hinfnorm(G)

コメント

タイトルとURLをコピーしました