Estimation of complex motions in video signals. Erhardt Barth. Institute for Neuro- and Bioinformatics University of Lübeck, Germany

Size: px
Start display at page:

Download "Estimation of complex motions in video signals. Erhardt Barth. Institute for Neuro- and Bioinformatics University of Lübeck, Germany"

Transcription

1 Esimaion o complex moions in video signals Erhard Barh Insiue or Neuro- and Bioinormaics Universiy o Lübec, Germany Nonlinear analysis o mulidimensional signals: LOCOMOTOR LOcal adapive esimaion o COmplex MOTion and ORienaion paerns LOCOMOTOR Franur Heidelberg Im DFG Schwerpun 4: Mahemaical mehods or ime series analysis and digial image processing hp:// Projec Ba-76/7 Lübec

2 Inrinsic dimension in 3D i0d: consan in all direcions: id: consan in direcions: id: consan in one direcion: i3d: no consan direcion: (x, y,) = cons. (x, y,) = g(ξ) (x, y,) = g(ξ,ζ) ( x,y,) = g(ξ,ζ,τ) Inrinsic dimension id and i3d regions are he leas requen bu mos signiican. i0d id id i3d

3 The classical moion model (x,y,) image-sequence inensiy v = (v,v ) T moion vecor α (v) = v + v + = V = 0 x y r = ( v,v, ) Why complex moions? Basic problems wih naural image sequences: muliple moions, occlusions, noisy daa, eaures a dieren scales,... 3

4 Transparen moions Movie (op) and x plo (below). Moion model or ransparen moions Opical low or one moion (x,y,) image sequence, v = (v,v ) T moion vecor α (v) = v + v + = V = 0 x y wih he derivaive operaor: α (v) = v + v + x y Opical low or N moions: = g (x-v )+...+g N (x-v N ) α(v ) α(v ) α (v ) = 0 N 4

5 Mixed-moion parameers Example or wo moions u, v : (x,)=g (x-v )+g (x-u ) α(v) α(u) = u v + ( u v xx + u v + u v ) + ( u + v ) x yy xy + ( u + v ) We deine he mixed-moion parameers as: y + c = vu c = vu c = uv +uv xx yy xy c = u+v c = u+v c = x y Solving or he mixed-moion parameers Wih he mixed-moion parameers we obain: α(v) α (u) = c + c + c xx xx yy yy xy xy + c + c + c x x y y = V L= 0 For one moion we had V = 0 The above consrain can be used in a number o ways o derive he mixed moion parameers in V rom, e.g. by deining he generalized srucure ensor e.g. JV= N 0 and solving he sysem m i im i(m ) i J N = T h LL V (M, M,...,( ) M ) M ij, i =,...,m are he minors o J N 5

6 Separaion o he mixed-moion parameers We inerpre he moion vecors u and v as complex numbers (v = v +i v, u = u +i u ) and observe ha u v = A0 = cxx cyy + icxy u+ v= A = c + ic x y A 0 and A are homogeneous symmerical uncions o he coordinaes u and v and hereore by Viea s rule he coeicien o he complex polynomial Q(z) = (z u)(z v) = z Az + A0 ha has he complex roos u and v. In case o N moions Q(z) = z A z + + ( ) A N N N N 0 Fourier analysis o muliple ransparen moions Coninuous ime ( x, ) = g α( u) α( v) ( x u) + g = 0 F = G δ( u ω+ ω ) + G ( x v) δ( v ω+ ω ) Discree ime ( x) = g F ( x u) + g ( ω) = φ G ( ω) + φ ( x v) G ( ω) α (u ) = u + φ -ju ω = e 6

