図7.2のプロット

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

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

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

#時間変数の定義
t = arange(0, 50, 0.01) #0から50まで0.01刻み

#ステップ応答の計算
y, t = matlab.step(sys,t) #ステップ応答

#図7．2のプロット
plt.plot(t, y) #ステップ応答をプロット
plt.xlim([0,40]) #横軸（時間軸）の範囲の指定
plt.ylim([-50,40]) #縦軸（出力)の範囲の指定）
plt.xticks([0,5,10,15,20,25,30,35,40]) #横軸の目盛りの値の設定
plt.yticks([-50,-40,-30,-20,-10,0,10,20,30,40,50]) #縦軸の目盛りの値の設定
plt.grid(color='gray') #罫線を灰色で表示
plt.xlabel("time   t[s]") #横軸のラベル表示
plt.ylabel("y(t)") #縦軸のラベル表示
plt.show() #グラフの表示