exm07417_1.m

精通matlab

%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