0% found this document useful (0 votes)
28 views3 pages

Advanced Control Systems Assignment 4

The document outlines an assignment for an advanced control systems course, detailing several problems related to state-space representation, eigenvalues, time-domain solutions, transfer functions, and system matrices. Each problem requires the application of control theory concepts to analyze and compute various aspects of dynamical systems. The assignment is due on February 17, 2025.

Uploaded by

diliputukuri76
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
28 views3 pages

Advanced Control Systems Assignment 4

The document outlines an assignment for an advanced control systems course, detailing several problems related to state-space representation, eigenvalues, time-domain solutions, transfer functions, and system matrices. Each problem requires the application of control theory concepts to analyze and compute various aspects of dynamical systems. The assignment is due on February 17, 2025.

Uploaded by

diliputukuri76
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

EE3302 - Advanced Control Systems

Name:.................................. Due Date of Submission:17.02.2025


Class Roll No.:.................... Assignment: Number 04

Problem 1 The state-space representation of a dynamical system is given as follows:

ẋ(t) = Ax(t) + Bu(t) (1)


y(t) = Cx(t)
     
0 1 1   −2
with A = ,B= , C = 2 1 , x0 = , D = 0.
0 −2 −1 3

(a) Find the eigenvalues and eigenvectors of the matrix A.

(b) Compute eAt using Laplace transform approach and C-H theorem.

(c) Assume that the control input is u(t) = 0, compute x(t) and y(t).

(d) Assume that the input is u(t) = 1 + 2e−2t , compute x(t) and y(t).
 
x1 (t)
(e) Given your answer to the previous question, compute x(t) = as t → ∞.
x2 (t)
Which state blows up? Also, find y(∞).

Problem 2

(a) For the following system, explicitly determine the time-domain solution x(t)
   
4 1 1
ẋ(t) = x(t) + u(t)
0 10 0
 
y(t) = 1 0 x(t)
 
1
where x(0) = and u(t) is a unit step function.
1

(b) Determine the transfer function G(s) for the zero initial state response, given the
system in part (a).

(c) Does the degree of state space model and transfer function match? Why or why
not?

1
Problem 3 For a system represented by the state equation ẋ(t) = Ax(t). The response
 −2t   
e 1
x(t) = when x(0) =
−2e−2t −2
and
e−t
   
1
x(t) = when x(0) =
−e−t −1

Determine the system matrix A and the state transition matrix(STM).

Problem 4 The block diagram of a linear time invariant system is given in fig.
(a) Write down the state variable equations for the system in matrix form assuming
 >
state vector to be x1 (t) x2 (t) .

(b) Find out the state transition matrix.

(c) Determine y(t), t ≥ 0, when the initial values of the state at time t = 0 are x1 (0) = 1
and x2 (0) = 1 with unit step input.

Problem 5 Given the following LTI dynamical system

ẋ(t) = Ax(t) + Bu(t); x(0) = xt0 (2)


y(t) = Cx(t) + Du(t)
   
1 −1 0 1 0  
0 1 0  T
where A =  1 −1 0 , B =  0 0 , C = , x(t0 ) = x(1) = 0 1 1
0 0 1
0 0 0 0 1
(a) Is A nilpotent of order 2?

(b) Determine eAt , eA(t−τ ) , eA(t−t0 ) , t0 = 1.

(c) If u(t) = 0, determine x(t) (or the zero input response) with provided initial condi-
tion.

2
 
1
(d) If x(t0 ) = x(1) = 0 and u(t) = 1+ u(t), determine x(t) (or the zero state
1
response).
 
1  T
(e) Determine y(t) if u(t) = 1+ u(t) and x(t0 ) = x(1) = 0 1 1 .
1

Problem 6 Obtain the controller form realization of the transfer function matrix
−1
 1 
G(s) = s+1 s+2
2 1
s+2 s2 +s+3

Problem 7 Consider the following system


     
0 1 1 0
ẋ = x+ u x(0) =
3 2 3 0
 
y= 0 1 x

1. Compute the eigenvalues of A.

2. Determine the associated eigenvectors of each eigenvalue you found.

3. Diagonalize the above system. Write T and T −1 matrix used in the diagonalization
process.

4. Calculate eAt using diagonalize approach.

5. Compute the output y(t) for a unit step u(t). Use the results from part (4) to
determine a solution for x(t), then find y(t).

6. With the matrices A, B, C determine the transfer function G(s) for the state space
representation. How do the poles of the transfer function compare to the eigenvalues
found in part (1)?

Problem 8 Consider the following system


   
3 2 −1 0
ẋ = −2 1 0
  + 0 u

4 3 1 1

Represent the system with respect to basis B = {B, AB, A2 B}

Common questions

Powered by AI

