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