hw8_2011_05_26_01_solutions

hw8_2011_05_26_01_solutions - EE263 S. Lall 2011.05.26.01...

Info iconThis preview shows pages 1–3. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: EE263 S. Lall 2011.05.26.01 Homework 8 Solutions Due Thursday 6/2 at 5 PM. No extensions permitted. 1. Scalar time-varying linear dynamical system. 50110 Show that the solution of x ( t ) = a ( t ) x ( t ), where x ( t ) R , is given by x ( t ) = exp parenleftbiggintegraldisplay t a ( ) d parenrightbigg x (0) . (You can just differentiate this expression, and show that it satisfies x ( t ) = a ( t ) x ( t ).) Find a specific example showing that the analogous formula does not hold when x ( t ) R n , with n > 1. Solution. Differentiating the given expression, we obtain x ( t ) = parenleftbigg d dt integraldisplay t a ( ) d parenrightbigg exp parenleftbiggintegraldisplay t a ( ) d parenrightbigg x (0) = a ( t ) exp parenleftbiggintegraldisplay t a ( ) d parenrightbigg x (0) = a ( t ) x ( t ) . For the second part, we look for a counterexample with x ( t ) R 2 . We let A ( t ) = braceleftbigg A 1 t < 1 A 2 t 1 , where A 1 = bracketleftbigg 1 bracketrightbigg , A 2 = bracketleftbigg 1 bracketrightbigg . Then we have x (2) = (exp A 2 )(exp A 1 ) x (0) = bracketleftbigg 1 1 1 bracketrightbiggbracketleftbigg 1 1 1 bracketrightbigg x (0) = bracketleftbigg 1 1 1 2 bracketrightbigg x (0) . The formula above gives x (2) = exp( A 1 + A 2 ) x (0) = bracketleftbigg 1 1 bracketrightbigg x (0) = bracketleftbigg 1 . 5431 1 . 1752 1 . 1752 1 . 5431 bracketrightbigg x (0) . Choosing almost any x (0) ( e.g. , x (0) = e 1 ) will give us a contradiction. 2. Toeplitz Matrices 0710 (a) Write a routine sys_toeplitz.m which computes the Toeplitz matrix corresponding to a discrete-time linear dynamical system. In particular, recall that, if x ( t + 1) = Ax ( t ) + Bu ( t ) y ( t ) = Cx ( t ) + Du ( t ) 1 EE263 S. Lall 2011.05.26.01 and the initial condition x (0) = 0, then y (0) y (1) y (2) . . . y ( N ) = T u (0) u (1) u (2) . . . u ( N ) where T is the Toeplitz matrix T = D CB D CAB CB D . . . . . . CA N- 1 B CA N- 2 B ... CB D We will use this matrix to find minimum-norm controllers for linear systems. Since we often are only concerned about the trajectory at particular times, it is useful to construct only those rows and columns of T necessary. The function should take arguments of the form T=sys_toeplitz(A,B,C,D,out_times,in_times); where out_times and in_times are row vectors. For example, if out_times=[1,2,4]; in_times=[0:10] then sys_toeplitz should return just the 3 11 block matrix T so that y (1) y (2) y (4) = T u (0) u (1) ....
View Full Document

This note was uploaded on 08/06/2011 for the course EE 263 taught by Professor Boyd,s during the Summer '08 term at Stanford.

Page1 / 13

hw8_2011_05_26_01_solutions - EE263 S. Lall 2011.05.26.01...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online