GROUP HOMEWORK 5 AERE573 Fall 2018 Due 11/16(F) NAMES The homework addresses the estimation of the 1-D position and velocity of a flexible robot arm
|
|
- Cory Parsons
- 5 years ago
- Views:
Transcription
1 1 GROUP HOMEWORK 5 AERE573 Fall 18 Due 11/16(F NAMES The homewor aresses he esimaion o he 1-D posiion an velociy o a lexible robo arm m or which a moel is shown a righ The lui has a mean viscosiy ha is inclue as a par o ( he viscous amping coeicien c However, he lui also has urbulence ha resuls in zeromean lucuaing loaing o he mass The ynamical moel or Figure 1 was evelope in he lecure noes enile Kalman Filering wih Known Forcing Funcion This moel is: r c ( where I mr n n n n r I ( / (, n / I, an c / I ( Figure 1 Arm moel The ynamical moel or Figure 1 was evelope in he lecure noes enile Kalman Filering wih Known Forcing Funcion The associae anser uncions are: ( s ns n ( s ( r / I H1( s an H( s ( s s s F s s s n n ( n n Numerical Inormaion: 5 r/ s, 5, n r 1 m, an I 1 g m A Moel or (: Clearly he mass will experience a mean rearing momen when in a uniorm low iel I will also experience perurbaions ue o he unseay naure o he iel In he presence o low Reynols number (RN (smooh low, he orce on he mass can be assume o have low specal banwih As he RN increases, his banwih will increase Furhermore, as he op plo a righ shows, he rag orce on he mass will begin o experience low ampliue oscillaions A high RN he small perioic oscillaions will evolve ino a sizeable sochasic sucure, as shown in he lower plo a righ We can moel he rag orce in each o hese Reynols number regimes as ollow: R1: (Low RN ( ~ C Gauss Marov(, where STR BW STR is he vorex sheing requency assoociae wih he Souhal number, which or a sphere is STR STRD / ( U [hps://enwiipeiaorg/wii/souhal_number ] Here, D is he sphere iameer, an U is he mean low velociy n R: (Me RN ( ~ C Orer Gauss Marov( res STR R3: (High RN ( ~ C Gauss Marov( BW STR Figure Drag coeicien C ( [hp://wwwicmil/ic//ullex/u/a494935p ] c c Hence, in regions (R1 an (R3 we have Fs ( In (R we have Fs ( wih s s s BW res res NOTE: These anser uncions o no inclue C This is a saic rag ha will be balance ou by a saic orque associae wih Hence, we are conucing a perurbaion analysis o he sysem in Figure 1 Numerical Inormaion: In his homewor we will ocus on region (R, wih 11 an 1 res n
2 (a(5ps Use he c comman in relaion o ( H ( s 1 o arrive a relaing ( an Use Nyquis requency N ra / s Soluion: [See (a] H ( z 1 From his, recover he ierence equaion (b(5ps I he riving whie noise ha resuls in requencies he magniue o Soluion: ( ( has variance Show ha his requires c ( u Inres r Suppose ha we require ha a very low o o (c(5ps Se 1 ( /18 17 ra Obain he anser uncions H ( z ( z / F( z an H ( z F( z / U( z 3 Soluion: [See (c] ((5ps Deine he sae x 1( 1( 1 ( ( 1 ( ( 1 Use he numerical inormaion in ( o arrive a he explici escripion o x Ax Bu Then give he C an D maices associae wih ( Cx Du Soluion: [See (] u u u ( 1 ( ( 1 ( (e(ps Simulae he sae an measuremen process or a oal o 55 ime sep Then iscar he irs 5 poin Give plos o an Soluion: [See (] ( ( Figure (e Plos o ( (le an ( (righ Commens: The mos challenging par o Kalman ilering is o arrive a a real-worl moel or he processes o be involve In his lecure we have move owar such a moel By varying parameers we can simulae all manner o urbulence, an all manner o weighs, siness, amping, ineria, ec in relaion o he penulum Having such a moel is
3 a naural guie o eveloping he sae an measuremen moels neee in he KF We will now procee o emonsae his DEVELOPMENT OF THE KALMAN FILTER EQUATIONS In he above moel: x 1( 1( 1 ( ( 1 ( ( 1 u ( 1 ( u( 1 u( an 3 ((5ps Explain why he sae inpu Explanaion: u is no in a orm ha is suiable o be he KF sae inpu (g(5ps In relaion o he KF sae equaion x 1 Fx Gu, le x ( ( 1 ( ( 1 an le u u( 1 u( Give he numerical values or F an G Answer: (h(5ps Noice ha he process Gu is never acually compue All ha is require by he KF is Q Cov( Gu For Cov( u I4 4, show ha Q GG / T Then compue is numerical value Soluion: (i(5ps The KF sae iniial coniion will be se o x coniion Then procee o arrive a he expression or Explanaion: P Explain why his is he mos logical iniial in erms o cross- an auo-correlaions Soluion: (j(5ps To compue he value or P irecly is a nonivial as [These are correlae ARMA(, processes] As an alernaive, recall ha Rx( E( x 1x1 an ha x 1 Fx Gu Subsiue his laer expression ino each elemen o he ormer, in orer o arrive a he relaion R ( FR ( F Q [Hin: E( xu E( ux, an ae avanage o (h] Soluion: x x
4 4 ((5ps Use he Malab coe lyap o obain he numerical value or Soluion: See (] R x ( (l(ps A KF coe or wss processes is inclue in he Appenix Use i o arrive a wo plos; one wih overlai an ; he oher wih overlai ( an ( ( ( [See (l] Figure ( KF esimaes overlai agains ( (le an ( (righ (m(5ps You shoul have oun ha he plos o ( an ( are exacly he same Explain why his is o be expece Explanaion: (n(5ps You shoul have also oun ha, wih he excepion o he irs ew poins, he plos o ( an ( are also essenially he same Explain why his is also o be expece To his en, begin by compuing he noncausal Wiener iler H ( ha gives he bes ( H( ( Soluion:
5 CONCLUSION Given all o he above, i is reasonable o as why a KF approach is even use or such a simple problem Aer all, even i he eerminisic inpu is inclue, since we now i, as well as he anser uncion ha relaes i o, an since ( ( (, we can simply subac rom he measure in orer o arrive a a ( perec esimae o ( ( ( All o his is ue The value o he KF meho is ha i allows us o move closer o a real-worl seing where, or example, we have (i aiive measuremen noise, (ii inexac esimaes o he ue anser uncions, (iii sensors which, hemselves have anser uncions, an (iii slow-ime-varying anser uncions 5
6 A Noe on Compuing Auo- an Cross-Power Specal Densiies As simple as he Malab lyap soluion is, i lacs insigh ino he naure o he elemens o P R x ( From (1 i is clear ha hose elemens are auo- an cross-correlaions a lags zero an one So, one migh as: How can we arrive a hese correlaions via ps inormaion We now presen a meho ha recovers all auo/cross-correlaions base on pss 1 Recall: Y ( z H( z F( z an c1 1 c u gives F( z U( z H ( ( 1 u z U z 1c z c z 1 Hence, Y ( z H( z H ( z U( z H( z U( z From he Wiener-Kinchine Theorem we hereore have: ( u ( ( u( u S z E H z U z H z (i: (ii: u y ( ( ( ( ( ( ( u( u S z E Y z F z E H z F z F z H z H z y ( ( ( ( ( ( ( ( u( u S z E Y z Y z E H z F z H z F z H z H z (iii: For numerical values, we can compue he quaniies an hen recover he associae correlaion uncions Noe ha his is no necessary in relaion o (i, as we have a irec meho o solving or he auocorrelaion or an AR( process NOTE: To recover he coeicien o a anser uncion H, ype: [num,en]=aa(h, v 6 The coe below recovers correlaion uncions rom pss speciie via anser uncions %PROGRAM NAME: psxcorrm %This coe recovers Rxy rom Sxy in iscree-ime omain %============================================ %Examples: m=1; %TF1: H1=([1 ],[1,-9],1; %Transer Funcion [H1n,H1]=aa(H,'v'; H1nw=(H1n,m; H1w=(H1,m; H1w=H1nw/H1w; %TF: %H=H1; %For Auocorrelaion H=([1 ],[1,-4],1; %For Crosscorrelaion [Hn,H]=aa(H,'v'; Hnw=(Hn,m; Hw=(H,m; Hw=Hnw/Hw; %Cross-Specum % Cross-Correlaion S1=H1w*conj(Hw; R1=real(i(S1; igure(1 plo(r1,'*' ile('crosscorrelaion or Two Speciie TFs' gri Commen: I chece o mae sure ha he imaginary pars o he i s were zero They were
7 7 Appenix Malab Coe %PROGRAM NAME: hw5m %(a: I=1; r=1; wn=; %(b: %(c: %(: %(e: no=55; n=5; igure(1 plo(,h ile('plo o hea(' xlabel('time (sec' ylabel('degrees' gri igure( plo(, ile('plo o (' xlabel('time (sec' ylabel('newons' gri %======================================== % **** KALMAN FILTERING SECTION **** %======================================== %(g: %Consucion o F an H Maices: %(h: %Consucion o Q an R Maices: %Measuremen Array: z=h; %(: %Iniial Coniions: xha_ol=zeros(4,1; P_ol= % %(l: % PROGRAM NAME: wssm % This program compues he Kalman iler sae esimae % or wss sae an measuremen processes % MODEL: % x( = F x(-1 + w(-1 % z( = H x( + v( % REQUIRED INPUT: % F=nxn sae ansiion maix % Q=nxn cov maix or nx1 sae whie noise
8 % H=mxn sae/msmn maix % R=mxm cov maix or mx1 msmn noise % xha_ol(1=iniial sae esimae % P_ol(1=iniial err_ol nxn cov maix % max = max ime over which o compue xha( % {z(:=1,max}=m x max msmn vecor % ********************************* nx=lengh(f; I=eye(nx; nps=lengh(z; K=[]; xha=[]; max=nps; or =1:max K= P_ol*H'*(H*P_ol*H' + R^(-1; K=[K,K]; zha_ol = H*xha_ol; xha=xha_ol + K*(z(-zha_ol; xha=[xha,xha]; zha( = H*xha; P=(I-K*H*P_ol; xha_ol=f*xha; P_ol=F*P*F' + Q; en igure( vec=1:nps; vec=t*vec; plo(vec,z,vec,zha,'r' legen('h','hha' gri ha = xha(3,:; igure(3 plo(vec,,vec,ha,'r' legen('','ha' gri 8
Homework 6 AERE331 Spring 2019 Due 4/24(W) Name Sec. 1 / 2
Homework 6 AERE33 Spring 9 Due 4/4(W) Name Sec / PROBLEM (5p In PROBLEM 4 of HW4 we used he frequency domain o design a yaw/rudder feedback conrol sysem for a plan wih ransfer funcion 46 Gp () s The conroller
More informationCharacteristics of Linear System
Characerisics o Linear Sysem h g h : Impulse response F G : Frequency ranser uncion Represenaion o Sysem in ime an requency. Low-pass iler g h G F he requency ranser uncion is he Fourier ransorm o he impulse
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationwhere the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).
Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and
More informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More informationTheory of! Partial Differential Equations-I!
hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
More informationI. OBJECTIVE OF THE EXPERIMENT.
I. OBJECTIVE OF THE EXPERIMENT. Swissmero raels a high speeds hrough a unnel a low pressure. I will hereore undergo ricion ha can be due o: ) Viscosiy o gas (c. "Viscosiy o gas" eperimen) ) The air in
More informationTheory of! Partial Differential Equations!
hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More informationConvolution. Lecture #6 2CT.3 8. BME 333 Biomedical Signals and Systems - J.Schesser
Convoluion Lecure #6 C.3 8 Deiniion When we compue he ollowing inegral or τ and τ we say ha he we are convoluing wih g d his says: ae τ, lip i convolve in ime -τ, hen displace i in ime by seconds -τ, and
More informationChapter Three Systems of Linear Differential Equations
Chaper Three Sysems of Linear Differenial Equaions In his chaper we are going o consier sysems of firs orer orinary ifferenial equaions. These are sysems of he form x a x a x a n x n x a x a x a n x n
More informationA dynamic AS-AD Model
A ynamic AS-AD Moel (Lecure Noes, Thomas Seger, Universiy of Leipzig, winer erm 10/11) This file escribes a ynamic AS-AD moel. The moel can be employe o assess he ynamic consequences of macroeconomic shocks
More informationF This leads to an unstable mode which is not observable at the output thus cannot be controlled by feeding back.
Lecure 8 Las ime: Semi-free configuraion design This is equivalen o: Noe ns, ener he sysem a he same place. is fixed. We design C (and perhaps B. We mus sabilize if i is given as unsable. Cs ( H( s = +
More informationZürich. ETH Master Course: L Autonomous Mobile Robots Localization II
Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),
More informationHomework 2 Solutions
Mah 308 Differenial Equaions Fall 2002 & 2. See he las page. Hoework 2 Soluions 3a). Newon s secon law of oion says ha a = F, an we know a =, so we have = F. One par of he force is graviy, g. However,
More informationME 391 Mechanical Engineering Analysis
Fall 04 ME 39 Mechanical Engineering Analsis Eam # Soluions Direcions: Open noes (including course web posings). No books, compuers, or phones. An calculaor is fair game. Problem Deermine he posiion of
More informationLecture 16 (Momentum and Impulse, Collisions and Conservation of Momentum) Physics Spring 2017 Douglas Fields
Lecure 16 (Momenum and Impulse, Collisions and Conservaion o Momenum) Physics 160-02 Spring 2017 Douglas Fields Newon s Laws & Energy The work-energy heorem is relaed o Newon s 2 nd Law W KE 1 2 1 2 F
More informationL07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms
L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationdt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.
Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies
More informationY 0.4Y 0.45Y Y to a proper ARMA specification.
HG Jan 04 ECON 50 Exercises II - 0 Feb 04 (wih answers Exercise. Read secion 8 in lecure noes 3 (LN3 on he common facor problem in ARMA-processes. Consider he following process Y 0.4Y 0.45Y 0.5 ( where
More informationLecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
More informationTopic Astable Circuits. Recall that an astable circuit has two unstable states;
Topic 2.2. Asable Circuis. Learning Objecives: A he end o his opic you will be able o; Recall ha an asable circui has wo unsable saes; Explain he operaion o a circui based on a Schmi inverer, and esimae
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationExam #3 Review. Skills #6 21.3, 24.2, 24.6, 25.1
Exam #3 Review Wha is on an exam? Read hrough he homework and class examples 4 mulipar quesions. Poins assigned based on complexiy. (4Q, 40-60 ps. Sp. 06/Fa. 06) Skills #5.3, 4., 4.6, 5. Skills #6.3, 4.,
More informationInstitute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler
MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion
More informationPhysics Notes - Ch. 2 Motion in One Dimension
Physics Noes - Ch. Moion in One Dimension I. The naure o physical quaniies: scalars and ecors A. Scalar quaniy ha describes only magniude (how much), NOT including direcion; e. mass, emperaure, ime, olume,
More informationAugmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004
Augmened Realiy II Kalman Filers Gudrun Klinker May 25, 2004 Ouline Moivaion Discree Kalman Filer Modeled Process Compuing Model Parameers Algorihm Exended Kalman Filer Kalman Filer for Sensor Fusion Lieraure
More informationState-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter
Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when
More informationLecture 10: Wave equation, solution by spherical means
Lecure : Wave equaion, soluion by spherical means Physical modeling eample: Elasodynamics u (; ) displacemen vecor in elasic body occupying a domain U R n, U, The posiion of he maerial poin siing a U in
More informationLecture #6: Continuous-Time Signals
EEL5: Discree-Time Signals and Sysems Lecure #6: Coninuous-Time Signals Lecure #6: Coninuous-Time Signals. Inroducion In his lecure, we discussed he ollowing opics:. Mahemaical represenaion and ransormaions
More information6.003 Homework #8 Solutions
6.003 Homework #8 Soluions Problems. Fourier Series Deermine he Fourier series coefficiens a k for x () shown below. x ()= x ( + 0) 0 a 0 = 0 a k = e /0 sin(/0) for k 0 a k = π x()e k d = 0 0 π e 0 k d
More informationLecture #7. EECS490: Digital Image Processing. Image Processing Example Fuzzy logic. Fourier Transform. Basics Image processing examples
Lecure #7 Image Processing Example Fuzzy logic Basics Image processing examples Fourier Transorm Inner produc, basis uncions Fourier series Image Processing Example original image Laplacian o image (c)
More informationThe Quantum Theory of Atoms and Molecules: The Schrodinger equation. Hilary Term 2008 Dr Grant Ritchie
e Quanum eory of Aoms and Molecules: e Scrodinger equaion Hilary erm 008 Dr Gran Ricie An equaion for maer waves? De Broglie posulaed a every paricles as an associaed wave of waveleng: / p Wave naure of
More informationModule 4: Time Response of discrete time systems Lecture Note 2
Module 4: Time Response of discree ime sysems Lecure Noe 2 1 Prooype second order sysem The sudy of a second order sysem is imporan because many higher order sysem can be approimaed by a second order model
More informationR t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t
Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,
More informationSUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL
HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN
More informationExcel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand
Excel-Based Soluion Mehod For The Opimal Policy Of The Hadley And Whiin s Exac Model Wih Arma Demand Kal Nami School of Business and Economics Winson Salem Sae Universiy Winson Salem, NC 27110 Phone: (336)750-2338
More informationAnnouncements: Warm-up Exercise:
Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple
More informationOBJECTIVES OF TIME SERIES ANALYSIS
OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More information2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006
2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)
More informationStochastic Structural Dynamics. Lecture-6
Sochasic Srucural Dynamics Lecure-6 Random processes- Dr C S Manohar Deparmen of Civil Engineering Professor of Srucural Engineering Indian Insiue of Science Bangalore 560 0 India manohar@civil.iisc.erne.in
More informationVanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law
Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing
More information( ) = Q 0. ( ) R = R dq. ( t) = I t
ircuis onceps The addiion of a simple capacior o a circui of resisors allows wo relaed phenomena o occur The observaion ha he ime-dependence of a complex waveform is alered by he circui is referred o as
More informationUNIVERSITY OF CALIFORNIA AT BERKELEY
Homework #10 Soluions EECS 40, Fall 2006 Prof. Chang-Hasnain Due a 6 pm in 240 Cory on Wednesday, 04/18/07 oal Poins: 100 Pu (1) your name and (2) discussion secion number on your homework. You need o
More informationLinear Gaussian State Space Models
Linear Gaussian Sae Space Models Srucural Time Series Models Level and Trend Models Basic Srucural Model (BSM Dynamic Linear Models Sae Space Model Represenaion Level, Trend, and Seasonal Models Time Varying
More informationProblem Set #1. i z. the complex propagation constant. For the characteristic impedance:
Problem Se # Problem : a) Using phasor noaion, calculae he volage and curren waves on a ransmission line by solving he wave equaion Assume ha R, L,, G are all non-zero and independen of frequency From
More informationHamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t
M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n
More informationMatlab and Python programming: how to get started
Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,
More informationHamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:
M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno
More informationSPH3U: Projectiles. Recorder: Manager: Speaker:
SPH3U: Projeciles Now i s ime o use our new skills o analyze he moion of a golf ball ha was ossed hrough he air. Le s find ou wha is special abou he moion of a projecile. Recorder: Manager: Speaker: 0
More informationSensors, Signals and Noise
Sensors, Signals and Noise COURSE OUTLINE Inroducion Signals and Noise: 1) Descripion Filering Sensors and associaed elecronics rv 2017/02/08 1 Noise Descripion Noise Waveforms and Samples Saisics of Noise
More informationAn recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes
WHAT IS A KALMAN FILTER An recursive analyical echnique o esimae ime dependen physical parameers in he presence of noise processes Example of a ime and frequency applicaion: Offse beween wo clocks PREDICTORS,
More informationData Fusion using Kalman Filter. Ioannis Rekleitis
Daa Fusion using Kalman Filer Ioannis Rekleiis Eample of a arameerized Baesian Filer: Kalman Filer Kalman filers (KF represen poserior belief b a Gaussian (normal disribuion A -d Gaussian disribuion is
More informationTwo-stage Benchmarking of Time-Series Models for. Small Area Estimation. Danny Pfeffermann, Richard Tiller
Two-sage Benchmarking of Time-Series Moels for Small Area Esimaion anny Pfeffermann, Souhampon Universiy, UK & Hebrew universiy, Israel Richar Tiller Bureau of Labor Saisics, U.S.A. Small Area Conference,
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationA Note on Fractional Electrodynamics. Abstract
Commun Nonlinear Sci Numer Simula 8 (3 589 593 A Noe on Fracional lecrodynamics Hosein Nasrolahpour Absrac We invesigae he ime evoluion o he racional elecromagneic waves by using he ime racional Maxwell's
More informationDifferential Equations
Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding
More information72 Calculus and Structures
72 Calculus and Srucures CHAPTER 5 DISTANCE AND ACCUMULATED CHANGE Calculus and Srucures 73 Copyrigh Chaper 5 DISTANCE AND ACCUMULATED CHANGE 5. DISTANCE a. Consan velociy Le s ake anoher look a Mary s
More informationLinear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation
Chaper 4 Linear Surface Graviy Waves 3., Dispersion, Group Velociy, and Energy Propagaion 4. Descripion In many aspecs of wave evoluion, he concep of group velociy plays a cenral role. Mos people now i
More informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More informationModeling Economic Time Series with Stochastic Linear Difference Equations
A. Thiemer, SLDG.mcd, 6..6 FH-Kiel Universiy of Applied Sciences Prof. Dr. Andreas Thiemer e-mail: andreas.hiemer@fh-kiel.de Modeling Economic Time Series wih Sochasic Linear Difference Equaions Summary:
More information4.5 Constant Acceleration
4.5 Consan Acceleraion v() v() = v 0 + a a() a a() = a v 0 Area = a (a) (b) Figure 4.8 Consan acceleraion: (a) velociy, (b) acceleraion When he x -componen of he velociy is a linear funcion (Figure 4.8(a)),
More informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
More informationRecursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems
8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationCOMPETITIVE GROWTH MODEL
COMPETITIVE GROWTH MODEL I Assumpions We are going o now solve he compeiive version of he opimal growh moel. Alhough he allocaions are he same as in he social planning problem, i will be useful o compare
More informationCointegration in Frequency Domain*
Coinegraion in Frequenc Domain* Daniel Lev Deparmen o Economics Bar-Ilan Universi Rama-Gan 59 ISRAEL Tel: 97-3-53-833 Fax: 97-3-535-38 LEVDA@MAIL.BIU.AC.IL and Deparmen o Economics Emor Universi Alana,
More informationHint: There's a table of particular solutions at the end of today's notes.
Mah 8- Fri Apr 4, Finish Wednesday's noes firs Then 94 Forced oscillaion problems via Fourier Series Today we will revisi he forced oscillaion problems of las Friday, where we prediced wheher or no resonance
More informationLecture 10 Estimating Nonlinear Regression Models
Lecure 0 Esimaing Nonlinear Regression Models References: Greene, Economeric Analysis, Chaper 0 Consider he following regression model: y = f(x, β) + ε =,, x is kx for each, β is an rxconsan vecor, ε is
More informationTesting the Random Walk Model. i.i.d. ( ) r
he random walk heory saes: esing he Random Walk Model µ ε () np = + np + Momen Condiions where where ε ~ i.i.d he idea here is o es direcly he resricions imposed by momen condiions. lnp lnp µ ( lnp lnp
More informationElectrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit
V() R L C 513 Elecrical Circuis Tools Used in Lab 13 Series Circuis Damped Vibraions: Energy Van der Pol Circui A series circui wih an inducor, resisor, and capacior can be represened by Lq + Rq + 1, a
More informationPosition, Velocity, and Acceleration
rev 06/2017 Posiion, Velociy, and Acceleraion Equipmen Qy Equipmen Par Number 1 Dynamic Track ME-9493 1 Car ME-9454 1 Fan Accessory ME-9491 1 Moion Sensor II CI-6742A 1 Track Barrier Purpose The purpose
More informationThe Brock-Mirman Stochastic Growth Model
c December 3, 208, Chrisopher D. Carroll BrockMirman The Brock-Mirman Sochasic Growh Model Brock and Mirman (972) provided he firs opimizing growh model wih unpredicable (sochasic) shocks. The social planner
More informationKinematics and kinematic functions
Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables and vice versa Direc Posiion
More informationHW6: MRI Imaging Pulse Sequences (7 Problems for 100 pts)
HW6: MRI Imaging Pulse Sequences (7 Problems for 100 ps) GOAL The overall goal of HW6 is o beer undersand pulse sequences for MRI image reconsrucion. OBJECTIVES 1) Design a spin echo pulse sequence o image
More informationProgram: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro. Category: Isotropic Linear Elasticity, Dynamics, Member
Verificaion Example Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro Caegory: Isoropic Linear Elasiciy, Dynamics, Member Verificaion Example: 0104 Canilever Beam wih Periodic Exciaion 0104 Canilever Beam
More informationAnno accademico 2006/2007. Davide Migliore
Roboica Anno accademico 2006/2007 Davide Migliore migliore@ele.polimi.i Today Eercise session: An Off-side roblem Robo Vision Task Measuring NBA layers erformance robabilisic Roboics Inroducion The Bayesian
More informationModule II, Part C. More Insight into Fiber Dispersion
Moule II Par C More Insigh ino Fiber Dispersion . Polariaion Moe Dispersion Fiber Birefringence: Imperfec cylinrical symmery leas o wha is known as birefringence. Recall he HE moe an is E x componen which
More informationA New Macro-Model for Power Diodes Reverse Recovery
Proceeings o he 7h WSEAS Inernaional Conerence on Power Sysems, Beijing, China, Sepember 15-17, 7 48 A New Macro-Moel or Power Dioes Reverse Recovery ALI DASTFAN Depaen o Elecrical & Roboic Engineering
More informationForecasting optimally
I) ile: Forecas Evaluaion II) Conens: Evaluaing forecass, properies of opimal forecass, esing properies of opimal forecass, saisical comparison of forecas accuracy III) Documenaion: - Diebold, Francis
More information12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =
1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of
More informationTemporal probability models
Temporal probabiliy models CS194-10 Fall 2011 Lecure 25 CS194-10 Fall 2011 Lecure 25 1 Ouline Hidden variables Inerence: ilering, predicion, smoohing Hidden Markov models Kalman ilers (a brie menion) Dynamic
More informationSequential Importance Resampling (SIR) Particle Filter
Paricle Filers++ Pieer Abbeel UC Berkeley EECS Many slides adaped from Thrun, Burgard and Fox, Probabilisic Roboics 1. Algorihm paricle_filer( S -1, u, z ): 2. Sequenial Imporance Resampling (SIR) Paricle
More informationNEWTON S SECOND LAW OF MOTION
Course and Secion Dae Names NEWTON S SECOND LAW OF MOTION The acceleraion of an objec is defined as he rae of change of elociy. If he elociy changes by an amoun in a ime, hen he aerage acceleraion during
More informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More informationEE363 homework 1 solutions
EE363 Prof. S. Boyd EE363 homework 1 soluions 1. LQR for a riple accumulaor. We consider he sysem x +1 = Ax + Bu, y = Cx, wih 1 1 A = 1 1, B =, C = [ 1 ]. 1 1 This sysem has ransfer funcion H(z) = (z 1)
More informationDesign of Controller for Robot Position Control
eign of Conroller for Robo oiion Conrol Two imporan goal of conrol: 1. Reference inpu racking: The oupu mu follow he reference inpu rajecory a quickly a poible. Se-poin racking: Tracking when he reference
More informationWeek 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)
Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =
More informationTutorial Sheet #2 discrete vs. continuous functions, periodicity, sampling
2.39 Tuorial Shee #2 discree vs. coninuous uncions, periodiciy, sampling We will encouner wo classes o signals in his class, coninuous-signals and discree-signals. The disinc mahemaical properies o each,
More informationDynamic Effects of Feedback Control!
Dynamic Effecs of Feedback Conrol! Rober Sengel! Roboics and Inelligen Sysems MAE 345, Princeon Universiy, 2017 Inner, Middle, and Ouer Feedback Conrol Loops Sep Response of Linear, Time- Invarian (LTI)
More informationContinuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.
Time-Domain Sysem Analysis Coninuous Time. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 1. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 2 Le a sysem be described by a 2 y ( ) + a 1
More informationLectures 29 and 30 BIQUADRATICS AND STATE SPACE OP AMP REALIZATIONS. I. Introduction
EE-202/445, 3/18/18 9-1 R. A. DeCarlo Lecures 29 and 30 BIQUADRATICS AND STATE SPACE OP AMP REALIZATIONS I. Inroducion 1. The biquadraic ransfer funcion has boh a 2nd order numeraor and a 2nd order denominaor:
More informationMath 527 Lecture 6: Hamilton-Jacobi Equation: Explicit Formulas
Mah 527 Lecure 6: Hamilon-Jacobi Equaion: Explici Formulas Sep. 23, 2 Mehod of characerisics. We r o appl he mehod of characerisics o he Hamilon-Jacobi equaion: u +Hx, Du = in R n, u = g on R n =. 2 To
More informationLecture #31, 32: The Ornstein-Uhlenbeck Process as a Model of Volatility
Saisics 441 (Fall 214) November 19, 21, 214 Prof Michael Kozdron Lecure #31, 32: The Ornsein-Uhlenbeck Process as a Model of Volailiy The Ornsein-Uhlenbeck process is a di usion process ha was inroduced
More informationPhysics 1402: Lecture 22 Today s Agenda
Physics 142: ecure 22 Today s Agenda Announcemens: R - RV - R circuis Homework 6: due nex Wednesday Inducion / A curren Inducion Self-Inducance, R ircuis X X X X X X X X X long solenoid Energy and energy
More informationAssignment 6. Tyler Shendruk December 6, 2010
Assignmen 6 Tyler Shendruk December 6, 1 1 Harden Problem 1 Le K be he coupling and h he exernal field in a 1D Ising model. From he lecures hese can be ransformed ino effecive coupling and fields K and
More information