Re: Disturbance Responce for a simple control system

[A. Scottedward Hodel]

# num_p, denom_p: numerator, denominator coefficients of plant
# num_k, denom_k: similar for controller

G = tf2sys(num_p, denom_p,0,"plant input","plant output");
K = tf2sys(num_k,denom_k,0,"contr input","contr output");
gk = sysgroup(G,K);    # parallel arrangement of G, K, no interconnections yet
sysout(gk)        # just to verify that input/output indices are what I think 
they are
gk = sysconnect(gk,[2 1],[1 2]);
# above is eqivalent to:
#    gk = sysconnect(gk,2,1);   # control output to plant input
#    gk = sysconnect(gk,1,2);   # plant output to control input
#    The above 2-step procedure avoids the warning about an algebraic loop
step(gk)

# or, if you just want the step response from the disturbance to
# the output:
T_out_dist = sysprune(gk,1,1);    # keep plant output, plant input
step(T_out_dist);

#  I should probably modify sysconnect (and sysprune) so that if
#    desired users can specify input/output names instead of indices

A S Hodel Assoc. Prof. Dept Elect Eng, Auburn Univ,AL  36849-5201
On leave at NASA Marshall Space Flight Center (256) 544-1426

----------
>From: "Hiroshi Shinozawa" <address@hidden>
>To: <address@hidden>
>Subject: Disturbance Responce for a simple control system
>Date: Fri, 25 Sep 1998 13:47:54 +0900
>

>I think about the following Feedback Loop
>
>DISTURBANCE -----------------------
>                                   | 
>STEP INPUT --->O---->|Controller|--O--|Process|--O---> OUTPUT
>               |                                 |
>               -----------------------------------
>I know Step Input Responce is easily calculated by OCST 
>"step" function.
>But I don't know how to calculate the Step Disturbance
>Responce.
>
>I use "lsode" function to simulate a step disturbance
>responce as well as a step input reponce.
>I am searching a simple function like OCST "step" to 
>simulate step disturbance respponce.
>Advise me please. 
>**********************************************
>address@hidden
>Hiroshi_Shinozawa
>YHC TC Phone 731-22535
>**********************************************
>