Systems of Differential Equations — Question 7

PDF ↗

Question 7

A controlled third derivative moves a state from rest to rest in unit time: y‴=u(t),(y,y′,y″)(0)=(0,0,0),(y,y′,y″)(1)=(1,0,0).y'''=u(t),\qquad (y,y',y'')(0)=(0,0,0),\qquad (y,y',y'')(1)=(1,0,0). For square-integrable controls uu on [0,1][0,1], interpret the derivative chain through its integral solution. The effort is J(u)=∫01u(t)2dtJ(u)=\int_0^1u(t)^2\,dt.

Tasks

  1. Write the first-order state system and its transition matrix. Derive the three terminal constraints as integrals of uu.

  2. Find the unique quadratic control satisfying the terminal conditions. Integrate it to obtain the displacement.

  3. Prove this quadratic is the unique minimum-effort control among all square-integrable admissible controls, up to equality almost everywhere. Compute the minimum effort.

  4. Verify both endpoint states and explain why the minimization proof covers controls that are not polynomials. Sketch the displacement and u/60u/60, explicitly labeling the scaling.

Original worksheet page 1: question and worked solution for 7-6-007
Show solutionHide solution

Question 7 – Solution

Strategy. Express the terminal state as three moments, then use orthogonality to prove optimality over the full control class.

Step 1: Derive the terminal moments. For X=(y,y′,y″)TX=(y,y',y'')^T, X′=NX+e3uX'=NX+e_3u, where NN has ones just above the diagonal and zeros elsewhere. Since N3=0N^3=0, eNr=(1rr2/201r001),X(1)=∫01((1−t)2/21−t1)u(t)dt.e^{Nr}=\begin{pmatrix}1&r&r^2/2\\0&1&r\\0&0&1\end{pmatrix},\qquad X(1)=\int_0^1\begin{pmatrix}(1-t)^2/2\\1-t\\1\end{pmatrix}u(t)\,dt. The required three moments are (1,0,0)T(1,0,0)^T.

Step 2: Solve within the quadratic space. For u=a+bt+ct2u=a+bt+ct^2, the constraints are a/6+b/24+c/60=1,a/2+b/6+c/12=0,a+b/2+c/3=0.a/6+b/24+c/60=1,\quad a/2+b/6+c/12=0,\quad a+b/2+c/3=0. Their solution is u*=60−360t+360t2,y*=10t3−15t4+6t5.\boxed{u_*=60-360t+360t^2,\qquad y_*=10t^3-15t^4+6t^5.} The moment kernels span all quadratics. A quadratic in the homogeneous null space is orthogonal to itself and hence zero, proving uniqueness here.

Step 3: Prove global minimum effort. For any admissible u=u*+vu=u_*+v, all three terminal moments of vv vanish. Because u*u_* is in their span, ∫01u*v=0\int_0^1u_*v=0. Therefore J(u)=J(u*)+∫01v2dt≥J(u*)=720.\boxed{J(u)=J(u_*)+\int_0^1v^2\,dt\geq J(u_*)=720.} Equality holds exactly when v=0v=0 almost everywhere. This argument only uses square integrability and moment constraints, not polynomial form for vv.

Step 4: Verify the endpoint motion. We have y*′=30t2(1−t)2y_*'=30t^2(1-t)^2 and y*″=60t(1−t)(1−2t)y_*''=60t(1-t)(1-2t), so the two endpoint states are exactly (0,0,0)(0,0,0) and (1,0,0)(1,0,0); differentiating again gives u*u_*. The graph uses u*/60u_*/60 so the control and displacement are both readable.

See the diagram in the original worksheet below.

Original worksheet page 2: question and worked solution for 7-6-007

Original worksheet layout. Use Enlarge or open the PDF for a closer view.