%%% インパルス応答 %%%
clear;

t=0:0.001:5;

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

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

%状態空間表現
A=[0 1
   -k/m -c/m];
B=[0
   1/m];
C=[1 0];
D=0;

%インパルス応答
yi=impulse(A,B,C,D,1,t);

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

