図12.18のプロット

#計算とグラフプロットに必要なモジュールの読み込み
import numpy as np
from control import matlab
from matplotlib import pyplot as plt
from scipy import arange 
from scipy import signal

#伝達関数の分子・分母多項式を与える
num = [0, 1] #分子多項式
den = [1, 1] #分母多項式

#伝達関数表現を与える
sys = matlab.tf(num,den) #伝達関数表現

w = np.logspace(-2, 2, 1000) #対数的に等間隔なベクトルの生成(10^{-2}から10^{2}で1000点)

matlab.nyquist(sys,w) #ナイキスト線図の計算とプロット
plt.xlim([-0.1,1.0]) #横軸（実軸）の範囲の指定
plt.ylim([-0.8,0.8]) #縦軸（虚軸）の範囲の指定
plt.grid(color='gray') #罫線を灰色で表示
plt.xlabel('Re') #横軸のラベル表示
plt.ylabel('Im') #縦軸のラベル表示
plt.show() #グラフの表示