Lecture 21: Bezier Approximation and de Casteljau s Algorithm. and thou shalt be near unto me Genesis 45:10

Size: px
Start display at page:

Download "Lecture 21: Bezier Approximation and de Casteljau s Algorithm. and thou shalt be near unto me Genesis 45:10"

Transcription

1 Lecure 2: Bezier Approximaion and de Caeljau Algorihm and hou hal be near uno me Genei 45:0. Inroducion In Lecure 20, we ued inerpolaion o pecify hape. Bu inerpolaion i no alway a good way o decribe he conour of a curve or urface. To accuraely reproduce complicaed hape, we may need o inerpolae lo of daa. Polynomial inerpolaion for many poin i impracical becaue he degree of he inerpolan can ge exremely high leading o low and numerically unable compuaion. Alo polynomial inerpolan may ocillae unnecearily and fail o reproduce he deired hape (ee Figure ). Thu, even if we were o pecify more and more poin, here i no guaranee ha he polynomial inerpolan would converge o he curve or urface we wih o repreen. Daa Polynomial Inerpolan Figure : Lagrange inerpolaion. Noice he ocillaion in he inerpolaing polynomial curve, even hough here i no ocillaion in he original daa poin. Spline inerpolaion -- ha i, inerpolaion by piecewie polynomial funcion -- i beer compuaionally becaue pline allow u o keep he degree low. Bu inerpolaing pline may ill ocillae unnecearily and fail o reproduce he deired hape. Our approach here i raher o abandon inerpolaion alogeher and o ake a very differen approach o decribing he hape of a curve or urface. Given a relaively mall collecion of poin in affine pace, we are going o inveigae mehod for generaing polynomial curve and urface ha approximae he hape decribed by hee poin. We hall no ini ha our curve and urface go hrough hee poin, bu we hall ini ha hee curve and urface come near o he poin and capure in ome mahemaically precie way he hape defined by hee poin. A uual we begin wih cheme for curve and laer exend our echnique o urface.

2 De Caeljau Algorihm Le u reurn for a momen o where we began our inveigaion of polynomial curve and urface: Lagrange inerpolaion and Neville algorihm. Recall ha Neville algorihm (Figure 2) i a dynamic programming procedure for compuing poin along a polynomial inerpolan. We are going o ar our inveigaion of approximaion cheme by uing he ame baic riangular rucure bu implifying he compuaion along he edge. The imple hing -- one migh almo ay he only hing -- we know how o do i linear inerpolaion. All our inerpolaion procedure, and epecially Neville algorihm, are baed on hi imple idea or ome varian hereof. Wha make Neville algorihm he lea bi complicaed i ha we perform a differen linear inerpolaion a each node of he diagram. To ake he ame riangular rucure and make he evaluaion algorihm a eay a poible, we will perform he ame linear inerpolaion a each node. Thi idea generae he algorihm repreened in Figure 3. 3 P 023 () 0 P 02 () P 23 () P 0 () P 2 () P 23 () P 0 P P 2 P 3 Figure 2: Neville algorihm (unnormalized) for cubic polynomial inerpolaion. b b B() a b a b a a b a b a P 0 P P 2 P 3 Figure 3: The de Caeljau algorihm for a cubic Bezier curve B() in he inerval [a, b]. The label on every edge mu be normalized by dividing by b a, o ha he label along arrow enering each node um o one. 2

3 b Thi algorihm repreened in Figure 3 i called de Caeljau evaluaion algorihm, and he curve ha emerge a he apex of hi diagram are called Bezier curve. Inermediae node marked and alo repreen Bezier curve, bu of lower degree. Thu he de Caeljau algorihm i a dynamic programming algorihm for compuing poin on a Bezier curve. Typically, for reaon ha will become clear in he nex ecion, Bezier curve are rericed o he inerval [a, b]. Uually, for impliciy, we ake a = 0 and b =, bu here are cae, a we hall ee laer on, where i i ueful o allow a and b o be arbirary a long a b > a. Noice ha when a = 0 and b =, no normalizaion i required. The de Caeljau algorihm ha he following elegan geomeric inerpreaion. Since each node repreen a linear inerpolaion, each node ymbolize a poin on he line egmen joining he wo poin whoe arrow poin ino he node. Drawing all hee line egmen generae he relli in Figure 4. a b P 0 P a P 2 b a b b b a a B() a a Figure 4: Geomeric conrucion algorihm for a poin on a cubic Bezier curve baed on a geomeric inerpreaion of he de Caeljau evaluaion algorihm. A he parameer, each line egmen in he relli i pli in he raio ( a)/(b ). b P 3 We are going o udy he geomeric characeriic of curve generaed by de Caeljau algorihm. We begin wih ome imple feaure; in ubequen lecure we will develop ome more advanced mahemaical properie of hi approximaion cheme. 3. Elemenary Properie of Bezier Curve Bezier curve have he following elemenary properie:. Polynomial Paramerizaion 2. Affine Invariance 3. Convex Hull 4. Symmery 5. Inerpolaion of End Poin 3

4 Below we briefly dicu and derive each of hee properie in urn, and we explain a well why hee feaure are imporan for compuer graphic. Mo of hee properie can be proved by direc obervaion or by eay inducive argumen uing he de Caeljau algorihm. To e up hee inducive argumen, le B[P 0,..., P n ]() denoe he Bezier curve over he inerval [a, b] wih affine conrol poin P 0,..., P n. Then he la age of he de Caeljau algorihm can be wrien a B[P 0,..., P n ]() = b b a B[P 0,..., P n ]() + a b a B[P,..., P n ](). (3.) We are now ready o proceed wih our derivaion.. Polynomial Paramerizaion In he de Caeljau algorihm, he only operaion we perform involving he funcion along he edge are addiion and muliplicaion (ee Figure 3). Since he funcion along he edge are linear polynomial, i follow ha a Bezier curve wih n + conrol poin i a polynomial curve of degree n becaue here are n level from he conrol poin a he bae o he curve a he apex of he riangle. (Thi reul alo follow by an eay inducion from Equaion (3.)). Since Bezier curve are polynomial curve, all he ool we know for polynomial apply. 2. Affine Invariance A curve i aid o be affinely invarian if i coni of a collecion of poin in affine pace. Equivalenly, a curve cheme i aid o be affinely invarian if applying an affine ranformaion o he conrol poin ranform every poin on he curve by he ame affine ranformaion. Affine invariance i an eay conequence of he de Caeljau algorihm. ince by linear inerpolaion every node in he algorihm i affinely invarian (ee Figure 5,6). (Thi reul alo follow by an eay inducion from Equaion (3.)). Affine invariance i a crucial feaure for any curve cheme becaue i aer ha he curve i independen of he choice of he coordinae yem. Thi propery i eenial for a good approximaion cheme, ince in a ypical geomeric model many differen coordinae yem are preen. Affine invariance guaranee ha he curve will be he ame no maer which coordinae yem i invoked. 4

