final

\(\def\Real{\mathbb{R}}
\def\Comp{\mathbb{C}}
\def\Rat{\mathbb{Q}}
\def\Field{\mathbb{F}}
\def\Fun{\mathbf{Fun}}
\def\e{\mathbf{e}}
\def\f{\mathbf{f}}
\def\bv{\mathbf{v}}
\)
Due by noon 12.18.

Numbers in brackets indicate the points for each (sub)problem.

Questions? Misprints suspected? Post as a comment here!

  1. Consider the vector field \(\bv\) given by
    \[
    \dot{x}=-xy^2-x^5;\\
    \dot{y}=-x^2y-y^5.
    \]

    1. [3] Find a Lyapunov function and prove thus that the origin is Lyapunov stable. Is it asymptotically stable? Is it exponentially stable?
    2. [1] Show that one can find an arbitrarily small linear perturbation of the vector field \(\bv\) (that is add arbitrarily small linear terms to \(\bv\) so that the origin becomes an unstable equilibrium.
    3. [3] What is the maximal degree of perturbation of the vector field that can make the origin unstable?
  2. Consider the system given by
    \[
    \begin{array}{ccl}
    \dot{x}_1&=&a_1(t)x_2,\\
    \dot{x}_2&=&u,\\
    \dot{x}_3&=&a_3(t)x_2,
    \end{array}
    \]
    where
    \[
    a_1=\left\{\begin{array}{ccc}
    0 &\mathrm{for}& 0\leq t<1;\\
    1 &\mathrm{for}& 1\leq t\leq 2;
    \end{array}\right. \mathrm{ and }\quad
    a_3=\left\{\begin{array}{ccc}
    1 &\mathrm{for}& 0\leq t<1;\\
    0 &\mathrm{for}& 1\leq t\leq 2;
    \end{array}\right.
    \]

    1. [3] Find the controllability Grammian \(W(0,t), 0\leq t\leq 2\).
    2. [1] Is this system controllable on \([0,1]\)? on \([1,2]\)?
    3. [2] Is this system controllable on \([0,2]\)?
  3. [5] Find the rank of the minimal realization of the following transfer function (in this \(101\times 101\)-matrix, the entries on the $k$-th anti-diagonal are \(1/(s+k-1)\)):
    \[
    P(s)=\left(
    \begin{array}{ccccc}
    \frac{1}{s}&\frac{1}{s+1}&\frac{1}{s+2}&\cdots&\frac{1}{s+100}\\
    \frac{1}{s+1}&\frac{1}{s+2}&\frac{1}{s+3}&\cdots&\frac{1}{s+101}\\
    \frac{1}{s+2}&\frac{1}{s+3}&\frac{1}{s+4}&\cdots&\frac{1}{s+102}\\
    \vdots&\vdots&\vdots&\ddots&\vdots\\
    \frac{1}{s+100}&\frac{1}{s+101}&\frac{1}{s+102}&\cdots&\frac{1}{s+200}
    \end{array}\right)
    \]
  4. Consider the optimal control problem
    \[
    \dot{x}_1=x_2;\\
    \dot{x}_2=x_1+u;\\
    |u|\leq 1; L\equiv 1; (x_1(0),x_2(0))=(x_o,v_o).
    \]

    1. [4] For \((x_o,v_o)=(1,-1)\), find the trajectory stopping at \((0,0)\) at the shortest time. (This is the shortest time problem, so the terminal time is not fixed, hence the Hamiltonian is \(\equiv 0\).)
    2. [3]Find for which \((x_o,v_o)\) it is possible at all to stop at \((0,0)\).
  5. [5] Consider the infinite horizon LQR
    \[
    \begin{array}{ccl}
    \dot{x}_1&=&x_2+b_1u;\\
    \dot{x}_2&=&-x_1+b_2u;
    \end{array}\\
    \int_0^\infty u^2+x_1^2 dt\to\min.
    \]
    Assume that the control vector \(b\) has norm one, i.e.
    \[
    b_1^2+b_2^2=1.
    \]
    Among such vectors, find one providing the speediest convergence of the close-loop feedback LQR regulator, that is minimizing the (negative) maximal real part of \(\dot{x}=(A-bk)x\).

18 Responses to final

  1. Joshua Aurich December 14, 2014 at 3:48 pm #

    On question 2(a),
    When you say from (0,t), do you mean from (0,2) or do you mean to some undecided t?
    If this is for some undecided t, can we assume a1(t > 2) = 1 and a3(t > 2) = 0?

    • yuliy December 14, 2014 at 4:29 pm #

      Only for t between 0 and 2.

  2. Qian Zhang December 15, 2014 at 8:25 pm #

    For the question 5, is it okay to just get a numerical value of b1 and b2? It seems a little bit complex to calculate in analytical way.

    • yuliy December 16, 2014 at 8:05 am #

      Sure.

  3. Zeyu Zhou December 16, 2014 at 2:48 pm #

    I don’t know if it is ok to ask but when you say a linear perturbation to the vector field v, do you mean: \dot(v) = f(v) + \epsilon * v, where epsilon is some small positive number ?

    • yuliy December 17, 2014 at 12:00 am #

      Arbitrarily small (in absolute value) number.

  4. Vishwa Shah December 16, 2014 at 3:35 pm #

    For question 5, I’m not quite sure what you mean when you say ‘ minimizing the (negative) maximal real part of x_dot = (A-bk)x’
    Do you mean, the largest negative value and take it to zero?
    or take the largest negative value and take it further away from zero?

    • yuliy December 16, 2014 at 11:59 pm #

      Further away from the origin (as stable as possible).

  5. Zeyu Zhou December 16, 2014 at 9:29 pm #

    In Problem 3, you asked for “the rank of the minimal realization”. Does it mean “the degree of the minimal realization” or “the rank of matrix A of the minimal realization” ?

    • yuliy December 16, 2014 at 11:59 pm #

      Latter (i.e. dimension of the state space of the minimal realization).

  6. Shankar Deka December 16, 2014 at 10:02 pm #

    Q4. (b) Do want to reach the origin only using the optimal control or do you want to just get there by ANY control (irrespective of whether or not it is optimal, as long as we get there) ?

    • yuliy December 16, 2014 at 11:58 pm #

      Reach the origin in minimal time – that is optimally for L=1.

  7. SEUNG WHAN CHUNG December 16, 2014 at 10:16 pm #

    Does L mean the Lagrangian of the system?

    • yuliy December 16, 2014 at 11:57 pm #

      Yes.

  8. Yuting Ng December 16, 2014 at 11:51 pm #

    Is it possible for there be a hint for Q4a? I can’t seem to solve it using the method from Page 222-224 of the class notes… Thank you!

    • yuliy December 16, 2014 at 11:56 pm #

      Form long Hamiltonian, then the short one (maximizing in the control). Write down the differential equations for the covectors. Figure out how many switches can there be.

  9. James December 20, 2014 at 11:18 pm #

    Will solutions be posted?

    • yuliy December 21, 2014 at 1:03 am #

      Yeah, hopefully next week.

Leave a Reply