%exm07417_1.m
clear,clf,shg
zeta2=[0.1 0.2 0.3 0.4 0.5 0.6 0.8 1.0];						%<2>
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:30)';				%時間采樣點
[Y,x]=step(S,t);				%單輸入多輸出系統的響應
tt=t*ones(size(zeta2));		%為畫彩帶圖，生成與函數值Y維數相同的時間矩陣。
ribbon(tt,Y,0.4)				%畫彩帶圖
%至此彩帶圖已經生成。以下指令都是為了使圖形效果更好、標識更清楚而用。
view([150,50]),shading interp,colormap(jet)%設置視角、明暗、色圖
light,lighting phong,box on					%設置光源、照射模式、坐標框
for k=1:n;str_lgd{k,1}=num2str(zeta2(k));end,legend(str_lgd)%圖例設置
str1='\itG = (s^{2} + 2\zetas + 1)^{-1}';
str2='\fontsize{16}\fontname{隸書}取不同';
str3='{\fontsize{10}\it\zeta}';
str4='\fontsize{16}\fontname{隸書}時的階躍響應';
title([str1,str2,str3,str4]),zlabel('\ity(\zeta,t) \rightarrow') 
disp('按任意鍵，將連續改變觀察視角,產生動的效果')
pause
for ii=1:24;view([150+15*ii,50]),drawnow,end