5 b B() + v a b P 0 + v b + v +v a b +v a + v + v a b a b a P + v P 2 + v P 3 + v Figure 5: Tranlaion invariance. Tranlaing each conrol poin by a vecor v ranlae every poin on he curve by he ame vecor v. b ( B(),) M a b (P 0,) M (,) M (,) M b a b a (,) M (,) M (,) M a b a b (P,) M (P 2,) M a (P 3,) M Figure 6: Affine invariance. Tranforming each conrol poin by an affine map M ranform every poin on he curve by he ame affine map M. 3. Convex Hull Propery A e S of poin in affine pace i aid o be convex if whenever P and Q are poin in S he enire line egmen from P o Q lie in S (ee Figure 7). The inerecion S of a collecion of convex e {S i } i a convex e becaue if P and Q are poin in S, hey mu alo be poin in each of he e S i. Since, by aumpion, he e S i are convex, he enire line egmen from P o Q lie in each e S i. Hence he enire line egmen from P o Q lie in he inerecion S, o S oo i convex. The convex hull of a collecion of poin in affine pace i he inerecion of all he convex e conaining he poin. Since he inerecion of convex e i a convex e, he convex hull i he malle convex e conaining he poin. For wo poin, he convex hull i he line egmen 5

6 joining he poin. For hree non-collinear poin, he convex hull i he riangle whoe verice are he hree poin. The convex hull of a finie collecion of poin in he plane can be found mechanically by placing a nail a each poin, reching a rubber band o ha i inerior conain all he nail, and hen releaing he rubber band. When he rubber band come o re on he nail, he inerior of he rubber band i he convex hull of he poin. P Q P Q (a) Convex Se (b) Non-Convex Se Figure 7: In a convex e (a) he line egmen joining any wo poin in he e lie enirely wihin he e. In a non-convex e (b) par of he line egmen joining wo poin in he e may lie ouide he e. Since he convex hull of wo poin i he line egmen joining he wo poin, ConvexHull{P 0, P } = {c 0 P 0 + c P c 0 + c = and c 0,c 0}. More generally i follow by a imple inducive argumen (ee Exercie ) ha n n ConvexHull{P 0,..., P n } = { c k P k c k = and c k 0}. k= 0 k =0 Bezier curve alway lie in he convex hull of heir conrol poin. Tha i, B[P 0,..., P n ]() ConvexHull{P 0,..., P n }. Again we can prove hi aerion by a imple inducive argumen. Fir recall ha, by convenion, we alway reric he Bezier curve B[P 0,..., P n ]() in (3.) o he parameer inerval a b. Wih hi rericion, he convex hull propery i cerainly rue for a Bezier curve wih only wo conrol poin ince, by conrucion, hi curve i he line egmen joining he wo conrol poin. More generally uppoe ha hi reul i valid for Bezier curve wih n conrol poin. By (3.), B[P 0,..., P n ]() lie on he line egmen joining he poin B[P 0,..., P n ]() and B[P,..., P n ](), and by he inducive hypohei B[P 0,..., P n ]() and B[P,..., P n ]() boh lie in he convex hull of he poin P 0,..., P n. Bu if wo poin lie in a convex e, he enire line egmen joining hem alo lie in he e; hu he enire Bezier curve B[P 0,..., P n ](), a b, mu lie in ConvexHull{P 0,..., P n } 6

7 0 2. The convex hull propery i imporan becaue i conrain Bezier curve o lie in he proximiy of heir conrol poin. Thi propery i a vial feaure for an approximaion cheme. Scieni and engineer no only require curve ha approximae he hape defined by heir conrol poin, hey alo demand curve ha lie in he ame region of pace a heir conrol poin. To be ueful in compuer graphic, he curve mu be viible on he graphic erminal. The convex hull propery guaranee ha if all he conrol poin are viible on he graphic erminal, hen he enire curve i viible a well. The rericion a b on he parameer i here preciely o guaranee he convex hull propery. 4. Symmery Replacing by a + b revere he order of he parameer domain. A he parameer varie from a o b, he curve B[P 0,..., P n ](a + b ) ravere he ame poin a B[P 0,..., P n ]() bu in he direcion from b o a raher han from a o b. Thu B[P 0,..., P n ](a + b ) i eenially he ame curve a B[P 0,..., P n ]() bu wih oppoie orienaion. Similarly, revering he order of he conrol poin of a Bezier curve generae he ame Bezier curve bu wih oppoie orienaion. Analyically hi mean ha, B[P n,..., P 0 ]() = B[P 0,..., P n ](a + b ) a b. (3.2) Thi ymmery propery i wha profeional engineer and mo naive uer would naurally expec of a imple approximaion cheme, o i i graifying o ee ha i hold for all Bezier curve. To prove (3.2), imply replace by a + b in he de Caeljau diagram and oberve ha he new diagram i he mirror image of he de Caeljau diagram for B[P n,..., P 0 ]() (ee Figure 8). B() a a b a b b a b a b a b P 0 P P 2 P 3 Figure 8: The de Caeljau algorihm for B[P 0,..., P n ](a + b ). Compare o Figure 3 wih he conrol poin in revere order. 7

8 5. Inerpolaion of End Poin Unlike Lagrange polynomial, Bezier curve generally do no inerpolae all heir conrol poin. Bu Bezier curve alway inerpolae heir fir and la conrol poin. In fac, B[P 0,..., P n ](a) = P 0 B[P 0,..., P n ](b) = P n. The fir reul follow eaily from eing = a in de Caeljau algorihm and oberving ha all he label on lef poining arrow become zero while all he label on righ poining arrow become one. If k 0, hen any pah from P k o he apex of he riangle mu ravere a lea one lef poining arrow, o here i no conribuion from P k o he value of he curve a = a. On he oher hand, when = a all he label on he ingle pah from P 0 o he apex of he riangle are one. Hence B[P 0,..., P n ](a) = P 0. A imilar argumen for = b how ha B[P 0,..., P n ](b) = P n. Again an eay inducive argumen baed on (3.) yield he ame reul. Inerpolaing end poin i imporan becaue we ofen wan o connec wo curve. To aure ha wo Bezier curve join a heir end poin, all we need o do i o make ure ha he fir conrol poin of he econd curve i he ame a he la conrol poin of he fir curve. Thi device inure coninuiy. In he nex ecion, we hall develop echnique for guaraneeing higher order moohne beween adjacen Bezier curve, bu before we can perform hi analyi we need o know how o differeniae Bezier curve. 4. Differeniaion Algorihm for Bezier Curve To compue he derivaive of a Bezier curve, we need o differeniae he de Caeljau algorihm. Differeniaing he de Caeljau algorihm for a degree n Bezier curve i acually quie eay: imply differeniae he label -- ha i, replace and -- on any level of he de Caeljau algorihm and muliply he oupu by n (ee Figure 9). We hall defer he proof of hi differeniaion algorihm for Bezier curve ill Lecure 23. In hi ecion we will udy he conequence of hi differeniaion algorihm. There are everal hing o noice abou Figure 9. Fir, by running only he lowe level of he differeniaion algorihm, we ee ha, up o a conan muliple, we can hink of he derivaive of a cubic Bezier curve wih conrol poin {P k } a a quadraic Bezier polynomial wih conrol vecor {P k+ P k } (ee Figure 0). 8

