%%% ステップ応答 %%%
clear;

t=0:0.001:5;

%プラントの定義
w=4.7;		%固有角振動数
tuta1=0.3;	%減衰比
tuta2=0.5;	%減衰比
tuta3=0.7;	%減衰比

m=0.8;		%質量
k=m*w^2;	%ばね定数
c=2*tuta2*w;	%減衰係数

%伝達関数表現
num=[0 0 1];
den=[m c k];

%ステップ応答
ys=step(num,den,t);

%グラフのプロット
plot(t,ys)				%グラフを描かせる
grid					%グラフの背景に等間隔に基準線を引く
xlabel('Time(sec)'),ylabel('x')		%x軸、y軸のラベル
title('Step response')			%グラフのタイトル

