Distance Between Two Ellipses in 3D
|
|
- Aubrie Terry
- 5 years ago
- Views:
Transcription
1 Disance Beween Two Ellipses in 3D David Eberly Magic Sofware 6006 Meadow Run Cour Chapel Hill, NC Inroducion An ellipse in 3D is represened by a cener C, uni lengh axes U and V wih corresponding axis lenghs a and b, and a plane conaining he ellipse, N ( X C) = 0 where N is a uni lengh normal o he plane. The vecors U, V, and N form a righ handed orhonormal coordinae sysem (he marix wih hese vecors as columns is orhonormal wih deerminan 1). The ellipse is parameerized as X = C + a cos(θ) U + b sin(θ) V for angles θ [0, 2π). The ellipse is also defined by he wo polynomial equaions N ( X C) = 0 ( X C) ( U U a 2 + V V b 2 ) ( X C) = 1 where he las equaion is wrien as a quadraic form. The firs equaion defines a plane and he second equaion defines an ellipsoid. The inersecion of plane and ellipsoid is an ellipse. 2 Soluion as Polynomial Sysem The wo ellipses are N 0 ( X C 0 ) = 0 and ( X C 0 ) A 0 ( X C 0 ) = 1 where A 0 = U 0U 0 /a V 0V 0 /b 2 0 and N 1 ( Y C 1 ) = 0 and ( Y C 1 ) A 1 ( Y C 1 ) = 1 where A 1 = U 1U 1 /a V 1V 1 /b 2 1. The problem is o minimize he squared disance Y Y 2 subjec o he four consrains menioned above. The problem can be solved wih he mehod of Lagrange mulipliers. Inroduce four new parameers, α, β, γ, and δ and minimize F ( X, Y ; α, β, γ, δ) = X Y 2 + α(( X C 0 ) A 0 ( X C 0 ) 1) + β( N 0 ( X C 0 ) 0) +γ(( Y C 1 ) A 1 ( Y C 1 ) 1) + δ( N 1 ( Y C 1 ) 0). 1
2 Taking derivaives yields F X = 2( X Y ) + 2αA 0 ( X C 0 ) + β N 0 F Y = 2( X Y ) + 2γA 1 ( Y C 1 ) + δ N 1 F α = ( X C 0 ) A 0 ( X C 0 ) 1 F β = N 0 ( X C 0 ) F γ = ( Y C 1 ) A 1 ( Y C 1 ) 1 F δ = N 1 ( Y C 1 ) Seing he las four equaions o zero yields he four original consrains. Seing he firs equaion o he zero vecor and muliplying by ( X C 0 ) yields α = 2( X C 0 ) ( X Y ). Seing he firs equaion o he zero vecor and muliplying by N 0 yields Similar manipulaions of he second equaion yield and The firs wo derivaive equaions become β = 2 N 0( X Y ). γ = 2( Y C 1 ) ( X Y ) δ = 2 N 1( X Y ). M 0 ( X Y ( ) = N0 N 0 + A 0 ( X C 0 )( X C ) 0 ) I ( X Y ) = 0 M 1 ( X Y ( ) = N1 N 1 + A 1 ( Y C 1 )( Y C ) 1 ) I ( X Y ) = 0 Observe ha M 0 N0 = 0, M 0 A 0 ( X C 0 ) = 0, and M 0 ( N 0 ( X C 0 )) = N 0 ( X C 0 ). Therefore, M 0 = W 0 W 0 / W 0 2 where W 0 = N 0 ( X C 0 ). Similarly, M 1 = W 1 W 1 / W 1 2 where W 1 = N 1 ( Y C 1 ). The previous displayed equaions are equivalen o W 0( X Y ) = 0 and W 1( X Y ) = 0. The poins X = (x 0, x 1, x 2 ) and Y = (y 0, y 1, y 2 ) ha aain minimum disance beween he wo ellipses are 2
3 soluions o he six quadraic equaions in six unknowns, p 0 (x 0, x 1, x 2 ) = N 0 ( X C 0 ) = 0, p 1 (x 0, x 1, x 2 ) = ( X C 0 ) A 0 ( X C 0 ) = 1, p 2 (x 0, x 1, x 2, y 0, y 1, y 2 ) = ( X Y ) N 0 ( X C 0 ) = 0, q 0 (y 0, y 1, y 2 ) = N 1 ( Y C 1 ) = 0, q 1 (y 0, y 1, y 2 ) = ( Y C 1 ) A 1 ( Y C 1 ) = 1, q 2 (x 0, x 1, x 2, y 0, y 1, y 2 ) = ( X Y ) N 1 ( Y C 1 ) = 0. On a compuer algebra sysem ha suppors he resulan operaion for eliminaing polynomial variables, he following se of operaions leads o a polynomial in one variable. Le Resulan[P, Q, z] denoe he resulan of polynomials P and Q where he variable z is eliminaed, r 0 (x 0, x 1, y 0, y 1, y 2 ) = Resulan[p 0, p 2, x 2 ] r 1 (x 0, x 1 ) = Resulan[p 1, p 2, x 2 ] r 2 (x 0, x 1, y 0, y 1 ) = Resulan[r 0, q 2, y 2 ] s 0 (x 0, x 1, x 2, y 0, y 1 ) = Resulan[q 0, q 2, y 2 ] s 1 (y 0, y 1 ) = Resulan[q 1, q 2, y 2 ] s 2 (x 0, x 1, y 0, y 1 ) = Resulan[s 0, p 2, x 2 ] r 3 (x 0, y 0, x 1 ) = Resulan[r 2, r 1, x 1 ] r 4 (x 0, y 0 ) = Resulan[r 3, s 1, y 1 ] s 3 (x 0, x 1, y 0 ) = Resulan[s 2, s 1, y 1 ] s 4 (x 0, y 0 ) = Resulan[s 3, r 1, x 1 ] φ(x 0 ) = Resulan[r 4, s 4, y 0 ] For wo circles, he degree of φ is 8. For a circle and an ellipse, he degree of φ is 12. For wo ellipses, he degree of φ is 16. 3
4 3 Soluion using Trigonomeric Approach Le he wo ellipses be X = C 0 + a 0 cos(θ) U 0 + b 0 sin(θ) V 0 Y = C 1 + a 1 cos(φ) U 1 + b 1 sin(φ) V 1 for θ [0, 2π) and φ [0, 2π). The squared disance beween any wo poins on he ellipses is F (θ, φ) = X(θ) Y (φ) 2. The problem is o minimize F (θ, φ). Define c 0 = cos(θ), s 0 = sin(θ), c 1 = cos(φ), and s 1 = sin(φ). Compue derivaives, F θ = ( X(θ) Y (φ)) X (θ) and F φ = ( X(θ) Y (φ)) Y (φ). Seing hese equal o zero leads o he wo polynomial equaions in c 0, s 0, c 1, and s 1. The wo polynomial consrains for he sines and cosines are also lised. p 0 = (a 2 0 b 2 0)s 0 c 0 + a 0 (α 00 + α 01 s 1 + α 02 c 1 )s 0 + b 0 (β 00 + β 01 s 1 + β 02 c 1 )c 0 = 0 p 1 = (a 2 1 b 2 1)s 1 c 1 + a 1 (α 10 + α 11 s 0 + α 12 c 0 )s 1 + b 1 (β 10 + β 11 s 0 + β 12 c 0 )c 1 = 0 q 0 = s c = 0 q 1 = s c = 0 This is a sysem of four quadraic polynomial equaions in four unknowns and can be solved wih resulans: r 0 (s 0, s 1, c 1 ) = Resulan[p 0, q 0, c 0 ] r 1 (s 0, s 1, c 0 ) = Resulan[p 1, q 1, c 1 ] r 2 (s 0, s 1 ) = Resulan[r 0, q 1, c 1 ] r 3 (s 0, s 1 ) = Resulan[r 1, q 0, c 0 ] φ(s 0 ) = Resulan[r 2, r 3, s 1 ] Alernaively, we can use he simple naure of q 0 and q 1 o do some of he eliminaion. Le p 0 = α 0 s 0 + β 0 c 0 + γ 0 s 0 c 0 where α 0 and β 0 are linear polynomials in s 1 and c 1. Similarly, p 1 = α 1 s 1 + β 1 c 1 + γ 1 s 1 c 1 where α 1 and β 1 are linear polynomials in s 0 and c 0. Solving for c 0 in p 0 = 0 and c 1 in p 1 = 0, squaring, and using he q i consrains leads o r 0 = (1 s 2 0)(γ 0 s 0 + β 0 ) 2 α 2 0s 2 0 = 0 r 1 = (1 s 2 1)(γ 1 s 1 + β 1 ) 2 α 2 1s 2 1 = 0 Using he q i consrains, wrie r i = r i 0 + r i 1c 1 i, i = 0, 1, where he r ij are polynomials in s 0 and s 1. The erms r i 0 are degree 4 and he erms r i1 is degree 3. Solving for c 0 in r 1 = 0 and c 1 in r 0 = 0, squaring, and 4
5 using he q i consrains leads o w 0 = (1 s 2 1)r 2 01 r 2 00 = 8 j=0 w 0js j 0 = 0 w 1 = (1 s 2 0)r 2 11 r 2 10 = 4 j=0 w 1js j 0 = 0 The coefficiens w ij are polynomials in s 1. The degrees of w 00 hrough w 08 respecively are 4, 3, 4, 3, 4, 3, 2, 1, and 0. The degree of w 1j is 8 j. Toal degree for each of w i is 8. The final eliminaion can be compued using a Bezou deerminan, φ(s 1 ) = de[e ij ], where he underlying marix is 8x8 and he enry is e ij = min(8,17 i j) k=max(9 i,9 j) v k,17 i j k where v i,j = w 0i w 1j w 0j w 1i. If he i or j index is ou of range in he w erms, hen he erm is assumed o be zero. The soluions o φ = 0 are he candidae poins for s 1. For each s 1, wo c 1 values are compued using s c 2 1 = 1. For each s 1, he roos of he polynomial w 1 (s 0 ) are compued. For each s 0, wo c 0 values are compued using s c 2 0 = 1. Ou of all such candidaes, X Y 2 can be compued and he minimum value is seleced. 4 Numerical Soluion Neiher algebraic mehod above seems reasonable. Each looks very slow o compue and you have he usual numerical problems wih polynomials of large degree. I have no implemened he following, bu my guess is i is an alernaive o consider. Implemen a disance calculaor for poin o ellipse (in 3D). This is a funcion of a single parameer, say F (θ) for θ [0, 2π]. Use a numerical minimizer ha does no require derivaive calculaion (Powell s mehod for example) and minimize F on he inerval [0, 2π]. The scheme is ieraive and hopefully converges rapidly o he soluion. 5
dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page
Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,
More informationAP Calculus BC Chapter 10 Part 1 AP Exam Problems
AP Calculus BC Chaper Par AP Eam Problems All problems are NO CALCULATOR unless oherwise indicaed Parameric Curves and Derivaives In he y plane, he graph of he parameric equaions = 5 + and y= for, is a
More informationMATH 128A, SUMMER 2009, FINAL EXAM SOLUTION
MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange
More information!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)
"#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5
More informationTHE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI
THE 2-BODY PROBLEM ROBERT J. VANDERBEI ABSTRACT. In his shor noe, we show ha a pair of ellipses wih a common focus is a soluion o he 2-body problem. INTRODUCTION. Solving he 2-body problem from scrach
More information3, 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 informationSolutions 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 informationMath 334 Fall 2011 Homework 11 Solutions
Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then
More informationChapter 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 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 informationParametrics and Vectors (BC Only)
Paramerics and Vecors (BC Only) The following relaionships should be learned and memorized. The paricle s posiion vecor is r() x(), y(). The velociy vecor is v(),. The speed is he magniude of he velociy
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 information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
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 informationTHE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2.
THE BERNOULLI NUMBERS The Bernoulli numbers are defined here by he exponenial generaing funcion ( e The firs one is easy o compue: (2 and (3 B 0 lim 0 e lim, 0 e ( d B lim 0 d e +e e lim 0 (e 2 lim 0 2(e
More informationa. Show that these lines intersect by finding the point of intersection. b. Find an equation for the plane containing these lines.
Mah A Final Eam Problems for onsideraion. Show all work for credi. Be sure o show wha you know. Given poins A(,,, B(,,, (,, 4 and (,,, find he volume of he parallelepiped wih adjacen edges AB, A, and A.
More informationChapter 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 information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
More informationSolutionbank Edexcel AS and A Level Modular Mathematics
Page of 4 Soluionbank Edexcel AS and A Level Modular Mahemaics Exercise A, Quesion Quesion: Skech he graphs of (a) y = e x + (b) y = 4e x (c) y = e x 3 (d) y = 4 e x (e) y = 6 + 0e x (f) y = 00e x + 0
More informationSystem 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 informationSolutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore
Soluions of Sample Problems for Third In-Class Exam Mah 6, Spring, Professor David Levermore Compue he Laplace ransform of f e from is definiion Soluion The definiion of he Laplace ransform gives L[f]s
More informationMath 116 Practice for Exam 2
Mah 6 Pracice for Exam Generaed Ocober 3, 7 Name: SOLUTIONS Insrucor: Secion Number:. This exam has 5 quesions. Noe ha he problems are no of equal difficuly, so you may wan o skip over and reurn o a problem
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 informationMath 315: Linear Algebra Solutions to Assignment 6
Mah 35: Linear Algebra s o Assignmen 6 # Which of he following ses of vecors are bases for R 2? {2,, 3, }, {4,, 7, 8}, {,,, 3}, {3, 9, 4, 2}. Explain your answer. To generae he whole R 2, wo linearly independen
More informationThe 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 informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
More information10.1 EXERCISES. y 2 t 2. y 1 t y t 3. y e
66 CHAPTER PARAMETRIC EQUATINS AND PLAR CRDINATES SLUTIN We use a graphing device o produce he graphs for he cases a,,.5,.,,.5,, and shown in Figure 7. Noice ha all of hese curves (ecep he case a ) have
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More informationGENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT
Inerna J Mah & Mah Sci Vol 4, No 7 000) 48 49 S0670000970 Hindawi Publishing Corp GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT RUMEN L MISHKOV Received
More informationChapter 7: Solving Trig Equations
Haberman MTH Secion I: The Trigonomeric Funcions Chaper 7: Solving Trig Equaions Le s sar by solving a couple of equaions ha involve he sine funcion EXAMPLE a: Solve he equaion sin( ) The inverse funcions
More informationConcourse 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 informationMATH 2050 Assignment 9 Winter Do not need to hand in. 1. Find the determinant by reducing to triangular form for the following matrices.
MATH 2050 Assignmen 9 Winer 206 Do no need o hand in Noe ha he final exam also covers maerial afer HW8, including, for insance, calculaing deerminan by row operaions, eigenvalues and eigenvecors, similariy
More informationt 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 informationChapters 6 & 7: Trigonometric Functions of Angles and Real Numbers. Divide both Sides by 180
Algebra Chapers & : Trigonomeric Funcions of Angles and Real Numbers Chapers & : Trigonomeric Funcions of Angles and Real Numbers - Angle Measures Radians: - a uni (rad o measure he size of an angle. rad
More informationFrom 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 information1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.
. Use Simpson s rule wih n 4 o esimae an () +. Soluion: Since we are using 4 seps, 4 Thus we have [ ( ) f() + 4f + f() + 4f 3 [ + 4 4 6 5 + + 4 4 3 + ] 5 [ + 6 6 5 + + 6 3 + ]. 5. Our funcion is f() +.
More informationBasilio 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 informationInnova Junior College H2 Mathematics JC2 Preliminary Examinations Paper 2 Solutions 0 (*)
Soluion 3 x 4x3 x 3 x 0 4x3 x 4x3 x 4x3 x 4x3 x x 3x 3 4x3 x Innova Junior College H Mahemaics JC Preliminary Examinaions Paper Soluions 3x 3 4x 3x 0 4x 3 4x 3 0 (*) 0 0 + + + - 3 3 4 3 3 3 3 Hence x or
More informationThe equation to any straight line can be expressed in the form:
Sring Graphs Par 1 Answers 1 TI-Nspire Invesigaion Suden min Aims Deermine a series of equaions of sraigh lines o form a paern similar o ha formed by he cables on he Jerusalem Chords Bridge. Deermine he
More informationCERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS
SARAJEVO JOURNAL OF MATHEMATICS Vol.10 (22 (2014, 67 76 DOI: 10.5644/SJM.10.1.09 CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS ALMA OMERSPAHIĆ AND VAHIDIN HADŽIABDIĆ Absrac. This paper presens sufficien
More information2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes
Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion
More informationAP Calculus BC - Parametric equations and vectors Chapter 9- AP Exam Problems solutions
AP Calculus BC - Parameric equaions and vecors Chaper 9- AP Exam Problems soluions. A 5 and 5. B A, 4 + 8. C A, 4 + 4 8 ; he poin a is (,). y + ( x ) x + 4 4. e + e D A, slope.5 6 e e e 5. A d hus d d
More informationMath 1. Two-Hours Exam December 10, 2017.
Mah. Two-Hours Exam December, 7. JE/JKL..7 Problem resar;wih(linearalgebra): Given he inhomogeneous linear sysem of equaions lign:=x-*x+3*x3=a^+*a-3; lign x x 3 x3 = a a 3 lign:=x+*x-*x3=a^+3; lign x x
More information2.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 informationRoller-Coaster Coordinate System
Winer 200 MECH 220: Mechanics 2 Roller-Coaser Coordinae Sysem Imagine you are riding on a roller-coaer in which he rack goes up and down, wiss and urns. Your velociy and acceleraion will change (quie abruply),
More informationChapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis
Chaper EEE83 EEE3 Chaper # EEE83 EEE3 Linear Conroller Design and Sae Space Analysis Ordinary Differenial Equaions.... Inroducion.... Firs Order ODEs... 3. Second Order ODEs... 7 3. General Maerial...
More informationMATH 5720: Gradient Methods Hung Phan, UMass Lowell October 4, 2018
MATH 5720: Gradien Mehods Hung Phan, UMass Lowell Ocober 4, 208 Descen Direcion Mehods Consider he problem min { f(x) x R n}. The general descen direcions mehod is x k+ = x k + k d k where x k is he curren
More informationLinear Dynamic Models
Linear Dnamic Models and Forecasing Reference aricle: Ineracions beween he muliplier analsis and he principle of acceleraion Ouline. The sae space ssem as an approach o working wih ssems of difference
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 information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationSolutions to Assignment 1
MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we
More informationAn Introduction to Malliavin calculus and its applications
An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214
More informationOscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson
PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,
More information( ) 2. Review Exercise 2. cos θ 2 3 = = 2 tan. cos. 2 x = = x a. Since π π, = 2. sin = = 2+ = = cotx. 2 sin θ 2+
Review Eercise sin 5 cos sin an cos 5 5 an 5 9 co 0 a sinθ 6 + 4 6 + sin θ 4 6+ + 6 + 4 cos θ sin θ + 4 4 sin θ + an θ cos θ ( ) + + + + Since π π, < θ < anθ should be negaive. anθ ( + ) Pearson Educaion
More informationLaplace Transforms. Examples. Is this equation differential? y 2 2y + 1 = 0, y 2 2y + 1 = 0, (y ) 2 2y + 1 = cos x,
Laplace Transforms Definiion. An ordinary differenial equaion is an equaion ha conains one or several derivaives of an unknown funcion which we call y and which we wan o deermine from he equaion. The equaion
More informationSimulation-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 informationON THE DEGREES OF RATIONAL KNOTS
ON THE DEGREES OF RATIONAL KNOTS DONOVAN MCFERON, ALEXANDRA ZUSER Absrac. In his paper, we explore he issue of minimizing he degrees on raional knos. We se a bound on hese degrees using Bézou s heorem,
More informationBEng (Hons) Telecommunications. Examinations for / Semester 2
BEng (Hons) Telecommunicaions Cohor: BTEL/14/FT Examinaions for 2015-2016 / Semeser 2 MODULE: ELECTROMAGNETIC THEORY MODULE CODE: ASE2103 Duraion: 2 ½ Hours Insrucions o Candidaes: 1. Answer ALL 4 (FOUR)
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 informationOperators related to the Jacobi setting, for all admissible parameter values
Operaors relaed o he Jacobi seing, for all admissible parameer values Peer Sjögren Universiy of Gohenburg Join work wih A. Nowak and T. Szarek Alba, June 2013 () 1 / 18 Le Pn α,β be he classical Jacobi
More informationMEI Mechanics 1 General motion. Section 1: Using calculus
Soluions o Exercise MEI Mechanics General moion Secion : Using calculus. s 4 v a 6 4 4 When =, v 4 a 6 4 6. (i) When = 0, s = -, so he iniial displacemen = - m. s v 4 When = 0, v = so he iniial velociy
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationLecture 23 Damped Motion
Differenial Equaions (MTH40) Lecure Daped Moion In he previous lecure, we discussed he free haronic oion ha assues no rearding forces acing on he oving ass. However No rearding forces acing on he oving
More informationADDITIONAL MATHEMATICS PAPER 1
000-CE A MATH PAPER HONG KONG EXAMINATIONS AUTHORITY HONG KONG CERTIFICATE OF EDUCATION EXAMINATION 000 ADDITIONAL MATHEMATICS PAPER 8.0 am 0.0 am ( hours This paper mus be answered in English. Answer
More informationAnswers to 1 Homework
Answers o Homework. x + and y x 5 y To eliminae he parameer, solve for x. Subsiue ino y s equaion o ge y x.. x and y, x y x To eliminae he parameer, solve for. Subsiue ino y s equaion o ge x y, x. (Noe:
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 informationLAPLACE TRANSFORM AND TRANSFER FUNCTION
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More informationWeek #13 - Integration by Parts & Numerical Integration Section 7.2
Week #3 - Inegraion by Pars & Numerical Inegraion Secion 7. From Calculus, Single Variable by Hughes-Halle, Gleason, McCallum e. al. Copyrigh 5 by John Wiley & Sons, Inc. This maerial is used by permission
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 informationThe expectation value of the field operator.
The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining
More informationk 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series
Final Review A Puzzle... Consider wo massless springs wih spring consans k 1 and k and he same equilibrium lengh. 1. If hese springs ac on a mass m in parallel, hey would be equivalen o a single spring
More informationApplication 5.4 Defective Eigenvalues and Generalized Eigenvectors
Applicaion 5.4 Defecive Eigenvalues and Generalized Eigenvecors The goal of his applicaion is he soluion of he linear sysems like where he coefficien marix is he exoic 5-by-5 marix x = Ax, (1) 9 11 21
More informationMath 116 Second Midterm March 21, 2016
Mah 6 Second Miderm March, 06 UMID: EXAM SOLUTIONS Iniials: Insrucor: Secion:. Do no open his exam unil you are old o do so.. Do no wrie your name anywhere on his exam. 3. This exam has pages including
More informationUnsteady Flow Problems
School of Mechanical Aerospace and Civil Engineering Unseady Flow Problems T. J. Craf George Begg Building, C41 TPFE MSc CFD-1 Reading: J. Ferziger, M. Peric, Compuaional Mehods for Fluid Dynamics H.K.
More informationThe average rate of change between two points on a function is d t
SM Dae: Secion: Objecive: The average rae of change beween wo poins on a funcion is d. For example, if he funcion ( ) represens he disance in miles ha a car has raveled afer hours, hen finding he slope
More informationElements of Computer Graphics
CS580: Compuer Graphics Min H. Kim KAIST School of Compuing Elemens of Compuer Graphics Geomery Maerial model Ligh Rendering Virual phoography 2 Foundaions of Compuer Graphics A PINHOLE CAMERA IN 3D 3
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 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 informationTracking Adversarial Targets
A. Proofs Proof of Lemma 3. Consider he Bellman equaion λ + V π,l x, a lx, a + V π,l Ax + Ba, πax + Ba. We prove he lemma by showing ha he given quadraic form is he unique soluion of he Bellman equaion.
More informationMathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation
Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion
More informationd 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3
and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or
More informationtranslational component of a rigid motion appear to originate.
!"#$%&'()(!"##$%&'(&')*+&'*',*-&.*'/)*/'0%'"1&.3$013'*'#".&45'/.*1%4*/0$1*4'-$/0$16'*1'/)*/' /)&',*-&.*6'(0/)'*'7$,*4'4&13/)'$7'8996'0%'-$+013'(0/)'+&4$,0/5'' :;6'?':@96'A96'B9>C' *>'D&1&.*4456'()*/'0%'/)&'7$,"%'$7'&E#*1%0$1'7$.'*'-$+013',*-&.*F'
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More information(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)
Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion
More informationDirectional Tubular Surfaces
Inernaional Journal of Algebra, Vol. 9, 015, no. 1, 57-535 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/10.1988/ija.015.5174 Direcional Tubular Surfaces Musafa Dede Deparmen of Mahemaics, Faculy of Ars
More informationAnalyze patterns and relationships. 3. Generate two numerical patterns using AC
envision ah 2.0 5h Grade ah Curriculum Quarer 1 Quarer 2 Quarer 3 Quarer 4 andards: =ajor =upporing =Addiional Firs 30 Day 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 andards: Operaions and Algebraic Thinking
More informationSINUSOIDAL WAVEFORMS
SINUSOIDAL WAVEFORMS The sinusoidal waveform is he only waveform whose shape is no affeced by he response characerisics of R, L, and C elemens. Enzo Paerno CIRCUIT ELEMENTS R [ Ω ] Resisance: Ω: Ohms Georg
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 informationMTH Feburary 2012 Final term PAPER SOLVED TODAY s Paper
MTH401 7 Feburary 01 Final erm PAPER SOLVED TODAY s Paper Toal Quesion: 5 Mcqz: 40 Subjecive quesion: 1 4 q of 5 marks 4 q of 3 marks 4 q of marks Guidelines: Prepare his file as I included all pas papers
More informationFractional 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 informationLet ( α, β 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 informationWavelet Methods for Time Series Analysis. What is a Wavelet? Part I: Introduction to Wavelets and Wavelet Transforms. sines & cosines are big waves
Wavele Mehods for Time Series Analysis Par I: Inroducion o Waveles and Wavele Transforms waveles are analysis ools for ime series and images as a subjec, waveles are relaively new (983 o presen) a synhesis
More informationEE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?
EE 35 Noes Gürdal Arslan CLASS (Secions.-.2) Wha is a signal? In his class, a signal is some funcion of ime and i represens how some physical quaniy changes over some window of ime. Examples: velociy of
More informationLie Derivatives operator vector field flow push back Lie derivative of
Lie Derivaives The Lie derivaive is a mehod of compuing he direcional derivaive of a vecor field wih respec o anoher vecor field We already know how o make sense of a direcional derivaive of real valued
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 information