PHYS 705: Classical Mechanics. Calculus of Variations II


 Alison Fletcher
 4 years ago
 Views:
Transcription
1 1 PHYS 705: Classcal Mechancs Calculus of Varatons II
2 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary value F F y ( x), y '( x); x 1,2,, N y x 1 ( ) N B I F y( x), y'( x); x dx We need to consder N varatons n N dmensons. y (, x) y (0, x) ( x) 1,2,, N A y x 1 ( ) N Agan, we have, as a parameter x C x 2 ( ) are smooth varaton at ( x ) ( x ) 0 A A
3 3 Varatonal Calculus: Generalzaton (no constrant) Everythng proceeds as before, and we get: x B di F d F ( xdx ) 0 d y dxy' xa ( x) If all the varatons are ndependent, ths equaton requres that each coeffcents n the ntegrant to vansh ndependently. Then, we have, y (note: sum over comes from chan rule n takng the dervatve: ) df ( y, y '; x) d F d F 0 1,2,, N y dxy' Ths s the same result (EulerLagrange Eq) as we have before. However
4 4 Varatonal Calculus: Wth Constrant (2D) When we have constrants, the coordnates and hence the varatons wll NOT be ndependent! e.g., the system s constrant to move on a lower dmensonal surface (or ( x) curve) and allowed varatons are coupled n order to satsfy the constrant. To extend our analyss to systems wth constrants, let frst consder the 2D case wth ( y, y ) y ( x) as our dependent varables wth one constrant. To specfy the Holonomc constrant, we have g( y1, y2; x) 0 lnkng the two dependent varables.
5 5 Varatonal Calculus: Wth Constrant (2D) Wth two dep varables, our varatonal equaton can be wrtten as, x B di F d F F d F 1( x) 2( x) dx0 (*) d y1 dx y1' y2 dx y2' x A Smlarly wth, y1(, x) y1(0, x) 1( x) y2(, x) y2(0, x) 2( x)
6 6 Varatonal Calculus: Wth Constrant (2D) The constrant lnks so that these two varatons are NOT ndependent! ( x) to ( x) di d 0 Thus, we can t say that (Eq. *) requres the ndvdual coeffcents wthn the ntegrant to go to zero ndependently as before. x B di F d F F d F 1( x) 2( x) dx0 (*) d y1 dxy1 ' y2 dxy2' x A
7 7 Varatonal Calculus: Wth Constrant (2D) The constrant lnks so that these two varatons are NOT ndependent! ( x) to ( x) di d 0 Thus, we can t say that (Eq. *) requres the ndvdual coeffcents wthn the ntegrant to go to zero ndependently as before. To get around ths, we consder the total varaton of g (total dervatve) wrt : dg g y g y d y d y d 0 (**) NOTE: Ths expresson needs to be dentcally zero snce the NET varaton must vansh and be consstent wth the constrant (everythng must stay on the surface defned by g( y, y ; x) 0.
8 8 Varatonal Calculus: Wth Constrant (2D) From the defnton of the varatons, we have y1 y2 1( x) and 2( x) Puttng these nto Eq. (**), dg g g d y y Rearrangng terms, ths then gves, g y g y2 (Note: the varatons are lnked through the equaton of constrant g.)
9 9 Varatonal Calculus: Wth Constrant (2D) 2 Substtutng nto Eq. (*), we have, di F d F F d F d y dx y y dx y g y A xb g y 1 1dx 1 1' 2 2' x Snce s an ndependent varaton, we can now argue that the coeffcents { } n front of t must vansh. F d F F d F g y y1 dxy1' y2 dxy2' g y2 1 0
10 10 Varatonal Calculus: Wth Constrant (2D) Rewrtng the expresson wth separately on the LHS and RHS, we have the followng equaton: y ( x)and y ( x) Note that these two terms are dentcal wth y y F d F 1 F d F 1 y1 dxy1' g y1 y2 dxy2' g y2 ( x) And, snce these two terms EQUAL to each other and depend on ndependently they can only be a functon of x ONLY! y1( x)and y2( x) And, we wll call ths as yet undetermned functon (multpler): ( x)
11 11 Varatonal Calculus: Wth Constrant (2D) Wth ths multpler, we can rewrte the LHS & RHS as two decoupled equatons: ( x) F d F g ( x) 0 1,2 (1) y dxy' y g( y, y ; x) 0 Together wth, these 3 equatons can be solved for the three unknowns: y ( x), y ( x), and ( x) ( x) s call the Lagrange undetermned multpler
12 12 Varatonal Calculus: Geometrc Vew Pont In our prevous 2D problem, the ntegral B I[ y Fy ( x), y '( x); xdx A s a functonal whch takes n two functons and gves back out a real number and the frst two terms on the LHS of EL equatons I F d F, 1,2 y y dxy' y ( x) and y ( x) s n fact a functonal dervatve of I [y wth respect to the functon y. Smlar to regular dervatves n vector calculus, t gves the rate of change of I wth respect to a change n the functon y. F d F g ( x) 0 y dx y' y
13 13 Varatonal Calculus: Geometrc Vew Pont So, one can thnk of the followng combned vector object as a functonal gradent smlar to a gradent vector n multvarable calculus. Ths gradent wll pont along the drecton of the largest ncrease of I n the space of functons correspondng to y ( x) and y ( x) and normal to the level curves of I. Smlarly, we can make a same assocaton wth the functonal dervatve of g wth respect to I I I[ y, y1 y2 y1( x) and y2( x) g g g g gy [,, y y y y (Prevously, we ddn t make ths dstncton but g s n fact a functonal.)
14 14 Varatonal Calculus: Geometrc Vew Pont In the space of, g y, y 0 y ( x), y ( x) [ I y gy [, y 0  defnes a curve n the space y ( x), y ( x) y 2 [ g y g[ y ponts n the drecton normal to the constrant curve defned by gy [ 1, y2 0 y 1
15 15 Varatonal Calculus: Geometrc Vew Pont In the space of, y ( x), y ( x) [ I y 1, 2 I y y c defne sets of level curves n the space of y ( x), y ( x) y 2 [ I y ponts n the drecton normal to the level curves defned by 1, 2 I y y c y 1
16 16 Varatonal Calculus: Geometrc Vew Pont Then, the EulerLagrange Equaton has a specfc geometrc meanng: g y, y 0 * I[ y F d F g ( x) y dxy' y y 2 * g[ y I[ y ( x) g[ y The soluton at are chosen such that y, y * * y 1 * I[ y * & g[ y are collnear!
17 17 Varatonal Calculus: Geometrc Vew Pont g y, y 0 * I[ y * g[ y Recall that y, * *  The soluton at to the EL s also a statonary pont to. I y y 1, 2  AND, now we are also requrng that y mn, max, or nflecton pt I y, y y 2 the soluton to be consstent wth the gy [, y 0 constrant. y 1 It s a problem n constraned optmzaton.
18 18 Varatonal Calculus: Geometrc Vew Pont y 2 g y, y 0 y 1 * I[ y * g[ y As we vary along the constrant gy [, y 0 curve (red), we wll typcally traverse dfferent level curves of I (blue curves). Only at the locaton where the two curves touch tangentally, the I value wll not change as we move along the constraned soluton. Ths geometrc condton can be compactly wrtten as: I[ y g[ y y * *, y
19 19 Varatonal Calculus: Generalzaton If the system has more than one constrant, let say M: g ( y; x) 0 1,2,, N k 1,2, M k (# dep vars) (# constrants) Each of these M equatons descrbes a (N1) dmensonal surface and the ntersecton of them wll be a (NM) dmensonal surface on whch the dynamcs must lve. There wll also be M normal (vectors) to each of these M (N1)dmensonal surfaces gven by: g [ y k 1,2,, M k
20 20 Varatonal Calculus: Generalzaton g [ y Each of these gradents wll pont n a drecton AWAY from k (normal to) ts prospectve constrant surface. g k ( y ; x) 0 0 g y k g [ y k I[ y In a smlar geometrc fashon, ponts n the drecton away from (normal I[ y to) the level set of and n the drecton of ts greatest rate of change. I[ y I[ y c c 1 2 level set of I[ y [ I y
21 21 Varatonal Calculus: Generalzaton I[ y The geometrc condton for to have a statonary value wth respect to the constrants wll be smlar to our 2D prevous case,.e., the optmal set of s the one wth, Or, equvalently y,, 1 yn I[ y span g [ y I[ y k Ths means that, can be wrtten as a lnear combnaton of the M normal vectors from the M constrant equatons. I[ y g [ y 0 k k (ponts normal to ALL the constran surfaces.) M I[ y g [ y k k k 1 Instead of one, we have M ndependent Lagrange Multplers
22 22 Varatonal Calculus: Generalzaton Puttng ths back nto the orgnal dfferental equaton format, we can wrte ths condton as, M F d F gk k ( x) 0 1, 2,, N y dxy' k 1 y g ( ; ) 0 k 1, 2,, M k y x (for statonarty) (to be on constrants) y ( x ) k ( x) Here, we have N+M unknowns: and And, we have N+M equatons: top (N) and bottom (M) Ths problem n general can be solved!
23 23 Hamlton s Prncple Now, back to mechancs and we wll defne the Acton as, where L = T U s the Lagrangan of the system I 2 1 L q, q, t dt Confguraton Space: The space of the q s. Each pont gves the full confguraton of the system. The moton of the system s a path through ths confguraton space. (Ths s not necessary the real space.) (Note: Here we don t requre all the generalzed coordnates to be necessarly q ' s ndependent. The can be lnked through constrants.) q Hamlton s Prncple: The moton of a system from t 1 to t 2 the acton evaluated along the actual path s statonary. s such that
24 24 Lagrange Equaton of Moton So, we apply our varatonal calculus results to the acton ntegral. We also further assume that the system s monogenc,.e., all forces except forces of constrant are dervable from a potental functon whch can be a functon of q, q,and t The resultant equaton s the Lagrange Equaton of Moton wth N generalzed coordnates (not necessary proper) and M Holonomc constrants: U q, q, t M L d L gk k () t 0 1, 2,, N q dtq k 1 q k 1, 2,, M gk( q; t) 0
25 25 Forces of Constrant Comments: 1. The Lagrange EOM can be formally wrtten as: M L d L gk k() t Q 1,2,, N q dtq k 1 q where the Q are the generalzed forces whch gve the magntudes of the forces needed to produce the ndvdual constrants. q  the generalzed coordnates are NOT necessary ndependent and they are lnked through constrants.  snce the choce of the sgn for k s arbtrary, the drecton of the forces of constrant forces cannot be determned.
26 26 Proper Generalzed Coordnates Comments: 2. If one chooses a set of proper (ndependent) generalzed coordnates n whch the (NM) are no longer lnked through the constrants, then the Lagrange EOM reduces to: q j ' s L d L 0 j 1,2,, N qj dt q j M  In practce, one typcally wll explctly use the constrant equatons to reduce the number of varables to the (NM) proper set of generalzed coordnates.  However, one CAN T solve for the forces of constrant here
Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationLagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013
Lagrange Multplers Monday, 5 September 013 Sometmes t s convenent to use redundant coordnates, and to effect the varaton of the acton consstent wth the constrants va the method of Lagrange undetermned
More informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationMechanics Physics 151
Mechancs Physcs 151 Lecture 3 Lagrange s Equatons (Goldsten Chapter 1) Hamlton s Prncple (Chapter 2) What We Dd Last Tme! Dscussed multpartcle systems! Internal and external forces! Laws of acton and
More information12. The HamiltonJacobi Equation Michael Fowler
1. The HamltonJacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
More informationPoisson brackets and canonical transformations
rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: mgk F 1 x O
More informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
More informationn α 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 informationThree views of mechanics
Three vews of mechancs John Hubbard, n L. Gross s course February 1, 211 1 Introducton A mechancal system s manfold wth a Remannan metrc K : T M R called knetc energy and a functon V : M R called potental
More information10. Canonical Transformations Michael Fowler
10. Canoncal Transformatons Mchael Fowler Pont Transformatons It s clear that Lagrange s equatons are correct for any reasonable choce of parameters labelng the system confguraton. Let s call our frst
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationPHYS 705: Classical Mechanics. HamiltonJacobi Equation
1 PHYS 705: Classcal Mechancs HamltonJacob Equaton HamltonJacob Equaton There s also a very elegant relaton between the Hamltonan Formulaton of Mechancs and Quantum Mechancs. To do that, we need to derve
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationClassical Mechanics Virtual Work & d Alembert s Principle
Classcal Mechancs Vrtual Work & d Alembert s Prncple Dpan Kumar Ghosh UMDAE Centre for Excellence n Basc Scences Kalna, Mumba 400098 August 15, 2016 1 Constrants Moton of a system of partcles s often
More informationLecture 20: Noether s Theorem
Lecture 20: Noether s Theorem In our revew of Newtonan Mechancs, we were remnded that some quanttes (energy, lnear momentum, and angular momentum) are conserved That s, they are constant f no external
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fntedmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationBezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0
Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the
More informationMechanics Physics 151
Mechancs Physcs 5 Lecture 0 Canoncal Transformatons (Chapter 9) What We Dd Last Tme Hamlton s Prncple n the Hamltonan formalsm Dervaton was smple δi δ p H(, p, t) = 0 Adonal endpont constrants δ t ( )
More informationLAGRANGIAN MECHANICS
LAGRANGIAN MECHANICS Generalzed Coordnates State of system of N partcles (Newtonan vew): PE, KE, Momentum, L calculated from m, r, ṙ Subscrpt covers: 1) partcles N 2) dmensons 2, 3, etc. PE U r = U x 1,
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the faketest data; fxed
More informationNMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582
NMT EE 589 & UNM ME 48/58 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 48/58 7. Robot Dynamcs 7.5 The Equatons of Moton Gven that we wsh to fnd the path q(t (n jont space) whch mnmzes the energy
More informationIntegrals and Invariants of EulerLagrange Equations
Lecture 16 Integrals and Invarants of EulerLagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3d rotaton around the xaxs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3d rotaton around the yaxs:
More informationSolutions to exam in SF1811 Optimization, Jan 14, 2015
Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 OO 4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 OO 5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable
More informationLagrangian Field Theory
Lagrangan Feld Theory Adam Lott PHY 391 Aprl 6, 017 1 Introducton Ths paper s a summary of Chapter of Mandl and Shaw s Quantum Feld Theory [1]. The frst thng to do s to fx the notaton. For the most part,
More informationCSci 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 informationprinceton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg
prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there
More informationTHEOREMS OF QUANTUM MECHANICS
THEOREMS OF QUANTUM MECHANICS In order to develop methods to treat manyelectron systems (atoms & molecules), many of the theorems of quantum mechancs are useful. Useful Notaton The matrx element A mn
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationLectures  Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures  Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1s tme nterval. The velocty of the partcle
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationClassical Field Theory
Classcal Feld Theory Before we embark on quantzng an nteractng theory, we wll take a dverson nto classcal feld theory and classcal perturbaton theory and see how far we can get. The reader s expected to
More informationLecture 10 Support Vector Machines. Oct
Lecture 10 Support Vector Machnes Oct  202008 Lnear Separators Whch of the lnear separators s optmal? Concept of Margn Recall that n Perceptron, we learned that the convergence rate of the Perceptron
More informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
More informationSalmon: Lectures on partial differential equations. Consider the general linear, secondorder PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of secondorder equatons There are general methods for classfyng hgherorder partal dfferental equatons. One s very general (applyng even to
More informationAssortment Optimization under MNL
Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenuemaxmzng assortment of products to offer when the prces of products are fxed.
More informationChapter Eight. Review and Summary. Two methods in solid mechanics  vectorial methods and energy methods or variational methods
Chapter Eght Energy Method 8. Introducton 8. Stran energy expressons 8.3 Prncpal of statonary potental energy; several degrees of freedom  Castglano s frst theorem  Examples 8.4 Prncpal of statonary
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationPhysics 106a, Caltech 11 October, Lecture 4: Constraints, Virtual Work, etc. Constraints
Physcs 106a, Caltech 11 October, 2018 Lecture 4: Constrants, Vrtual Work, etc. Many, f not all, dynamcal problems we want to solve are constraned: not all of the possble 3 coordnates for M partcles (or
More informationIterative General Dynamic Model for SerialLink Manipulators
EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for SeralLnk Manpulators In ths set of notes, we are gong to develop a method for computng a general
More informationWhich Separator? Spring 1
Whch Separator? 6.034  Sprng 1 Whch Separator? Mamze the margn to closest ponts 6.034  Sprng Whch Separator? Mamze the margn to closest ponts 6.034  Sprng 3 Margn of a pont " # y (w $ + b) proportonal
More informationFrom BiotSavart Law to Divergence of B (1)
From BotSavart Law to Dvergence of B (1) Let s prove that BotSavart gves us B (r ) = 0 for an arbtrary current densty. Frst take the dvergence of both sdes of BotSavart. The dervatve s wth respect to
More information14 Lagrange Multipliers
Lagrange Multplers 14 Lagrange Multplers The Method of Lagrange Multplers s a powerful technque for constraned optmzaton. Whle t has applcatons far beyond machne learnng t was orgnally developed to solve
More information3.1 Expectation of Functions of Several Random Variables. )' be a kdimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationwhere the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the snglepartcle potental and V nt
Physcs 543 Quantum Mechancs II Fall 998 HartreeFock and the Selfconsstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 16 () () () (v) (v) Overall Mass Balance Momentum
More informationIndeterminate pinjointed frames (trusses)
Indetermnate pnjonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationSupport Vector Machines CS434
Support Vector Machnes CS434 Lnear Separators Many lnear separators exst that perfectly classfy all tranng examples Whch of the lnear separators s the best? + + + + + + + + + Intuton of Margn Consder ponts
More informationQuantum Mechanics for Scientists and Engineers. David Miller
Quantum Mechancs for Scentsts and Engneers Davd Mller Types of lnear operators Types of lnear operators Blnear expanson of operators Blnear expanson of lnear operators We know that we can expand functons
More information9 Characteristic classes
THEODORE VORONOV DIFFERENTIAL GEOMETRY. Sprng 2009 [under constructon] 9 Characterstc classes 9.1 The frst Chern class of a lne bundle Consder a complex vector bundle E B of rank p. We shall construct
More informationFeature Selection: Part 1
CSE 546: Machne Learnng Lecture 5 Feature Selecton: Part 1 Instructor: Sham Kakade 1 Regresson n the hgh dmensonal settng How do we learn when the number of features d s greater than the sample sze n?
More informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
More informationLECTURE 9 CANONICAL CORRELATION ANALYSIS
LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of
More informationTensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q
For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals
More informationAnalytical classical dynamics
Analytcal classcal ynamcs by Youun Hu Insttute of plasma physcs, Chnese Acaemy of Scences Emal: yhu@pp.cas.cn Abstract These notes were ntally wrtten when I rea tzpatrck s book[] an were later revse to
More informationAPPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14
APPROXIMAE PRICES OF BASKE AND ASIAN OPIONS DUPON OLIVIER Prema 14 Contents Introducton 1 1. Framewor 1 1.1. Baset optons 1.. Asan optons. Computng the prce 3. Lower bound 3.1. Closed formula for the prce
More informationCalculus of Variations Basics
Chapter 1 Calculus of Varatons Bascs 1.1 Varaton of a General Functonal In ths chapter, we derve the general formula for the varaton of a functonal of the form J [y 1,y 2,,y n ] F x,y 1,y 2,,y n,y 1,y
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationThe Geometry of Logit and Probit
The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationMechanics Physics 151
Mechancs Physcs 5 Lecture 7 Specal Relatvty (Chapter 7) What We Dd Last Tme Worked on relatvstc knematcs Essental tool for epermental physcs Basc technques are easy: Defne all 4 vectors Calculate com
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationCHAPTER 5: Lie Differentiation and Angular Momentum
CHAPTER 5: Le Dfferentaton and Angular Momentum Jose G. Vargas 1 Le dfferentaton Kähler s theory of angular momentum s a specalzaton of hs approach to Le dfferentaton. We could deal wth the former drectly,
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton
More information2.3 Nilpotent endomorphisms
s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms
More information(2mn, m 2 n 2, m 2 + n 2 )
MATH 16T Homewk Solutons 1. Recall that a natural number n N s a perfect square f n = m f some m N. a) Let n = p α even f = 1,,..., k. be the prme factzaton of some n. Prove that n s a perfect square f
More informationRigid body simulation
Rgd bod smulaton Rgd bod smulaton Once we consder an object wth spacal etent, partcle sstem smulaton s no longer suffcent Problems Problems Unconstraned sstem rotatonal moton torques and angular momentum
More informationThe Second AntiMathima on Game Theory
The Second AntMathma on Game Theory Ath. Kehagas December 1 2006 1 Introducton In ths note we wll examne the noton of game equlbrum for three types of games 1. 2player 2acton zerosum games 2. 2player
More information6. Hamilton s Equations
6. Hamlton s Equatons Mchael Fowler A Dynamcal System s Path n Confguraton Sace and n State Sace The story so far: For a mechancal system wth n degrees of freedom, the satal confguraton at some nstant
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5*  NOTE (Summary) ECON 5*  NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationWeek 5: Neural Networks
Week 5: Neural Networks Instructor: Sergey Levne Neural Networks Summary In the prevous lecture, we saw how we can construct neural networks by extendng logstc regresson. Neural networks consst of multple
More informationCHEMICAL REACTIONS AND DIFFUSION
CHEMICAL REACTIONS AND DIFFUSION A.K.A. NETWORK THERMODYNAMICS BACKGROUND Classcal thermodynamcs descrbes equlbrum states. Nonequlbrum thermodynamcs descrbes steady states. Network thermodynamcs descrbes
More information4DVAR, according to the name, is a fourdimensional variational method.
4DVaratonal Data Assmlaton (4DVar) 4DVAR, accordng to the name, s a fourdmensonal varatonal method. 4DVar s actually a drect generalzaton of 3DVar to handle observatons that are dstrbuted n tme. The
More informationRobert Eisberg Second edition CH 09 Multielectron atoms ground states and xray excitations
Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and xray exctatons 901 By gong through the procedure ndcated n the text, develop the tmendependent Schroednger equaton
More informationADVANCED MACHINE LEARNING ADVANCED MACHINE LEARNING
1 ADVANCED ACHINE LEARNING ADVANCED ACHINE LEARNING Nonlnear regresson technques 2 ADVANCED ACHINE LEARNING Regresson: Prncple N ap Ndm. nput x to a contnuous output y. Learn a functon of the type: N
More informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationSolution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method
Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another wellknown teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222
More informationAERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTINGLINE THEORY
LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTINGLINE THEORY The BotSavart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot Savart formula
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationk t+1 + c t A t k t, t=0
Macro II (UC3M, MA/PhD Econ) Professor: Matthas Kredler Fnal Exam 6 May 208 You have 50 mnutes to complete the exam There are 80 ponts n total The exam has 4 pages If somethng n the queston s unclear,
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationAffine transformations and convexity
Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/
More informationDECOUPLING THEORY HW2
8.8 DECOUPLIG THEORY HW2 DOGHAO WAG DATE:OCT. 3 207 Problem We shall start by reformulatng the problem. Denote by δ S n the delta functon that s evenly dstrbuted at the n ) dmensonal unt sphere. As a temporal
More informationA particle in a state of uniform motion remain in that state of motion unless acted upon by external force.
The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,
More informationSpinrotation coupling of the angularly accelerated rigid body
Spnrotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 Emal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More information8.592J: Solutions for Assignment 7 Spring 2005
8.59J: Solutons for Assgnment 7 Sprng 5 Problem 1 (a) A flament of length l can be created by addton of a monomer to one of length l 1 (at rate a) or removal of a monomer from a flament of length l + 1
More information