reinf3_12.m

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

%%%%%%%%%% Reinforcement Problem 3.12 %%%%%%%%%%%
%   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                 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   ---- Feedback connection ----
%
clear
disp('Reinforcement Problem 3.12')
Ts = 1;
G = tf(0.5*[3 -1],[13 -5 8],Ts)        % forward path in TF form
H = tf([7 4],conv([1 0.5],[4 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:40;
y = step(T,dtime);                     % compute step response
figure
stem(dtime,y,':','filled');grid        % plot step response
title('Step response for Reinforcement Problem 3.12')
xlabel('Discrete time k')
%%%%%%%%%%