%図7.1
%
clear all; close all;


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

%伝達関数表現を与える
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,50]) %横軸(時間軸)の範囲の指定
ylim([0.0 3.5]) %縦軸の範囲の指定
xticks([0 10 20 30 40 50]) %横軸の目盛りの値の設定
yticks([0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5]) %縦軸の目盛りの値の設定
grid; %罫線を表示
xlabel('time  t [s]'); %横軸のラベル表示
ylabel('y(t)') %縦軸のラベル表示