9 0 2 B () = 3 _ P 0 P P 2 P 3 Figure 9: Compuing he fir derivaive of a cubic Bezier curve wih conrol poin P 0,P, P 2,P 3 by differeniaing he label on he boom level of he de Caeljau algorihm. To ge he correc derivaive B (), we mu muliply he oupu of hi algorihm by n = 3. B () = 3 _ P P 0 P2 P P 3 P 2 Figure 0: The derivaive of a cubic Bezier curve wih conrol poin {P k } i, up o a conan muliple, a quadraic Bezier polynomial wih conrol vecor {P k+ P k }. In general, up o a conan muliple, he derivaive of a Bezier curve wih conrol poin {P k } i a Bezier polynomial of one lower degree wih conrol vecor {P k+ P k }. Therefore, by inducion, we can expre he kh derivaive of a degree n Bezier curve a a Bezier polynomial of degree n k. Moreover, o compue he kh derivaive of a Bezier curve, we can differeniae he label on any k n! level he de Caeljau algorihm and muliply he oupu by (ee Figure ). (n k)! 9

10 0 2 B () = 6 _ P 0 P P 2 P 3 Figure : Compuing he econd derivaive of a cubic Bezier curve wih conrol poin P 0,P, P 2,P 3 by differeniaing he label on he boom wo level of he de Caeljau algorihm. Here o ge he correc econd derivaive B (), we mu muliply he oupu of hi algorihm by 6. We prefaced our dicuion of differeniaion by aying ha we waned o be able o join wo Bezier curve ogeher moohly. To do o, we need o calculae derivaive a heir end poin. Le B() be a Bezier curve wih conrol poin P 0,K,P n. Subiuing = 0 or = ino he differeniaion algorihm, we ee ha: B(0) = P 0 B() = P n B (0) = n(p P 0 ) B () = n(p n P n ) B (0) = n(n )(P 2 2P + P 0 ) B () = n(n )(P n 2P n + P n 2 ) In general, i follow by inducion on k ha he kh derivaive a = 0 depend only on he fir k + conrol poin, and he kh derivaive a = depend only on he la k + conrol poin. Suppoe hen ha we are given a Bezier curve wih conrol poin P 0,K,P n and we wan o conruc anoher Bezier curve wih conrol poin Q 0,K,Q n ha mee he fir curve and mache i fir r derivaive a i end poin. From he reul in he previou paragraph, we find ha he conrol poin Q 0,K,Q n mu aify he following conrain: r = 0 : Q 0 = P n r = : Q Q 0 = P n P n Q = P n + (P n P n ) r = 2 : Q 2 2Q + Q 0 = P n 2P n + P n 2 Q 2 = P n 2 + 4(P n P n ) and o on for higher value of r. Each addiional derivaive allow u o olve for one addiional conrol poin. Clearly we could go on in hi manner olving for one poin a a ime. An alernaive approach ha avoid all hi ediou compuaion will be preened in Lecure 22. 0

11 5. Tenor Produc Bezier Surface A enor produc Bezier pach i a recangular urface pach generaed by a recangular array of conrol poin. Connecing conrol poin wih adjacen indice by raigh line generae a conrol polyhedron ha conrol he hape of he Bezier pach in much he ame way ha he Bezier conrol polygon conrol he hape of a Bezier curve (ee Figure 2 and 3). In paricular, dragging a conrol poin pull he urface pach in he ame general direcion a he conrol poin. = d P 03 P 3 P 23 P 33 = a = b P 02 P 2 P 22 P 32 P 0 P P 2 P 3 = c P 00 P 0 P 20 P 30 (a) Domain -- Recangle (b) Range -- Recangular Array of Poin Figure 2: Daa for a bicubic enor produc Bezier pach. Noice ha he domain i imply a recangle and ha, unlike Lagrange inerpolaion, he domain ha no node and no grid. Compare o Figure 8 in Lecure 20. Figure 3: A bicubic enor produc Bezier urface wih i conrol polyhedron, formed by connecing conrol poin wih adjacen indice. The conrol poin are he ame a hoe for he Lagrange urface in Figure 0 of Lecure 20. To conruc a enor produc Bezier pach B(,) from a recangular array of conrol poin {P ij }, we proceed in a manner imilar o he conrucion of a enor produc Lagrange urface,

12 excep ha Lagrange curve and Neville algorihm are replaced by Bezier curve and de Caeljau algorihm. Le P i () = B[P i,0,k,p i,n ]() i = 0,K,m, (4.) be he Bezier curve for he conrol poin P i,0,k,p i,n. For each fixed value of, le B(,) = B[P 0 (),K, P m ()]() (4.2) be he Bezier curve for he conrol poin P 0 (),K, P m (). Then a varie from a o b and varie from c o d, he curve B(,) weep ou a urface (ee Figure 4,5). Thi urface i called a enor produc Bezier pach. Equaion (4.), (4.2) and Figure 4,5 ugge he following evaluaion algorihm for enor produc Bezier pache: fir ue de Caeljau algorihm m + ime o compue he poin a he parameer along he degree n Bezier curve P 0 (),..., P m () ; hen ue de Caeljau algorihm one more ime o compue he poin a he parameer along he degree m Bezier curve wih conrol poin P 0 (),..., P m () (ee Figure 6). P 3 P 23 P 33 P 03 P 2 P 22 P 02 P 0 P 0 () P 00 P P 0 B(, ) P 2 () P ( ) P 2 P 20 P 30 P 3 ( ) Figure 4: Conrucion of poin on a bicubic enor produc Bezier urface B(, ). Fir he Bezier curve P i (), i = 0,..., 3 are conruced from he conrol poin P i0, P i, P i2, P i3. Then for a fixed value of, he Bezier curve B(, ) i conruced uing he poin P 0 (), P (), P 2 (), P 3 () a conrol poin. In general, he Bezier urface B(, ) doe no inerpolae i conrol poin. Compare o Figure 9 of Lecure 20 for bicubic Lagrange inerpolaion. P 32 P 3 2

