Multiple view geometry

Size: px
Start display at page:

Download "Multiple view geometry"

Transcription

1 EECS 442 Computer vson Multple vew geometry Perspectve Structure from Moton - Perspectve structure from moton prolem - mgutes - lgerc methods - Fctorzton methods - Bundle djustment - Self-clrton Redng: [HZ] Chpters:,8,9 [FP] Chpter: 3 Some sldes of ths lectures re courtesy of prof. S. Lzenk

2 Structure from moton prolem X j j M m j 2j M 2 M m From the mn correspondences j, estmte: m projecton mtrces M n 3D ponts X j moton structure

3 Structure from moton mguty - Poston mguty: t s mpossle sed on the mges lone to estmte the solute locton nd pose of the scene w.r.t. 3D world coordnte frme

4 Structure from moton mguty -Scle mguty: t s mpossle sed on the mges lone to estmte the solute scle of the scene (.e. house heght)

5 Structure from moton mguty -he scene s determned y the mges only up smlrty trnsformton (rotton, trnslton nd sclng) H s s s R t

6 Structure from moton mguty he mguty ests even for clrted cmers For clrted cmers, the smlrty mguty s the only mguty [Longuet-Hggns 8]

7 Structure from moton mguty In the generl cse (nothng s known) the mguty s epressed y n rtrry ffne or projectve trnsformton M j X j M R H X M H j j j M X j M H - H X j

8 Projectve mguty

9 ffne mguty

10 Structure from moton prolem X j j M m j 2j M 2 M m Gven m mges of n fed 3D ponts j = M X j, =,, m, j =,, n

11 Structure from moton prolem X j j M m j 2j M 2 M m m cmers M M m M

12 he Structure-from-Moton Prolem Gven m mges of n fed ponts X j we cn wrte M j X j Prolem: estmte the m 34 mtrces M nd the n postons X j from the mn correspondences j. Wth no clrton nfo, cmers nd ponts cn only e recovered up to 44 projectve Gven two cmers, how mny ponts re needed? How mny equtons nd how mny unknown? 2m n equtons n m+3n 5 unknowns So 7 ponts! [227 = 28; = 28]

13 Structure-from-Moton lgorthms lgerc pproch (y fundmentl mtr) Fctorzton method (y SVD) Bundle djustment

14 lgerc pproch (2-vew cse) j M 2j M j X j M 2 pply projectve trnsformton H such tht: M H I M H Cnoncl perspectve cmers 2

