AE470F23Homework6

pdf

School

University of Michigan *

*We aren’t endorsed by this school

Course

470

Subject

Mechanical Engineering

Date

Dec 6, 2023

Type

pdf

Pages

3

Uploaded by MasterIron5790

Report
AE470 Fall 2023 Homework #6 Be sure to check for the latest version of the course notes in case I post an update. Due Friday, October 20, by 11:00 PM uploaded to Canvas as a single PDF file . NOTE: There is a 45-minute grace period for late uploads. If you try to upload your homework past 11:45 PM, you will not be able to submit it. However, I will drop the lowest HW score, so if you miss just one homework, it will not affect your grade. Instructions: i ) Be sure to follow the honor code guidelines in the course information posted on the Canvas site. In particular, conceptual discussion is allowed, but all detailed work must be your own. ii ) Reminder: No use of solutions from prior offerings of this course is allowed. iii ) Symbolic computing is not allowed. iv ) Homework must be neat and professional in appearance. You may type it in Word or Latex if you wish (Professor Bernstein uses Latex.) Use a ruler to draw all lines and diagrams. Messy homework will not be graded. No crossouts of any kind may appear anywhere. v ) Put a box around your final answer to help the grader. vi ) Label your file as: LastNameAE470F23HW6.pdf vii ) Please upload a single pdf file that contains the following: (a) HW solutions (b) Figures with axes labels, captions and legends (in the case of multiple curves) (c) PDF published from Matlab script 1
Problem 1 [20 points]. Consider the URB with position output, no input, and disturbance and sensor noise ˙ x = 0 1 0 0 x + 1 2 × 1 w 1 , y = 1 0 x + w 2 . i ) Manually compute the observability matrix O ( A, C ). Is this system observable? ii ) Given the covariance matrices V 1 and V 2 V 1 = 0 . 1 0 0 0 . 1 , V 2 = 1 , determine the solution Q to the estimator Riccati equation. Hint: You can still use the MATLAB function care , but recall the difference between the Riccati equation for LQR and the Riccati equation used for estimation. iii ) Using Q from part ( ii ), calculate the observer gain matrix F associated with this observer problem. iv ) Now given the observer dynamics ˙ ˆ x = A ˆ x + F ( y ˆ y ) , ˆ y = C ˆ x, determine the numerical values of the matrices ˜ A and ˜ C such that ˙ x ˆ ˙ x = ˜ A x ˆ x + 1 2 × 1 0 2 × 1 w 1 + 0 2 × 1 F w 2 , y ˆ y = ˜ C x ˆ x + 1 0 w 2 . What are the eigenvalues of ˜ A ? v ) Initializing the system dynamics with x (0) = [1 1] T and the observer dynamics with ˆ x (0) = [0 0] T , simulate the system for t [0 , 15] with no disturbance or sensor noise ( w 1 = w 2 = 0). On one plot include the actual trajectory x 1 and the estimated trajectory ˆ x 1 . Do these trajectories converge? On another plot include the actual trajectory x 2 and the estimated trajectory ˆ x 2 . Do these trajectories converge? vi ) Adjust V 1 and V 2 such that the trajectories of ˆ x converge to the actual trajectories of x within 5 seconds. Using all the same conditions as in part ( v ), simulate the system with your chosen V 1 and V 2 and provide similar plots to those in part ( v ). Problem 2 [25 points]. Consider the UO with position output, no input, and disturbance and sensor noise ˙ x = 0 1 1 0 x + 1 2 × 1 w 1 , (1) y = 1 0 x + w 2 . (2) i ) Manually compute the observability matrix O ( A, C ). Is this system observable? ii ) Now given the observer dynamics ˙ ˆ x = A ˆ x + F ( y ˆ y ) , ˆ y = C ˆ x, determine covariance matrices V 1 and V 2 such that ˆ x converges to x within 5 seconds. Initialize the system dynamics with x (0) = [2 1] T and the observer dynamics with ˆ x (0) = x (0), and assume w 1 = w 2 = 0. Plot the trajectories of x 1 and ˆ x 1 on one figure and x 2 and ˆ x 2 on another figure. 2
iii ) Now consider the same UO model but with velocity output y = x 2 and the V 1 and V 2 covariance matrices V 1 = 0 . 1 0 . 1 0 . 1 0 . 1 , V 2 = 0 . 01 . Consider the same observer described in part ( ii ) with the same initial conditions for both the system and observer but with only w 1 = 0. Let w 2 be modeled as white noise with variance V 2 . Plot the trajectories of x 1 and ˆ x 1 on one figure and x 2 and ˆ x 2 on another figure. Hint: In MATLAB, σ × randn(1,k) generates a normally distributed vector of length k with variance σ 2 . iv ) Repeat part ( iii ) with the same covariance matrices and initial conditions. However, consider 2 separate cases for w 2 . One case where w 2 has a variance >> V 2 and one where w 2 has a variance << V 2 . Based on these results and those from part ( iii ), discuss how the variance of the sensor noise w 2 affects the estimates ˆ x . Hint: Recall what V 2 represents. Problem 3 [30 points]. Consider the UO again with position output, but with an input, disturbance, and sensor noise ˙ x = 0 1 1 0 x + 0 1 u + 1 2 × 1 w 1 , y = 1 0 x + w 2 . i ) Manually compute the controllability matrix C ( A, B ). Is this system controllable? ii ) Given initial conditions x (0) = [2 1] T , design R 1 and R 2 weighting matrices such that the output y converges within 20 seconds and the input u does not exceed ± 1. Plot the input u and output y . What are the eigenvalues of your closed loop A + BK ? Hint: Assume no disturbance or measurement noise. iii ) Given the observer ˙ ˆ x = A ˆ x + Bu + F ( y ˆ y ) , ˆ y = C ˆ x, design V 1 and V 2 weighting matrices such that the estimate ˆ y converges to y within 2.5 seconds given initial conditions x (0) = [2 1] T and ˆ x (0) = [0 0] T . For simplicity, assume all elements of V 1 are the same. Assume no input, disturbance, or sensor noise when testing your simulation. On a single figure, plot y and ˆ y . What are the eigenvalues of your closed observer loop A FC ? iv ) Provide the numerical values for the matrices ˜ A , ˜ B , ˜ C , E , and G for the augmented system ˙ x ˙ ˆ x = ˜ A x ˆ x + ˜ Bu + E w 1 w 2 , (3) y ˆ y = ˜ C x ˆ x + G w 1 w 2 . (4) v ) Using your gain matrix K from part ( ii ) where u = K ˆ x , reduce (3) to closed-loop form ˙ x ˙ ˆ x = ˜ A CL x ˆ x + E w 1 w 2 , (5) y ˆ y = ˜ C x ˆ x + G w 1 w 2 . (6) What are the eigenvalues of ˜ A CL ? Are they the same as the eigenvalues of A + BK and A FC ? vi ) Simulate your closed loop system (5) and (6) given the same initial conditions in part ( iii ). Include disturbance and sensor noise w 1 and w 2 where w 1 and w 2 can be modeled as white noise both with variance σ 2 = 0 . 001. On one figure plot y and ˆ y . On a separate figure, plot the input u . How do the combined results compare to the individual design considerations? 3
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help