13 Figure 5: The bicubic Bezier pach in Figure 3 along wih i cubic Bezier conrol curve. Noice ha only he boundary conrol curve are inerpolaed by he urface. Compare o he enor produc Lagrange urface in Figure 0 of Lecure 20. B(,) P 0 () P () P 2 () P 00 P 0 P 02 P 0 P P 2 P 20 P 2 P 22 Figure 6: De Caeljau evaluaion algorihm for a biquadraic Bezier pach. The hree lower riangle repreen Bezier curve in he direcion, and he upper riangle blend hee curve in he direcion. Compare o Figure of Lecure Neville algorihm for a biquadraic Lagrange inerpolaing pach. 3

14 Alernaively, inead of aring wih he Bezier curve P i () = B[P i,0,k,p i,n ]() i = 0,K,m, we could begin wih he Bezier curve P j () = B[P0, j,k,p m, j ]() j = 0,K,n, and for each fixed value of le B(,) = B[P 0 (),K,Pn ()]() be he Bezier curve for he conrol poin P 0 (),K,Pn (). Again a varie from a o b and varie from c o d, he curve B(,) weep ou a urface. We hall how in a laer lecure ha hi urface i idenical o he enor produc Bezier pach in our previou conrucion. Thu we have he following alernaive evaluaion algorihm for enor produc Bezier pache: fir ue de Caeljau algorihm n + ime o compue he poin a he parameer along he degree m Bezier curve P 0 (),...,Pn (); hen ue de Caeljau algorihm one more ime o compue he poin a he parameer along he degree n Bezier curve wih conrol poin P 0 (),...,Pn () (ee Figure 7). When m < n, hen hi algorihm i more efficien han he previou algorihm. B(,) P 0 () P () P 2 () P 00 P 0 P 20 P 0 P P 2 P 02 P 2 P 22 Figure 7: An alernaive verion of he de Caeljau evaluaion algorihm for a biquadraic Bezier pach. The hree lower riangle repreen Bezier curve in he direcion, and he upper riangle blend hee curve in he direcion. Compare o Figure 6. 4

15 Tenor produc Bezier pache inheri many of he characeriic properie of Bezier curve: hey are affine invarian, lie in he convex hull of heir conrol poin, and inerpolae heir corner conrol poin (ee Exercie 6-8). Thee properie follow eaily from Figure 4,6 and he correponding properie of Bezier curve. Moreover he boundarie of a enor produc Bezier pach are he Bezier curve deermined by heir boundary conrol poin, ince B(0,) = P 0 () = B[P 0,0,K,P 0,n ]() and B(,) = P m () = B[P m,0,k,p m,n ]() B(,0) = P 0 () = B[P0,0,K,P m,0 ]() B(,) = P n () = B[P0,n,K,P m,n ](). I follow ha alhough enor produc Bezier pache do no generally inerpolae heir conrol poin, hey alway inerpolae he four corner poin P 00,P m0, P 0n,P mn. To compue he parial derivaive of a Bezier pach, we can apply our procedure for differeniaing he de Caeljau algorihm for Bezier curve. Conider Figure 6. We can compue B / imply by differeniaing he de Caeljau algorihm for each of he Bezier curve P 0 (),K, P m () a he bae of he diagram. Similarly, we can compue B / by differeniaing any one of he upper level of he diagram and muliplying he reul by he degree in (ee Figure 8). Symmeric reul hold for differeniaing he algorihm in Figure 7: imply revere he role of and. The normal vecor N o a Bezier pach i given by N = B B. B B = 2 _ P 0 () P () P 2 () P 0 () P () P 2 () Figure 8: Compuing he parial derivaive of a biquadraic Bezier pach by applying he procedure for differeniaing he de Caeljau algorihm for Bezier curve. To find B / (lef) imply differeniae he de Caeljau algorihm for each of he Bezier curve P 0 (), P (), P 2 (). To find B / (righ) imply differeniae one level of he de Caeljau algorihm i and muliply he reul by he degree in. 5

16 Exercie:. Le P 0,..., P n be a collecion of poin in an affine pace. Prove ha n n ConvexHull{P 0,..., P n } = { c k P k c k = and c k 0}. k= 0 k =0 2. Prove ha B[P 0,..., P n ](b) = P n. 3. Show ha Lagrange inerpolaing polynomial curve are affine invarian. 4. Give an example o how ha a degree n Lagrange inerpolaing polynomial wih node 0 < < L < n doe no necearily aify he convex hull propery on he inerval [ 0, n ]. 5. Give an example o how ha Bezier curve do no necearily reproduce polynomial curve. Tha i, give an example o how ha if here i a degree n polynomial P() uch ha P k = P( k ), hen i doe no necearily follow ha B[P 0,K,P n ]() = P(). 6. Show ha enor produc Bezier pache are affine invarian. 7. Show ha every enor produc Bezier pach lie in he convex hull of i conrol poin. 8. Show ha every enor produc Bezier pach inerpolae i corner conrol poin. 9. Le P i, j, 0 i m, 0 j n, be a recangular array of conrol poin, and le B[P i, j ](,), 0,, denoe he correponding enor produc Bezier pach of bidegree (m, n). Show ha he urface pach B[P i, j ](,) ha he following ymmery properie: a. B[P j,i ](,) = B[P i, j ](,) ; b. B[P i,n j ](,) = B[P i, j ](, ) ; c. B[P m i, j ](,) = B[P i, j ](,). 0. Le P i, j, 0 i m, 0 j n, be a recangular array of conrol poin, and le P i () = B[P i,0,k,p i,n ]() i = 0,K,m, be a equence of Bezier curve. Show how o combine Neville algorihm and de Caeljau algorihm o generae a urface C(,) ha inerpolae he curve P 0 (),...,P m () a he parameer value 0,..., m. 6

17 . Given poin and derivaive daa (R 0,v 0 ),...,(R n,v n ), explain how o place Bezier conrol poin o generae a C piecewie cubic curve o inerpolae hi daa. 2. The formula for he uni angen U(), he curvaure K(), and he orion T() of a parameric curve P() are given by i. U() = P () P () ii. K() = P () P () P () 3 iii. T() = P () ( P () P ()) P () P () 2 a. Compue he uni angen, curvaure, and orion of a Bezier curve a = 0,. b. Find condiion on he conrol poin of a Bezier curve C() o ha i mache a given Bezier curve B() wih coninuou uni angen, curvaure, and orion. 7

Problem Set If all directed edges in a network have distinct capacities, then there is a unique maximum flow.

Problem Set If all directed edges in a network have distinct capacities, then there is a unique maximum flow. CSE 202: Deign and Analyi of Algorihm Winer 2013 Problem Se 3 Inrucor: Kamalika Chaudhuri Due on: Tue. Feb 26, 2013 Inrucion For your proof, you may ue any lower bound, algorihm or daa rucure from he ex