15 X H X ~ X I X ~ ] [ ~ H M lgerc pproch (Fundmentl mtr) ~ ] [I X ~ X ~ X ~ ] [I X ~ X ~ X ~ ] [ ~ ] [ X X X X ~ ] [ ~ H M 2 ) ( ) (? ) (

16 ] [ z y y z y z Cross product s mtr multplcton

17 X H X ~ ) ( X I X ~ ] [ ~ H M lgerc pproch (Fundmentl mtr) ~ ] [I X ~ X ~ X ~ ] [I X ~ X ~ X ~ ] [ ~ ] [ X X X ~ ] [ H M 2? ] [ F ] [ s ths fmlr? F ) ( ) (

18 Eppolr Constrnt [lecture 6] X 2 e e 2 O O 2 F 2 s the eppolr lne ssocted wth 2 (l = F 2 ) F s the eppolr lne ssocted wth (l 2 = F ) F s sngulr (rnk two) F e 2 = nd F e = F s 33 mtr; 7 DOF

19 lgerc pproch (Fundmentl mtr) F F [ ] F [ ] Snce F =, =e Compute s lest sq. soluton of F = det(f)=; = = [ ] - F = [ ] F M p I M p [e ]F e 2 Perspectve cmers re known

20 Structure-from-Moton lgorthms lgerc pproch (y fundmentl mtr) Fctorzton method (y SVD) Bundle djustment

21 Projectve fctorzton D z z z 2 m 2 m z z z 2 22 m m2 z z z n 2n mn n 2n mn M M M 2 m cmers (3 m 4) X X 2 ponts (4 n) X n D = MS hs rnk 4 If we knew the depths z, we could fctorze D to estmte M nd S If we knew M nd S, we could solve for z Soluton: tertve pproch (lternte etween ove two steps)

22 Structure-from-Moton lgorthms lgerc pproch (y fundmentl mtr) Fctorzton method (y SVD) Bundle djustment

23 Bundle djustment Non-lner method for refnng structure nd moton Mnmzng re-projecton error E(M, X) m n j D 2, X j M j X j? M X j j 3j P M 2 X j 2j M 3 X j P 2 P 3

24 Bundle djustment Non-lner method for refnng structure nd moton Mnmzng re-projecton error E(M, X) m n j dvntges Hndle lrge numer of vews Hndle mssng dt D 2, X j M Lmttons Lrge mnmzton prolem (prmeters grow wth numer of vews) requres good ntl condton j Used s the fnl step of SFM

25 Removng the mgutes: the Strtfed reconstructon up grde reconstructon from perspectve to ffne [y mesurng the plne t nfnty] up grde reconstructon from ffne to metrc [y mesurng the solute conc] Recoverng the metrc reconstructon from the perspectve one s clled self-clrton

26 Self-clrton Process of determnng ntrnsc cmer prmeters drectly from un-clrted mges Suppose we hve projectve reconstructon M,X } { j GOL: fnd rectfyng (non-sngulr) homogrphy H such tht { M H,H X j} M M M X H j s metrc reconstructon m M [R If world ref. system = cmer ref. system: [I ] M If the perspectve cmer s cnoncl: M [I ] ]

27 Self-clrton M M H M [I M [I ] ] [ ] [I ] H H v t k t We cn set k= (ths fes the scle of the reconstructon) H v

28 Plnes t nfnty (lecture 5) z y In the metrc (Euclden) world coordntes 2 plnes re prllel ff ther ntersectons s lne tht elongs to he projectve trnsformton of plne t nfnty cn e epressed s H p

29 Self-clrton H v v v H p v p

30 Self-clrton GOL: fnd rectfyng homogrphy H such tht { M,X j} { M H,H X j} s metrc reconstructon H p cot sn u v o o = clrton mtr of frst cmer [p ] = plne t nfnty n the projectve reconstructon 5 unknowns 3 unknowns

31 Self-clrton sc equton ] [ M = perspectve reconstructon of the cmer (known) p H ] [R M p R m 2 H M M R p p p R = metrc reconstructon of the cmer = rectfyng homogrphy (unknown)

32 Self-clrton sc equton p R p R I R R I p p p p?

33 solute conc ny [From lecture 5] stsfes: s C Projectve trnsformton of 2 3 ( ) * Dul mge of the solute conc

34 ) ( It s not functon of R, symmetrc (5 unknowns) [From lecture 5] Propertes of

35 Self-clrton sc equton p p * * p p =2 m How mny unknowns? 3 from p 5 from How mny equtons? [ nd re known] [per vew] 5 ndependent equtons [per vew] rt of self-clrton: use constrnts on () to generte enough equtons on the unknowns

36 Self-clrton dentcl s * * p p * * p p For m vews, 5(m-) constrnts Numer of unknowns: 8 m>=3 provdes enough constrnts o solve the self-clrton prolem wth dentcl cmers we need t lest 3 vews

37 Propertes of [From lecture 5] ( ) zero-skew squre pel zero-offset

38 Self-clrton other constrnts * * p p zero-offset m lner constrnts zero-skew etc 2 m lner constrnts

39 Self-clrton - summry Condton N. Vews Constnt nternl prmeters 3 spect rto nd skew known Focl length nd offset vry spect rto nd skew known Focl length nd offset vry 4* 5* skew =, ll other prmeters vry 8* Issue: the lrger s the numer of vew, the hrder s the correspondence prolem Bundle djustment helps!

40 SFM prolem - summry. Estmte structure nd moton up perspectve trnsformton. lgerc 2. fctorzton method 3. undle djustment 2. Convert from perspectve to metrc (self-clrton) 3. Bundle djustment ** or **. Bundle djustment wth self-clrton constrnts

41 Self-clrton - summry Constrnts on cmer moton cn e ncorported Lnerly trnsltng cmer Sngle s of rotton: turntle moton

42 pplctons Courtesy of Oford Vsul Geometry Group

43 pplctons D. Nstér, PhD thess

44 pplctons M. Pollefeys et l 98---

45 pplctons M. Brown nd D. G. Lowe. Unsupervsed 3D Oject Recognton nd Reconstructon n Unordered Dtsets. (3DIM25)

46 Photo synth Noh Snvely, Steven M. Setz, Rchrd Szelsk, "Photo toursm: Eplorng photo collectons n 3D," CM rnsctons on Grphcs (SIGGRPH Proceedngs),26,

47

Lecture 3 Camera Models 2 & Camera Calibration. Professor Silvio Savarese Computational Vision and Geometry Lab

Lecture 3 Camera Models 2 & Camera Calibration. Professor Silvio Savarese Computational Vision and Geometry Lab Lecture Cer Models Cer Clbrton rofessor Slvo Svrese Coputtonl Vson nd Geoetry Lb Slvo Svrese Lecture - Jn 7 th, 8 Lecture Cer Models Cer Clbrton Recp of cer odels Cer clbrton proble Cer clbrton wth rdl

More information

Announcements. Image Formation: Outline. The course. Image Formation and Cameras (cont.)

Announcements. Image Formation: Outline. The course. Image Formation and Cameras (cont.) nnouncements Imge Formton nd Cmers (cont.) ssgnment : Cmer & Lenses, gd Trnsformtons, nd Homogrph wll be posted lter tod. CSE 5 Lecture 5 CS5, Fll CS5, Fll CS5, Fll The course rt : The phscs of mgng rt

More information

Lecture 8: Camera Calibration

Lecture 8: Camera Calibration Lecture 8: Cer Clbrton rofessor Fe-Fe L Stnford Vson Lb Fe-Fe L 9-Oct- Wht we wll lern tody? Revew cer preters Affne cer odel (roble Set (Q4)) Cer clbrton Vnshng ponts nd lnes (roble Set (Q)) Redng: [F]

More information

Definition of Tracking

Definition of Tracking Trckng Defnton of Trckng Trckng: Generte some conclusons bout the moton of the scene, objects, or the cmer, gven sequence of mges. Knowng ths moton, predct where thngs re gong to project n the net mge,

More information

Least squares. Václav Hlaváč. Czech Technical University in Prague

Least squares. Václav Hlaváč. Czech Technical University in Prague Lest squres Václv Hlváč Czech echncl Unversty n Prgue hlvc@fel.cvut.cz http://cmp.felk.cvut.cz/~hlvc Courtesy: Fred Pghn nd J.P. Lews, SIGGRAPH 2007 Course; Outlne 2 Lner regresson Geometry of lest-squres

More information

Review of linear algebra. Nuno Vasconcelos UCSD

Review of linear algebra. Nuno Vasconcelos UCSD Revew of lner lgebr Nuno Vsconcelos UCSD Vector spces Defnton: vector spce s set H where ddton nd sclr multplcton re defned nd stsf: ) +( + ) (+ )+ 5) λ H 2) + + H 6) 3) H, + 7) λ(λ ) (λλ ) 4) H, - + 8)

More information

Lecture 8: Camera Calibra0on

