% exp4_15.m
clear
close all
clc
% 狀態空間系統描述
a=[-0.6 -1.04 0 0;1.04 0 0 0;0 0.96 -0.7 -0.32;
   0 0 0.32 0];
b=[1 0 0 0]';
c=[0 0 0 0.32];
d=0;
% 圖1繪制波特圖
figure(1)
bode(a,b,c,d);
% 圖2繪制幅相曲線
figure(2)
subplot(211)
nyquist(a,b,c,d);
[z,p,k]=ss2zp(a,b,c,d);
p
subplot(212)
[rm,im]=nyquist(a,b,c,d);
plot(rm,im)
% 圖3繪制帶裕度及相應頻率顯示的波特圖
figure(3)
margin(a,b,c,d);
% 圖4繪制沖激響應曲線
figure(4)
[ac,bc,cc,dc]=cloop(a,b,c,d);
impulse(ac,bc,cc,dc)