例5.2の図5.6のプロット

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

#例5.2

#パラメータの値を与える
R = 100 #抵抗の値
L = 10 #インダクタンスの値

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

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

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

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

#図5.6のプロット
plt.plot(t, y) #ステップ応答をプロット
plt.xlim([0,1]) #横軸（時間軸の範囲の指定）
plt.yticks([0,0.002,0.004,0.006,0.008,0.010,0.012]) #縦軸の目盛りの値の設定
plt.grid(color='gray') #罫線を灰色で表示
plt.xlabel("time   t[s]") #横軸のラベル表示
plt.ylabel("i(t) [A]") #縦軸のラベル表示
plt.show() #グラフの表示