imppro.m

华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 1. 安装MATLAB软件。安装时

function [mp,tm,ts]=imppro(tt,yy,tss)
% 计算脉冲响应曲线的性能指标

% Copyright 2001-2002 ECUST. 
%   $Revision: 1.2 $  $Date: 2001/05/12 07:15:33 $

if nargin < 3,tss=0.05;end;
if nargin==2 & length(tt)~=length(yy)
    error('The input and output must have same length');
end;
if nargin < 2, error('Not enough data');end;
tend=tt(length(tt));[mp,ind]=max(yy);tm=tt(ind);
i=length(tt);n=0;
while n==0,i=i-1;
    if i==1| yy(i)>=tss*mp,n=1;end;
end;t1=i;
i=length(tt);n=0;
while n==0,i=i-1;
    if i==1 | yy(i)<=-tss*mp,n=1;end;
end;t2=i;
if t2>t1,ts=tt(t2);j=t2;else ts=tt(t1);j=t1;end;
if nargout== 0,
plot(tt,yy,tm,yy(ind),'*k',ts,yy(j),'*k');line([0 tend],[0 0]);
line([0 tend],[tss*mp tss*mp],'lines','--');line([tm tm],[mp 0],'lines','--');
line([0 tend],[-tss*mp -tss*mp],'lines','--');
text(tm,mp,['\fontsize {10} \bf','  MP=',num2str(mp),';  TM=',num2str(tm)]);
text(ts,yy(j),['\fontsize {10}\bf',' TS=',num2str(ts)]);
end;
% imppro
% num=[0.0225];den=[1 0.24 0.0225];
% G=tf(num,den);
% tt=0:.01:50;yy=impulse(G,tt);