reinf3_13.m

来自「离散控制系统设计的MATLAB 代码」· M 代码 · 共 27 行

%%%%%%%%%% Reinforcement Problem 3.13 %%%%%%%%%%%
%   Discrete-Time Control Problems using        %
%       MATLAB and the Control System Toolbox   %
%   by J.H. Chow, D.K. Frederick, & N.W. Chbat  %
%         Brooks/Cole Publishing Company        %
%                September 2002                 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   ---- Another feedback connection ----
%
clear
disp('Reinforcement Problem 3.13')
Ts = 1;
denG = conv([1 0.5],[69 -5 5]);
G = tf(6*[3 1],denG,Ts)         % forward path in TF form
H = tf([15 1],[5 -1],Ts)        % feedback path in TF form
T = feedback(G,H)                   % closed-loop system
[zT,pT,kT] = zpkdata(T,'v')         % CL zeros,poles,gain
[magpT,thetapT] = xy2p(pT)          % polar form for CL poles

dtime = 0:Ts:40;
y = step(T,dtime);                      % compute step response
figure
stem(dtime,y,':','filled');grid         % plot step response
title('Step response for Reinforcement Problem 3.13')
xlabel('Discrete time k')
%%%%%%%%%%