More information

Randomized Perfect Bipartite Matching

Randomized Perfect Bipartite Matching Inenive Algorihm Lecure 24 Randomized Perfec Biparie Maching Lecurer: Daniel A. Spielman April 9, 208 24. Inroducion We explain a randomized algorihm by Ahih Goel, Michael Kapralov and Sanjeev Khanna for

More information

Algorithmic Discrete Mathematics 6. Exercise Sheet

Algorithmic Discrete Mathematics 6. Exercise Sheet Algorihmic Dicree Mahemaic. Exercie Shee Deparmen of Mahemaic SS 0 PD Dr. Ulf Lorenz 7. and 8. Juni 0 Dipl.-Mah. David Meffer Verion of June, 0 Groupwork Exercie G (Heap-Sor) Ue Heap-Sor wih a min-heap

More information

18 Extensions of Maximum Flow

18 Extensions of Maximum Flow Who are you?" aid Lunkwill, riing angrily from hi ea. Wha do you wan?" I am Majikhie!" announced he older one. And I demand ha I am Vroomfondel!" houed he younger one. Majikhie urned on Vroomfondel. I

More information

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship Laplace Tranform (Lin & DeCarlo: Ch 3) ENSC30 Elecric Circui II The Laplace ranform i an inegral ranformaion. I ranform: f ( ) F( ) ime variable complex variable From Euler > Lagrange > Laplace. Hence,

More information

CS 473G Lecture 15: Max-Flow Algorithms and Applications Fall 2005

CS 473G Lecture 15: Max-Flow Algorithms and Applications Fall 2005 CS 473G Lecure 1: Max-Flow Algorihm and Applicaion Fall 200 1 Max-Flow Algorihm and Applicaion (November 1) 1.1 Recap Fix a direced graph G = (V, E) ha doe no conain boh an edge u v and i reveral v u,

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 31 Signal & Syem Prof. Mark Fowler Noe Se #27 C-T Syem: Laplace Tranform Power Tool for yem analyi Reading Aignmen: Secion 6.1 6.3 of Kamen and Heck 1/18 Coure Flow Diagram The arrow here how concepual

More information

Linear Algebra Primer

Linear Algebra Primer Linear Algebra rimer And a video dicuion of linear algebra from EE263 i here (lecure 3 and 4): hp://ee.anford.edu/coure/ee263 lide from Sanford CS3 Ouline Vecor and marice Baic Mari Operaion Deerminan,

More information

Introduction to Congestion Games

Introduction to Congestion Games Algorihmic Game Theory, Summer 2017 Inroducion o Congeion Game Lecure 1 (5 page) Inrucor: Thoma Keelheim In hi lecure, we ge o know congeion game, which will be our running example for many concep in game

More information

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method ME 352 GRHICL VELCITY NLYSIS 52 GRHICL VELCITY NLYSIS olygon Mehod Velociy analyi form he hear of kinemaic and dynamic of mechanical yem Velociy analyi i uually performed following a poiion analyi; ie,

More information

Introduction to SLE Lecture Notes

Introduction to SLE Lecture Notes Inroducion o SLE Lecure Noe May 13, 16 - The goal of hi ecion i o find a ufficien condiion of λ for he hull K o be generaed by a imple cure. I urn ou if λ 1 < 4 hen K i generaed by a imple curve. We will

More information

16 Max-Flow Algorithms and Applications

16 Max-Flow Algorithms and Applications Algorihm A proce canno be underood by opping i. Underanding mu move wih he flow of he proce, mu join i and flow wih i. The Fir Law of Mena, in Frank Herber Dune (196) There a difference beween knowing

More information

1 Motivation and Basic Definitions

1 Motivation and Basic Definitions CSCE : Deign and Analyi of Algorihm Noe on Max Flow Fall 20 (Baed on he preenaion in Chaper 26 of Inroducion o Algorihm, 3rd Ed. by Cormen, Leieron, Rive and Sein.) Moivaion and Baic Definiion Conider

More information

Exponential Sawtooth

Exponential Sawtooth ECPE 36 HOMEWORK 3: PROPERTIES OF THE FOURIER TRANSFORM SOLUTION. Exponenial Sawooh: The eaie way o do hi problem i o look a he Fourier ranform of a ingle exponenial funcion, () = exp( )u(). From he able

More information

u(t) Figure 1. Open loop control system

u(t) Figure 1. Open loop control system Open loop conrol v cloed loop feedbac conrol The nex wo figure preen he rucure of open loop and feedbac conrol yem Figure how an open loop conrol yem whoe funcion i o caue he oupu y o follow he reference

More information

Chapter 7: Inverse-Response Systems

Chapter 7: Inverse-Response Systems Chaper 7: Invere-Repone Syem Normal Syem Invere-Repone Syem Baic Sar ou in he wrong direcion End up in he original eady-ae gain value Two or more yem wih differen magniude and cale in parallel Main yem

More information

6.8 Laplace Transform: General Formulas

6.8 Laplace Transform: General Formulas 48 HAP. 6 Laplace Tranform 6.8 Laplace Tranform: General Formula Formula Name, ommen Sec. F() l{ f ()} e f () d f () l {F()} Definiion of Tranform Invere Tranform 6. l{af () bg()} al{f ()} bl{g()} Lineariy

More information

2. VECTORS. R Vectors are denoted by bold-face characters such as R, V, etc. The magnitude of a vector, such as R, is denoted as R, R, V

2. VECTORS. R Vectors are denoted by bold-face characters such as R, V, etc. The magnitude of a vector, such as R, is denoted as R, R, V ME 352 VETS 2. VETS Vecor algebra form he mahemaical foundaion for kinemaic and dnamic. Geomer of moion i a he hear of boh he kinemaic and dnamic of mechanical em. Vecor anali i he imehonored ool for decribing

More information

Admin MAX FLOW APPLICATIONS. Flow graph/networks. Flow constraints 4/30/13. CS lunch today Grading. in-flow = out-flow for every vertex (except s, t)

Admin MAX FLOW APPLICATIONS. Flow graph/networks. Flow constraints 4/30/13. CS lunch today Grading. in-flow = out-flow for every vertex (except s, t) /0/ dmin lunch oday rading MX LOW PPLIION 0, pring avid Kauchak low graph/nework low nework direced, weighed graph (V, ) poiive edge weigh indicaing he capaciy (generally, aume ineger) conain a ingle ource

More information

Mon Apr 2: Laplace transform and initial value problems like we studied in Chapter 5

Mon Apr 2: Laplace transform and initial value problems like we studied in Chapter 5 Mah 225-4 Week 2 April 2-6 coninue.-.3; alo cover par of.4-.5, EP 7.6 Mon Apr 2:.-.3 Laplace ranform and iniial value problem like we udied in Chaper 5 Announcemen: Warm-up Exercie: Recall, The Laplace