Lecture 8: Camera Calibra0on Lecture 8: Cer Clbron rofessor Fe- Fe L Stnford Vson Lb Lecture 8 -! Wht we wll lern tody? Revew cer preters Affne cer odel (roble Set (Q4)) Cer clbron Vnshng ponts nd lnes (roble Set (Q)) Redng: [F] Chpter

More information

CISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting

CISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting CISE 3: umercl Methods Lecture 5 Topc 4 Lest Squres Curve Fttng Dr. Amr Khouh Term Red Chpter 7 of the tetoo c Khouh CISE3_Topc4_Lest Squre Motvton Gven set of epermentl dt 3 5. 5.9 6.3 The reltonshp etween

More information

Lecture 4: Piecewise Cubic Interpolation

Lecture 4: Piecewise Cubic Interpolation Lecture notes on Vrtonl nd Approxmte Methods n Appled Mthemtcs - A Perce UBC Lecture 4: Pecewse Cubc Interpolton Compled 6 August 7 In ths lecture we consder pecewse cubc nterpolton n whch cubc polynoml

More information

Variable time amplitude amplification and quantum algorithms for linear algebra. Andris Ambainis University of Latvia

Variable time amplitude amplification and quantum algorithms for linear algebra. Andris Ambainis University of Latvia Vrble tme mpltude mplfcton nd quntum lgorthms for lner lgebr Andrs Ambns Unversty of Ltv Tlk outlne. ew verson of mpltude mplfcton;. Quntum lgorthm for testng f A s sngulr; 3. Quntum lgorthm for solvng

More information

7.2 Volume. A cross section is the shape we get when cutting straight through an object.

