%exm095_1.m
clear
clf
%圖形數據的產生
zeta2=[0.2 0.4 0.6 0.8 1.0];
n=length(zeta2);
for k=1:n;
    Num{k,1}=1;
    Den{k,1}=[1 2*zeta2(k) 1];
end
S=tf(Num,Den);					%產生單輸入多輸出系統
t=(0:0.4:20)';					%時間采樣點
[Y,x]=step(S,t);				%單輸入多輸出系統的響應
tt=t*ones(size(zeta2));		%為畫彩帶圖，生成與函數值Y維數相同的時間矩陣。
%產生全窗軸位框，并隱去坐標軸
clf reset,H=axes('Position', [0, 0, 1, 1], 'Visible', 'off');
%產生包含多行字符串的元胞數組
str{1}='\fontname{隸書}二階系統階躍響應';									%<11>
str{2}='y(t) = 1 - \beta^{-1}e^{-\zetat}sin(\betat + \theta)';
str{3}=' ';str{4}='\fontname{隸書}其中：';
str{5}='\beta = (1 - \zeta^{2})^{0.5}';
str{6}='\theta = arctg(\beta/\zeta)';
str{7}='\zeta = .2, .4, .6, .8, 1';										%<15>
%使H句柄軸對象成為當前軸，然后注釋多行文字。
set(gcf, 'CurrentAxes',H)													%<18>
text(0.01, 0.73, str, 'FontSize', 12)									%<19>
h1=axes('Position', [0.45, 0.45, 0.5, 0.5]);	%產生右半窗的軸位框
ribbon(tt,Y,0.4)										%在h1軸位框中畫彩帶圖
%對X軸、Z軸重標刻度值
set(h1,'XTickLabelMode','manual','XTickLabel','0|0.4|0.8|1.2');%<23>
set(h1,'ZTickLabel','0|1.0|2.0');										%<24>
%低層指令標識軸名
set(get(h1,'XLabel'),'String','\zeta \rightarrow','Rotation',17.5)
set(get(h1,'YLabel'),'String','\leftarrow t','Rotation',-25)	%<27>
set(get(h1,'Zlabel'),'String','y \rightarrow')
h2=axes('Position',[0.03, 0.08, 0.27, 0.27]);%在左下角，產生小的軸位框。
plot(tt,Y)												%在h2軸對象上畫二維圖
%在右下方畫系統方塊框圖
h3=axes('Position',[0.37,0.04,0.63,0.32]);	%設置畫框圖的軸位框
set(h3,'Xlim',[0,1.2],'Ylim',[0,0.5])			%設置軸的刻度范圍
set(h3,'DataAspectRatio',[1 1 1])				%設置刻度比例
set(h3,'ColorOrder',[0,0,0])						%設置繪線的首選用色
set(h3,'Visible','off')								%隱去坐標軸
hh1=rectangle('Position',[0.5,0.2,0.4,0.2],'Curvature',[0,0]);
																%畫方框				<37>
hh2=rectangle('Position',[0.2,0.26,0.08,0.08],'Curvature',[1,1]);
																%畫圓框				<38>
xx1=0.05:0.01:0.2;xx2=0.28:0.02:0.5;
xx3=0.9:0.02:1.1;xx4=0.24:0.02:1;
yy5=0.1:0.02:0.26;yy6=0.1:0.02:0.3;
yy1=0.3*ones(size(xx1));yy2=0.3*ones(size(xx2));
yy3=0.3*ones(size(xx3));yy4=0.1*ones(size(xx4));
xx5=0.24*ones(size(yy5));xx6=ones(size(yy6));
line(xx1,yy1);line(xx2,yy2);line(xx3,yy3);line(xx4,yy4);
line(xx5,yy5);line(xx6,yy6)
line(0.17,0.3,'Marker','>','MarkerFaceColor','k')
line(0.47,0.3,'Marker','>','MarkerFaceColor','k')
line(1.1,0.3,'Marker','>','MarkerFaceColor','k')
line(0.24,0.23,'Marker','^','MarkerFaceColor','k')
line(0.17,0.35,'Marker','+')
text(0.27,0.23,'-')
text(0.05,0.35,'u(t)')
text(1,0.35,'y(t)')
text(0.6,0.26,'s{^2} + 2{\zeta}s');
xx7=0.56:0.02:0.84;yy7=0.3*ones(size(xx7));line(xx7,yy7)
text(0.68,0.35,'1')