7 Bloc-maching or ransparen moions F ( ω) = φ G( ω) + φg ( ω) Layer eliminaion F a ( ) ( ) 0( ) ω a F ω + a F ω = 0 = φ +φ a = φ φ Bloc-maching consrain Bac o space domain x) ( x u) ( x v) + ( x u v) = 0 ( 0 Synheic resuls or bloc-maching SNR : 35 db Bloc size : 5 x 5 pixels 7

8 Fourier analysis o muliple occluded moions ( x, ) = χ( x u) g( x u) + χ( x u) g ( x v) g ( x) g ( x) occluding objec occluded objec χ characerisic uncion o he occluding objec χ = χ How does he Fourier ransorm o (x,) deviae rom he wo moion planes? α ( u ) α ( v ) ( x, )? Occluded moion α( u) α( v) ( x, ) = q( x, u, v, ) δ( B( x u)) q( x, u, v, ) = ( u v) N( x u) ( u v) g ( x v) Foreground moion Bacground moion Normal N o he border B(x)=0 The equaion is he same as or ransparen moions excep a he occluding boundary. q( x, u, v, ) is maximal i u - v is orhogonal o he boundary is null i u - v is angen o he boundary 8

9 Occluded moion in he Fourier domain F( ω, ω ) = A( ω) δ( u ω+ω ) + B( ω) δ( v ω+ω ) + C( ω, ω ) C( ω, ω ) d /( v ω+ω ) along he plane u ω+ω = c The disorion C has hyperbolic decay Log(+F) Log(+F) The hyperbolic decay o C has been irs recognized by Beauchemin a al. or he paricular case o sraigh line boundary Reerences C. Moa, I. Sue, T. Aach, and E. Barh: Spaio-emporal moion esimaion or ransparency and occlusions. Signal Processing: Image Communicaion. Elsevier Science, 005. C. Moa, I. Sue, T. Aach, and E. Barh: Divide-and-Conquer Sraegies or Esimaing Muliple Transparen Moions. In: Proceedings o he s Inernaional Worshop on Complex Moion, Schloss Reisensburg, Germany. Lecure Noes on Compuer Science, LNCS Vol. 347, 005. C. Moa, I. Sue, T. Aach, and E. Barh: Esimaion o muliple orienaions a corners and juncions. In: 6h Paern Recogniion Symposium (DAGM'04), Tübingen, 63-70, 004. C. Moa, T. Aach, I. Sue, and E. Barh: Esimaion o muliple orienaions in muli-dimensional signals. In: IEEE Inernaional Conerence on Image Processing, , 004. I. Sue, T. Aach, E. Barh, and C. Moa: Muliple-moion-esimaion by bloc-maching using MRF. Inernaional Journal o Compuer & Inormaion Science, Vol. 5, No., pp. 4-5, 004. I. Sue, T. Aach, E. Barh, and C. Moa: Muliple-moion esimaion by bloc-maching using Marov random ields. In: Inernaional Journal o Compuer & Inormaion Science, Vol. 5, No., pp. 4-5, June 004. E. Barh, I. Sue, T. Aach, and C. Moa: Spaio-Temporal Moion Esimaion or Transparency and Occlusions.In: Proceedings o IEEE Inernaional Conerence on Image Processing (ICIP 03), Barcelona, Spain, Sepember 4-7, 003, Vol. III, pp E. Barh, I. Sue, and C. Moa: Analysis o moion and curvaure in image sequences (invied paper). In: Proc. o he 5h IEEE Souhwes Symposium on Image Analysis and Inerpreaion, 06-0, 00. C. Moa, I. Sue, and E. Barh: Analyic Soluions or Muliple Moions. In: Proceedings o he IEEE Inernaional Conerence on Image Processing (ICIP 0), Thessalonii, Greece, Ocober 7-0, 00, Vol. II, pp

Theory of! Partial Differential Equations!

Theory 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 information

Theory of! Partial Differential Equations-I!

Theory 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 information

Lecture #6: Continuous-Time Signals

Lecture #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 information

Math 221: Mathematical Notation

Math 221: Mathematical Notation Mah 221: Mahemaical Noaion Purpose: One goal in any course is o properly use he language o ha subjec. These noaions summarize some o he major conceps and more diicul opics o he uni. Typing hem helps you

More information

Rotation, Scale and Translation Resilient Public Watermarking for Images

Rotation, Scale and Translation Resilient Public Watermarking for Images Roaion Scale and Translaion Resilien Public Waermarking or Images Ching-Yung Lin a Min Wu b Jerey A. Bloom c Ma L. Miller c Ingemar J. Cox c Yui Man Lui d a Columbia Universiy New York NY b Princeon Universiy

More information

Tracking. Announcements

Tracking. Announcements Tracking Tuesday, Nov 24 Krisen Grauman UT Ausin Announcemens Pse 5 ou onigh, due 12/4 Shorer assignmen Auo exension il 12/8 I will no hold office hours omorrow 5 6 pm due o Thanksgiving 1 Las ime: Moion

More information

Solutions for homework 12

Solutions for homework 12 y Soluions for homework Secion Nonlinear sysems: The linearizaion of a nonlinear sysem Consider he sysem y y y y y (i) Skech he nullclines Use a disincive marking for each nullcline so hey can be disinguished

More information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

3, so θ = arccos

3, so θ = arccos Mahemaics 210 Professor Alan H Sein Monday, Ocober 1, 2007 SOLUTIONS This problem se is worh 50 poins 1 Find he angle beween he vecors (2, 7, 3) and (5, 2, 4) Soluion: Le θ be he angle (2, 7, 3) (5, 2,

More information

Chapter 1 Fundamental Concepts

Chapter 1 Fundamental Concepts Chaper 1 Fundamenal Conceps 1 Signals A signal is a paern of variaion of a physical quaniy, ofen as a funcion of ime (bu also space, disance, posiion, ec). These quaniies are usually he independen variables

More information

Some Basic Information about M-S-D Systems

Some 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 information

EELE Lecture 8 Example of Fourier Series for a Triangle from the Fourier Transform. Homework password is: 14445

EELE Lecture 8 Example of Fourier Series for a Triangle from the Fourier Transform. Homework password is: 14445 EELE445-4 Lecure 8 Eample o Fourier Series or a riangle rom he Fourier ransorm Homework password is: 4445 3 4 EELE445-4 Lecure 8 LI Sysems and Filers 5 LI Sysem 6 3 Linear ime-invarian Sysem Deiniion o

More information

Scientific Research of the Institute of Mathematics and Computer Science DIFFERENT VARIANTS OF THE BOUNDARY ELEMENT METHOD FOR PARABOLIC EQUATIONS

Scientific Research of the Institute of Mathematics and Computer Science DIFFERENT VARIANTS OF THE BOUNDARY ELEMENT METHOD FOR PARABOLIC EQUATIONS Scieniic Research o he Insiue o Mahemaics and Compuer Science DIERENT VARIANTS O THE BOUNDARY ELEMENT METHOD OR PARABOLIC EQUATIONS Ewa Majchrzak,, Ewa Ładyga Jerzy Mendakiewicz, Alicja Piasecka Belkhaya

More information

CS376 Computer Vision Lecture 6: Optical Flow

CS376 Computer Vision Lecture 6: Optical Flow CS376 Compuer Vision Lecure 6: Opical Flow Qiing Huang Feb. 11 h 2019 Slides Credi: Krisen Grauman and Sebasian Thrun, Michael Black, Marc Pollefeys Opical Flow mage racking 3D compuaion mage sequence

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

Mean Square Projection Error Gradient-based Variable Forgetting Factor FAPI

Mean Square Projection Error Gradient-based Variable Forgetting Factor FAPI 3rd Inernaional Conference on Advances in Elecrical and Elecronics Engineering (ICAEE'4) Feb. -, 4 Singapore Mean Square Projecion Error Gradien-based Variable Forgeing Facor FAPI Young-Kwang Seo, Jong-Woo

More information

UNIVERSITY OF TRENTO MEASUREMENTS OF TRANSIENT PHENOMENA WITH DIGITAL OSCILLOSCOPES. Antonio Moschitta, Fabrizio Stefani, Dario Petri.

UNIVERSITY OF TRENTO MEASUREMENTS OF TRANSIENT PHENOMENA WITH DIGITAL OSCILLOSCOPES. Antonio Moschitta, Fabrizio Stefani, Dario Petri. UNIVERSIY OF RENO DEPARMEN OF INFORMAION AND COMMUNICAION ECHNOLOGY 385 Povo reno Ialy Via Sommarive 4 hp://www.di.unin.i MEASUREMENS OF RANSIEN PHENOMENA WIH DIGIAL OSCILLOSCOPES Anonio Moschia Fabrizio

More information

Announcements: Warm-up Exercise:

Announcements: 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 information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = 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 information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

Kinematics and kinematic functions

Kinematics 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 information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

More information

Lecture 9: Advanced DFT concepts: The Exchange-correlation functional and time-dependent DFT

Lecture 9: Advanced DFT concepts: The Exchange-correlation functional and time-dependent DFT Lecure 9: Advanced DFT conceps: The Exchange-correlaion funcional and ime-dependen DFT Marie Curie Tuorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dep. of Chemisry and Couran Insiue

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing. MECHANICS APPLICATIONS OF SECOND-ORDER ODES 7 Mechanics applicaions of second-order ODEs Second-order linear ODEs wih consan coefficiens arise in many physical applicaions. One physical sysems whose behaviour

More information

A Limit Symmetry of Modified KdV Equation and Its Applications

A Limit Symmetry of Modified KdV Equation and Its Applications Commun. Theor. Phys. 55 011 960 964 Vol. 55 No. 6 June 15 011 A Limi Symmery o Modiied KdV Equaion and Is Applicaions ZHANG Jian-Bing Ï 1 JI Jie SHEN Qing ã 3 and ZHANG Da-Jun 3 1 School o Mahemaical Sciences

More information

Tutorial Sheet #2 discrete vs. continuous functions, periodicity, sampling

Tutorial 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 information

SOUND AND VIBRATION SIGNAL ANALYSIS USING IMPROVED SHORT-TIME FOURIER REPRESENTATION. June-Yule Lee

SOUND AND VIBRATION SIGNAL ANALYSIS USING IMPROVED SHORT-TIME FOURIER REPRESENTATION. June-Yule Lee Inernaional Journal o Auomoive and Mechanical Engineering IJAME ISSN: 9-8649 Prin; ISSN: 80-606 Online; Volume 7 pp. 8-89 January-June 03 Universii Malaysia Pahang DOI: hp://d.doi.org/0.58/ijame.7.0..0066

More information

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves Rapid Terminaion Evaluaion for Recursive Subdivision of Bezier Curves Thomas F. Hain School of Compuer and Informaion Sciences, Universiy of Souh Alabama, Mobile, AL, U.S.A. Absrac Bézier curve flaening

More information

On the Separation Theorem of Stochastic Systems in the Case Of Continuous Observation Channels with Memory

On the Separation Theorem of Stochastic Systems in the Case Of Continuous Observation Channels with Memory Journal of Physics: Conference eries PAPER OPEN ACCE On he eparaion heorem of ochasic ysems in he Case Of Coninuous Observaion Channels wih Memory o cie his aricle: V Rozhova e al 15 J. Phys.: Conf. er.

More information

Fractional Method of Characteristics for Fractional Partial Differential Equations

Fractional Method of Characteristics for Fractional Partial Differential Equations Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics

More information

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

L07. 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 information

Optical Flow I. Guido Gerig CS 6320, Spring 2015

Optical Flow I. Guido Gerig CS 6320, Spring 2015 Opical Flow Guido Gerig CS 6320, Spring 2015 (credis: Marc Pollefeys UNC Chapel Hill, Comp 256 / K.H. Shafique, UCSF, CAP5415 / S. Narasimhan, CMU / Bahadir K. Gunurk, EE 7730 / Bradski&Thrun, Sanford

More information

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems. di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x

More information

Let ( α, β be the eigenvector associated with the eigenvalue λ i

Let ( α, β be the eigenvector associated with the eigenvalue λ i ENGI 940 4.05 - Sabiliy Analysis (Linear) Page 4.5 Le ( α, be he eigenvecor associaed wih he eigenvalue λ i of he coefficien i i) marix A Le c, c be arbirary consans. a b c d Case of real, disinc, negaive

More information

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

SUFFICIENT 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 information

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that ODEs, Homework #4 Soluions. Check ha y ( = is a soluion of he second-order ODE ( cos sin y + y sin y sin = 0 and hen use his fac o find all soluions of he ODE. When y =, we have y = and also y = 0, so

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

More information

From Particles to Rigid Bodies

From Particles to Rigid Bodies Rigid Body Dynamics From Paricles o Rigid Bodies Paricles No roaions Linear velociy v only Rigid bodies Body roaions Linear velociy v Angular velociy ω Rigid Bodies Rigid bodies have boh a posiion and

More information

Extracting a Fluid Dynamic Texture and the Background from Video

Extracting a Fluid Dynamic Texture and the Background from Video Exracing a Fluid Dynamic Texure and he Background from Video Bernard Ghanem and Narendra Ahuja Beckman nsiue for Advanced Science and Technology Deparmen of Elecrical and Compuer Engineering Universiy

More information

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler

Institute 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 information

Dirac Particles Emission from Reissner-Nordstrom-Vaidya Black Hole

Dirac Particles Emission from Reissner-Nordstrom-Vaidya Black Hole Journal o Physics: Conerence Series PAPER OPEN ACCESS Dirac Paricles Emission rom Reissner-Nordsrom-Vaidya Black Hole To cie his aricle: Yuan Tiandho and Triyana 016 J. Phys.: Con. Ser. 739 01146 View

More information

4.5 Constant Acceleration

4.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 information

Single and Double Pendulum Models

Single and Double Pendulum Models Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double

More information

ME 391 Mechanical Engineering Analysis

ME 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 information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

Brock University Physics 1P21/1P91 Fall 2013 Dr. D Agostino. Solutions for Tutorial 3: Chapter 2, Motion in One Dimension

Brock University Physics 1P21/1P91 Fall 2013 Dr. D Agostino. Solutions for Tutorial 3: Chapter 2, Motion in One Dimension Brock Uniersiy Physics 1P21/1P91 Fall 2013 Dr. D Agosino Soluions for Tuorial 3: Chaper 2, Moion in One Dimension The goals of his uorial are: undersand posiion-ime graphs, elociy-ime graphs, and heir

More information

Improved Video Coding Using Long-term Global Motion Compensation

Improved Video Coding Using Long-term Global Motion Compensation Please veriy ha () all pages are presen, (2) all igures are accepable, (3) all ons and special characers are correc, and (4) all e and igures i wihin he Improved Video Coding Using Long-erm Global Moion

More information

Augmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004

Augmented 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 information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Trajectory planning in Cartesian space

Trajectory planning in Cartesian space Roboics 1 Trajecory planning in Caresian space Prof. Alessandro De Luca Roboics 1 1 Trajecories in Caresian space in general, he rajecory planning mehods proposed in he join space can be applied also in

More information

Plasma Astrophysics Chapter 3: Kinetic Theory. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University

Plasma Astrophysics Chapter 3: Kinetic Theory. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University Plasma Asrophysics Chaper 3: Kineic Theory Yosuke Mizuno Insiue o Asronomy Naional Tsing-Hua Universiy Kineic Theory Single paricle descripion: enuous plasma wih srong exernal ields, imporan or gaining

More information

Fitting Parameters in a Differential Equation With a Non-Analytical Solution

Fitting Parameters in a Differential Equation With a Non-Analytical Solution NDSolve u.nb 1 Fiing Parameers in a Dierenial Equaion Wih a Non-Analyical Soluion Mahew Smih 9-Jan-005 In[]:= Remove "Global` " O General::"spell" O General::"spell1" Graphics` Oen dierenial equaions do

More information

Lecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples

Lecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples EE263 Auumn 27-8 Sephen Boyd Lecure 1 Overview course mechanics ouline & opics wha is a linear dynamical sysem? why sudy linear sysems? some examples 1 1 Course mechanics all class info, lecures, homeworks,

More information

Basilio Bona ROBOTICA 03CFIOR 1

Basilio Bona ROBOTICA 03CFIOR 1 Indusrial Robos Kinemaics 1 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

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

Welcome Back to Physics 215!

Welcome Back to Physics 215! Welcome Back o Physics 215! (General Physics I) Thurs. Jan 19 h, 2017 Lecure01-2 1 Las ime: Syllabus Unis and dimensional analysis Today: Displacemen, velociy, acceleraion graphs Nex ime: More acceleraion

More information

Optimal trajectory generation for wheeled mobile robot

Optimal trajectory generation for wheeled mobile robot Opimal rajecory generaion or wheeled mobile robo Bin Qin Yeng Chai Soh Ming Xie Danwei Wang Nanyang echnological Universiy, Nangyang Avenue, Singapore mbqin@nu.edu.sg, eycsoh@nu.edu.sg Absrac. he problem

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math 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 information

k The function Ψ(x) is called wavelet function and shows band-pass behavior. The wavelet coefficients d a,b

k The function Ψ(x) is called wavelet function and shows band-pass behavior. The wavelet coefficients d a,b Wavele Transform Wavele Transform The wavele ransform corresponds o he decomposiion of a quadraic inegrable funcion sx ε L 2 R in a family of scaled and ranslaed funcions Ψ,l, ψ, l 1/2 = ψ l The funcion

More information

1. Kinematics I: Position and Velocity

1. Kinematics I: Position and Velocity 1. Kinemaics I: Posiion and Velociy Inroducion The purpose of his eperimen is o undersand and describe moion. We describe he moion of an objec by specifying is posiion, velociy, and acceleraion. In his

More information

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

EECS 2602 Winter Laboratory 3 Fourier series, Fourier transform and Bode Plots in MATLAB

EECS 2602 Winter Laboratory 3 Fourier series, Fourier transform and Bode Plots in MATLAB EECS 6 Winer 7 Laboraory 3 Fourier series, Fourier ransform and Bode Plos in MATLAB Inroducion: The objecives of his lab are o use MATLAB:. To plo periodic signals wih Fourier series represenaion. To obain

More information

A Note on Fractional Electrodynamics. Abstract

A 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 information

Using the Kalman filter Extended Kalman filter

Using the Kalman filter Extended Kalman filter Using he Kalman filer Eended Kalman filer Doz. G. Bleser Prof. Sricker Compuer Vision: Objec and People Tracking SA- Ouline Recap: Kalman filer algorihm Using Kalman filers Eended Kalman filer algorihm

More information

The physics most important to 20th Century. The most useful interactions - Electricity & PHYSICS. » Physics on the atomic scale.

The physics most important to 20th Century. The most useful interactions - Electricity & PHYSICS. » Physics on the atomic scale. Physic301 Where o in Tomorrow s Assignmen Class Web Page hp://www.physics.umn.edu» Class normaion» 1301.100 Syllabus can be downloaded rom web» Schedule» Grading Lecure noes Te problem soluion oulines

More information

6.003: Signal Processing

6.003: Signal Processing 6.003: Signal Processing Working wih Signals Overview of Subjec Signals: Definiions, Examples, and Operaions Basis Funcions and Transforms Sepember 6, 2018 Welcome o 6.003 Piloing a new version of 6.003

More information

Ordinary Differential Equations

Ordinary Differential Equations Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described

More information

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations Concourse Mah 80 Spring 0 Worked Examples: Marix Mehods for Solving Sysems of s Order Linear Differenial Equaions The Main Idea: Given a sysem of s order linear differenial equaions d x d Ax wih iniial

More information

Physics 20 Lesson 5 Graphical Analysis Acceleration

Physics 20 Lesson 5 Graphical Analysis Acceleration Physics 2 Lesson 5 Graphical Analysis Acceleraion I. Insananeous Velociy From our previous work wih consan speed and consan velociy, we know ha he slope of a posiion-ime graph is equal o he velociy of

More information

Bianchi Type-II, VIII & IX Universe Filled with Wet Dark Fluid in f ( R,

Bianchi Type-II, VIII & IX Universe Filled with Wet Dark Fluid in f ( R, Inernaional Journal o heoreical Mahemaical Physics 0, (5): 9-6 DOI: 059/jmp00050 ianchi ype-ii, VIII & IX Universe Filled wih We Dark Fluid in ( R, K S dhav *, M V Dawe, R G Deshmukh Deparmen o Mahemaics,

More information

Human-Robot Cooperative Manipulation with Motion Estimation

Human-Robot Cooperative Manipulation with Motion Estimation Human-Robo Cooperaive Manipulaion wih Moion Esimaion Yusuke MAEDA, Takayuki HARA and Tamio ARAI (The Universiy o Tokyo) 1. Inroducion 2. Virual Compliance Conrol 3. Esimaion o Human Moion 4. Eperimens

More information

11/16/99 (T.F. Weiss) Lecture #18: Continuous time periodic signals. Outline: Fourier series of periodic functions. Motivation:

11/16/99 (T.F. Weiss) Lecture #18: Continuous time periodic signals. Outline: Fourier series of periodic functions. Motivation: /6/99 F Weiss) Lecure #8: Coninuous ime periodic signals Moivaion: Represenaion o coninuous ime, periodic signals in he requency domain Periodic signals occur requenly moion o planes and heir saellies,

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

More information

Mean-square Stability Control for Networked Systems with Stochastic Time Delay

Mean-square Stability Control for Networked Systems with Stochastic Time Delay JOURNAL OF SIMULAION VOL. 5 NO. May 7 Mean-square Sabiliy Conrol for Newored Sysems wih Sochasic ime Delay YAO Hejun YUAN Fushun School of Mahemaics and Saisics Anyang Normal Universiy Anyang Henan. 455

More information

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

More information

I. OBJECTIVE OF THE EXPERIMENT.

I. 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 information

Ch.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

Ch.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC Ch.. Group Work Unis Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Uni 2 Jusify wheher he following saemens are rue or false: a) Two sreamlines, corresponding o a same insan of ime, can never inersec

More information

STAD57 Time Series Analysis. Lecture 5

STAD57 Time Series Analysis. Lecture 5 STAD57 Time Series Analysis Lecure 5 1 Exploraory Daa Analysis Check if given TS is saionary: µ is consan σ 2 is consan γ(s,) is funcion of h= s If no, ry o make i saionary using some of he mehods below:

More information

Two Coupled Oscillators / Normal Modes

Two 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 information

ENGINEERING FOR RURAL DEVELOPMENT Jelgava, MECHANISM MOTION STUDIES WITH COLLISIONS AT SEVERAL POINTS

ENGINEERING FOR RURAL DEVELOPMENT Jelgava, MECHANISM MOTION STUDIES WITH COLLISIONS AT SEVERAL POINTS EGIEERIG FOR RURAL DEVELOPMET Jelgava, 5.-7.5.16. MECHAISM MOTIO STUDIES WITH COLLISIOS AT SEVERAL POITS Edgars Kovals 1, Janis Viba 1, Lauris Sals 1, Svelana Sokolova 1, Vialy Krupenin 1 Riga Technical

More information

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance Ch: Moion along a sraigh line Moion Posiion and Displacemen Average Velociy and Average Speed Insananeous Velociy and Speed Acceleraion Consan Acceleraion: A Special Case Anoher Look a Consan Acceleraion

More information

Production Analysis Reserves Estimation in Unconventional Reservoirs New Rate-Time Relations

Production Analysis Reserves Estimation in Unconventional Reservoirs New Rate-Time Relations Naural Gas Engineering Producion Analysis Reserves Esimaion in Unconvenional Reservoirs New Rae-Time Relaions T.A. Blasingame, Texas A&M U. Deparmen of Peroleum Engineering Texas A&M Universiy College

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

Physics for Scientists and Engineers I

Physics for Scientists and Engineers I Physics for Scieniss and Engineers I PHY 48, Secion 4 Dr. Beariz Roldán Cuenya Universiy of Cenral Florida, Physics Deparmen, Orlando, FL Chaper - Inroducion I. General II. Inernaional Sysem of Unis III.

More information

Signals and Systems Prof. Brian L. Evans Dept. of Electrical and Computer Engineering The University of Texas at Austin

Signals and Systems Prof. Brian L. Evans Dept. of Electrical and Computer Engineering The University of Texas at Austin EE 345S Real-Time Digial Signal Processing Lab Spring 26 Signals and Sysems Prof. Brian L. Evans Dep. of Elecrical and Compuer Engineering The Universiy of Texas a Ausin Review Signals As Funcions of Time

More information

On Safe Navigation Support System using GPS Compass

On Safe Navigation Support System using GPS Compass Proceedings o he 6 h Inernaional Ship Sabiliy Workshop, 5-7 June 27, Belgrade, Serbia On Sae Navigaion Suppor Sysem using GPS Compass Daisuke erada, Naional Deense Academy, derada@nda.ac.jp asahiro amashima,

More information

Reduced-order Observer Design for Structure and Motion Estimation

Reduced-order Observer Design for Structure and Motion Estimation Reduced-order Observer Design for Srucure and Moion Esimaion Downloaded /8/7 o 46.3.4.58. Redisribuion subjec o SIAM license or copyrigh; see hp://www.siam.org/journals/ojsa.php Absrac Hui Xie Romeo Tasambon

More information

Accurate RMS Calculations for Periodic Signals by. Trapezoidal Rule with the Least Data Amount

Accurate RMS Calculations for Periodic Signals by. Trapezoidal Rule with the Least Data Amount Adv. Sudies Theor. Phys., Vol. 7, 3, no., 3-33 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/.988/asp.3.3999 Accurae RS Calculaions for Periodic Signals by Trapezoidal Rule wih he Leas Daa Amoun Sompop Poomjan,

More information

ln 2 1 ln y x c y C x

ln 2 1 ln y x c y C x Lecure 14 Appendi B: Some sample problems from Boas Here are some soluions o he sample problems assigned for Chaper 8 8: 6 Soluion: We wan o find he soluion o he following firs order equaion using separaion

More information

EMS SCM joint meeting. On stochastic partial differential equations of parabolic type

EMS SCM joint meeting. On stochastic partial differential equations of parabolic type EMS SCM join meeing Barcelona, May 28-30, 2015 On sochasic parial differenial equaions of parabolic ype Isván Gyöngy School of Mahemaics and Maxwell Insiue Edinburgh Universiy 1 I. Filering problem II.

More information

Lecture 10: Wave equation, solution by spherical means

Lecture 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 information

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC CH.7. PLANE LINEAR ELASTICITY Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Overview Plane Linear Elasici Theor Plane Sress Simplifing Hpohesis Srain Field Consiuive Equaion Displacemen Field The Linear

More information

IMPROVED HYPERBOLIC FUNCTION METHOD AND EXACT SOLUTIONS FOR VARIABLE COEFFICIENT BENJAMIN-BONA-MAHONY-BURGERS EQUATION

IMPROVED HYPERBOLIC FUNCTION METHOD AND EXACT SOLUTIONS FOR VARIABLE COEFFICIENT BENJAMIN-BONA-MAHONY-BURGERS EQUATION THERMAL SCIENCE, Year 015, Vol. 19, No. 4, pp. 1183-1187 1183 IMPROVED HYPERBOLIC FUNCTION METHOD AND EXACT SOLUTIONS FOR VARIABLE COEFFICIENT BENJAMIN-BONA-MAHONY-BURGERS EQUATION by Hong-Cai MA a,b*,

More information

Anti-Disturbance Control for Multiple Disturbances

Anti-Disturbance Control for Multiple Disturbances Workshop a 3 ACC Ani-Disurbance Conrol for Muliple Disurbances Lei Guo (lguo@buaa.edu.cn) Naional Key Laboraory on Science and Technology on Aircraf Conrol, Beihang Universiy, Beijing, 9, P.R. China. Presened

More information

Lecture 10 Estimating Nonlinear Regression Models

Lecture 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 information

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary

More information