The zero-state response accounts for system behavior due only to an external input. For a linear system described by ˙x(t) = Ax(t) + Bu(t), set initial conditions x(0) = 0, then solve for x(zero-state) by integrating state equations with respect to u(t). Separately, the zero-input response considers A = matrix dynamics with x(0) as the non-zero initial condition and u(t) = 0 leading to evolution based solely on initial states. These dual analysis responses are superimposed for complete system behavior, often during responses to a sudden change like a step input .

A matrix is nilpotent if a positive integer k exists such that A^k = 0. For A to be nilpotent of order 2, A^2 must result in the zero matrix. Implications include that the system has no structural modes that persist indefinitely; dynamics linked to each mode decay to zero, reflecting in unobservable or unreachable states. In practical terms, this reveals transient-only dynamics, crucial for control schemes where lingering state presence skews desired system response, impacting controller/observer design .

Poles of the transfer function G(s) reflect the system's natural response modes and stability, directly connected to eigenvalues of the system matrix A. Poles, obtained from denominator roots of G(s), equate to eigenvalues for SISO systems, analyzing their real/imaginary components reveals stability (negative real poles) and oscillatory behavior (complex poles). A significant deviation from expected relation suggests latent states or modes affecting dynamics beyond observable outputs, pivotal in designing stable and responsive control systems. A mismatch often points to model reduction or estimation errors .

The transfer function G(s) is derived by taking the Laplace transform of the state-space representation to find Y(s)/U(s), which involves the matrices A, B, and C in G(s) = C(sI - A)^-1B + D. The degree of the transfer function corresponds to the highest power of the denominator when written in polynomial form, which is the number of poles of the system. However, this degree may differ from the number of states if the system has unattainable or unobservable states that do not influence output dynamics. This discrepancy involves controllability and observability theory and structural characteristics of matrix connections .

When u(t) = 1 + 2e^-2t, the state equation ˙x(t) = Ax(t) + Bu(t) results in a solution where the transient component disappears as t approaches infinity, leaving a steady state or causing parts of the state to grow unbounded depending on the eigenvalues' real parts. For x(t) = [x1(t) x2(t)] as t → ∞, x1 approaches a constant, while x2 diverges if the system is unstable. The output y(t) also reflects this behavior, becoming predominantly influenced by the steady-state values of the states, computed as y(∞). The specific dynamic of each variable is tied to how A's eigenvalues dictate stability .

The state transition matrix Φ(t), significant in determining system state x(t) over time, is defined as e^At where A is the system matrix. It transforms initial conditions into state behavior x(t) = Φ(t)x(0). To compute e^At, methods like matrix diagonalization (eigen-decomposition), Laplace transform for Φ(s) = (sI-A)^-1, or series expansion, are applied, depending on A's properties. This matrix's accuracy is crucial for forecasting x(t) development and resolving control inputs' effects, impacting observability and stability assessments .

Stability of the state-space model can be determined by analyzing the eigenvalues of the system matrix A. For A = [4 1; 0 10], solve λI - A = 0 to find eigenvalues λ = 4 and λ = 10. Since stable systems require all eigenvalues to have negative real parts (or zero real parts for marginal stability), this system with positive eigenvalues indicates instability. Thus, no eigenvalue has a negative real part, confirming that even a small disturbance can cause the state variables to grow indefinitely, showcasing instability .

Using different basis sets like {B, AB, A^2B} for system representation facilitates simplifying control and observability analyses by transforming the state space. This form highlights system dynamic properties: controllable subspace identification aids design by focusing on directly influenceable states, and in observable forms, states distinctly affect output derivation, streamlining observer design. This structured representation reveals clearer geometric insight into state transformations relative to input-output relation, essential in advanced control strategies for efficacy in complex dynamical systems modeling .

To compute e^At using the Laplace transform, observe that L{e^At} = (sI - A)^(-1). For matrix A = [0 1; 0 -2], (sI - A) = [s 1; 0 s+2]. The inverse (sI - A)^(-1) leads to the expression involving a partial fraction decomposition. Alternatively, the Cayley-Hamilton (C-H) Theorem states that matrix A satisfies its own characteristic polynomial, which here is λ^2 + 2λ = 0. Using this, express e^At as a polynomial in A based on the characteristic polynomial, giving a solution derived from power series expansion terms based on e^0 and e^-2t factors .

To compute the eigenvalues of matrix A = [0 1; 0 -2], we solve the characteristic equation det(A - λI) = 0, where I is the identity matrix. This results in the equation (0-λ)(-2-λ) - (1)(0) = λ^2 + 2λ. Solving λ^2+2λ = 0 gives λ = 0, -2. To find the eigenvectors, for λ = 0, solve (A - 0I)x = 0, leading to the eigenvector associated with λ = 0 as any scalar multiple of vector [1; 0]. For λ = -2, solve (A + 2I)x = 0 leading to the eigenvector as any scalar multiple of vector [1; -2].

You might also like