More information

Network Flows: Introduction & Maximum Flow

Network Flows: Introduction & Maximum Flow CSC 373 - lgorihm Deign, nalyi, and Complexiy Summer 2016 Lalla Mouaadid Nework Flow: Inroducion & Maximum Flow We now urn our aenion o anoher powerful algorihmic echnique: Local Search. In a local earch

More information

, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max

, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max ecure 8 7. Sabiliy Analyi For an n dimenional vecor R n, he and he vecor norm are defined a: = T = i n i (7.) I i eay o how ha hee wo norm aify he following relaion: n (7.) If a vecor i ime-dependen, hen

More information

CHAPTER 7: SECOND-ORDER CIRCUITS

CHAPTER 7: SECOND-ORDER CIRCUITS EEE5: CI RCUI T THEORY CHAPTER 7: SECOND-ORDER CIRCUITS 7. Inroducion Thi chaper conider circui wih wo orage elemen. Known a econd-order circui becaue heir repone are decribed by differenial equaion ha

More information

18.03SC Unit 3 Practice Exam and Solutions

18.03SC Unit 3 Practice Exam and Solutions Sudy Guide on Sep, Dela, Convoluion, Laplace You can hink of he ep funcion u() a any nice mooh funcion which i for < a and for > a, where a i a poiive number which i much maller han any ime cale we care

More information

Geometric Path Problems with Violations

Geometric Path Problems with Violations Click here o view linked Reference 1 1 1 0 1 0 1 0 1 0 1 Geomeric Pah Problem wih Violaion Anil Mahehwari 1, Subha C. Nandy, Drimi Paanayak, Saanka Roy and Michiel Smid 1 1 School of Compuer Science, Carleon

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

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

More information

arxiv: v1 [cs.cg] 21 Mar 2013

arxiv: v1 [cs.cg] 21 Mar 2013 On he rech facor of he Thea-4 graph Lui Barba Proenji Boe Jean-Lou De Carufel André van Renen Sander Verdoncho arxiv:1303.5473v1 [c.cg] 21 Mar 2013 Abrac In hi paper we how ha he θ-graph wih 4 cone ha

More information

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves

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

More information

EE Control Systems LECTURE 2

EE Control Systems LECTURE 2 Copyrigh F.L. Lewi 999 All righ reerved EE 434 - Conrol Syem LECTURE REVIEW OF LAPLACE TRANSFORM LAPLACE TRANSFORM The Laplace ranform i very ueful in analyi and deign for yem ha are linear and ime-invarian

More information

CSC 364S Notes University of Toronto, Spring, The networks we will consider are directed graphs, where each edge has associated with it

CSC 364S Notes University of Toronto, Spring, The networks we will consider are directed graphs, where each edge has associated with it CSC 36S Noe Univeriy of Torono, Spring, 2003 Flow Algorihm The nework we will conider are direced graph, where each edge ha aociaed wih i a nonnegaive capaciy. The inuiion i ha if edge (u; v) ha capaciy

More information

Flow networks. Flow Networks. A flow on a network. Flow networks. The maximum-flow problem. Introduction to Algorithms, Lecture 22 December 5, 2001

Flow networks. Flow Networks. A flow on a network. Flow networks. The maximum-flow problem. Introduction to Algorithms, Lecture 22 December 5, 2001 CS 545 Flow Nework lon Efra Slide courey of Charle Leieron wih mall change by Carola Wenk Flow nework Definiion. flow nework i a direced graph G = (V, E) wih wo diinguihed verice: a ource and a ink. Each

More information

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2)

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2) Laplace Tranform Maoud Malek The Laplace ranform i an inegral ranform named in honor of mahemaician and aronomer Pierre-Simon Laplace, who ued he ranform in hi work on probabiliy heory. I i a powerful

More information

Algorithms and Data Structures 2011/12 Week 9 Solutions (Tues 15th - Fri 18th Nov)

Algorithms and Data Structures 2011/12 Week 9 Solutions (Tues 15th - Fri 18th Nov) Algorihm and Daa Srucure 2011/ Week Soluion (Tue 15h - Fri 18h No) 1. Queion: e are gien 11/16 / 15/20 8/13 0/ 1/ / 11/1 / / To queion: (a) Find a pair of ube X, Y V uch ha f(x, Y) = f(v X, Y). (b) Find

More information

Piecewise-Defined Functions and Periodic Functions

Piecewise-Defined Functions and Periodic Functions 28 Piecewie-Defined Funcion and Periodic Funcion A he ar of our udy of he Laplace ranform, i wa claimed ha he Laplace ranform i paricularly ueful when dealing wih nonhomogeneou equaion in which he forcing

More information

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

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

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

Introduction to Probability and Statistics Slides 4 Chapter 4 Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random

More information

Let. x y. denote a bivariate time series with zero mean.

Let. x y. denote a bivariate time series with zero mean. Linear Filer Le x y : T denoe a bivariae ime erie wih zero mean. Suppoe ha he ime erie {y : T} i conruced a follow: y a x The ime erie {y : T} i aid o be conruced from {x : T} by mean of a Linear Filer.

More information

CHAPTER 7. Definition and Properties. of Laplace Transforms

