Problem Set 4 for MAE280A Linear Systems Theory, Fall 2017: due Tuesday November 28, 2017, in class
|
|
- Sylvia Nichols
- 6 years ago
- Views:
Transcription
1 Problem Set 4 for MAE8A Linear Systems Theory, Fall 7: due Tuesday November 8, 7, in class Problem Controllability and reachability Consider the C-circuit depicted in class on Slide 3 of the Controllability lecture sequence. u + _ C + C + v v - - C-circuit example state is the set of capacitor voltage! apple state equation ẋ(t) = C C The state-space system is described by Transition function/matrix exponential 3 / C! 4 C " # 5t ẋ t = / C x e C t / C t + u /! e C t. C = e C t Solution " # e C t Z " For our particular choice of resistors and capacitors, we rewrite the state equations in units of seconds t x x(t) = () 4 as! e C t + C e C (t x () C e C (t ẋ t = x t + u t. Find a control signal, u If C=C t, for t [, ] that transfers the state fromand applex()= then the state is always a sca So [t,t ]= : 5 x =, to x Likewise, =. to transfer 7 x to the origin =e C t x() +! apple [Bonus] Find the minimum-energy control signal. Compute the energy. So C[t,t ]= as well : Hint for both parts: look at the Slides 4 and 5 of the Controllability notes. Solution: The solution for LTI system is When C C we will find that t x(t) = e A(t t) x(t ) + e A(t τ) Bu(τ)dτ, t we substitute the matrix A, B and state of the systems at two different time. 5 x() = = e A( ) + e A( τ) u(τ)dτ, 7 since matrix A is diagonal then we have e At e t = e t and also, e A( τ) = Example e ( τ) e (τ ) e ( τ) = e (τ ). [t,t ]=C[t,t ]=
2 Now let us check if the system is reachable, C = [ B AB ] = rank of matrix C is, hence the system is reachable. 4 Use the solution equation above and see that we want to drive the zero-state solution x zs (t) = e A( τ) u(τ)dτ, 5 from x zs () = to x zs () = e 7 A. Now we are allow to apply results from Lecture 4, W (, ) = e A( τ) [ ] [ ] e A( τ) T dτ = e (τ ) [ e (τ ) ] e (τ ) e (τ ) dτ W (, ) = [ e (τ ) e (τ ) ] [e (τ ) e (τ )] [ e (τ ) e dτ = 3(τ ) ] dτ e 3(τ ) 4e 4(τ ) W (, ) = [ ( e )/ ( e 3 ] )/3 ( e 3 )/3 ( e 4 ) According to Lecture 4, we need to find η such that 5 x zs () = e A = W 7 (, )η ImW (, ), then, for t [, ] control function, U(t) = B(t) T e A( t)t η = [ ] [ e t ] e (t ) η. The minimum energy is E = U(τ) dτ = x zs () T W (, ) x zs (). Problem Observers with real aeroengine data The file signal.dat contains data from an experiment at United Technologies esearch Center examining combustion instabilities in a jet engine single-nozzle rig. The combustion instability manifests as an audible single tone, sometimes called hooting or howling coming from the engine. The instability is undesirable and causes NO x pollution. The data file contains three columns of numbers: the current time, the pressure signal inside the combustor, and the heat release rate measured with a thermal sensor. You task is to build an observer to estimate the time delay between the peaks in the pressure signal and the peaks in the heat release rate signal. You need to build an observer to do this.
3 (i) Plot the data versus time to see the main oscillation in the data is around Hz. You can use MATLAB s play function to listen to it, since the data is sampled at 5kHz. Here is the zoomed plot of pressure versus time. Note that there are about.5 cycles in.6 seconds or 8 Hz approximately pressure (units) Figure : Plot for task (i). (ii) Plot the absolute value of the FFT of the pressure data to show that the main frequency is at.573 Hz. Here a zoomed in view using stem plot of the abs(fft(pressure)). You see that bin 69 is the max. This represents 69 cycles in the data of 6384 samples at 5 Hz or f = =.573Hz. (iii) Convert the heat release rate, which is necessarily positive, to a more symmetric signal by taking its logarithm. Call this lheat. (iv) Show that, the signals x i (t) = A cos(ωt + φ) and x q (t) = A sin(ωt + φ) satisfy the state equation xi (t + ) cos ω sin ω xi (t) =. x q (t + ) sin ω cos ω x q (t) x i (t + ) = A cos(ω(t + ) + φ) = A cos(ωt + φ + ω) = A[cos(ωt + φ) cos(ω) sin(ωt + φ) sin(ω)] = x i (t) cos(ω) x q (t) sin(ω) and x q (t + ) = A sin(ω(t + ) + φ) = A sin(ωt + φ + ω) = A[sin(ωt + φ) cos(ω) + cos(ωt + φ) sin(ω)] = x q (t) sin(ω) + x i (t) cos(ω)
4 amplitude bin number + Figure : Plot for task (ii). so we have, xi (t + ) = x q (t + ) [ cos ω sin ω sin ω cos ω ] [ xi (t) x q (t) ] Note: since we are incrementing time, not by second, but by. seconds, we need to take ω = π.573. =.646 c /s.
5 (v) Accordingly, argue for the following being a state-space realization for pressure and for lheat. x p i (t + ) cos ω sin ω x p x p i (t) q(t + ) x h i (t + ) = sin ω cos ω x p q(t) cos ω sin ω x h x h i (t), q (t + ) sin ω cos ω x h q (t) x p i (t) pressure(t) = x p q(t) lheat(t) x h i (t). x h q (t) That is, this model describes pressure= A cos(ωt + φ p ) and lheat= B cos(ωt + φ h ) with the same frequency ω but differing amplitudes A and B and differing phases φ p and φ h. Note that this model has no input signal, since the audible oscillation is self-supporting. This is fairly obvious, since both pressure and lheat are described by separate signals with the same frequency. (vi) If you had the full dimension 4 state x(t) available, determine an approach to compute A, B, ωt + φ p and ωt + φ h at that time from the state elements. From part (iv) we see that A = (x p i (t)) + (x p q(t)) and ωt + φ p = atan(x p q(t), x p i (t)). and also B = (x h i (t)) + (x h q (t)) and ωt + φ h = atan(x h q (t), x h i (t)). (vii) Work a bit harder and derive a formula for estimating φ h φ p as a function of the state. Work even harder still and express this phase difference at this ω as a time delay. Start with the following trigonometric formulæ. So, AB cos(φ h φ p ) = AB cos[(ωt + φ h ) (ωt + φ p )], = AB cos(ωt + φ h ) cos(ωt + φ p ) + AB sin(ωt + φ h ) sin(ωt + φ p ), = x h i (t) x p i (t) + xh q (t) x p q(t). AB sin(φ h φ p ) = AB sin[(ωt + φ h ) (ωt + φ p )], = AB sin(ωt + φ h ) cos(ωt + φ p ) AB cos(ωt + φ h ) sin(ωt + φ p ), = x h q (t) x p i (t) xh i (t) x p q(t). φ h φ p = atan(ab sin(φ h φ p ), AB cos(φ h φ p )), [ ] = atan x h q (t) x p i (t) xh i (t) x p q(t), x h i (t) x p i (t) + xh q (t) x p q(t). Since this is a radian phase difference of a signal at radian frequency ω radians/sec, the time difference is τ = φh φ p. ω
6 (viii) Using the value ω = π.573 rad/s, write this model in MATLAB. Call the matrices A and C as usual. Do not forget to include the sampling time of T s =. s. (ix) Show that the entire state is observable from the two outputs for this model. We can check the rank of the observability matrix, which is the controllability matrix of A and C. rank(ctrb(a,c ))=4 (x) Select an appropriate place for observer poles and determine using place an observer output-injection gain. Note that this is all in discrete time. I tried two selections:.98 eig(a) and.975 eig(a). Both were found using commands like L=place(A,C,.98*eig(A)). (xi) Write the MATLAB code to do the following. (a) Start from an initial state estimate of, say, zero and run this observer, taking the two measurements, pressure(t) and lheat(t) as inputs and yielding state estimate ˆx(t) as its output. In general we have observer In this problem we do have u(t) = so ˆx(t + ) = (A LC)ˆx(t) + Bu(t) + Ly(t). ˆx(t + ) = (A LC)ˆx(t) + Ly(t) Since we want all elements of the state estimate, we have Ā = (A LC), B = L, C = eye(4), D = and our observer is defined as the following system, observ=ss(a-lc, L,eye(4),zeros(4,),/5) Here are the zoomed plots of state estimates with each set of poles. (b) Apply your formula from step?? to compute the time delay estimate between pressure and heat release rate over the entire data set. (xii) Plot your estimate over time. Let s take things slowly. Here are the sin term plot(time,xhq.*xpi-xhi.*xpq) and the cos term plot(time,xhi.*xpi+xhq.*xpq), where xpi = xhat(:,); xpq = xhat(:,); xhi = xhat(:,3); xhq = xhat(:,4); Now consider the computed phase difference: atan(xhq.*xpi-xhi.*xpq,xhi.*xpi+xhq.*xpq). We see that the computed phase difference is jumping over π a lot. Let s unwrap it. This suggests that the computed phase difference via atan(.,.) is unreliable, since it is drifting over time. Take a backward step and use the mean sin and cos terms to compute
7 5 poles.98*eig(a) 4 3 xhat atan(mean(xhq.*xpi-xhi.*xpq),mean(xhi.*xpi+xhq.*xpq)) = Now add the phase difference from the differentiator, π, and the negatively sloped memoryless nonlinearity, π, to get a total phase difference of 4.73 radians or a time delay of τ = 4.73/(ω.) =.36 seconds. Hint: I get something like delay τ=3.7 ms on average. Welcome to the world of state estimation. The time delay is a central modeling parameter describing the propagation time of the fuel exiting the nozzle to the flame front.
8 4 poles.975*eig(a) 3 xhat ABsin( h - p ) plot ABsin( h - p ) Figure 3: AB sin plot: mean e-35.
9 .5-3 ABcos( h - p ) plot.5 h - p ) ABcos( Figure 4: AB cos plot: mean e computed phase difference plot 3 computed phase difference Figure 5: φ h φ p plot: mean
10 Unwrapped computed phase difference plot 9 8 unwrapped phase difference Figure 6: Unwrappped φ h φ p plot.
MAE143A Signals & Systems, Final Exam - Wednesday March 16, 2005
MAE13A Signals & Systems, Final Exam - Wednesday March 16, 5 Instructions This quiz is open book. You may use whatever written materials you choose including your class notes and the textbook. You may
More informationHOMEWORK 4: MATH 265: SOLUTIONS. y p = cos(ω 0t) 9 ω 2 0
HOMEWORK 4: MATH 265: SOLUTIONS. Find the solution to the initial value problems y + 9y = cos(ωt) with y(0) = 0, y (0) = 0 (account for all ω > 0). Draw a plot of the solution when ω = and when ω = 3.
More informationBasics of Electric Circuits
António Dente Célia de Jesus February 2014 1 Alternating Current Circuits 1.1 Using Phasors There are practical and economic reasons justifying that electrical generators produce emf with alternating and
More informationEM Oscillations. David J. Starling Penn State Hazleton PHYS 212
I ve got an oscillating fan at my house. The fan goes back and forth. It looks like the fan is saying No. So I like to ask it questions that a fan would say no to. Do you keep my hair in place? Do you
More informationName:... Section:... Physics 208 Quiz 8. April 11, 2008; due April 18, 2008
Name:... Section:... Problem 1 (6 Points) Physics 8 Quiz 8 April 11, 8; due April 18, 8 Consider the AC circuit consisting of an AC voltage in series with a coil of self-inductance,, and a capacitor of
More informationChapter 10: Sinusoidal Steady-State Analysis
Chapter 10: Sinusoidal Steady-State Analysis 1 Objectives : sinusoidal functions Impedance use phasors to determine the forced response of a circuit subjected to sinusoidal excitation Apply techniques
More informationLecture 35: FRI 17 APR Electrical Oscillations, LC Circuits, Alternating Current I
Physics 3 Jonathan Dowling Lecture 35: FRI 7 APR Electrical Oscillations, LC Circuits, Alternating Current I Nikolai Tesla What are we going to learn? A road map Electric charge è Electric force on other
More informationModule 25: Outline Resonance & Resonance Driven & LRC Circuits Circuits 2
Module 25: Driven RLC Circuits 1 Module 25: Outline Resonance & Driven LRC Circuits 2 Driven Oscillations: Resonance 3 Mass on a Spring: Simple Harmonic Motion A Second Look 4 Mass on a Spring (1) (2)
More informationChapter 3: Capacitors, Inductors, and Complex Impedance
hapter 3: apacitors, Inductors, and omplex Impedance In this chapter we introduce the concept of complex resistance, or impedance, by studying two reactive circuit elements, the capacitor and the inductor.
More informationCircuit Analysis-III. Circuit Analysis-II Lecture # 3 Friday 06 th April, 18
Circuit Analysis-III Sinusoids Example #1 ü Find the amplitude, phase, period and frequency of the sinusoid: v (t ) =12cos(50t +10 ) Signal Conversion ü From sine to cosine and vice versa. ü sin (A ± B)
More informationThis is the number of cycles per unit time, and its units are, for example,
16 4. Sinusoidal solutions Many things in nature are periodic, even sinusoidal. We will begin by reviewing terms surrounding periodic functions. If an LTI system is fed a periodic input signal, we have
More informationModule 03 Linear Systems Theory: Necessary Background
Module 03 Linear Systems Theory: Necessary Background Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html September
More informationLecture 7. Please note. Additional tutorial. Please note that there is no lecture on Tuesday, 15 November 2011.
Lecture 7 3 Ordinary differential equations (ODEs) (continued) 6 Linear equations of second order 7 Systems of differential equations Please note Please note that there is no lecture on Tuesday, 15 November
More informationPhysics 11b Lecture #15
Physics 11b ecture #15 and ircuits A ircuits S&J hapter 3 & 33 Administravia Midterm # is Thursday If you can t take midterm, you MUST let us (me, arol and Shaun) know in writing before Wednesday noon
More informationPhysics 4B Chapter 31: Electromagnetic Oscillations and Alternating Current
Physics 4B Chapter 31: Electromagnetic Oscillations and Alternating Current People of mediocre ability sometimes achieve outstanding success because they don't know when to quit. Most men succeed because
More informationto have roots with negative real parts, the necessary and sufficient conditions are that:
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 543 LINEAR SYSTEMS AND CONTROL H O M E W O R K # 7 Sebastian A. Nugroho November 6, 7 Due date of the homework is: Sunday, November 6th @ :59pm.. The following
More informatione jωt y (t) = ω 2 Ke jωt K =
BME 171, Sec 2: Homework 2 Solutions due Tue, Sep 16 by 5pm 1. Consider a system governed by the second-order differential equation a d2 y(t) + b dy(t) where a, b and c are nonnegative real numbers. (a)
More informationAC analysis. EE 201 AC analysis 1
AC analysis Now we turn to circuits with sinusoidal sources. Earlier, we had a brief look at sinusoids, but now we will add in capacitors and inductors, making the story much more interesting. What are
More informationControl Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli
Control Systems I Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)
More informationE2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)
E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,
More informationLecture XXVI. Morris Swartz Dept. of Physics and Astronomy Johns Hopkins University November 5, 2003
Lecture XXVI Morris Swartz Dept. of Physics and Astronomy Johns Hopins University morris@jhu.edu November 5, 2003 Lecture XXVI: Oscillations Oscillations are periodic motions. There are many examples of
More informationMAE143A Signals & Systems - Homework 5, Winter 2013 due by the end of class Tuesday February 12, 2013.
MAE43A Signals & Systems - Homework 5, Winter 23 due by the end of class Tuesday February 2, 23. If left under my door, then straight to the recycling bin with it. This week s homework will be a refresher
More informationECE 201 Fall 2009 Final Exam
ECE 01 Fall 009 Final Exam December 16, 009 Division 0101: Tan (11:30am) Division 001: Clark (7:30 am) Division 0301: Elliott (1:30 pm) Instructions 1. DO NOT START UNTIL TOLD TO DO SO.. Write your Name,
More information2.3 Oscillation. The harmonic oscillator equation is the differential equation. d 2 y dt 2 r y (r > 0). Its solutions have the form
2. Oscillation So far, we have used differential equations to describe functions that grow or decay over time. The next most common behavior for a function is to oscillate, meaning that it increases and
More informationProf. Shayla Sawyer CP08 solution
What does the time constant represent in an exponential function? How do you define a sinusoid? What is impedance? How is a capacitor affected by an input signal that changes over time? How is an inductor
More informationModule 08 Observability and State Estimator Design of Dynamical LTI Systems
Module 08 Observability and State Estimator Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha November
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics 8.02 Spring 2003 Experiment 17: RLC Circuit (modified 4/15/2003) OBJECTIVES
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics 8. Spring 3 Experiment 7: R Circuit (modified 4/5/3) OBJECTIVES. To observe electrical oscillations, measure their frequencies, and verify energy
More informationDriven RLC Circuits Challenge Problem Solutions
Driven LC Circuits Challenge Problem Solutions Problem : Using the same circuit as in problem 6, only this time leaving the function generator on and driving below resonance, which in the following pairs
More informationDynamic circuits: Frequency domain analysis
Electronic Circuits 1 Dynamic circuits: Contents Free oscillation and natural frequency Transfer functions Frequency response Bode plots 1 System behaviour: overview 2 System behaviour : review solution
More informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationFirst and Second Order Circuits. Claudio Talarico, Gonzaga University Spring 2015
First and Second Order Circuits Claudio Talarico, Gonzaga University Spring 2015 Capacitors and Inductors intuition: bucket of charge q = Cv i = C dv dt Resist change of voltage DC open circuit Store voltage
More informationNotes on the Periodically Forced Harmonic Oscillator
Notes on the Periodically orced Harmonic Oscillator Warren Weckesser Math 38 - Differential Equations 1 The Periodically orced Harmonic Oscillator. By periodically forced harmonic oscillator, we mean the
More informationA system that is both linear and time-invariant is called linear time-invariant (LTI).
The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped
More informationChapter 33. Alternating Current Circuits
Chapter 33 Alternating Current Circuits 1 Capacitor Resistor + Q = C V = I R R I + + Inductance d I Vab = L dt AC power source The AC power source provides an alternative voltage, Notation - Lower case
More information4. Sinusoidal solutions
16 4. Sinusoidal solutions Many things in nature are periodic, even sinusoidal. We will begin by reviewing terms surrounding periodic functions. If an LTI system is fed a periodic input signal, we have
More information1 (30 pts) Dominant Pole
EECS C8/ME C34 Fall Problem Set 9 Solutions (3 pts) Dominant Pole For the following transfer function: Y (s) U(s) = (s + )(s + ) a) Give state space description of the system in parallel form (ẋ = Ax +
More informationLaboratory handout 5 Mode shapes and resonance
laboratory handouts, me 34 82 Laboratory handout 5 Mode shapes and resonance In this handout, material and assignments marked as optional can be skipped when preparing for the lab, but may provide a useful
More information(Refer Slide Time: 00:32)
Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 12 Scilab simulation of Lotka Volterra
More informationUnstable Oscillations!
Unstable Oscillations X( t ) = [ A 0 + A( t ) ] sin( ω t + Φ 0 + Φ( t ) ) Amplitude modulation: A( t ) Phase modulation: Φ( t ) S(ω) S(ω) Special case: C(ω) Unstable oscillation has a broader periodogram
More informationChapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
More informationPHYSICS 110A : CLASSICAL MECHANICS HW 2 SOLUTIONS. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see
PHYSICS 11A : CLASSICAL MECHANICS HW SOLUTIONS (1) Taylor 5. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see 1.5 1 U(r).5.5 1 4 6 8 1 r Figure 1: Plot for problem
More informationPhysics for Scientists & Engineers 2
Electromagnetic Oscillations Physics for Scientists & Engineers Spring Semester 005 Lecture 8! We have been working with circuits that have a constant current a current that increases to a constant current
More informationMAT292 - Calculus III - Fall Solution for Term Test 2 - November 6, 2014 DO NOT WRITE ON THE QR CODE AT THE TOP OF THE PAGES.
MAT9 - Calculus III - Fall 4 Solution for Term Test - November 6, 4 Time allotted: 9 minutes. Aids permitted: None. Full Name: Last First Student ID: Email: @mail.utoronto.ca Instructions DO NOT WRITE
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationHomework 6 Solutions
8-290 Signals and Systems Profs. Byron Yu and Pulkit Grover Fall 208 Homework 6 Solutions. Part One. (2 points) Consider an LTI system with impulse response h(t) e αt u(t), (a) Compute the frequency response
More informationWave Phenomena Physics 15c
Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory
More informationWave Phenomena Physics 15c
Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory
More informationSprings: Part I Modeling the Action The Mass/Spring System
17 Springs: Part I Second-order differential equations arise in a number of applications We saw one involving a falling object at the beginning of this text (the falling frozen duck example in section
More informationOrdinary differential equations
Class 7 Today s topics The nonhomogeneous equation Resonance u + pu + qu = g(t). The nonhomogeneous second order linear equation This is the nonhomogeneous second order linear equation u + pu + qu = g(t).
More informationALTERNATING CURRENT. with X C = 0.34 A. SET UP: The specified value is the root-mean-square current; I. EXECUTE: (a) V = (0.34 A) = 0.12 A.
ATENATING UENT 3 3 IDENTIFY: i Icosωt and I I/ SET UP: The specified value is the root-mean-square current; I 34 A EXEUTE: (a) I 34 A (b) I I (34 A) 48 A (c) Since the current is positive half of the time
More information2.4 Harmonic Oscillator Models
2.4 Harmonic Oscillator Models In this section we give three important examples from physics of harmonic oscillator models. Such models are ubiquitous in physics, but are also used in chemistry, biology,
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More information12 Chapter Driven RLC Circuits
hapter Driven ircuits. A Sources... -. A ircuits with a Source and One ircuit Element... -3.. Purely esistive oad... -3.. Purely Inductive oad... -6..3 Purely apacitive oad... -8.3 The Series ircuit...
More informationLecture 15. LC Circuit. LC Oscillation - Qualitative. LC Oscillator
Lecture 5 Phys. 07: Waves and Light Physics Department Yarmouk University 63 Irbid Jordan &KDSWHUElectromagnetic Oscillations and Alternating urrent L ircuit In this chapter you will see how the electric
More informationTransform Solutions to LTI Systems Part 3
Transform Solutions to LTI Systems Part 3 Example of second order system solution: Same example with increased damping: k=5 N/m, b=6 Ns/m, F=2 N, m=1 Kg Given x(0) = 0, x (0) = 0, find x(t). The revised
More informationSinusoids. Amplitude and Magnitude. Phase and Period. CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation
Sinusoids CMPT 889: Lecture Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 6, 005 Sinusoids are
More informationChapter 21: RLC Circuits. PHY2054: Chapter 21 1
Chapter 21: RC Circuits PHY2054: Chapter 21 1 Voltage and Current in RC Circuits AC emf source: driving frequency f ε = ε sinωt ω = 2π f m If circuit contains only R + emf source, current is simple ε ε
More informationCMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation
CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 26, 2005 1 Sinusoids Sinusoids
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More information2.4 Models of Oscillation
2.4 Models of Oscillation In this section we give three examples of oscillating physical systems that can be modeled by the harmonic oscillator equation. Such models are ubiquitous in physics, but are
More informationUniversity of TN Chattanooga Physics 1040L 8/18/2012 PHYSICS 1040L LAB LAB 4: R.C. TIME CONSTANT LAB
PHYSICS 1040L LAB LAB 4: R.C. TIME CONSTANT LAB OBJECT: To study the discharging of a capacitor and determine the time constant for a simple circuit. APPARATUS: Capacitor (about 24 μf), two resistors (about
More informationObservability. It was the property in Lyapunov stability which allowed us to resolve that
Observability We have seen observability twice already It was the property which permitted us to retrieve the initial state from the initial data {u(0),y(0),u(1),y(1),...,u(n 1),y(n 1)} It was the property
More informationLecture 2 and 3: Controllability of DT-LTI systems
1 Lecture 2 and 3: Controllability of DT-LTI systems Spring 2013 - EE 194, Advanced Control (Prof Khan) January 23 (Wed) and 28 (Mon), 2013 I LTI SYSTEMS Recall that continuous-time LTI systems can be
More informationOscillations. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring of constant k is
Dr. Alain Brizard College Physics I (PY 10) Oscillations Textbook Reference: Chapter 14 sections 1-8. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring
More information2 Growth, Decay, and Oscillation
2 Growth, Decay, and Oscillation b The city of Suzhou in Jiangsu Province, China. Suzhou is the fastest growing city in the world, with an annual population growth of 6.5% between the years 2000 and 2014.1
More informationTLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall TLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall Problem 1.
TLT-5/56 COMMUNICATION THEORY, Exercise 3, Fall Problem. The "random walk" was modelled as a random sequence [ n] where W[i] are binary i.i.d. random variables with P[W[i] = s] = p (orward step with probability
More informationNetworks and Systems Prof. V. G. K. Murti Department of Electrical Engineering Indian Institution of Technology, Madras
Networks and Systems Prof. V. G. K. Murti Department of Electrical Engineering Indian Institution of Technology, Madras Lecture - 32 Network Function (3) 2-port networks: Symmetry Equivalent networks Examples
More information2B30 Formal Report Simon Hearn Dr Doel
DEPARTMENT OF PHYSICS & ASTRONOMY SECOND YEAR LAB REPORT DECEMBER 2001 EXPERIMENT E7: STUDY OF AN OSCILLATING SYSTEM DRIVEN INTO RESONANCE PERFORMED BY SIMON HEARN, LAB PARTNER CAROLINE BRIDGES Abstract
More information1. (10 points) Find the general solution to the following second-order differential equation:
Math 307A, Winter 014 Midterm Solutions Page 1 of 8 1. (10 points) Find the general solution to the following second-order differential equation: 4y 1y + 9y = 9t. To find the general solution to this nonhomogeneous
More informationKeble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II
Tomi Johnson 1 Keble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II Prepare full solutions to the problems with a self assessment of your progress on a cover
More informationInductance, Inductors, RL Circuits & RC Circuits, LC, and RLC Circuits
Inductance, Inductors, RL Circuits & RC Circuits, LC, and RLC Circuits Self-inductance A time-varying current in a circuit produces an induced emf opposing the emf that initially set up the timevarying
More informationLinear Algebra and ODEs review
Linear Algebra and ODEs review Ania A Baetica September 9, 015 1 Linear Algebra 11 Eigenvalues and eigenvectors Consider the square matrix A R n n (v, λ are an (eigenvector, eigenvalue pair of matrix A
More informationThe Harmonic Oscillator
The Harmonic Oscillator Math 4: Ordinary Differential Equations Chris Meyer May 3, 008 Introduction The harmonic oscillator is a common model used in physics because of the wide range of problems it can
More informationAlternating Currents. The power is transmitted from a power house on high voltage ac because (a) Electric current travels faster at higher volts (b) It is more economical due to less power wastage (c)
More information2.1 Exponential Growth
2.1 Exponential Growth A mathematical model is a description of a real-world system using mathematical language and ideas. Differential equations are fundamental to modern science and engineering. Many
More informationSignals and Systems. Problem Set: The z-transform and DT Fourier Transform
Signals and Systems Problem Set: The z-transform and DT Fourier Transform Updated: October 9, 7 Problem Set Problem - Transfer functions in MATLAB A discrete-time, causal LTI system is described by the
More information/ \ ( )-----/\/\/\/ \ / In Lecture 3 we offered this as an example of a first order LTI system.
18.03 Class 17, March 12, 2010 Linearity and time invariance [1] RLC [2] Superposition III [3] Time invariance [4] Review of solution methods [1] We've spent a lot of time with mx" + bx' + cx = q(t). There
More informationELECTROMAGNETIC OSCILLATIONS AND ALTERNATING CURRENT
Chapter 31: ELECTROMAGNETIC OSCILLATIONS AND ALTERNATING CURRENT 1 A charged capacitor and an inductor are connected in series At time t = 0 the current is zero, but the capacitor is charged If T is the
More informationOscillations and Electromagnetic Waves. March 30, 2014 Chapter 31 1
Oscillations and Electromagnetic Waves March 30, 2014 Chapter 31 1 Three Polarizers! Consider the case of unpolarized light with intensity I 0 incident on three polarizers! The first polarizer has a polarizing
More informationEECE 301 Signals & Systems Prof. Mark Fowler
EECE 301 Signals & Systems Prof. Mark Fowler Note Set #15 C-T Systems: CT Filters & Frequency Response 1/14 Ideal Filters Often we have a scenario where part of the input signal s spectrum comprises what
More informationI. Impedance of an R-L circuit.
I. Impedance of an R-L circuit. [For inductor in an AC Circuit, see Chapter 31, pg. 1024] Consider the R-L circuit shown in Figure: 1. A current i(t) = I cos(ωt) is driven across the circuit using an AC
More informationAnalysis III for D-BAUG, Fall 2017 Lecture 10
Analysis III for D-BAUG, Fall 27 Lecture Lecturer: Alex Sisto (sisto@math.ethz.ch Convolution (Faltung We have already seen that the Laplace transform is not multiplicative, that is, L {f(tg(t} L {f(t}
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More information16.362: Signals and Systems: 1.0
16.362: Signals and Systems: 1.0 Prof. K. Chandra ECE, UMASS Lowell September 1, 2016 1 Background The pre-requisites for this course are Calculus II and Differential Equations. A basic understanding of
More informationEE102 Homework 2, 3, and 4 Solutions
EE12 Prof. S. Boyd EE12 Homework 2, 3, and 4 Solutions 7. Some convolution systems. Consider a convolution system, y(t) = + u(t τ)h(τ) dτ, where h is a function called the kernel or impulse response of
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems
CDS 101 1. Åström and Murray, Exercise 1.3 2. Åström and Murray, Exercise 1.4 3. Åström and Murray, Exercise 2.6, parts (a) and (b) CDS 110a 1. Åström and Murray, Exercise 1.4 2. Åström and Murray, Exercise
More informationFrequency Response of Linear Time Invariant Systems
ME 328, Spring 203, Prof. Rajamani, University of Minnesota Frequency Response of Linear Time Invariant Systems Complex Numbers: Recall that every complex number has a magnitude and a phase. Example: z
More informationEXAM 3: SOLUTIONS. B = B. A 2 = BA 2 cos 0 o = BA 2. =Φ(2) B A 2 = A 1 cos 60 o = A 1 2 =0.5m2
EXAM : S Q.The normal to a certain m area makes an angle of 6 o with a uniform magnetic field. The magnetic flux through this area is the same as the flux through a second area that is perpendicular to
More informationAssignment 8. [η j, η k ] = J jk
Assignment 8 Goldstein 9.8 Prove directly that the transformation is canonical and find a generating function. Q 1 = q 1, P 1 = p 1 p Q = p, P = q 1 q We can establish that the transformation is canonical
More informationMET 487 Instrumentation and Automatic Control. Lecture 3
MET 487 Instrumentation and Automatic Control Lecture 3 Electrical Circuits and Components http://www.etcs.ipfw.edu/~lin Lecture 2 - By P. Lin 1 Electrical Circuits and Components Basic Electrical Elements
More informationMassachusetts Institute of Technology Physics 8.03 Fall 2004 Final Exam Thursday, December 16, 2004
You have 3 hours Do all eight problems You may use calculators Massachusetts Institute of Technology Physics 8.03 Fall 004 Final Exam Thursday, December 16, 004 This is a closed-book exam; no notes are
More informationTopic 1: Simple harmonic motion
Topic 1: Simple harmonic motion Introduction Why do we need to know about waves 1. Ubiquitous in science nature likes wave solutions to equations 2. They are an exemplar for some essential Physics skills:
More informationMATH 251 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam
MATH 51 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam A collection of previous exams could be found at the coordinator s web: http://www.math.psu.edu/tseng/class/m51samples.html
More information221B Lecture Notes on Resonances in Classical Mechanics
1B Lecture Notes on Resonances in Classical Mechanics 1 Harmonic Oscillators Harmonic oscillators appear in many different contexts in classical mechanics. Examples include: spring, pendulum (with a small
More informationLecture 2: Controllability of nonlinear systems
DISC Systems and Control Theory of Nonlinear Systems 1 Lecture 2: Controllability of nonlinear systems Nonlinear Dynamical Control Systems, Chapter 3 See www.math.rug.nl/ arjan (under teaching) for info
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More informationChap 4. State-Space Solutions and
Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations
More informationPhysics 4 Spring 1989 Lab 5 - AC Circuits
Physics 4 Spring 1989 Lab 5 - AC Circuits Theory Consider the series inductor-resistor-capacitor circuit shown in figure 1. When an alternating voltage is applied to this circuit, the current and voltage
More informationVibrations and waves: revision. Martin Dove Queen Mary University of London
Vibrations and waves: revision Martin Dove Queen Mary University of London Form of the examination Part A = 50%, 10 short questions, no options Part B = 50%, Answer questions from a choice of 4 Total exam
More information