reinf3_14.m

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

%%%%%%%%%% Reinforcement Problem 3.14 %%%%%%%%%%%
%   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                 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   ---- One more feedback connection ----
%
clear
disp('Reinforcement Problem 3.14')
Ts = 1;
denG = conv([10 -0.5],[3 0.25]);
G = tf([12 -5],denG,Ts)          % forward path in TF form
H = tf(6,[3 -2 -0.3],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.14')
xlabel('Discrete time k')
%%%%%%%%%%