Simplified Robotics Joint-Space Trajectory Generation with Via Point using a Single Polynomial
|
|
- Leonard Greer
- 5 years ago
- Views:
Transcription
1 R.L. Willims II, 13, Simplified Robotics Joint-Spce Trjectory Genertion with i Point using Single Polynomil, Hindwi Journl of Robotics, Article ID , 6 pges, Mrch, doi: /13/ Simplified Robotics Joint-Spce Trjectory Genertion with i Point using Single Polynomil Robert L. Willims II Mechnicl Engineering, Ohio University Athens, Ohio, USA willir4@ohio.edu ABSTRACT This rticle presents novel fourth- nd sixth-order polynomils to solve the problem of jointspce trjectory genertion with vi point. These new polynomils use single function rther thn two polynomil functions mtched t the vi point s in previous methods. The problem of infinite spikes in jerk is lso ddressed. KEYWORDS Robot joint-spce trjectory genertion, polynomil, vi point, finite jerk.
2 1. INTRODUCTION Joint-spce trjectory genertion is in common usge in robotics to provide smooth, continuous motion from one set of n joint ngles to nother, for instnce for moving between two distinct Crtesin poses for which the inverse pose solution hs yielded two distinct sets of n joint ngles. The joint-spce trjectory genertion occurs t runtime for ll n joints independently but simultneously. There is n entire body of literture devoted to trjectory genertion (k motion plnning nd pth plnning) t the joint level. Pul nd Zong (1984) were mong the first to suggest the use of polynomils for robot trjectory genertion. A common joint-spce trjectory genertion method (linerly-chnging joint velocity using strting nd ending prbolic blends) is identiclly presented mny uthors (Koivo 1989, Prkin 1991, Scivicco nd Sicilino 1996, Crig 5, Spong et l. 5, nd Jzr 1). Despite the widespred populrity of this method, it suffers from infinite spikes in jerk (the derivtive of ccelertion) nd requires three seprte functions insted of one. An identicl thirdorder polynomil joint-spce trjectory genertion pproch is lso presented by mny uthors (Koivo 1989, Scivicco nd Sicilino 1996, Crig 5, Spong et l. 5, nd Jzr 1). urther, n identicl fifth-order polynomil joint-spce trjectory genertion pproch is presented by mny uthors (Koivo 1989, Scivicco nd Sicilino 1996, Crig 5, nd Jzr 1). rom these uthor lists it my pper tht Koivo ws the first to present these methods, when in fct they were lredy presented in Crig s first edition in u et l. (1987) deprt from these stndrd methods, suggesting initil, intermedite, nd finl polynomils of order 4-3-4, 3-5-3, or 5 third-order polynomils, for single joint motion. These re by fr the most (unnecessrily) complicted methods, nd re presented without justifiction or comprison with simpler methods. or deling with vi point in which the robot need not stop t the vi point (such s for obstcle voidnce), Crig (5) suggests mtching two thirdorder polynomils. Apprently lone mongst ll of the mjor robotics textbook uthors, Angeles (7, Section 6.5) derives seventh-order polynomil to fit two vi points nd ensure finite joint jerk t the strt
3 3 nd end of motion. Angeles (7, Section 6.6) lso presents pproximtion of this seventh-order polynomil with cubic spline nd discusses the ssocited errors. An erly version of the current rticle ws presented by the uthor t conference (Willims, 11). The originl contribution of this rticle is single sixth-order polynomil to provide smooth, continuous joint motion through vi point. Previous pproches hve mtched two third-order polynomils t the vi point. We lso present single fourth-order polynomil to ccomplish the sme purpose. However, this cse should not be used since it leds to discontinuous ccelertions, leding to infinite spikes in jerk, which is uncceptble for relible, smooth, long-life robotic systems. This infinite jerk problem, prevlent in ll of the robotics textbooks, hs lso been pointed out by Mcfrlne nd Croft (3), who present jerk-bounded trjectories nd Gosselin nd Hdj-Messoud (1993) nd Petrinec nd Kovcic (7) who go one step further to ensure continuous, not just bounded, jerk.
4 4. JOINT-SPACE TRAJECTORY GENERATION WITH A IA POINT Stndrd joint-spce trjectory genertion ssumes two sets of n discrete joint prmeters re known (ngles for revolute joints, lengths for prismtic joints) nd it is required to move smoothly in joint spce from one set to the next. or instnce, given required Crtesin poses X S nd X, inverse pose kinemtics clcultes the required n joint vlue sets for chieving ech of these; cll them S nd. Note subscript S stnds for strt (or initil) nd subscript stnds for finish (or finl). The stndrd joint-spce trjectory method then moves smoothly from S nd on ll n joints independently but simultneously. Polynomils re nturl choices for providing smooth, continuous motion, with some level of continuous derivtives. In mny robotics motion plnning problems, the robot must pss through intermedite point(s) between the strt nd finish poses, such s for obstcle voidnce. The joint rtes nd ccelertions need not go to zero t these so-clled vi points, but they must be mtched between functions meeting t the vi point(s)..1 Two Third-Order Polynomils Crig (5) suggested the use of two third-order polynomils meeting t the vi point. This entire Section.1 is summry of Crig s work (with n originl exmple by the current uthor). Two third-order polynomils will provide smooth motion with continuous position nd velocity nd zero velocity t the strt nd end. The two third-order polynomils re: () t t t t 1() t 313t 1t11 () t 6 t () t t t t () t 33t t1 () t 6 t (1) time points: These two polynomils require eight constrints. our constrints come from the initil nd finl
5 1() S () 1 ( t) ( t ) 5 () It is convenient to shift the time xis so tht the (t) polynomil hs zero time strting t the vi time. t 1 is the end time of the first rnge (t 1 = t ) nd t is the reltive end time of the second rnge (t = t t ). We need four more constrints. We force the first polynomil to end t the vi ngle nd the second polynomil to strt t the vi ngle. We lso ensure the velocities nd ccelertions mtch t the vi point (they need not go to zero). This will ensure tht the jerk stys finite during this vi point trnsition between the two polynomils. The four dditionl constrints re: ( t ) 1 1 () 1( t1) () ( t ) () 1 1 (3) The constrints yield eight liner equtions in the eight unknowns (two third-order polynomils, four unknown coefficients ech). Three of the unknown polynomil coefficients re found immeditely, from the initil nd vi time constrints: 1 = S, 11 =, nd =. The simplified 5x5 mtrix/vector eqution to solve for the remining three unknowns is: 1 t1 1 S t 1 1 t t 13 t 1 t 3t 1 t1 3t1 1 6t 1 3 (4) Solve (4) for the remining five unknowns for use in (1). or the specil cse of t 1 t t t T, the nlyticl solution is 1 = S, 11 =, nd =, nd:
6 6 S 1 S T 3 5 4T 3 S 4T 3 3 T 5 3 4T S S 3 3 (5) Exmple Two Third-Order Polynomils with i Point Use two third-order polynomils for smooth joint spce trjectory genertion, plus motion through vi point, for one joint. Given 3, 18, 1, t 1.5, t 3 ( t 1 t 1.5 ) sec, find 1 (t) nd (t): S t t t 1( t) 6.67t 34t 1( t) t34 ( t) ( ) t t t t ( t) t 8t45 ( t) t8 ( t) ( ) (6) deg units re used throughout this exmple. These results re plotted in igure d dd ddd time (sec) igure 1. Two Third-Order Polynomils with i Point
7 Note tht the joint ngle psses through the vi point 18 nd keeps going for brief time this is becuse the velocity is still positive beyond the first 1.5 sec. The mximum ngle is 7 MAX 183.9, occurring t t = 1.68 sec. At MAX the ngulr velocity goes to zero (the slope of the ngle is zero t tht point since the ngle is chnging direction). Also note tht the velocity term in (t), the t coefficient, is non-zero since the velocity does not need to go to zero t this vi point. urther note tht the velocity nd ccelertion hve been successfully mtched t the 1.5 sec vi point trnsition s required. Note tht the jerk is not mtched t the trnsition (the slope of the two polynomil ccelertions re the sme mgnitude but different signs), but jerk remins finite since the ccelertion is mtched. The jerk hs infinite spikes t the strt nd end, which is common in the robotics books, but is uncceptble.. Two ifth-order Polynomils The two third-order polynomils presented bove, mtched t vi point, pper in mny robotics textbooks (e.g. Crig, 5). However, this pproch suffers from discontinuous ccelertion functions t the initil nd finl time, leding to infinite spikes in jerk t these points in time. This is uncceptble for the repeted motions of mechnicl systems, due to uncceptble wer, noise, nd dynmic excittion. To fix this, we could mtch two fifth-order polynomils t the vi point, specifying zero ccelertion t the initil nd finl times (not t the vi time insted mtching velocities nd ccelertions between the two fifth-order polynomil functions t the vi time). However, this will not be presented becuse we cme up with better method, single polynomil vs. two mtched polynomils this originl work is presented in the next two subsections.
8 8.3 Single ourth-order Polynomil The first originl contribution of this rticle is presented in this subsection. We cn chieve the sme gols s the two third-order polynomils meeting t vi point much more simply: let us use only one polynomil, forced to go through the vi point. Here re the constrints for meeting the required ngles with smooth motion (since we use single continuous polynomil, the velocity nd ccelertions re gurnteed to mtch nd be continuous t the vi point): () S () ( ) t ( t ) f ( t ) f (7) Note in this cse since there is only one time rnge, it is convenient to tret ll times s bsolute, rther thn reltive s when we did mtching of two third-order polynomils. With five constrints single fourth-order polynomil is required: () t t t t t ( t) 4 t 3 t t (8) ive liner equtions in the five unknown polynomil coefficients i, i =,1,,3,4 result from the five constrints (7). Two of the unknown polynomil coefficients re found immeditely, from the initil time constrints: = S nd 1 =. The simplified 3x3 mtrix/vector eqution to solve for the remining three unknowns is: 1 t t S t 1 t t 3 S t 3t 4t 4 (9) The single fourth-order polynomil solution for ech joint is:
9 t t t t S S t t t t t St t S 3 3 t t t t 3t t 1 S S 4 3 t t t t (1) or the specil cse of t t, the solution is: t 3 3 t 4 4 t S S 16 7 S S 8 S S (11) Exmple Single ourth-order Polynomil with i Point Use single fourth-order polynomil for smooth joint spce trjectory genertion, plus motion through vi point, for one joint. Given 3, 18, 1, t 1.5, t 3 sec, find (t): S t t t t 3 ( t) 8.96t t t ( t) 48.89t t ( t) t ( ) (1) Agin, deg units re used throughout this exmple. These results re plotted in igure.
10 d dd ddd time (sec) igure. Single ourth-order Polynomil with i Point The (t) shpe for the single fourth-order polynomil is very similr to tht for the two thirdorder polynomils. Agin the joint ngle psses through the vi point 18 nd keeps going briefly, due to the fct tht the velocity is still positive beyond the first 1.5 sec. The pek for this fourth-order cse is slightly greter nd occurs in time slightly fter the pek for the two third-order polynomils exmple. The mximum ngle is 185.4, occurring t t = 1.74 sec. or the fourth-order MAX polynomil, side benefit hs risen: the jerk is now continuous t the vi time, where it ws discontinuous for the two mtched third-order polynomils. The jerk still hs n infinite spike t the strt nd end s we sw with the two third-order polynomils exmple, which is uncceptble. We now improve upon this with single sixth-order polynomil in the next subsection gin, this next subsection is originl work.
11 11.4 Single Sixth-Order Polynomil This section presents nother originl contribution. We cn chieve the sme gols s the single fourth-order polynomil with vi point nd eliminte the infinite spikes in jerk t the strt nd end s follows. To the previous five constrints, dd two more, for zero ccelertion t the strt nd end of the single motion rnge. Here re the seven constrints: () S () () ( ) t ( t ) f ( t ) f ( t ) f (13) With seven constrints single sixth-order polynomil is required: () t t t t t t t ( t) 6 t 5 t 4 t 3 t t ( t) 3 t t 1 t 6 t (14) Seven liner equtions in the seven unknown polynomil coefficients i, i =,1,,3,4,5,6 result from the seven constrints (13). Three of the unknown polynomil coefficients re found immeditely, from the initil time constrints: = S, 1 =, nd =. The simplified 4x4 mtrix/vector eqution to solve for the remining four unknowns is: t t t 3 S t t t t 4 S t 3 3 4t 5t 6t t t 3t 6 (15) The single sixth-order polynomil solution for ech joint is:
12 3 t t t t t t S S t t t t 3 3 t S5t 9tt 5t 3 S t t t t 3 3 t S8t 9tt t 3 S t t t t t t t S S 1t 15t t 5 t 6t 1 (16) or the specil cse of t t, the nlyticl solution is: t 4 4 t 5 5 t 6 6 t S S S S S S 3 S S (17) Exmple Single Sixth-Order Polynomil with i Point Use single sixth-order polynomil for smooth joint spce trjectory genertion, plus motion through vi point, for one joint. Given 3, 18, 1, t 1.5, t 3 sec, find (t): S t t t t t ( t) 55.3t 45.9t 16.6t 846.7t 4 3 ( t) 76.5t 173.7t t t 3 ( t) 116.t t t ( ) (18) Agin, deg units re used throughout this exmple. These results re plotted in igure 3.
13 d dd ddd time (sec) igure 3. Single Sixth-Order Polynomil with i Point Agin the (t) shpe is similr to the previous two exmples. The mximum ngle is MAX 185.6, occurring t t = 1.7 sec. Note the mgnitude of the ngulr velocity, ccelertion, nd jerk re not gretly different for the three cses we hve presented with vi point. Generlly they re lowest for the two third-order polynomils. This single sixth-order polynomil pproch hs eliminted the problem of infinite spikes in jerk t the strt nd end of ech motion. Here the jerk is discontinuous but it remins finite, which obeys the rule of thumb for mechnicl design/motion.
14 14 3. DISCUSSION All joint-spce trjectory genertion methods presented by previous uthors nd the new ones proposed in this rticle cn esily be pplied to multiple joints, simultneously but independently. This min originl contribution of this rticle is to present new joint spce trjectory genertion pproch using single polynomil for motion through given vi point. The robot does not need to stop t the vi point, but the motion must be smooth nd continuous through the vi point. A single fourth-order polynomil ws developed to chieve this smooth motion through vi point, to replce the two three-order polynomils mtched t the vi point in common usge tody. With the new pproch this vi point mtching comes utomticlly nd there is no need for two polynomil functions. An importnt secondry contribution of this rticle is to expose bd prctice in common usge in joint-spce trjectory genertion in robotics tody. We re not the first to notice this (Mcfrlne nd Croft 3, Gosselin nd Hdj-Messoud 1993, nd Petrinec nd Kovcic 7), but judging from the four robotics textbooks in the reference list published within the pst five yers, the messge hs not been widely understood (only Angeles, 7, is concerned with preventing infinite jerk). or ny functions in which the ccelertion is discontinuous, the ssocited jerk (time derivtive of ccelertion) function will hve infinite spikes t those ccelertion discontinuities. Cm design teches tht these infinite spikes in jerk re uncceptble nd must be voided by the designed/motion controller. By extension, in robot joint spce trjectory genertion, ny method which llows discontinuous ccelertion functions re uncceptble nd must not be used. Unfortuntely, the liner velocity with prbolic blends, the third-order polynomil, nd the two third-order polynomils mtched to meet t vi point ll suffer from discontinuity in ccelertion nd thus infinite spikes in jerk. These methods re still very much in use in the stndrd robotics textbooks tody these methods must be discrded. Otherwise, uncceptble noise, wer, stress, reduced life, nd the introduction of bd dynmics my result.
15 15 The novel single fourth-order polynomil introduced in this rticle to go through vi point lso suffers from this discontinuous ccelertion. Therefore, to remedy this we introduced novel sixthorder polynomil to move smoothly nd continuously through vi point, with finite jerk throughout the entire motion, which is cceptble. 4. CONCLUSION or joint-spce trjectory genertion in robotics, the recommendtions of this rticle re simple for ll joint motions without vi point, use the stndrd fifth-order polynomil from Crig (5). or ll joint motions with vi point, use the new single sixth-order polynomil introduced in this rticle. Not only does this new sixth-order polynomil keep the jerk finite, it lso chieves the vi point motion with only one function.
16 16 REERENCES J. Angeles, 7, undmentls of Robotic Mechnicl Systems, 3 rd edition, Springer, New York, NY. J.J. Crig, 5, Introduction to Robotics: Mechnics nd Control, 3 rd edition, Person Prentice Hll, Upper Sddle River, NJ. K.S. u, R.C. Gonzlez, nd C.S.G. Lee, 1987, Robotics: Control, Sensing, ision, nd Intelligence, McGrw-Hill, New York. C.M. Gosselin nd A. Hdj-Messoud, 1993, Automtic plnning of smooth trjectories for pick-nd-plce opertions, Journl of Mechnicl Design. Trns ASME, 115(3): R.N. Jzr, 1, Theory of Applied Robotics: Kinemtics, Dynmics, nd Control, nd edition, Springer, New York, NY. York. A.J. Koivo, 1989, undmentls for Control of Robotic Mnipultors, John Wiley & Sons, New S. Mcfrlne nd E.A. Croft, 3, "Jerk-bounded mnipultor trjectory plnning: Design for rel-time pplictions", IEEE Trnsctions on Robotics nd Automtion, 19(1): 4 5. R.E. Prkin, 1991, Applied Robotic Anlysis, Prentice Hll, Englewood Cliffs, NJ. R.P. Pul nd H. Zhong, 1984, Robot Motion Trjectory Specifiction nd Genertion, Second Interntionl Symposium on Robotics Reserch, Kyoto, Jpn, August. K. Petrinec nd Z. Kovcic, 7, "Trjectory Plnning Algorithm Bsed on the Continuity of Jerk, Proceedings of the 15 th Mediterrnen Conference on Control nd Automtion, IEEE, Athens, Greece, July. L. Scivicco nd B. Sicilino, 8, Robotics: Modeling, Plnning, nd Control, McGrw-Hill, New York. M.W. Spong, S. Hutchinson, nd M. idysgr, 5, Robot Modeling nd Control, John Wiley & Sons, New York. R.L. Willims II, 11, Improved Robotics Joint-Spce Trjectory Genertion with i Point, Proceedings of the ASME IDETC/CIE, Wshington DC, August 8-31, 11, DETC
Research Article Simplified Robotics Joint-Space Trajectory Generation with a via Point Using a Single Polynomial
Robotics Volume, Article ID 75958, 6 pages http://dx.doi.org/.55//75958 Research Article Simplified Robotics Joint-Space Trajectory Generation with a via Point Using a Single Polynomial Robert L. Williams
More informationDepartment of Mechanical Engineering MECE 551 Final examination Winter 2008 April 16, 9:00 11:30. Question Value Mark
Deprtment of Mechnicl Engineering MECE 55 Finl exmintion Winter 8 April 6, 9: :3 Notes: You my hve your text book nd one pge formul sheet Electronic devices re not llowed except n pproved clcultor NAME:
More informationMath Lecture 23
Mth 8 - Lecture 3 Dyln Zwick Fll 3 In our lst lecture we delt with solutions to the system: x = Ax where A is n n n mtrix with n distinct eigenvlues. As promised, tody we will del with the question of
More informationUnit #9 : Definite Integral Properties; Fundamental Theorem of Calculus
Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl
More informationI1 = 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 informationTaylor Polynomial Inequalities
Tylor Polynomil Inequlities Ben Glin September 17, 24 Abstrct There re instnces where we my wish to pproximte the vlue of complicted function round given point by constructing simpler function such s polynomil
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationPolynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230
Polynomil Approimtions for the Nturl Logrithm nd Arctngent Functions Mth 23 You recll from first semester clculus how one cn use the derivtive to find n eqution for the tngent line to function t given
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationAQA Further Pure 1. Complex Numbers. Section 1: Introduction to Complex Numbers. The number system
Complex Numbers Section 1: Introduction to Complex Numbers Notes nd Exmples These notes contin subsections on The number system Adding nd subtrcting complex numbers Multiplying complex numbers Complex
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
More informationSeptember 13 Homework Solutions
College of Engineering nd Computer Science Mechnicl Engineering Deprtment Mechnicl Engineering 5A Seminr in Engineering Anlysis Fll Ticket: 5966 Instructor: Lrry Cretto Septemer Homework Solutions. Are
More informationP 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)
1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this
More informationChapter 0. What is the Lebesgue integral about?
Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More informationHandout: Natural deduction for first order logic
MATH 457 Introduction to Mthemticl Logic Spring 2016 Dr Json Rute Hndout: Nturl deduction for first order logic We will extend our nturl deduction rules for sententil logic to first order logic These notes
More informationGeneralized Fano and non-fano networks
Generlized Fno nd non-fno networks Nildri Ds nd Brijesh Kumr Ri Deprtment of Electronics nd Electricl Engineering Indin Institute of Technology Guwhti, Guwhti, Assm, Indi Emil: {d.nildri, bkri}@iitg.ernet.in
More informationLooking for All Palindromes in a String
Looking or All Plindromes in String Shih Jng Pn nd R C T Lee Deprtment o Computer Science nd Inormtion Engineering, Ntionl Chi-Nn University, Puli, Nntou Hsien,, Tiwn, ROC sjpn@lgdoccsiencnuedutw, rctlee@ncnuedutw
More informationp-adic Egyptian Fractions
p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction
More informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
More informationReview 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 informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk bout solving systems of liner equtions. These re problems tht give couple of equtions with couple of unknowns, like: 6 2 3 7 4
More information1 Probability Density Functions
Lis Yn CS 9 Continuous Distributions Lecture Notes #9 July 6, 28 Bsed on chpter by Chris Piech So fr, ll rndom vribles we hve seen hve been discrete. In ll the cses we hve seen in CS 9, this ment tht our
More informationWe will see what is meant by standard form very shortly
THEOREM: For fesible liner progrm in its stndrd form, the optimum vlue of the objective over its nonempty fesible region is () either unbounded or (b) is chievble t lest t one extreme point of the fesible
More informationZ b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...
Chpter 7 Numericl Methods 7. Introduction In mny cses the integrl f(x)dx cn be found by finding function F (x) such tht F 0 (x) =f(x), nd using f(x)dx = F (b) F () which is known s the nlyticl (exct) solution.
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 7.. Introduction to Integrtion. 7. Integrl Clculus As ws the cse with the chpter on differentil
More informationHere we study square linear systems and properties of their coefficient matrices as they relate to the solution set of the linear system.
Section 24 Nonsingulr Liner Systems Here we study squre liner systems nd properties of their coefficient mtrices s they relte to the solution set of the liner system Let A be n n Then we know from previous
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationENGI 3424 Engineering Mathematics Five Tutorial Examples of Partial Fractions
ENGI 44 Engineering Mthemtics Five Tutoril Exmples o Prtil Frctions 1. Express x in prtil rctions: x 4 x 4 x 4 b x x x x Both denomintors re liner non-repeted ctors. The cover-up rule my be used: 4 4 4
More informationChapter 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 information13.3 CLASSICAL STRAIGHTEDGE AND COMPASS CONSTRUCTIONS
33 CLASSICAL STRAIGHTEDGE AND COMPASS CONSTRUCTIONS As simple ppliction of the results we hve obtined on lgebric extensions, nd in prticulr on the multiplictivity of extension degrees, we cn nswer (in
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationa a a a a a a a a a a a a a a a a a a a a a a a In this section, we introduce a general formula for computing determinants.
Section 9 The Lplce Expnsion In the lst section, we defined the determinnt of (3 3) mtrix A 12 to be 22 12 21 22 2231 22 12 21. In this section, we introduce generl formul for computing determinnts. Rewriting
More informationDuality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.
Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we
More informationLecture Solution of a System of Linear Equation
ChE Lecture Notes, Dept. of Chemicl Engineering, Univ. of TN, Knoville - D. Keffer, 5/9/98 (updted /) Lecture 8- - Solution of System of Liner Eqution 8. Why is it importnt to e le to solve system of liner
More informationThomas Whitham Sixth Form
Thoms Whithm Sith Form Pure Mthemtics Unit C Alger Trigonometry Geometry Clculus Vectors Trigonometry Compound ngle formule sin sin cos cos Pge A B sin Acos B cos Asin B A B sin Acos B cos Asin B A B cos
More informationMath 360: A primitive integral and elementary functions
Mth 360: A primitive integrl nd elementry functions D. DeTurck University of Pennsylvni October 16, 2017 D. DeTurck Mth 360 001 2017C: Integrl/functions 1 / 32 Setup for the integrl prtitions Definition:
More informationReinforcement learning II
CS 1675 Introduction to Mchine Lerning Lecture 26 Reinforcement lerning II Milos Huskrecht milos@cs.pitt.edu 5329 Sennott Squre Reinforcement lerning Bsics: Input x Lerner Output Reinforcement r Critic
More informationCS 275 Automata and Formal Language Theory
CS 275 Automt nd Forml Lnguge Theory Course Notes Prt II: The Recognition Problem (II) Chpter II.6.: Push Down Automt Remrk: This mteril is no longer tught nd not directly exm relevnt Anton Setzer (Bsed
More informationEnergy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon
Phys463.nb 49 7 Energy Bnds Ref: textbook, Chpter 7 Q: Why re there insultors nd conductors? Q: Wht will hppen when n electron moves in crystl? In the previous chpter, we discussed free electron gses,
More informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More informationTHREE-DIMENSIONAL KINEMATICS OF RIGID BODIES
THREE-DIMENSIONAL KINEMATICS OF RIGID BODIES 1. TRANSLATION Figure shows rigid body trnslting in three-dimensionl spce. Any two points in the body, such s A nd B, will move long prllel stright lines if
More informationand that at t = 0 the object is at position 5. Find the position of the object at t = 2.
7.2 The Fundmentl Theorem of Clculus 49 re mny, mny problems tht pper much different on the surfce but tht turn out to be the sme s these problems, in the sense tht when we try to pproimte solutions we
More informationThe Wave Equation I. MA 436 Kurt Bryan
1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string
More informationLecture XVII. Vector functions, vector and scalar fields Definition 1 A vector-valued function is a map associating vectors to real numbers, that is
Lecture XVII Abstrct We introduce the concepts of vector functions, sclr nd vector fields nd stress their relevnce in pplied sciences. We study curves in three-dimensionl Eucliden spce nd introduce the
More informationJim Lambers MAT 169 Fall Semester Lecture 4 Notes
Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of
More informationWe partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.
Mth 255 - Vector lculus II Notes 4.2 Pth nd Line Integrls We begin with discussion of pth integrls (the book clls them sclr line integrls). We will do this for function of two vribles, but these ides cn
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationAQA Further Pure 2. Hyperbolic Functions. Section 2: The inverse hyperbolic functions
Hperbolic Functions Section : The inverse hperbolic functions Notes nd Emples These notes contin subsections on The inverse hperbolic functions Integrtion using the inverse hperbolic functions Logrithmic
More information1 Bending of a beam with a rectangular section
1 Bending of bem with rectngulr section x3 Episseur b M x 2 x x 1 2h M Figure 1 : Geometry of the bem nd pplied lod The bem in figure 1 hs rectngur section (thickness 2h, width b. The pplied lod is pure
More informationParse trees, ambiguity, and Chomsky normal form
Prse trees, miguity, nd Chomsky norml form In this lecture we will discuss few importnt notions connected with contextfree grmmrs, including prse trees, miguity, nd specil form for context-free grmmrs
More information1.9 C 2 inner variations
46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for
More informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
More informationLecture 6. Notes. Notes. Notes. Representations Z A B and A B R. BTE Electronics Fundamentals August Bern University of Applied Sciences
Lecture 6 epresenttions epresenttions TE52 - Electronics Fundmentls ugust 24 ern University of pplied ciences ev. c2d5c88 6. Integers () sign-nd-mgnitude representtion The set of integers contins the Nturl
More informationApplicable Analysis and Discrete Mathematics available online at
Applicble Anlysis nd Discrete Mthemtics vilble online t http://pefmth.etf.rs Appl. Anl. Discrete Mth. 4 (2010), 23 31. doi:10.2298/aadm100201012k NUMERICAL ANALYSIS MEETS NUMBER THEORY: USING ROOTFINDING
More informationPART 1 MULTIPLE CHOICE Circle the appropriate response to each of the questions below. Each question has a value of 1 point.
PART MULTIPLE CHOICE Circle the pproprite response to ech of the questions below. Ech question hs vlue of point.. If in sequence the second level difference is constnt, thn the sequence is:. rithmetic
More informationVyacheslav Telnin. Search for New Numbers.
Vycheslv Telnin Serch for New Numbers. 1 CHAPTER I 2 I.1 Introduction. In 1984, in the first issue for tht yer of the Science nd Life mgzine, I red the rticle "Non-Stndrd Anlysis" by V. Uspensky, in which
More informationCS 275 Automata and Formal Language Theory
CS 275 Automt nd Forml Lnguge Theory Course Notes Prt II: The Recognition Problem (II) Chpter II.5.: Properties of Context Free Grmmrs (14) Anton Setzer (Bsed on book drft by J. V. Tucker nd K. Stephenson)
More informationIn-Class Problems 2 and 3: Projectile Motion Solutions. In-Class Problem 2: Throwing a Stone Down a Hill
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Deprtment of Physics Physics 8T Fll Term 4 In-Clss Problems nd 3: Projectile Motion Solutions We would like ech group to pply the problem solving strtegy with the
More informationarxiv:math/ v2 [math.ho] 16 Dec 2003
rxiv:mth/0312293v2 [mth.ho] 16 Dec 2003 Clssicl Lebesgue Integrtion Theorems for the Riemnn Integrl Josh Isrlowitz 244 Ridge Rd. Rutherford, NJ 07070 jbi2@njit.edu Februry 1, 2008 Abstrct In this pper,
More informationSOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı
Mthemticl nd Computtionl Applictions, Vol., No., pp. 37-5, 5 http://dx.doi.org/.99/mc-5- SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD Mehmet Pkdemirli nd Gözde
More informationConvert the NFA into DFA
Convert the NF into F For ech NF we cn find F ccepting the sme lnguge. The numer of sttes of the F could e exponentil in the numer of sttes of the NF, ut in prctice this worst cse occurs rrely. lgorithm:
More informationTime Optimal Control of the Brockett Integrator
Milno (Itly) August 8 - September, 011 Time Optiml Control of the Brockett Integrtor S. Sinh Deprtment of Mthemtics, IIT Bomby, Mumbi, Indi (emil : sunnysphs4891@gmil.com) Abstrct: The Brockett integrtor
More information19 Optimal behavior: Game theory
Intro. to Artificil Intelligence: Dle Schuurmns, Relu Ptrscu 1 19 Optiml behvior: Gme theory Adversril stte dynmics hve to ccount for worst cse Compute policy π : S A tht mximizes minimum rewrd Let S (,
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk out solving systems of liner equtions. These re prolems tht give couple of equtions with couple of unknowns, like: 6= x + x 7=
More informationEngineering Analysis ENG 3420 Fall Dan C. Marinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00
Engineering Anlysis ENG 3420 Fll 2009 Dn C. Mrinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00 Lecture 13 Lst time: Problem solving in preprtion for the quiz Liner Algebr Concepts Vector Spces,
More information10 Vector Integral Calculus
Vector Integrl lculus Vector integrl clculus extends integrls s known from clculus to integrls over curves ("line integrls"), surfces ("surfce integrls") nd solids ("volume integrls"). These integrls hve
More informationState space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies
Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response
More informationChapter 3 Solving Nonlinear Equations
Chpter 3 Solving Nonliner Equtions 3.1 Introduction The nonliner function of unknown vrible x is in the form of where n could be non-integer. Root is the numericl vlue of x tht stisfies f ( x) 0. Grphiclly,
More informationLinear Motion. Kinematics Quantities
Liner Motion Physics 101 Eyres Kinemtics Quntities Time Instnt t Fundmentl Time Interl Defined Position x Fundmentl Displcement Defined Aerge Velocity g Defined Aerge Accelertion g Defined 1 Kinemtics
More informationA. Limits - L Hopital s Rule ( ) How to find it: Try and find limits by traditional methods (plugging in). If you get 0 0 or!!, apply C.! 1 6 C.
A. Limits - L Hopitl s Rule Wht you re finding: L Hopitl s Rule is used to find limits of the form f ( x) lim where lim f x x! c g x ( ) = or lim f ( x) = limg( x) = ". ( ) x! c limg( x) = 0 x! c x! c
More information4.4 Areas, Integrals and Antiderivatives
. res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order
More informationTABLE OF CONTENTS 3 CHAPTER 1
TABLE OF CONTENTS 3 CHAPTER 1 Set Lnguge & Nottion 3 CHAPTER 2 Functions 3 CHAPTER 3 Qudrtic Functions 4 CHAPTER 4 Indices & Surds 4 CHAPTER 5 Fctors of Polynomils 4 CHAPTER 6 Simultneous Equtions 4 CHAPTER
More informationNondeterminism and Nodeterministic Automata
Nondeterminism nd Nodeterministic Automt 61 Nondeterminism nd Nondeterministic Automt The computtionl mchine models tht we lerned in the clss re deterministic in the sense tht the next move is uniquely
More information1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.
1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt
More informationFig. 1. Open-Loop and Closed-Loop Systems with Plant Variations
ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses
More informationPredict Global Earth Temperature using Linier Regression
Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id
More information13: Diffusion in 2 Energy Groups
3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups
More informationExploring parametric representation with the TI-84 Plus CE graphing calculator
Exploring prmetric representtion with the TI-84 Plus CE grphing clcultor Richrd Prr Executive Director Rice University School Mthemtics Project rprr@rice.edu Alice Fisher Director of Director of Technology
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More informationIntegral points on the rational curve
Integrl points on the rtionl curve y x bx c x ;, b, c integers. Konstntine Zeltor Mthemtics University of Wisconsin - Mrinette 750 W. Byshore Street Mrinette, WI 5443-453 Also: Konstntine Zeltor P.O. Box
More informationLecture 20: Numerical Integration III
cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed
More informationf(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral
Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one
More informationLesson 1: Quadratic Equations
Lesson 1: Qudrtic Equtions Qudrtic Eqution: The qudrtic eqution in form is. In this section, we will review 4 methods of qudrtic equtions, nd when it is most to use ech method. 1. 3.. 4. Method 1: Fctoring
More informationIntermediate Math Circles Wednesday, November 14, 2018 Finite Automata II. Nickolas Rollick a b b. a b 4
Intermedite Mth Circles Wednesdy, Novemer 14, 2018 Finite Automt II Nickols Rollick nrollick@uwterloo.c Regulr Lnguges Lst time, we were introduced to the ide of DFA (deterministic finite utomton), one
More informationMath 3B Final Review
Mth 3B Finl Review Written by Victori Kl vtkl@mth.ucsb.edu SH 6432u Office Hours: R 9:45-10:45m SH 1607 Mth Lb Hours: TR 1-2pm Lst updted: 12/06/14 This is continution of the midterm review. Prctice problems
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More informationLecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.
Lecture 3 3 Solving liner equtions In this lecture we will discuss lgorithms for solving systems of liner equtions Multiplictive identity Let us restrict ourselves to considering squre mtrices since one
More information20 MATHEMATICS POLYNOMIALS
0 MATHEMATICS POLYNOMIALS.1 Introduction In Clss IX, you hve studied polynomils in one vrible nd their degrees. Recll tht if p(x) is polynomil in x, the highest power of x in p(x) is clled the degree of
More informationThe Properties of Stars
10/11/010 The Properties of Strs sses Using Newton s Lw of Grvity to Determine the ss of Celestil ody ny two prticles in the universe ttrct ech other with force tht is directly proportionl to the product
More informationalong the vector 5 a) Find the plane s coordinate after 1 hour. b) Find the plane s coordinate after 2 hours. c) Find the plane s coordinate
L8 VECTOR EQUATIONS OF LINES HL Mth - Sntowski Vector eqution of line 1 A plne strts journey t the point (4,1) moves ech hour long the vector. ) Find the plne s coordinte fter 1 hour. b) Find the plne
More informationImproper Integrals. Type I Improper Integrals How do we evaluate an integral such as
Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph
More informationMath 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of
More informationINTRODUCTION. The three general approaches to the solution of kinetics problems are:
INTRODUCTION According to Newton s lw, prticle will ccelerte when it is subjected to unblnced forces. Kinetics is the study of the reltions between unblnced forces nd the resulting chnges in motion. The
More informationGenetic Programming. Outline. Evolutionary Strategies. Evolutionary strategies Genetic programming Summary
Outline Genetic Progrmming Evolutionry strtegies Genetic progrmming Summry Bsed on the mteril provided y Professor Michel Negnevitsky Evolutionry Strtegies An pproch simulting nturl evolution ws proposed
More informationCalculus of Variations
Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function
More information