7.2 Volume. A cross section is the shape we get when cutting straight through an object. 7. Volume Let s revew the volume of smple sold, cylnder frst. Cylnder s volume=se re heght. As llustrted n Fgure (). Fgure ( nd (c) re specl cylnders. Fgure () s rght crculr cylnder. Fgure (c) s ox. A

More information

Lecture 3. Camera Models 2 & Camera Calibration. Professor Silvio Savarese Computational Vision and Geometry Lab. 13- Jan- 15.

Lecture 3. Camera Models 2 & Camera Calibration. Professor Silvio Savarese Computational Vision and Geometry Lab. 13- Jan- 15. Lecture Caera Models Caera Calbraton rofessor Slvo Savarese Coputatonal Vson and Geoetry Lab Slvo Savarese Lecture - - Jan- 5 Lecture Caera Models Caera Calbraton Recap of caera odels Caera calbraton proble

More information

Principle Component Analysis

Principle Component Analysis Prncple Component Anlyss Jng Go SUNY Bufflo Why Dmensonlty Reducton? We hve too mny dmensons o reson bout or obtn nsghts from o vsulze oo much nose n the dt Need to reduce them to smller set of fctors

More information

Model Fitting and Robust Regression Methods

Model Fitting and Robust Regression Methods Dertment o Comuter Engneerng Unverst o Clorn t Snt Cruz Model Fttng nd Robust Regresson Methods CMPE 64: Imge Anlss nd Comuter Vson H o Fttng lnes nd ellses to mge dt Dertment o Comuter Engneerng Unverst

More information

Rank One Update And the Google Matrix by Al Bernstein Signal Science, LLC

Rank One Update And the Google Matrix by Al Bernstein Signal Science, LLC Introducton Rnk One Updte And the Google Mtrx y Al Bernsten Sgnl Scence, LLC www.sgnlscence.net here re two dfferent wys to perform mtrx multplctons. he frst uses dot product formulton nd the second uses

More information

6 Roots of Equations: Open Methods

6 Roots of Equations: Open Methods HK Km Slghtly modfed 3//9, /8/6 Frstly wrtten t Mrch 5 6 Roots of Equtons: Open Methods Smple Fed-Pont Iterton Newton-Rphson Secnt Methods MATLAB Functon: fzero Polynomls Cse Study: Ppe Frcton Brcketng

More information

UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II

UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II Mcroeconomc Theory I UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS MSc n Economcs MICROECONOMIC THEORY I Techng: A Lptns (Note: The number of ndctes exercse s dffculty level) ()True or flse? If V( y )

More information

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Chapter Newton-Raphson Method of Solving a Nonlinear Equation Chpter.4 Newton-Rphson Method of Solvng Nonlner Equton After redng ths chpter, you should be ble to:. derve the Newton-Rphson method formul,. develop the lgorthm of the Newton-Rphson method,. use the Newton-Rphson

More information

VECTORS AND TENSORS IV.1.1. INTRODUCTION

VECTORS AND TENSORS IV.1.1. INTRODUCTION Chpter IV Vector nd Tensor Anlyss IV. Vectors nd Tensors Septemer 5, 08 05 IV. VECTORS AND TENSORS IV... INTRODUCTION In mthemtcs nd mechncs, we he to operte wth qunttes whch requre dfferent mthemtcl ojects

More information

Structure from Motion. Forsyth&Ponce: Chap. 12 and 13 Szeliski: Chap. 7

Structure from Motion. Forsyth&Ponce: Chap. 12 and 13 Szeliski: Chap. 7 Structure from Moton Forsyth&once: Chap. 2 and 3 Szelsk: Chap. 7 Introducton to Structure from Moton Forsyth&once: Chap. 2 Szelsk: Chap. 7 Structure from Moton Intro he Reconstructon roblem p 3?? p p 2

More information

International Journal of Pure and Applied Sciences and Technology

International Journal of Pure and Applied Sciences and Technology Int. J. Pure Appl. Sc. Technol., () (), pp. 44-49 Interntonl Journl of Pure nd Appled Scences nd Technolog ISSN 9-67 Avlle onlne t www.jopst.n Reserch Pper Numercl Soluton for Non-Lner Fredholm Integrl

More information

VECTORS VECTORS VECTORS VECTORS. 2. Vector Representation. 1. Definition. 3. Types of Vectors. 5. Vector Operations I. 4. Equal and Opposite Vectors

VECTORS VECTORS VECTORS VECTORS. 2. Vector Representation. 1. Definition. 3. Types of Vectors. 5. Vector Operations I. 4. Equal and Opposite Vectors 1. Defnton A vetor s n entt tht m represent phsl quntt tht hs mgntude nd dreton s opposed to slr tht ls dreton.. Vetor Representton A vetor n e represented grphll n rrow. The length of the rrow s the mgntude

More information

Jens Siebel (University of Applied Sciences Kaiserslautern) An Interactive Introduction to Complex Numbers

Jens Siebel (University of Applied Sciences Kaiserslautern) An Interactive Introduction to Complex Numbers Jens Sebel (Unversty of Appled Scences Kserslutern) An Interctve Introducton to Complex Numbers 1. Introducton We know tht some polynoml equtons do not hve ny solutons on R/. Exmple 1.1: Solve x + 1= for

More information

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Chapter Newton-Raphson Method of Solving a Nonlinear Equation Chpter 0.04 Newton-Rphson Method o Solvng Nonlner Equton Ater redng ths chpter, you should be ble to:. derve the Newton-Rphson method ormul,. develop the lgorthm o the Newton-Rphson method,. use the Newton-Rphson

More information

Math 497C Sep 17, Curves and Surfaces Fall 2004, PSU

Math 497C Sep 17, Curves and Surfaces Fall 2004, PSU Mth 497C Sep 17, 004 1 Curves nd Surfces Fll 004, PSU Lecture Notes 3 1.8 The generl defnton of curvture; Fox-Mlnor s Theorem Let α: [, b] R n be curve nd P = {t 0,...,t n } be prtton of [, b], then the

More information

An Introduction to Support Vector Machines

An Introduction to Support Vector Machines An Introducton to Support Vector Mchnes Wht s good Decson Boundry? Consder two-clss, lnerly seprble clssfcton problem Clss How to fnd the lne (or hyperplne n n-dmensons, n>)? Any de? Clss Per Lug Mrtell

More information

Partially Observable Systems. 1 Partially Observable Markov Decision Process (POMDP) Formalism

Partially Observable Systems. 1 Partially Observable Markov Decision Process (POMDP) Formalism CS294-40 Lernng for Rootcs nd Control Lecture 10-9/30/2008 Lecturer: Peter Aeel Prtlly Oservle Systems Scre: Dvd Nchum Lecture outlne POMDP formlsm Pont-sed vlue terton Glol methods: polytree, enumerton,

More information

FUNDAMENTALS ON ALGEBRA MATRICES AND DETERMINANTS

FUNDAMENTALS ON ALGEBRA MATRICES AND DETERMINANTS Dol Bgyoko (0 FUNDAMENTALS ON ALGEBRA MATRICES AND DETERMINANTS Introducton Expressons of the form P(x o + x + x + + n x n re clled polynomls The coeffcents o,, n re ndependent of x nd the exponents 0,,,

More information

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x) DCDM BUSINESS SCHOOL NUMEICAL METHODS (COS -8) Solutons to Assgnment Queston Consder the followng dt: 5 f() 8 7 5 () Set up dfference tble through fourth dfferences. (b) Wht s the mnmum degree tht n nterpoltng

More information

Haddow s Experiment:

Haddow s Experiment: schemtc drwng of Hddow's expermentl set-up movng pston non-contctng moton sensor bems of sprng steel poston vres to djust frequences blocks of sold steel shker Hddow s Experment: terr frm Theoretcl nd

More information

Effects of polarization on the reflected wave

Effects of polarization on the reflected wave Lecture Notes. L Ros PPLIED OPTICS Effects of polrzton on the reflected wve Ref: The Feynmn Lectures on Physcs, Vol-I, Secton 33-6 Plne of ncdence Z Plne of nterfce Fg. 1 Y Y r 1 Glss r 1 Glss Fg. Reflecton

More information

4. Eccentric axial loading, cross-section core

4. Eccentric axial loading, cross-section core . Eccentrc xl lodng, cross-secton core Introducton We re strtng to consder more generl cse when the xl force nd bxl bendng ct smultneousl n the cross-secton of the br. B vrtue of Snt-Vennt s prncple we

More information

Proof that if Voting is Perfect in One Dimension, then the First. Eigenvector Extracted from the Double-Centered Transformed

Proof that if Voting is Perfect in One Dimension, then the First. Eigenvector Extracted from the Double-Centered Transformed Proof tht f Votng s Perfect n One Dmenson, then the Frst Egenvector Extrcted from the Doule-Centered Trnsformed Agreement Score Mtrx hs the Sme Rn Orderng s the True Dt Keth T Poole Unversty of Houston

More information

The Schur-Cohn Algorithm

The Schur-Cohn Algorithm Modelng, Estmton nd Otml Flterng n Sgnl Processng Mohmed Njm Coyrght 8, ISTE Ltd. Aendx F The Schur-Cohn Algorthm In ths endx, our m s to resent the Schur-Cohn lgorthm [] whch s often used s crteron for

More information

COMPLEX NUMBERS INDEX

COMPLEX NUMBERS INDEX COMPLEX NUMBERS INDEX. The hstory of the complex numers;. The mgnry unt I ;. The Algerc form;. The Guss plne; 5. The trgonometrc form;. The exponentl form; 7. The pplctons of the complex numers. School

More information

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede. with respect to λ. 1. χ λ χ λ ( ) λ, and thus:

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede. with respect to λ. 1. χ λ χ λ ( ) λ, and thus: More on χ nd errors : uppose tht we re fttng for sngle -prmeter, mnmzng: If we epnd The vlue χ ( ( ( ; ( wth respect to. χ n Tlor seres n the vcnt of ts mnmum vlue χ ( mn χ χ χ χ + + + mn mnmzes χ, nd

More information

Quiz: Experimental Physics Lab-I

Quiz: Experimental Physics Lab-I Mxmum Mrks: 18 Totl tme llowed: 35 mn Quz: Expermentl Physcs Lb-I Nme: Roll no: Attempt ll questons. 1. In n experment, bll of mss 100 g s dropped from heght of 65 cm nto the snd contner, the mpct s clled

More information

Estimating the Fundamental Matrix by Transforming Image Points in Projective Space 1

Estimating the Fundamental Matrix by Transforming Image Points in Projective Space 1 Estmatng the Fundamental Matrx by Transformng Image Ponts n Projectve Space 1 Zhengyou Zhang and Charles Loop Mcrosoft Research, One Mcrosoft Way, Redmond, WA 98052, USA E-mal: fzhang,cloopg@mcrosoft.com

More information

Pyramid Algorithms for Barycentric Rational Interpolation

Pyramid Algorithms for Barycentric Rational Interpolation Pyrmd Algorthms for Brycentrc Rtonl Interpolton K Hormnn Scott Schefer Astrct We present new perspectve on the Floter Hormnn nterpolnt. Ths nterpolnt s rtonl of degree (n, d), reproduces polynomls of degree

More information

consider in the case of 1) internal resonance ω 2ω and 2) external resonance Ω ω and small damping

consider in the case of 1) internal resonance ω 2ω and 2) external resonance Ω ω and small damping consder n the cse o nternl resonnce nd externl resonnce Ω nd smll dmpng recll rom "Two_Degs_Frdm_.ppt" tht θ + μ θ + θ = θφ + cos Ω t + τ where = k α α nd φ + μ φ + φ = θ + cos Ω t where = α τ s constnt

More information

Advanced Machine Learning. An Ising model on 2-D image

Advanced Machine Learning. An Ising model on 2-D image Advnced Mchne Lernng Vrtonl Inference Erc ng Lecture 12, August 12, 2009 Redng: Erc ng Erc ng @ CMU, 2006-2009 1 An Isng model on 2-D mge odes encode hdden nformton ptchdentty. They receve locl nformton

More information

CS4495/6495 Introduction to Computer Vision. 3C-L3 Calibrating cameras

CS4495/6495 Introduction to Computer Vision. 3C-L3 Calibrating cameras CS4495/6495 Introducton to Computer Vson 3C-L3 Calbratng cameras Fnally (last tme): Camera parameters Projecton equaton the cumulatve effect of all parameters: M (3x4) f s x ' 1 0 0 0 c R 0 I T 3 3 3 x1

More information

(x',y') (x', y') (x, y) (x, y) Review: Rendering Pipeline. Review: Geometry Pipeline. Review: Graphics State. Transformations

(x',y') (x', y') (x, y) (x, y) Review: Rendering Pipeline. Review: Geometry Pipeline. Review: Graphics State. Transformations Unverst of rtsh olum PS 34 omputer Grphcs Jn-pr 23 Tmr Munzner Trnsformtons Geometr Dtse Scn onverson Revew: Renderng Ppelne Model/Vew Trnsform. Teturng Lghtng Depth Test Perspectve Trnsform. lendng lppng

More information

A New Algorithm Linear Programming

A New Algorithm Linear Programming A New Algorthm ner Progrmmng Dhnnjy P. ehendle Sr Prshurmhu College, Tlk Rod, Pune-400, Ind dhnnjy.p.mehendle@gml.com Astrct In ths pper we propose two types of new lgorthms for lner progrmmng. The frst

More information

Abhilasha Classes Class- XII Date: SOLUTION (Chap - 9,10,12) MM 50 Mob no

Abhilasha Classes Class- XII Date: SOLUTION (Chap - 9,10,12) MM 50 Mob no hlsh Clsses Clss- XII Dte: 0- - SOLUTION Chp - 9,0, MM 50 Mo no-996 If nd re poston vets of nd B respetvel, fnd the poston vet of pont C n B produed suh tht C B vet r C B = where = hs length nd dreton

More information

COMPLEX NUMBER & QUADRATIC EQUATION

COMPLEX NUMBER & QUADRATIC EQUATION MCQ COMPLEX NUMBER & QUADRATIC EQUATION Syllus : Comple numers s ordered prs of rels, Representton of comple numers n the form + nd ther representton n plne, Argnd dgrm, lger of comple numers, modulus

More information

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 More oundr-vlue Prolems nd genvlue Prolems n Os ovemer 9, 7 More oundr-vlue Prolems nd genvlue Prolems n Os Lrr retto Menl ngneerng 5 Semnr n ngneerng nlss ovemer 9, 7 Outlne Revew oundr-vlue prolems Soot

More information

GAUSS ELIMINATION. Consider the following system of algebraic linear equations

GAUSS ELIMINATION. Consider the following system of algebraic linear equations Numercl Anlyss for Engneers Germn Jordnn Unversty GAUSS ELIMINATION Consder the followng system of lgebrc lner equtons To solve the bove system usng clsscl methods, equton () s subtrcted from equton ()

More information

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 9

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 9 CS434/541: Pttern Recognton Prof. Olg Veksler Lecture 9 Announcements Fnl project proposl due Nov. 1 1-2 prgrph descrpton Lte Penlt: s 1 pont off for ech d lte Assgnment 3 due November 10 Dt for fnl project

More information

Electrochemical Thermodynamics. Interfaces and Energy Conversion

Electrochemical Thermodynamics. Interfaces and Energy Conversion CHE465/865, 2006-3, Lecture 6, 18 th Sep., 2006 Electrochemcl Thermodynmcs Interfces nd Energy Converson Where does the energy contrbuton F zϕ dn come from? Frst lw of thermodynmcs (conservton of energy):

More information

Katholieke Universiteit Leuven Department of Computer Science

Katholieke Universiteit Leuven Department of Computer Science Updte Rules for Weghted Non-negtve FH*G Fctorzton Peter Peers Phlp Dutré Report CW 440, Aprl 006 Ktholeke Unverstet Leuven Deprtment of Computer Scence Celestjnenln 00A B-3001 Heverlee (Belgum) Updte Rules

More information

CONIC SECTIONS. Chapter 11

CONIC SECTIONS. Chapter 11 CONIC SECTIONS Chpter. Overview.. Sections of cone Let l e fied verticl line nd m e nother line intersecting it t fied point V nd inclined to it t n ngle α (Fig..). Fig.. Suppose we rotte the line m round

More information

perturbation theory and its applications

perturbation theory and its applications Second-order order guge-nvrnt perturton theory nd ts pplctons (Short revew of my poster presentton) Some detls cn e seen n my poster Kouj Nkmur (Grd. Unv. Adv. Stud. (NAOJ)) References : K.N. Prog. Theor.

More information

Introduction to Numerical Integration Part II

Introduction to Numerical Integration Part II Introducton to umercl Integrton Prt II CS 75/Mth 75 Brn T. Smth, UM, CS Dept. Sprng, 998 4/9/998 qud_ Intro to Gussn Qudrture s eore, the generl tretment chnges the ntegrton prolem to ndng the ntegrl w

More information

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3 2 The Prllel Circuit Electric Circuits: Figure 2- elow show ttery nd multiple resistors rrnged in prllel. Ech resistor receives portion of the current from the ttery sed on its resistnce. The split is

More information

Applied Statistics Qualifier Examination

Applied Statistics Qualifier Examination Appled Sttstcs Qulfer Exmnton Qul_june_8 Fll 8 Instructons: () The exmnton contns 4 Questons. You re to nswer 3 out of 4 of them. () You my use ny books nd clss notes tht you mght fnd helpful n solvng

More information

Geometric Correction or Georeferencing

Geometric Correction or Georeferencing Geoetrc Correcton or Georeferencng GEOREFERENCING: fro ge to p Coordntes on erth: (λ, φ) ge: (, ) p: (, ) rel nteger Trnsfortons (nvolvng deforton): erth-to-ge: χ erth-to-p: ψ (crtogrphc proecton) ge-to-p:

More information

Course Review Introduction to Computer Methods

Course Review Introduction to Computer Methods Course Revew Wht you hopefully hve lerned:. How to nvgte nsde MIT computer system: Athen, UNIX, emcs etc. (GCR). Generl des bout progrmmng (GCR): formultng the problem, codng n Englsh trnslton nto computer

More information

CS 532: 3D Computer Vision 2 nd Set of Notes

CS 532: 3D Computer Vision 2 nd Set of Notes CS 532: 3D Computer Vson 2 nd Set of Notes Instructor: Plppos Mordoa Webpage: www.cs.stevens.edu/~mordoa E-mal: Plppos.Mordoa@stevens.edu Offce: Leb 25 Lecture Outlne 2D projectve transformatons Homograpes

More information

PART 1: VECTOR & TENSOR ANALYSIS

PART 1: VECTOR & TENSOR ANALYSIS PART : VECTOR & TENSOR ANALYSIS wth LINEAR ALGEBRA Obectves Introduce the concepts, theores, nd opertonl mplementton of vectors, nd more generlly tensors, n dvnced engneerng nlyss. The emphss s on geometrc

More information

Dennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1

Dennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1 Denns Brcker, 2001 Dept of Industrl Engneerng The Unversty of Iow MDP: Tx pge 1 A tx serves three djcent towns: A, B, nd C. Ech tme the tx dschrges pssenger, the drver must choose from three possble ctons:

More information

Demand. Demand and Comparative Statics. Graphically. Marshallian Demand. ECON 370: Microeconomic Theory Summer 2004 Rice University Stanley Gilbert

Demand. Demand and Comparative Statics. Graphically. Marshallian Demand. ECON 370: Microeconomic Theory Summer 2004 Rice University Stanley Gilbert Demnd Demnd nd Comrtve Sttcs ECON 370: Mcroeconomc Theory Summer 004 Rce Unversty Stnley Glbert Usng the tools we hve develoed u to ths ont, we cn now determne demnd for n ndvdul consumer We seek demnd

More information

Basic projective geometry

Basic projective geometry Bsc projectve geoetr Augusto Srt Stefno ubro Mrco Mrcon Fbo Antoncc Dprtento d Eettronc Inforzone e Bongegner Potecnco d Mno Overvew Bsc projectve geoetr (P nd P ) rnsfortons n P Affne rectfcton n P Metrc

More information

Correction and rectification of light field

Correction and rectification of light field Correcton nd rectfcton of lgt feld Ke Deng, Lfeng Wng b, Zoucen Ln b, To Feng b, Zdong Deng, Bnng Guo b* Tsngu Unverst, b Mcrosoft Reserc As ABSTRACT Te lgt feld s well known mge-bsed renderng tecnolog.

More information

18.7 Artificial Neural Networks

18.7 Artificial Neural Networks 310 18.7 Artfcl Neurl Networks Neuroscence hs hypotheszed tht mentl ctvty conssts prmrly of electrochemcl ctvty n networks of brn cells clled neurons Ths led McCulloch nd Ptts to devse ther mthemtcl model

More information

Support vector machines for regression

Support vector machines for regression S 75 Mchne ernng ecture 5 Support vector mchnes for regresson Mos Huskrecht mos@cs.ptt.edu 539 Sennott Squre S 75 Mchne ernng he decson oundr: ˆ he decson: Support vector mchnes ˆ α SV ˆ sgn αˆ SV!!: Decson

More information

Lecture 36. Finite Element Methods

Lecture 36. Finite Element Methods CE 60: Numercl Methods Lecture 36 Fnte Element Methods Course Coordntor: Dr. Suresh A. Krth, Assocte Professor, Deprtment of Cvl Engneerng, IIT Guwht. In the lst clss, we dscussed on the ppromte methods

More information

CMSC 828D: Fundamentals of Computer Vision Homework 4

CMSC 828D: Fundamentals of Computer Vision Homework 4 CMSC 828D: Fundmentls of Computer Vision Homework 4 nstructors: Lrry Dvis, Rmni Duriswmi, Dniel DeMenthon, nd Yinnis loimonos Solution sed on homework sumitted y Hiying Liu. Show y direct clcultion for

More information

Matrix Eigenvalues and Eigenvectors September 13, 2017

Matrix Eigenvalues and Eigenvectors September 13, 2017 Mtri Eigenvlues nd Eigenvectors September, 7 Mtri Eigenvlues nd Eigenvectors Lrry Cretto Mechnicl Engineering 5A Seminr in Engineering Anlysis September, 7 Outline Review lst lecture Definition of eigenvlues

More information

CALIBRATION OF SMALL AREA ESTIMATES IN BUSINESS SURVEYS

CALIBRATION OF SMALL AREA ESTIMATES IN BUSINESS SURVEYS CALIBRATION OF SMALL AREA ESTIMATES IN BUSINESS SURVES Rodolphe Prm, Ntle Shlomo Southmpton Sttstcl Scences Reserch Insttute Unverst of Southmpton Unted Kngdom SAE, August 20 The BLUE-ETS Project s fnnced

More information

A Theoretical Study on the Rank of the Integral Operators for Large- Scale Electrodynamic Analysis

A Theoretical Study on the Rank of the Integral Operators for Large- Scale Electrodynamic Analysis Purdue Unversty Purdue e-pubs ECE Techncl Reports Electrcl nd Computer Engneerng -22-2 A Theoretcl Study on the Rnk of the Integrl Opertors for Lrge- Scle Electrodynmc Anlyss Wenwen Ch Purdue Unversty,

More information

Systematic Construction of examples for cycling in the simplex method

Systematic Construction of examples for cycling in the simplex method Systemtc Constructon of emples for cyclng n the smple method Prof. r. Peter Zörng Prof. Assocdo Unversdde de Brsíl UnB eprtmento de Esttístc e-ml: peter@un.r ) Lner progrmmng prolem cnoncl form: stndrd

More information

Work and Energy (Work Done by a Varying Force)

Work and Energy (Work Done by a Varying Force) Lecture 1 Chpter 7 Physcs I 3.5.14 ork nd Energy (ork Done y Vryng Force) Course weste: http://fculty.uml.edu/andry_dnylov/techng/physcsi Lecture Cpture: http://echo36.uml.edu/dnylov13/physcs1fll.html

More information

Remember: Project Proposals are due April 11.

Remember: Project Proposals are due April 11. Bonformtcs ecture Notes Announcements Remember: Project Proposls re due Aprl. Clss 22 Aprl 4, 2002 A. Hdden Mrov Models. Defntons Emple - Consder the emple we tled bout n clss lst tme wth the cons. However,

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

In this Chapter. Chap. 3 Markov chains and hidden Markov models. Probabilistic Models. Example: CpG Islands

In this Chapter. Chap. 3 Markov chains and hidden Markov models. Probabilistic Models. Example: CpG Islands In ths Chpter Chp. 3 Mrov chns nd hdden Mrov models Bontellgence bortory School of Computer Sc. & Eng. Seoul Ntonl Unversty Seoul 5-74, Kore The probblstc model for sequence nlyss HMM (hdden Mrov model)

More information

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz

More information

Linear and Nonlinear Optimization

Linear and Nonlinear Optimization Lner nd Nonlner Optmzton Ynyu Ye Deprtment of Mngement Scence nd Engneerng Stnford Unversty Stnford, CA 9430, U.S.A. http://www.stnford.edu/~yyye http://www.stnford.edu/clss/msnde/ Ynyu Ye, Stnford, MS&E

More information

Mathematics. Area under Curve.

Mathematics. Area under Curve. Mthemtics Are under Curve www.testprepkrt.com Tle of Content 1. Introduction.. Procedure of Curve Sketching. 3. Sketching of Some common Curves. 4. Are of Bounded Regions. 5. Sign convention for finding

More information

Linear Regression & Least Squares!

Linear Regression & Least Squares! Lner Regresson & Lest Squres Al Borj UWM CS 790 Slde credt: Aykut Erdem Ths&week Lner&regresson&prolem&& ' con0nuous&outputs& ' smple&model Introduce&key&concepts:&& ' loss&func0ons& ' generlz0on& ' op0mz0on&

More information

Many-Body Calculations of the Isotope Shift

Many-Body Calculations of the Isotope Shift Mny-Body Clcultons of the Isotope Shft W. R. Johnson Mrch 11, 1 1 Introducton Atomc energy levels re commonly evluted ssumng tht the nucler mss s nfnte. In ths report, we consder correctons to tomc levels

More information

ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS

ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS Numercl Alyss for Egeers Germ Jord Uversty ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS Numercl soluto of lrge systems of ler lgerc equtos usg drect methods such s Mtr Iverse, Guss

More information

? plate in A G in

? plate in A G in Proble (0 ponts): The plstc block shon s bonded to rgd support nd to vertcl plte to hch 0 kp lod P s ppled. Knong tht for the plstc used G = 50 ks, deterne the deflecton of the plte. Gven: G 50 ks, P 0

More information

Strong Gravity and the BKL Conjecture

Strong Gravity and the BKL Conjecture Introducton Strong Grvty nd the BKL Conecture Dvd Slon Penn Stte October 16, 2007 Dvd Slon Strong Grvty nd the BKL Conecture Introducton Outlne The BKL Conecture Ashtekr Vrbles Ksner Sngulrty 1 Introducton

More information

Lecture 4 Single View Metrology

Lecture 4 Single View Metrology Lecture 4 Single View Metrology Professor Silio Srese Computtionl Vision nd Geometry Lb Silio Srese Lecture 4-6-Jn-5 Lecture 4 Single View Metrology Reiew clibrtion nd 2D trnsformtions Vnishing points

More information

The practical version

The practical version Roerto s Notes on Integrl Clculus Chpter 4: Definite integrls nd the FTC Section 7 The Fundmentl Theorem of Clculus: The prcticl version Wht you need to know lredy: The theoreticl version of the FTC. Wht

More information

The area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O

The area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O 1 Section 5. The Definite Integrl Suppose tht function f is continuous nd positive over n intervl [, ]. y = f(x) x The re under the grph of f nd ove the x-xis etween nd is denoted y f(x) dx nd clled the

More information

Chapter 5 Supplemental Text Material R S T. ij i j ij ijk

Chapter 5 Supplemental Text Material R S T. ij i j ij ijk Chpter 5 Supplementl Text Mterl 5-. Expected Men Squres n the Two-fctor Fctorl Consder the two-fctor fxed effects model y = µ + τ + β + ( τβ) + ε k R S T =,,, =,,, k =,,, n gven s Equton (5-) n the textook.

More information

Chemical Reaction Engineering

Chemical Reaction Engineering Lecture 20 hemcl Recton Engneerng (RE) s the feld tht studes the rtes nd mechnsms of chemcl rectons nd the desgn of the rectors n whch they tke plce. Lst Lecture Energy Blnce Fundmentls F 0 E 0 F E Q W

More information

Nonabelian Dualization of Plane Wave Backgrounds

Nonabelian Dualization of Plane Wave Backgrounds Journl of Modern Physcs 88-95 http://dxdoorg/6/jmp9 Pulshed Onlne Septemer (http://wwwscrporg/journl/jmp) Noneln Dulzton of Plne Wve Bckgrounds Ldslv Hlvtý Mroslv Turek Fculty of Nucler Scences nd Physcl

More information

PHYS 2421 Fields and Waves

PHYS 2421 Fields and Waves PHYS 242 Felds nd Wves Instucto: Joge A. López Offce: PSCI 29 A, Phone: 747-7528 Textook: Unvesty Physcs e, Young nd Feedmn 23. Electc potentl enegy 23.2 Electc potentl 23.3 Clcultng electc potentl 23.4

More information

ESCI 342 Atmospheric Dynamics I Lesson 1 Vectors and Vector Calculus

ESCI 342 Atmospheric Dynamics I Lesson 1 Vectors and Vector Calculus ESI 34 tmospherc Dnmcs I Lesson 1 Vectors nd Vector lculus Reference: Schum s Outlne Seres: Mthemtcl Hndbook of Formuls nd Tbles Suggested Redng: Mrtn Secton 1 OORDINTE SYSTEMS n orthonorml coordnte sstem

More information

The Delta Parallel Robot: Kinematics Solutions Robert L. Williams II, Ph.D., Mechanical Engineering, Ohio University, April 2015

The Delta Parallel Robot: Kinematics Solutions Robert L. Williams II, Ph.D., Mechanical Engineering, Ohio University, April 2015 he Delt rllel Robot: Knemtcs Solutons Robert L. Wllms II, h.d., wllr4@oho.edu Mechncl Engneerng, Oho Unversty, Aprl 05 Clvel s Delt Robot s rgubly the most successful commercl prllel robot to dte. he left

More information

1B40 Practical Skills

1B40 Practical Skills B40 Prcticl Skills Comining uncertinties from severl quntities error propgtion We usully encounter situtions where the result of n experiment is given in terms of two (or more) quntities. We then need

More information

Review of Gaussian Quadrature method

Review of Gaussian Quadrature method Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge

More information

INTRODUCTION TO COMPLEX NUMBERS

INTRODUCTION TO COMPLEX NUMBERS INTRODUCTION TO COMPLEX NUMBERS The numers -4, -3, -, -1, 0, 1,, 3, 4 represent the negtve nd postve rel numers termed ntegers. As one frst lerns n mddle school they cn e thought of s unt dstnce spced

More information

Chemical Reaction Engineering

Chemical Reaction Engineering Lecture 20 hemcl Recton Engneerng (RE) s the feld tht studes the rtes nd mechnsms of chemcl rectons nd the desgn of the rectors n whch they tke plce. Lst Lecture Energy Blnce Fundmentls F E F E + Q! 0

More information

PROPERTIES OF AREAS In general, and for an irregular shape, the definition of the centroid at position ( x, y) is given by

PROPERTIES OF AREAS In general, and for an irregular shape, the definition of the centroid at position ( x, y) is given by PROPERTES OF RES Centroid The concept of the centroid is prol lred fmilir to ou For plne shpe with n ovious geometric centre, (rectngle, circle) the centroid is t the centre f n re hs n is of smmetr, the

More information

CENTROID (AĞIRLIK MERKEZİ )

CENTROID (AĞIRLIK MERKEZİ ) CENTOD (ĞLK MEKEZİ ) centrod s geometrcl concept rsng from prllel forces. Tus, onl prllel forces possess centrod. Centrod s tougt of s te pont were te wole wegt of pscl od or sstem of prtcles s lumped.

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information