%図7.2
%
clear all; close all;


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

%伝達関数表現を与える
sys = tf( num, den ); %伝達関数表現

%時間変数の定義
t = 0:0.01:50; %0から50まで0.01刻み

%ステップ応答の計算
[ y, t ] = step( sys, t ); %ステップ応答

%図7.1のプロット
plot(t,y); %ステップ応答をプロット
xlim([0,40]) %横軸(時間軸)の範囲の指定
ylim([-50 50]) %縦軸の範囲の指定
xticks([0 5 10 15 20 25 30 35 40]) %横軸の目盛りの値の設定
yticks([-50 -40 -30 -20 -10 0 10 20 30 40 50]) %縦軸の目盛りの値の設定
grid; %罫線を表示
xlabel('time  t [s]'); %横軸のラベル表示
ylabel('y(t)') %縦軸のラベル表示