Loading web-font TeX/Math/Italic

MATLABでナイキスト線図を書く

MATLAB/simulink

ナイキスト線図を使えば複雑なシステムの安定判別をより容易に行えるようになる場合がある。今回は細かいことを気にせずただナイキスト線図を描画する。正しく描けているかのチェックにはMATLABの関数「nyquist」を使用した。

いま開ループ伝達関数が

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

で与えられているとする。ナイキスト線図を書くにはs=j \omegaを代入し、\omegaを変化させればいい。つまり

\begin{align} G(j \omega) =\frac{2(j \omega)^2+5(j \omega)+1}{(j \omega)^2+2(j \omega)+3} \end{align}

これを計算してプロットした結果がこちら

「nyquist」で描画したときと同様の結果が得られた。

MATLABコードはこちら

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
s=tf('s');
 
G =(2*s^2+5*s+1)/(s^2+2*s+3);
nyquist(G)
 
j=sqrt(-1);
 
omega=-100:0.001:100;
G =(2*(j*omega).^2+5*(j*omega)+1)./((j*omega).^2+2.*(j*omega)+3);
 
figure;
plot(G)
xlabel('Re');
ylabel('Im');
grid on

コメント

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