CHAPTER 7. Definition and Properties. of Laplace Transforms SERIES OF CLSS NOTES FOR 5-6 TO INTRODUCE LINER ND NONLINER PROBLEMS TO ENGINEERS, SCIENTISTS, ND PPLIED MTHEMTICINS DE CLSS NOTES COLLECTION OF HNDOUTS ON SCLR LINER ORDINRY DIFFERENTIL EQUTIONS (ODE")

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

PHYSICS Solving Equations

PHYSICS Solving Equations Sepember 20, 2013 PHYSIS Solving Equaion Sepember 2013 www.njcl.org Solving for a Variable Our goal i o be able o olve any equaion for any variable ha appear in i. Le' look a a imple equaion fir. The variable

More information

Notes on cointegration of real interest rates and real exchange rates. ρ (2)

Notes on cointegration of real interest rates and real exchange rates. ρ (2) Noe on coinegraion of real inere rae and real exchange rae Charle ngel, Univeriy of Wiconin Le me ar wih he obervaion ha while he lieraure (mo prominenly Meee and Rogoff (988) and dion and Paul (993))

More information

NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY

NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY Ian Crawford THE INSTITUTE FOR FISCAL STUDIES DEPARTMENT OF ECONOMICS, UCL cemmap working paper CWP02/04 Neceary and Sufficien Condiion for Laen

More information

Flow Networks. Ma/CS 6a. Class 14: Flow Exercises

Flow Networks. Ma/CS 6a. Class 14: Flow Exercises 0/0/206 Ma/CS 6a Cla 4: Flow Exercie Flow Nework A flow nework i a digraph G = V, E, ogeher wih a ource verex V, a ink verex V, and a capaciy funcion c: E N. Capaciy Source 7 a b c d e Sink 0/0/206 Flow

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Chapter 9 - The Laplace Transform

Chapter 9 - The Laplace Transform Chaper 9 - The Laplace Tranform Selece Soluion. Skech he pole-zero plo an region of convergence (if i exi) for hee ignal. ω [] () 8 (a) x e u = 8 ROC σ ( ) 3 (b) x e co π u ω [] ( ) () (c) x e u e u ROC

More information

CS4445/9544 Analysis of Algorithms II Solution for Assignment 1

CS4445/9544 Analysis of Algorithms II Solution for Assignment 1 Conider he following flow nework CS444/944 Analyi of Algorihm II Soluion for Aignmen (0 mark) In he following nework a minimum cu ha capaciy 0 Eiher prove ha hi aemen i rue, or how ha i i fale Uing he

More information

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese

More information

From Complex Fourier Series to Fourier Transforms

From Complex Fourier Series to Fourier Transforms Topic From Complex Fourier Series o Fourier Transforms. Inroducion In he previous lecure you saw ha complex Fourier Series and is coeciens were dened by as f ( = n= C ne in! where C n = T T = T = f (e

More information

Analyze patterns and relationships. 3. Generate two numerical patterns using AC

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

Graphs III - Network Flow

Graphs III - Network Flow Graph III - Nework Flow Flow nework eup graph G=(V,E) edge capaciy w(u,v) 0 - if edge doe no exi, hen w(u,v)=0 pecial verice: ource verex ; ink verex - no edge ino and no edge ou of Aume every verex v

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

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

Lecture 33: November 29

Lecture 33: November 29 36-705: Inermediae Saisics Fall 2017 Lecurer: Siva Balakrishnan Lecure 33: November 29 Today we will coninue discussing he boosrap, and hen ry o undersand why i works in a simple case. In he las lecure

More information

SOLUTIONS TO ECE 3084

SOLUTIONS TO ECE 3084 SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

CHAPTER. Forced Equations and Systems { } ( ) ( ) 8.1 The Laplace Transform and Its Inverse. Transforms from the Definition.

CHAPTER. Forced Equations and Systems { } ( ) ( ) 8.1 The Laplace Transform and Its Inverse. Transforms from the Definition. CHAPTER 8 Forced Equaion and Syem 8 The aplace Tranform and I Invere Tranform from he Definiion 5 5 = b b {} 5 = 5e d = lim5 e = ( ) b {} = e d = lim e + e d b = (inegraion by par) = = = = b b ( ) ( )

More information

arxiv: v2 [cs.cg] 30 Apr 2015

arxiv: v2 [cs.cg] 30 Apr 2015 Diffue Reflecion Diameer in Simple Polygon arxiv:1302.2271v2 [c.cg] 30 Apr 2015 Gill Bareque Deparmen of Compuer Science Technion Haifa, Irael bareque@c.echnion.ac.il Eli Fox-Epein Deparmen of Compuer

More information

t )? How would you have tried to solve this problem in Chapter 3?

t )? How would you have tried to solve this problem in Chapter 3? Exercie 9) Ue Laplace ranform o wrie down he oluion o 2 x x = F in x = x x = v. wha phenomena do oluion o hi DE illurae (even hough we're forcing wih in co )? How would you have ried o olve hi problem

More information

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015 Explaining Toal Facor Produciviy Ulrich Kohli Universiy of Geneva December 2015 Needed: A Theory of Toal Facor Produciviy Edward C. Presco (1998) 2 1. Inroducion Toal Facor Produciviy (TFP) has become

More information

Linear Algebra Primer

Linear Algebra Primer Linear Algebra Primer Juan Carlo Nieble and Ranja Krihna Sanford Viion and Learning Lab Anoher, ver in-deph linear algebra review from CS229 i available here: hp://c229.anford.edu/ecion/c229-linalg.pdf

More information

CONTROL SYSTEMS. Chapter 3 Mathematical Modelling of Physical Systems-Laplace Transforms. Prof.Dr. Fatih Mehmet Botsalı

CONTROL SYSTEMS. Chapter 3 Mathematical Modelling of Physical Systems-Laplace Transforms. Prof.Dr. Fatih Mehmet Botsalı CONTROL SYSTEMS Chaper Mahemaical Modelling of Phyical Syem-Laplace Tranform Prof.Dr. Faih Mehme Boalı Definiion Tranform -- a mahemaical converion from one way of hinking o anoher o make a problem eaier

More information

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

MATH 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

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

More information

Discussion Session 2 Constant Acceleration/Relative Motion Week 03

Discussion Session 2 Constant Acceleration/Relative Motion Week 03 PHYS 100 Dicuion Seion Conan Acceleraion/Relaive Moion Week 03 The Plan Today you will work wih your group explore he idea of reference frame (i.e. relaive moion) and moion wih conan acceleraion. You ll

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

y z P 3 P T P1 P 2. Werner Purgathofer. b a

y z P 3 P T P1 P 2. Werner Purgathofer. b a Einführung in Viual Compuing Einführung in Viual Compuing 86.822 in co T P 3 P co in T P P 2 co in Geomeric Tranformaion Geomeric Tranformaion W P h f Werner Purgahofer b a Tranformaion in he Rendering

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

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

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

Interpolation and Pulse Shaping

Interpolation and Pulse Shaping EE345S Real-Time Digial Signal Proceing Lab Spring 2006 Inerpolaion and Pule Shaping Prof. Brian L. Evan Dep. of Elecrical and Compuer Engineering The Univeriy of Texa a Auin Lecure 7 Dicree-o-Coninuou

More information

Chapter 7: Solving Trig Equations

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

CONTROL SYSTEMS. Chapter 10 : State Space Response

CONTROL SYSTEMS. Chapter 10 : State Space Response CONTROL SYSTEMS Chaper : Sae Space Repone GATE Objecive & Numerical Type Soluion Queion 5 [GATE EE 99 IIT-Bombay : Mark] Conider a econd order yem whoe ae pace repreenaion i of he form A Bu. If () (),

More information

Christos Papadimitriou & Luca Trevisan November 22, 2016

Christos Papadimitriou & Luca Trevisan November 22, 2016 U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream

More information

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation

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

Average Case Lower Bounds for Monotone Switching Networks

Average Case Lower Bounds for Monotone Switching Networks Average Cae Lower Bound for Monoone Swiching Nework Yuval Filmu, Toniann Piai, Rober Robere, Sephen Cook Deparmen of Compuer Science Univeriy of Torono Monoone Compuaion (Refreher) Monoone circui were

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

Effects of Coordinate Curvature on Integration

Effects of Coordinate Curvature on Integration Effecs of Coordinae Curvaure on Inegraion Chrisopher A. Lafore clafore@gmail.com Absrac In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure of he manifold

More information

10. State Space Methods

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

The Residual Graph. 11 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm

The Residual Graph. 11 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm Augmening Pah Algorihm Greedy-algorihm: ar wih f (e) = everywhere find an - pah wih f (e) < c(e) on every edge augmen flow along he pah repea a long a poible The Reidual Graph From he graph G = (V, E,

More information

Main Reference: Sections in CLRS.

Main Reference: Sections in CLRS. Maximum Flow Reied 09/09/200 Main Reference: Secion 26.-26. in CLRS. Inroducion Definiion Muli-Source Muli-Sink The Ford-Fulkeron Mehod Reidual Nework Augmening Pah The Max-Flow Min-Cu Theorem The Edmond-Karp

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Integration Over Manifolds with Variable Coordinate Density

Integration Over Manifolds with Variable Coordinate Density Inegraion Over Manifolds wih Variable Coordinae Densiy Absrac Chrisopher A. Lafore clafore@gmail.com In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure

More information

Price of Stability and Introduction to Mechanism Design

Price of Stability and Introduction to Mechanism Design Algorihmic Game Theory Summer 2017, Week 5 ETH Zürich Price of Sabiliy and Inroducion o Mechanim Deign Paolo Penna Thi i he lecure where we ar deigning yem which involve elfih player. Roughly peaking,

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

Computing Similarity between Piecewise-Linear Functions

Computing Similarity between Piecewise-Linear Functions Compuing Similariy beween Piecewie-Linear Funcion Pankaj K. Agarwal Dep. of Compuer Science, Duke Univeriy Maaren Löffler Compuer Science Dep. Univeriy of California, Irvine Bori Aronov Dep. of Compuer

More information

13.1 Accelerating Objects

13.1 Accelerating Objects 13.1 Acceleraing Objec A you learned in Chaper 12, when you are ravelling a a conan peed in a raigh line, you have uniform moion. However, mo objec do no ravel a conan peed in a raigh line o hey do no

More information

Lecture 10: The Poincaré Inequality in Euclidean space

Lecture 10: The Poincaré Inequality in Euclidean space Deparmens of Mahemaics Monana Sae Universiy Fall 215 Prof. Kevin Wildrick n inroducion o non-smooh analysis and geomery Lecure 1: The Poincaré Inequaliy in Euclidean space 1. Wha is he Poincaré inequaliy?

More information

Reminder: Flow Networks

Reminder: Flow Networks 0/0/204 Ma/CS 6a Cla 4: Variou (Flow) Execie Reminder: Flow Nework A flow nework i a digraph G = V, E, ogeher wih a ource verex V, a ink verex V, and a capaciy funcion c: E N. Capaciy Source 7 a b c d

More information

13.1 Circuit Elements in the s Domain Circuit Analysis in the s Domain The Transfer Function and Natural Response 13.

13.1 Circuit Elements in the s Domain Circuit Analysis in the s Domain The Transfer Function and Natural Response 13. Chaper 3 The Laplace Tranform in Circui Analyi 3. Circui Elemen in he Domain 3.-3 Circui Analyi in he Domain 3.4-5 The Tranfer Funcion and Naural Repone 3.6 The Tranfer Funcion and he Convoluion Inegral

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Network Flow. Data Structures and Algorithms Andrei Bulatov

Network Flow. Data Structures and Algorithms Andrei Bulatov Nework Flow Daa Srucure and Algorihm Andrei Bulao Algorihm Nework Flow 24-2 Flow Nework Think of a graph a yem of pipe We ue hi yem o pump waer from he ource o ink Eery pipe/edge ha limied capaciy Flow

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

How to Solve System Dynamic s Problems

How to Solve System Dynamic s Problems How o Solve Sye Dynaic Proble A ye dynaic proble involve wo or ore bodie (objec) under he influence of everal exernal force. The objec ay uliaely re, ove wih conan velociy, conan acceleraion or oe cobinaion

More information

Chapter 6. Laplace Transforms

Chapter 6. Laplace Transforms Chaper 6. Laplace Tranform Kreyzig by YHLee;45; 6- An ODE i reduced o an algebraic problem by operaional calculu. The equaion i olved by algebraic manipulaion. The reul i ranformed back for he oluion of

More information

Module 2: Analysis of Stress

Module 2: Analysis of Stress Module/Leon Module : Anali of Sre.. INTRODUCTION A bod under he acion of eernal force, undergoe diorion and he effec due o hi em of force i ranmied hroughou he bod developing inernal force in i. To eamine

More information

ESTIMATES FOR THE DERIVATIVE OF DIFFUSION SEMIGROUPS

ESTIMATES FOR THE DERIVATIVE OF DIFFUSION SEMIGROUPS Elec. Comm. in Probab. 3 (998) 65 74 ELECTRONIC COMMUNICATIONS in PROBABILITY ESTIMATES FOR THE DERIVATIVE OF DIFFUSION SEMIGROUPS L.A. RINCON Deparmen of Mahemaic Univeriy of Wale Swanea Singleon Par

More information

18 Biological models with discrete time

18 Biological models with discrete time 8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so

More information

Sample Final Exam (finals03) Covering Chapters 1-9 of Fundamentals of Signals & Systems

Sample Final Exam (finals03) Covering Chapters 1-9 of Fundamentals of Signals & Systems Sample Final Exam Covering Chaper 9 (final04) Sample Final Exam (final03) Covering Chaper 9 of Fundamenal of Signal & Syem Problem (0 mar) Conider he caual opamp circui iniially a re depiced below. I LI

More information

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,

More information

FLAT CYCLOTOMIC POLYNOMIALS OF ORDER FOUR AND HIGHER

FLAT CYCLOTOMIC POLYNOMIALS OF ORDER FOUR AND HIGHER #A30 INTEGERS 10 (010), 357-363 FLAT CYCLOTOMIC POLYNOMIALS OF ORDER FOUR AND HIGHER Nahan Kaplan Deparmen of Mahemaic, Harvard Univeriy, Cambridge, MA nkaplan@mah.harvard.edu Received: 7/15/09, Revied:

More information

Echocardiography Project and Finite Fourier Series

Echocardiography Project and Finite Fourier Series Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every

More information

Physics 240: Worksheet 16 Name

Physics 240: Worksheet 16 Name Phyic 4: Workhee 16 Nae Non-unifor circular oion Each of hee proble involve non-unifor circular oion wih a conan α. (1) Obain each of he equaion of oion for non-unifor circular oion under a conan acceleraion,

More information