Approximate merging of a pair of BeÂzier curves
|
|
- Esther Flynn
- 5 years ago
- Views:
Transcription
1 COMPUTER-AIDED DESIGN Computer-Aded Desgn 33 (1) 15±136 Approxmate mergng of a par of BeÂzer curves Sh-Mn Hu a,b, *, Rou-Feng Tong c, Tao Ju a,b, Ja-Guang Sun a,b a Natona CAD Engneerng Center, Tsnghua Unversty, Beng 184, Peope's Repubc of Chna b Department of Computer Scence and Technoogy, Tsnghua Unversty, Beng 184, Peope's Repubc of Chna c Department of Computer Scence, Zheang Unversty, Hangzhou 317, Peope's Repubc of Chna Receved 18 September 1999; receved n revsed form 4 Apr ; accepted 9 Apr Abstract Ths paper deas wth the mergng probem,.e. to approxmate two adacent BeÂzer curves by a snge BeÂzer curve. A nove approach for approxmate mergng s ntroduced n the paper by usng the constraned optmzaton method. The basc dea of ths method s to nd condtons for the precse mergng of BeÂzer curves rst, and then compute the constraned optmzaton souton by movng the contro ponts. ªDscreteº coef cent norm n L sense and ªsquared dfference ntegraº norm are used n our method. Contnuty at the endponts of curves are consdered n the mergng process, and approxmate mergng wth ponts constrants are aso dscussed. Further, t s shown that the degree eevaton of orgna BeÂzer curves w reduce the mergng error. q 1 Esever Scence Ltd. A rghts reserved. Keywords: BeÂzer curve; Approxmate converson; Mergng 1. Introducton and probem statement Parametrc poynoma representatons are wdey used n CAD systems whch mode free-form curves and surfaces. Bernsten±Bezer, Schoenberg-B-spne and Hermte± Coons type bass functons are frequenty used n dfferent systems. Wth the avaabty of a fast growng varety of modeng systems the demand has rsen for exchangng curve and surface descrptons between varous CAD systems [1]. The genera am when transferrng geometrc nformaton from one system to another s to ensure a hgh degree of accuracy, the east possbe oss of nformaton and a sma amount of geometrc data for communcaton. Converson from one poynoma base to another can be acheved by drect matrx mutpcaton. For reducng the amount of communcatng data, approxmate converson s consdered by Hoschek []. As mentoned by Hoschek, the approxmate converson ncudes the foowng. ² Degree reducton: to nd a parametrc curve of degree n to approxmate the gven curve of degree m n, m : ² Mergng: to merge as many as possbe curve segments of degree M to one curve segment of degree N M # N : Degree reducton methods of BeÂzer and Ba curves have * Correspondng author. Te.: ; fax: E-ma address: shmn@tsnghua.edu.cn (S.-M. Hu). been extensvey nvestgated by Watkns and Worsey [3], Lachance [4], Eck [5,6], Bogack et a. [7] and Hu et a. [8,9] etc., degree reducton of B-spne curves has been consdered by Peg et a. [1], degree reducton of parametrc surfaces aso has been consdered by Hu et a. [11,1]. However mergng of BeÂzer curves s st an nterestng and open probem. Probem. For two adacent BeÂzer curves P(u) and Q(v) # u; v # 1 wth correspondng contro ponts P and Q ˆ ; 1; ¼; n ; mergng of P(u) and Q(v) s a process that amounts to ndng an n degree BeÂzer curve R(t) wth contro ponts R ˆ ; 1; ¼n ; such that a sutabe dstance functon d R; R between R(t) and 8 >< R t ˆ >: P B n t Q B n t 1 # t # # t # 1 s mnmzed on the nterva [, 1], where s a subdvson parameter. The basc dea of ths artce s to nd condtons for precse mergng of BeÂzer curves. We modfy contro ponts of two BeÂzer curves such that the mod ed curves satsfy the precse mergng condtons, then new contro ponts of /1/$ - see front matter q 1 Esever Scence Ltd. A rghts reserved. PII: S1-4485()83-X
2 16 S.-M. Hu et a. / Computer-Aded Desgn 33 (1) 15±136 the merged curve can be obtaned by an extrapoaton agorthm. The rest of the paper s organzed as foows. Two constraned optmzaton methods wth dfferent optmzaton crtera are ntroduced, both of whch acheve mergng of BeÂzer curves by sovng near equaton systems. An exstence proof of the system of near equatons s gven n Secton 3. Condtons for matchng orgna endponts and dervatves are consdered n Secton 4, and mergng wth ponts constrants s deat wth n Secton 5. Fnay, we dscuss approxmate mergng wth degree eevaton n Secton 6.. Constraned optmzaton method for mergng.1. Condtons of precse mergng Consder two n degree BeÂzer curves P(u) and Q(v), we rst gve condtons of precse mergng. If P(u) and Q(v) can be merged precsey, ths s to say, there exsts an n degree BeÂzer curve R t ˆPn R B n t such that R t ˆ R t ( R t s de ned as n Eq. (1)), so we have P u u Note that uˆ1ˆ Q v v ˆ ; 1; ¼n: vˆ u ˆ t ; v ˆ t 1 ; by the dervatve formua of BeÂzer curves, Eq. () yeds 1 n n D P n ˆ 1 n 1 n D Q ; where D s the dfference operator de ned as DP ˆ P 11 P ; so we have the foowng theorem. Theorem 1. Two n degree BeÂzer curves P(u) and Q(v) can be merged precsey, f and ony f there exsts m m. ; such that D P n ˆ m D Q for ˆ ; 1; ¼; n: Proof. 1. Necessty. If P(u) and Q(v) can be merged precsey, from Eq. (3) we have D P n ˆ m D Q for ˆ ; 1; ¼; n; where m ˆ 1 : Therefore condton (4) s necessary.. Suf cency. Frst, we represent BeÂzer curves n the 3 4 A exponenta form. For any n degree BeÂzer curve S(t) wth contro ponts S ˆ ; 1; ¼; n we have S t ˆXn n t D S ˆ Xn n t 1 D S n : 5 Note that Eq. (5) hods true for any t t $ : Next we construct a new BeÂzer curve R(t), whose contro ponts R ˆ ; 1; ¼; n are de ned by R ˆ X P B 1 ; where ˆ m m 1 1 and B t s the degree BeÂzer base functon. We w see that P(u) and Q(v) can be precsey merged nto R(t). We cacuate D k R k ˆ ; 1; ¼n as " D k R ˆ Xk k P B 1 # ˆ ˆ Xk X 4 k P X k 1 B ˆ 1 k 3 1 k 5 ˆ Xk k P 1 k 1 k ˆ k Dk P : 6 Substtute Eq. (6) nto exponenta form (5), for t # t # ; we have R t ˆXn n t D R ˆ Xn t n D P ˆ P t Gven condton (4), for t # t # 1 ; we have R t ˆXn t n D P ˆ Xn ˆ Xn t 1 n D P n t 1 n m D Q ˆ Q t 1 Accordng to Eq. (1), P(u) and Q(v) can be precsey merged nto R(t) wth as the subdvson parameter. Ths competes the proof of Theorem 1... Constraned optmzaton method (1) We now consder the approxmate mergng probem
3 stated n the prevous Secton.1. In order to obtan the merged BeÂzer curve, we w expore two dfferent optmzaton crtera n the foowng dscusson. Frst, we use the ªdscreteº coef cent norm,.e. we set d E t ; F t ˆ P n E F ; where E(t) and F(t) are n degree BeÂzer curves wth contro ponts E and F, s the Eucdean norm. Mnmzaton of the norm w resut n the east dfference between correspondng contro ponts of the two curves. For BeÂzer curves P(u) and Q(v), perturbaton e ˆ e x ; e y ; ez Š T and d ˆ d x ; d y ; dz Š T can be obtaned for contro ponts P of P(u) and Q of Q(v), respectvey, such that the mod ed curves ^P u ˆXn ^P B n u ˆXn S.-M. Hu et a. / Computer-Aded Desgn 33 (1) 15± D d ˆ X 1 k d k ; kˆ k 1 by settng P 1 e B n u # u # 1 L e x ; L e y ; L e z ; L d x ; to be zero, we have e 1 Xn d Xn ˆn ˆ 1 n m 1 L d y ; L d z ; L L ; x y and L z ˆ ˆ ; 1; ¼; n; n ˆ ˆ ; 1; ¼; n; 11 1 ^Q v ˆXn ^Q B n v ˆXn Q 1 d B n v # v # 1 sats es the condtons of Theorem 1,.e., D ^P n ˆ m D ^Q : Then, ^P u and ^Q v can be precsey merged nto abeâzer curve R(t) of degree n, whch s consdered to be the approxmate merged curve of P(u) and Q(v). We determne e ; d ˆ ; 1; ¼; n by settng the optmzaton crteron O(e, d )as O e ; d ˆXn e 1 d ˆMn and Lagrange functon s de ned by L ˆ Xn 1 Xn e 1 d D P n 1 e n m D Q 1 d Š; where ˆ x ; y ; z Š are Lagrange mutpers. In order to obtan an approprate vaue of m, we de ne a sequence m ˆ 1; ; ¼; n such that D P n ˆ m D Q Accordng to Theorem 1, the average vaue of the sequence w be a reasonabe choce for parameter m. Thus we take m as m ˆ ˆ1 Note that D e n ˆ n m Xn kˆn ˆ ˆ1 1 nk 7 s 1 D P n D Q n C : 8 A e k ; n k 9 D P n 1 D e n m D Q m D d ˆ ˆ ; 1; ¼; n: 13 Substtutng Eqs. (11) and (1) nto Eq. (13), produces the system of near equatons " # 1 X m 1 mn ; ˆ k k kˆ ˆ D P n m D Q ˆ ; 1; ¼; n: 14 Eq. (14) can be further smp ed as " # m 1 ˆ D P n m D Q ˆ ˆ ; 1; ¼; n: (15) So ˆ ; 1; ¼; n can be obtaned by sovng the above system of near equatons, and e ; d ˆ ; 1; ¼; n can be obtaned from Eqs. (11) and (1). The exstng proof of souton of the near equaton system (15) w be gven n Secton Constraned optmzaton method () We can aso use the ªsquared dfference ntegraº norm, such that d E t ; F t ˆ Z1 E t F t dt; where E(t) and F(t) are n degree BeÂzer curves. The norm can be vewed geometrcay as the squared area between two BeÂzer curves. Therefore, mnmzaton of the norm w resut n the east encosng area between the two curves. We determne perturbatons e ; d ˆ ; 1; ¼; n ; de ned n Secton., by settng the optmzaton crteron O(e, d ) as " Z1 # O e ; d ˆ B n t e 1 Xn B n t d dt ˆ Mn
4 18 S.-M. Hu et a. / Computer-Aded Desgn 33 (1) 15±136 (a) Fg. 1. Compute new contro ponts by extrapoaton. and Lagrange functon s de ned by " Z1 # L ˆ B n t e 1 Xn B n t d dt 1 Xn D P n 1 e n m D Q 1 d Š 16 where ˆ x ; y ; z Š are Lagrange mutpers and m s de ned as n Secton.. By settng L e x ; L e y ; L e z ; L d x ; L ; d y L d z ; L L ; x y and L z (b) (c) to be zero, we have Xn e N 1 Xn 1 n ˆn Xn d N Xn m 1 ˆ ˆ ˆ ; 1; ¼; n; n 17 ˆ ˆ ; 1; ¼; n; 18 D P n 1 D e n m D Q m D d ˆ ˆ ; 1; ¼; n: 19 where N s de ned as N ˆ Z1 B n t B n t dt ; ˆ ; 1; ¼; n: Thus e, d and ˆ ; 1; ¼; n can be obtaned by sovng the near equaton system consstng of Eqs. (17)±(19). The exstence proof w be gven n Secton 3. (d).4. Cacuaton of new contro ponts R We are now n a poston to gve a recursve agorthm to compute new contro ponts of the merged curve R(t). Snce contro ponts of ^P u and ^Q v can be derved by subdvdng R(t) att ˆ ˆ m= 1 1 m (see Fg. 1), we have the foowng agorthm to compute new contro ponts by usng extrapoaton. Agorthm Set ˆ m= 1 1 m Set P ˆ P 1 e for ˆ ; 1; ¼; n For ˆ 1; ; ¼; n For ˆ ; 1 1; ¼; n Fg.. Mergng of two cubc BeÂzer curves. P ˆ 1 1= P 1 1 Next Next R k ˆ P k k for k ˆ ; 1; ¼; n: 1 1=P 1 Exampes. In Fg., we demonstrate mergng of cubc BeÂzer curves par usng the constraned optmzaton methods dscussed so far. Exampes (a) and (c) adopt the rst constraned optmzaton method, whe (b) and (d) use the
5 S.-M. Hu et a. / Computer-Aded Desgn 33 (1) 15± second optmzaton method n mergng. The orgna curves and contro poygons are shown as sod nes, whe the merged curves and assocated poygons are shown as dashed nes. We can see that the rst method works better n (a), whe the second method acheves ncer resuts n (d). 3. An exstence proof of souton In ths secton, we rst gve the exstence proof of system (15). By rewrtng Eqs. (7) and (13) n the matrx form, we get L ˆ h T h 1 T Ah 1 B and Ah 1 B ˆ 1 n whch hˆ e ; e 1 ; ¼;e n ;d ;d 1 ; ¼; d n T ;ˆ ; 1 ; ¼; n T ; A ˆ a n11 n1 de nte matrx, thus t's non-snguar. Ths competes the exstence proof of the system (15). Smary, by rewrtng Eq. (16) n matrx form, we get L ˆ h T Ch 1 T Ah 1 B : 5 here C s a symmetrc matrx de ned as C ˆ c n1 n1 n whch 8 ; for # # n and n, # n 1 1; >< or # # n and n, # n 1 1 c ˆ ; N ; for # ; # n >: N n1 n1 ; for n, ; # n 1 1 where N are de ned n Eq. (). Accordng to Eq. (16), h T Ch s aways postve for any h h ± ; therefore C s postve de nte. Dfferentatng Eq. (5) wth respect to h, we get (Eqs. (17) and (18) n matrx form) Ch 1 A T ˆ : 6 B ˆ b ; b 1 ; ¼; b n Š T ; Eq. (19) can aso be wrtten n matrx form where Ah 1 B ˆ : 8 ; for # 1, n or n 1 1 1, # n 1 1 >< 1 n ; for # ; # n; 1 $ n a ˆ n >: 1 n1 m ; for # # n; n 1 1 # # n n 1 7 b ˆ D P n m D Q : Dfferentatng Eq. (1) wth respect to h, we obtan (Eqs. (11) and (1) n the matrx form) h 1 A T ˆ 3 By substtutng Eq. (3) nto Eq. (), we obtan (Eq. (15) n matrx form) 1 AAT ˆ B: 4 Note that 1 ¼ p p ¼ ¼ p p p p ¼ ¼ p p p p p p ¼... A ˆ... ] ] p ¼ p p p p p p ¼ p p p ¼ p p p p p p ¼ p p C A p p p ¼ p p p p p p ¼ p p p n11 n1 s obvousy a row-fu-rank matrx. So 1 AAT s a postve Substtutng Eq. (6) nto Eq. (7), we obtan 1 AC1 A T ˆ B 8 because A s a row-fu-rank matrx and C 1 s symmetrc postve de nte, so 1 AC1 A T s non-snguar [13]. By substtutng Eq. (8) nto Eq. (6), we obtan h as h 1 C 1 A T AC 1 A T 1 B ˆ Ths competes the exstence proof of the near equaton system consstng of Eqs. (17)±(19). 4. Condtons for matchng orgna endponts and dervatves To match orgna endponts P and Q n, two constrants e ˆ ; d n ˆ shoud be added. Take the rst constraned optmzaton method for exampe, by dscusson smar to
6 13 S.-M. Hu et a. / Computer-Aded Desgn 33 (1) 15±136 (a) (b) Fg. 3. Mergng of two quartc BeÂzer curves wth endponts matchng. that n Secton., we obtan a system of near equatons 8 " # m 1 ˆ D P n m D Q ˆ >< ˆ ; 1; ¼; n 1 " # 1 X 1 1 m 1n mn ;n 1 n ˆ D n P m n D n Q >: ˆ kˆ k k Then, ˆ ; 1; ¼; n can be obtaned, and e ; d ˆ ; 1; ¼; n can be obtaned as foows 8 e ˆ e ˆ 1 >< ˆn d ˆ 1 1 ˆ >: d n ˆ 1 n ˆ 1; ¼; n n m ˆ ; 1; ¼; n 1 Further, constrants e 1 ˆ and d n1 ˆ coud be added to the above equaton system so that the merged curve woud match the orgna dervatves at the eft and rght ends of the par of BeÂzer curves. The system of near equatons woud become 8 " # m 1 ˆ D P n m D Q ˆ >< ˆ ; 1; ¼; n " # 1 X 1 1 m 1 mn ;n ˆ D P m D Q ˆ kˆ k k >: ˆ n 1; n Then e ; d ˆ ; 1; ¼; n can be obtaned as foows 8 e ˆ e 1 ˆ e ˆ 1 1 n ˆ ; ¼; n >< ˆn n d ˆ 1 1 m ˆ ; 1; ¼; n ˆ >: d n ˆ d n1 ˆ Smary, for the second constraned optmzaton method, matchng of orgna endponts and dervatves can be acheved by addng constrants e ˆ e 1 ˆ and d n ˆ d n1 ˆ to Lagrange functon (16). Perturbatons e, d can then be obtaned by sovng the correspondng system of near equatons. Exampes. In Fg. 3, we ustrate mergng resuts of two quartc BeÂzer curves usng the rst optmzaton method. Endponts matchng s mpemented n (a), and matchng of dervatves at endponts s acheved n (b). The orgna curves and contro poygons are shown as sod nes, and the merged curves and assocated poygons are shown as dashed nes. 5. Approxmate mergng wth ponts constrants A more genera probem of matchng endponts s for the merged curve to pass through target ponts, not ony endponts on the orgna BeÂzer curves. Ths can be acheved by addng ponts constraned condtons to the Lagrange functon. Gven that B n u P ˆ Xn B n v k Q ˆ Xn B n u P 1 e ˆ ; 1; ¼; ; B n v k Q 1 d k ˆ ; 1; ¼; m; where u ˆ ; 1; ¼; and v k ˆ ; 1; ¼; m are parameters of dfferent target ponts on P(u) and Q(v), respectvey, we have B n u e ˆ ˆ ; 1; ¼; ; 9
7 S.-M. Hu et a. / Computer-Aded Desgn 33 (1) 15± (a) (b) Note that unke the equaton system n Secton.3, the above near system s not sovabe when the tota number of constraned condtons n Eqs. (13), (9) and (3) exceeds the number of constraned varabes (e and d ). Generay speakng, the tota number of dfferent ponts constrants on P(u) and Q(v) 1 m 1 shoud not be arger than n 1 1 for a sovabe souton. Exampe. Ponts constrants can be used n preservng the geometrc detas of the orgna par of BeÂzer curves. In Fg. 4(a), two quntc BeÂzer curves are merged usng the second optmzaton method wth endponts matchng,.e. u ˆ :; v ˆ 1:; whe (b) ustrates the mergng resuts after three more ponts constrants at u ˆ :5; u 1 ˆ 1:; v ˆ :5 are added. 6. Constraned optmzaton wth degree eevaton Fg. 4. Mergng of two quntc BeÂzer curves wth ponts constrants. B n v k d ˆ k ˆ ; 1; ¼; m: 3 By addng the above-constraned condtons, we rewrte the Lagrange functon as L ˆ O e ; d 1 Xn 1 X n111 ˆ D P n 1 e n m D Q 1 d Š 1 Xm n11k1 kˆ B n u e B n v k d ; (31) where O(e, d ) s the approprate optmzaton crteron. Dfferentatng Eq. (31) wth respect to e ; d ˆ ; 1; ¼; n and ˆ ; 1; ¼; n 1 1 m 1 ; usng smar technques descrbed n Sectons. and.3, we obtan a near system of 3n 1 m equatons consstng of Eqs. (13), (9), (3) and O e ; d e 1 Xn ˆn 1 n n 1 X ˆ n111 B n u ˆ ˆ ; 1; ¼; n; (3) O e ; d d Xn ˆ m 1 1 Xm ˆ n111 B n v ˆ ˆ ; 1; ¼; n: 33 In ths secton, we w dscuss constraned optmzaton wth degree eevaton. Sometmes, two adacent BeÂzer curves are not sutabe to be approxmated by a new BeÂzer curve wth the same degree. For exampe, a cubc BeÂzer curve can have at most one n ecton pont, so mergng a par of BeÂzer curves that each has an n ecton pont nto a new cubc BeÂzer curve w not be a good dea. A reasonabe choce n ths case s to merge two cubc curves nto a curve of degree 4 or 5. In effect, better approxmaton can be acheved when two n degree BeÂzer curves P(u) and Q(v) are merged nto a new BeÂzer curve after the degree eevaton. We can thnk of the ncrease n approxmaton precson n ths way. Constraned optmzaton can be vewed as a two-step process: rsty, a set of BeÂzer curve pars ^P u and ^Q v are obtaned such that each par can be precsey merged nto a BeÂzer curve; secondy, the curves par wth the east optmzaton norm s ªchosenº and merged as the na approxmate merged curve. Snce each par of ^P u and ^Q v obtaned n the rst step can be eevated to a hgher degree and st satsfy the condton of precse mergng, the set of curve pars obtaned n the rst step of constraned optmzaton of degree n s ony a subset of those obtaned n the rst step for degree n 1 1 or above. In other words, optmzaton method w have more egbe curve pars to ªchoose fromº after degree eevaton. Consequenty better approxmate merged curve w be obtaned. We ustrate ths process by usng the second optmzaton method. Frst, we prove a theorem about degree eevaton. Theorem. P(u) s a n degree BeÂzer curve wth contro ponts P ˆ ; 1; ¼; n ; and D s the dfference operator, we have D P p n11 ˆ 1 D P n 1 1 n ˆ ; 1; ¼; n 1 1; 34
8 13 S.-M. Hu et a. / Computer-Aded Desgn 33 (1) 15±136 (a) (b) (c) (d) (e) (f) Fg. 5. Mergng of two quartc BeÂzer curves wth degree eevaton. D P p ˆ 1 D P n 1 1 ˆ ; 1; ¼; n where P p ˆ ; 1; ¼; n 1 1 are n 1 1 degree contro ponts after degree eevaton of P(u). Proof. Note that P p ˆ 1 P n n 1 1 P 1 ˆ ; 1; ¼; n 1 1; 36
9 S.-M. Hu et a. / Computer-Aded Desgn 33 (1) 15± (a) (b) (c) (d) (e) (f) Fg. 6. Mergng of two quartc BeÂzer curves wth degree eevaton.
10 134 S.-M. Hu et a. / Computer-Aded Desgn 33 (1) 15±136 Tabe 1 Comparson of mergng error after degree eevaton and D P p n11 ˆ X D P p ˆ X kˆ kˆ 1 k 1 k P p n111k ˆ ; 1; ¼; n 1 1; k 37 P p k ˆ ; 1; ¼; n 1 1: 38 k Substtutng Eq. (36) nto Eqs. (37) and (38) woud yed Eqs. (34) and (35). Ths competes the proof of Theorem. Theorem guarantees that the vaue of m (de ned n Eq. (8)) does not change after degree eevaton. Thus mergng methods dscussed n ths paper w proceed wth the same parameter m for BeÂzer curves par after degree eevaton. Gven two adacent BeÂzer curves of degree n, P(u) and Q(v), we rst consder the approxmate merged curve R(t) usng the second optmzaton method. Here, we de ne the mergng error E as n Secton.3, such that E R t ; m ˆ m E R t ; m Z1 Degree 4 Degree 5 Degree 6 Exampe (a) Exampe (d) ^P u P u Z1 du 1 ^Q v Q v dv; 39 where ^P u and ^Q v can be precsey merged nto R(t) wth subdvson parameter m= 1 1 m : After the degree eevaton of BeÂzer curves P(u) and Q(v), we obtan two new BeÂzer curves P p (u) and Q p (v) of degree n 1 1: Smary, we obtan n 1 1 degree BeÂzer curves ^P p u and ^Q p v after degree eevaton of ^P u and ^Q v : Therefore BeÂzer curves ^P p u and ^Q p v can aso be precsey merged nto R p (t), whch s the n 1 1 degree BeÂzer curve eevated from R(t). Examne the mergng Tabe Comparson of mergng error after degree eevaton m E R t ; m Degree 4 Degree 5 Degree 6 Exampe (a) Exampe (d) error of R p (t) wth Eq. (39), we have E R p t ; m ˆ ˆ Z1 Z1 ^P p u P p u Z1 du 1 ^Q p v Q p v dv ^P u P u Z1 du 1 ^Q v Q v dv ˆ E R t ; m On the other hand, for n 1 1 degree BeÂzer curves P p (u) and Q p (v), we compute an approxmate merged curve R t of degree n 1 1 wth parameter m usng the second optmzaton method. Accordng to the optmzaton crteron, we have E R t ; m ˆMn # E R p t ; m ˆE R t ; m 4 From Eq. (4) we can concude that R t resuts n ess mergng error E than ts n degree counterpart R(t), and thus acheves better approxmaton after degree eevaton. Smar resuts can be obtaned for the rst optmzaton method, wth the de nton of mergng error E as E R t ; m ˆXn e 1 d ; where e and d are perturbatons of contro ponts of ^P u and ^Q v from P(u) andq(v). Practca exampes of mergng wth both methods are provded n Fgs. 5 and 6. Exampes. Two mergng exampes of quartc cases are presented n Fg. 5 by usng the rst optmzaton method. Two orgna quartc BeÂzer curves are merged nto a new quartc BeÂzer curve each n (a) and (d), and are then merged nto BeÂzer curves of degree 5 and 6, respectvey, n (b), (c) and (e), (f) after degree eevaton. The orgna curves and contro poygons are shown as sod nes, and the merged curves and assocated poygons are shown as dashed nes. Tabe 1 contans the comparson resuts of mergng error of each exampe presented n Fg. 5. From Tabe 1 we can see that the mergng error E(R(t),m) decreases each tme after degree eevaton. Exampes. Two mergng exampes of quartc cases are presented n Fg. 6 by usng the second optmzaton method. Two orgna quartc BeÂzer curves are merged nto a new quartc BeÂzer curve each n (a) and (d), and are then merged nto BeÂzer curves of degree 5 and 6, respectvey, n (b), (c) and (e), (f) after degree eevaton. The orgna curves and contro poygons are shown as sod nes, and the merged curves and assocated poygons are shown as dashed nes. Tabe contans resuts of mergng error of each exampe
11 S.-M. Hu et a. / Computer-Aded Desgn 33 (1) 15± presented n Fg. 6. From Tabe we can see that mergng error E(R(t),m) decreases each tme after degree eevaton. 7. Remarks Remark 1. In genera, the dscusson n ths paper can appy to a pars of curves wth G onng. However, Theorem 1 prescrbes that the merged curve shoud have at east G 1 contnuty at the connectng pont. Therefore we recommend usng the methods presented n the paper to pars of BeÂzer curves wth at east G 1 onng n order to avod dscrepancy at connectng pont. Remark. Strcty speakng, the subdvson parameter m n Eq. (7) shoud be dependent on e ; d ˆ ; 1; ¼; n ; such that m ˆ P n 1 e n P n1 e n1 : Q 1 1 d 1 Q d Under such a choce, the Lagrange functon L w get too compcated. So we take m as n Eq. (8), where m s ndependent on e ; d ˆ ; 1; ¼; n : However, a better choce m can be obtaned from m so that the mergng error w get smaer. By usng m de ned n Eq. (8), the merged curve R(t) of P(u) and Q(v) can be computed rsty as n Secton. We hope to obtan better parameter m. Wth the unknown parameter m, we subdvde the merged curve R(t) nto two BeÂzer curves wth contro ponts S and T ˆ ; 1; ¼; n ; and de ne e ; d ˆ ; 1; ¼; n as e ˆ S P ˆ X ˆ d ˆ T Q ˆ Xn ˆ B m R P ˆ ; 1; ¼; n 41 B n m R 1 Q ˆ ; 1; ¼; n 4 Snce we wsh to obtan m such that optmzaton crteron O be mnmzed, we have do e ; d dm ˆ 43 where O e ; d s the approprate optmzaton crteron. Substtutng Eqs. (41) and (4) nto Eq. (43) w produce a non-near equaton of m. Practca agorthms can be derved to cacuate m by usng Newton's method. Remark 3. Mergng methods dscussed n ths paper can be further extended for cases of mutpe BeÂzer curves. Gven a sequence of consecutve BeÂzer curves, our obectve now s to obtan a new sequence of consecutve BeÂzer curves that best approxmates the orgna BeÂzer curves set. We ntend to mnmze the number of curves n the new sequence wth the condton that the mergng error stays wthn a requred error range E. It can proceed as foows. Frst we group the nta sequence nto adacent pars, eavng the snge BeÂzer curve at the end f the number of curves s odd. Then for each par, we appy optmzaton method wth endponts matchng and examne the mergng error. If the error exceeds a spec ed error bound e, the par w not be merged and becomes two curves n the resut sequence. Otherwse the par w be merged and a merged curves from these pars form a new sequence of BeÂzer curves (not necessary consecutve) for the next round of groupng. Repeat the process unt no two adacent BeÂzer curves can be found n the sequence, and we nay obtan the resut sequence of BeÂzer curves. Note that error bound of each step of mergng e must be carefuy spec ed such that accumuated mergng error of the process meets the requrement E. Acknowedgements Ths work was supported by the Natura Scence Foundaton of Chna (proect number 6994). We thank Prof. Tong-Guang Jn and Prof. Guo-Zhao Wang for encouragng ths work, and Dr Wen-Shen Xu for hs usefu suggestons. We especay thank the revewers for ther hepfu comments. References [1] Dannenberg L, Nowack H. Approxmate converson of surface representatons wth poynoma bases. Computer-Aded Geometrc Desgn 1985;:13±31. [] Hoschek J. Approxmate converson of spne curves. Computer- Aded Geometrc Desgn 1987;4:59±66. [3] Watkns MA, Worsey AJ. Degree reducton of BeÂzer curves. Computer-Aded Desgn 1988;:398±45. [4] Lachance MA. Chebyshev economzaton for parametrc surfaces. Computer-Aded Geometrc Desgn 1988;5:19±8. [5] Eck M. Degree reducton of BeÂzer curves. Computer-Aded Geometrc Desgn 1993;1:37±51. [6] Eck M. Least square degree reducton of BeÂzer curves. Computer- Aded Desgn 1995;7:845±51. [7] Bogack P, Wensten SE, Xu Y. Degree reducton of BeÂzer curves by unform approxmaton wth endpont nterpoaton. Computer-Aded Desgn 1995;7:651±61. [8] Hu Sh-Mn, Jn Tong-Guang. Approxmate degree reducton of BeÂzer curves. Proceedngs of Symposum on Computatona Geometry, Hed n Hangzhou, Chna, 199. p. 11±6 (aso n: Tsnghua Scence and Technoogy 1998;3:997±1). [9] Hu Sh-Mn, Wang Guo-Zhao, Jn Tong-Guang. Propertes of two types of generazed Ba curves. Computer-Aded Desgn 1996;8:18±. [1] Peg L. Agorthm for degree reducton of B-spne curves. Computer-Aded Desgn 1995;7:11±1. [11] Hu Sh-Mn, Zheng Guo-Qn, Sun Ja-Guang. Approxmate degree reducton of rectanguar Bezer surfaces. Chnese Journa of Software Research 1997;4:353±61. [1] Hu Sh-Mn, Zuo Zheng, Sun Ja-Guang. Approxmate degree reducton of tranguar Bezer surfaces. Tsnghua Scence and Technoogy 1998;3:11±4. [13] Zhhao Cao, Zhang Yude, et a. Matrx computaton and equaton sovng. Chna: Advanced Educaton Press, 1984.
12 136 S.-M. Hu et a. / Computer-Aded Desgn 33 (1) 15±136 Sh-Mn Hu s an assocate professor n the Department of Computer Scence and Technoogy at Tsnghua Unversty. He receved the BS n mathematcs from the Jn Unversty of Chna n 199, the MS and the PhD n CAGD and Computer Graphcs from Zheang Unversty of Chna n 1993 and 1996, he nshed hs postdoctora research at Tsnghua Unversty n Hs current research nterests are n computer-aded geometrc desgn, computer graphcs and content-based mages ndexng. Tao Ju s an Undergraduate at Tsnghua Unversty for hs degree both n Computer Scence and Engsh Language. Hs current research nterests are n computer-aded geometrc desgn. Ruo-Feng Tong s an assocate professor n the Department of Computer Scence at Zheang Unversty. He receved the BS n Apped Mathematcs from Zheang Unversty of Chna n 1991, and the PhD n CAGD and Computer Graphcs from the Zheang unversty of Chna n Hs current research nterests are n computer-aded geometrc desgn and computer graphcs. Ja-Guang Sun s a professor n the Department of Computer Scence and Technoogy at Tsnghua Unversty. He s aso the Drector of Natona CAD Engneerng Center at Tsnghua Unversty and Academcan of the Chnese Academy of Engneerng. He receved the BS n Computer Scence from Tsnghua Unversty n 197. Hs current research nterests are n computer-aded geometrc desgn, computer graphcs and Product Data Management.
DEGREE REDUCTION OF BÉZIER CURVES USING CONSTRAINED CHEBYSHEV POLYNOMIALS OF THE SECOND KIND
ANZIAM J. 45(003), 195 05 DEGREE REDUCTION OF BÉZIER CURVES USING CONSTRAINED CHEBYSHEV POLYNOMIALS OF THE SECOND KIND YOUNG JOON AHN 1 (Receved 3 August, 001; revsed 7 June, 00) Abstract In ths paper
More informationA finite difference method for heat equation in the unbounded domain
Internatona Conerence on Advanced ectronc Scence and Technoogy (AST 6) A nte derence method or heat equaton n the unbounded doman a Quan Zheng and Xn Zhao Coege o Scence North Chna nversty o Technoogy
More informationResearch on Complex Networks Control Based on Fuzzy Integral Sliding Theory
Advanced Scence and Technoogy Letters Vo.83 (ISA 205), pp.60-65 http://dx.do.org/0.4257/ast.205.83.2 Research on Compex etworks Contro Based on Fuzzy Integra Sdng Theory Dongsheng Yang, Bngqng L, 2, He
More informationNote 2. Ling fong Li. 1 Klein Gordon Equation Probablity interpretation Solutions to Klein-Gordon Equation... 2
Note 2 Lng fong L Contents Ken Gordon Equaton. Probabty nterpretaton......................................2 Soutons to Ken-Gordon Equaton............................... 2 2 Drac Equaton 3 2. Probabty nterpretaton.....................................
More informationThe Quadratic Trigonometric Bézier Curve with Single Shape Parameter
J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma
More informationMARKOV CHAIN AND HIDDEN MARKOV MODEL
MARKOV CHAIN AND HIDDEN MARKOV MODEL JIAN ZHANG JIANZHAN@STAT.PURDUE.EDU Markov chan and hdden Markov mode are probaby the smpest modes whch can be used to mode sequenta data,.e. data sampes whch are not
More informationCyclic Codes BCH Codes
Cycc Codes BCH Codes Gaos Feds GF m A Gaos fed of m eements can be obtaned usng the symbos 0,, á, and the eements beng 0,, á, á, á 3 m,... so that fed F* s cosed under mutpcaton wth m eements. The operator
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
More informationNeural network-based athletics performance prediction optimization model applied research
Avaabe onne www.jocpr.com Journa of Chemca and Pharmaceutca Research, 04, 6(6):8-5 Research Artce ISSN : 0975-784 CODEN(USA) : JCPRC5 Neura networ-based athetcs performance predcton optmzaton mode apped
More informationBezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0
Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the
More informationwe have E Y x t ( ( xl)) 1 ( xl), e a in I( Λ ) are as follows:
APPENDICES Aendx : the roof of Equaton (6 For j m n we have Smary from Equaton ( note that j '( ( ( j E Y x t ( ( x ( x a V ( ( x a ( ( x ( x b V ( ( x b V x e d ( abx ( ( x e a a bx ( x xe b a bx By usng
More informationL-Edge Chromatic Number Of A Graph
IJISET - Internatona Journa of Innovatve Scence Engneerng & Technoogy Vo. 3 Issue 3 March 06. ISSN 348 7968 L-Edge Chromatc Number Of A Graph Dr.R.B.Gnana Joth Assocate Professor of Mathematcs V.V.Vannaperuma
More informationNested case-control and case-cohort studies
Outne: Nested case-contro and case-cohort studes Ørnuf Borgan Department of Mathematcs Unversty of Oso NORBIS course Unversty of Oso 4-8 December 217 1 Radaton and breast cancer data Nested case contro
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationPredicting Model of Traffic Volume Based on Grey-Markov
Vo. No. Modern Apped Scence Predctng Mode of Traffc Voume Based on Grey-Marov Ynpeng Zhang Zhengzhou Muncpa Engneerng Desgn & Research Insttute Zhengzhou 5005 Chna Abstract Grey-marov forecastng mode of
More informationImage Classification Using EM And JE algorithms
Machne earnng project report Fa, 2 Xaojn Sh, jennfer@soe Image Cassfcaton Usng EM And JE agorthms Xaojn Sh Department of Computer Engneerng, Unversty of Caforna, Santa Cruz, CA, 9564 jennfer@soe.ucsc.edu
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationApplication of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationLECTURE 21 Mohr s Method for Calculation of General Displacements. 1 The Reciprocal Theorem
V. DEMENKO MECHANICS OF MATERIALS 05 LECTURE Mohr s Method for Cacuaton of Genera Dspacements The Recproca Theorem The recproca theorem s one of the genera theorems of strength of materas. It foows drect
More informationLower bounds for the Crossing Number of the Cartesian Product of a Vertex-transitive Graph with a Cycle
Lower bounds for the Crossng Number of the Cartesan Product of a Vertex-transtve Graph wth a Cyce Junho Won MIT-PRIMES December 4, 013 Abstract. The mnmum number of crossngs for a drawngs of a gven graph
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationNP-Completeness : Proofs
NP-Completeness : Proofs Proof Methods A method to show a decson problem Π NP-complete s as follows. (1) Show Π NP. (2) Choose an NP-complete problem Π. (3) Show Π Π. A method to show an optmzaton problem
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationPerron Vectors of an Irreducible Nonnegative Interval Matrix
Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationThe lower and upper bounds on Perron root of nonnegative irreducible matrices
Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College
More informationValuated Binary Tree: A New Approach in Study of Integers
Internatonal Journal of Scentfc Innovatve Mathematcal Research (IJSIMR) Volume 4, Issue 3, March 6, PP 63-67 ISS 347-37X (Prnt) & ISS 347-34 (Onlne) wwwarcournalsorg Valuated Bnary Tree: A ew Approach
More informationDeriving the Dual. Prof. Bennett Math of Data Science 1/13/06
Dervng the Dua Prof. Bennett Math of Data Scence /3/06 Outne Ntty Grtty for SVM Revew Rdge Regresson LS-SVM=KRR Dua Dervaton Bas Issue Summary Ntty Grtty Need Dua of w, b, z w 2 2 mn st. ( x w ) = C z
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationBoundary Value Problems. Lecture Objectives. Ch. 27
Boundar Vaue Probes Ch. 7 Lecture Obectves o understand the dfference between an nta vaue and boundar vaue ODE o be abe to understand when and how to app the shootng ethod and FD ethod. o understand what
More informationCubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two
World Academy of Scence Engneerng and echnology Internatonal Journal of Mathematcal and omputatonal Scences Vol: No:0 00 ubc rgonometrc B-Splne Appled to Lnear wo-pont Boundary Value Problems of Order
More informationAssociative Memories
Assocatve Memores We consder now modes for unsupervsed earnng probems, caed auto-assocaton probems. Assocaton s the task of mappng patterns to patterns. In an assocatve memory the stmuus of an ncompete
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationThe line method combined with spectral chebyshev for space-time fractional diffusion equation
Apped and Computatona Mathematcs 014; 3(6): 330-336 Pubshed onne December 31, 014 (http://www.scencepubshnggroup.com/j/acm) do: 10.1164/j.acm.0140306.17 ISS: 3-5605 (Prnt); ISS: 3-5613 (Onne) The ne method
More informationON AUTOMATIC CONTINUITY OF DERIVATIONS FOR BANACH ALGEBRAS WITH INVOLUTION
European Journa of Mathematcs and Computer Scence Vo. No. 1, 2017 ON AUTOMATC CONTNUTY OF DERVATONS FOR BANACH ALGEBRAS WTH NVOLUTON Mohamed BELAM & Youssef T DL MATC Laboratory Hassan Unversty MORO CCO
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationOn the Power Function of the Likelihood Ratio Test for MANOVA
Journa of Mutvarate Anayss 8, 416 41 (00) do:10.1006/jmva.001.036 On the Power Functon of the Lkehood Rato Test for MANOVA Dua Kumar Bhaumk Unversty of South Aabama and Unversty of Inos at Chcago and Sanat
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationXin Li Department of Information Systems, College of Business, City University of Hong Kong, Hong Kong, CHINA
RESEARCH ARTICLE MOELING FIXE OS BETTING FOR FUTURE EVENT PREICTION Weyun Chen eartment of Educatona Informaton Technoogy, Facuty of Educaton, East Chna Norma Unversty, Shangha, CHINA {weyun.chen@qq.com}
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationMultispectral Remote Sensing Image Classification Algorithm Based on Rough Set Theory
Proceedngs of the 2009 IEEE Internatona Conference on Systems Man and Cybernetcs San Antono TX USA - October 2009 Mutspectra Remote Sensng Image Cassfcaton Agorthm Based on Rough Set Theory Yng Wang Xaoyun
More informationCOXREG. Estimation (1)
COXREG Cox (972) frst suggested the modes n whch factors reated to fetme have a mutpcatve effect on the hazard functon. These modes are caed proportona hazards (PH) modes. Under the proportona hazards
More informationSupplementary Material: Learning Structured Weight Uncertainty in Bayesian Neural Networks
Shengyang Sun, Changyou Chen, Lawrence Carn Suppementary Matera: Learnng Structured Weght Uncertanty n Bayesan Neura Networks Shengyang Sun Changyou Chen Lawrence Carn Tsnghua Unversty Duke Unversty Duke
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationCubic Trigonometric Rational Wat Bezier Curves
Cubc Trgonometrc Ratonal Wat Bezer Curves Urvash Mshra Department of Mathematcs Mata Gujr Mahla Mahavdyalaya Jabalpur Madhya Pradesh Inda Abstract- A new knd of Ratonal cubc Bézer bass functon by the blendng
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationComparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method
Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method
More informationExample: Suppose we want to build a classifier that recognizes WebPages of graduate students.
Exampe: Suppose we want to bud a cassfer that recognzes WebPages of graduate students. How can we fnd tranng data? We can browse the web and coect a sampe of WebPages of graduate students of varous unverstes.
More informationVQ widely used in coding speech, image, and video
at Scalar quantzers are specal cases of vector quantzers (VQ): they are constraned to look at one sample at a tme (memoryless) VQ does not have such constrant better RD perfomance expected Source codng
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationCOMPLEX NUMBERS AND QUADRATIC EQUATIONS
COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not
More informationSupporting Information
Supportng Informaton The neural network f n Eq. 1 s gven by: f x l = ReLU W atom x l + b atom, 2 where ReLU s the element-wse rectfed lnear unt, 21.e., ReLUx = max0, x, W atom R d d s the weght matrx to
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationFormulas for the Determinant
page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use
More informationGreyworld White Balancing with Low Computation Cost for On- Board Video Capturing
reyword Whte aancng wth Low Computaton Cost for On- oard Vdeo Capturng Peng Wu Yuxn Zoe) Lu Hewett-Packard Laboratores Hewett-Packard Co. Pao Ato CA 94304 USA Abstract Whte baancng s a process commony
More informationQUARTERLY OF APPLIED MATHEMATICS
QUARTERLY OF APPLIED MATHEMATICS Voume XLI October 983 Number 3 DIAKOPTICS OR TEARING-A MATHEMATICAL APPROACH* By P. W. AITCHISON Unversty of Mantoba Abstract. The method of dakoptcs or tearng was ntroduced
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationA new construction of 3-separable matrices via an improved decoding of Macula s construction
Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationMaejo International Journal of Science and Technology
Maejo Int. J. Sc. Technol. () - Full Paper Maejo Internatonal Journal of Scence and Technology ISSN - Avalable onlne at www.mjst.mju.ac.th Fourth-order method for sngularly perturbed sngular boundary value
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationLecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem.
prnceton u. sp 02 cos 598B: algorthms and complexty Lecture 20: Lft and Project, SDP Dualty Lecturer: Sanjeev Arora Scrbe:Yury Makarychev Today we wll study the Lft and Project method. Then we wll prove
More informationOn the Multicriteria Integer Network Flow Problem
BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 5, No 2 Sofa 2005 On the Multcrtera Integer Network Flow Problem Vassl Vasslev, Marana Nkolova, Maryana Vassleva Insttute of
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
More informationCS 468 Lecture 16: Isometry Invariance and Spectral Techniques
CS 468 Lecture 16: Isometry Invarance and Spectral Technques Justn Solomon Scrbe: Evan Gawlk Introducton. In geometry processng, t s often desrable to characterze the shape of an object n a manner that
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationNumerical Investigation of Power Tunability in Two-Section QD Superluminescent Diodes
Numerca Investgaton of Power Tunabty n Two-Secton QD Superumnescent Dodes Matta Rossett Paoo Bardea Ivo Montrosset POLITECNICO DI TORINO DELEN Summary 1. A smpfed mode for QD Super Lumnescent Dodes (SLD)
More informationA parametric Linear Programming Model Describing Bandwidth Sharing Policies for ABR Traffic
parametrc Lnear Programmng Mode Descrbng Bandwdth Sharng Poces for BR Traffc I. Moschoos, M. Logothets and G. Kokknaks Wre ommuncatons Laboratory, Dept. of Eectrca & omputer Engneerng, Unversty of Patras,
More informationACTM State Calculus Competition Saturday April 30, 2011
ACTM State Calculus Competton Saturday Aprl 30, 2011 ACTM State Calculus Competton Sprng 2011 Page 1 Instructons: For questons 1 through 25, mark the best answer choce on the answer sheet provde Afterward
More informationSimulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests
Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth
More informationChapter 6. Rotations and Tensors
Vector Spaces n Physcs 8/6/5 Chapter 6. Rotatons and ensors here s a speca knd of near transformaton whch s used to transforms coordnates from one set of axes to another set of axes (wth the same orgn).
More informationANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)
Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of
More informationAn Effective Space Charge Solver. for DYNAMION Code
A. Orzhehovsaya W. Barth S. Yaramyshev GSI Hemhotzzentrum für Schweronenforschung (Darmstadt) An Effectve Space Charge Sover for DYNAMION Code Introducton Genera space charge agorthms based on the effectve
More informationDevelopment of whole CORe Thermal Hydraulic analysis code CORTH Pan JunJie, Tang QiFen, Chai XiaoMing, Lu Wei, Liu Dong
Deveopment of whoe CORe Therma Hydrauc anayss code CORTH Pan JunJe, Tang QFen, Cha XaoMng, Lu We, Lu Dong cence and technoogy on reactor system desgn technoogy, Nucear Power Insttute of Chna, Chengdu,
More informationSPATIAL KINEMATICS OF GEARS IN ABSOLUTE COORDINATES
SATIAL KINEMATICS OF GEARS IN ABSOLUTE COORDINATES Dmtry Vasenko and Roand Kasper Insttute of Mobe Systems (IMS) Otto-von-Guercke-Unversty Magdeburg D-39016, Magdeburg, Germany E-ma: Dmtr.Vasenko@ovgu.de
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationProblem Set 9 Solutions
Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem
More informationTR/95 February Splines G. H. BEHFOROOZ* & N. PAPAMICHAEL
TR/9 February 980 End Condtons for Interpolatory Quntc Splnes by G. H. BEHFOROOZ* & N. PAPAMICHAEL *Present address: Dept of Matematcs Unversty of Tabrz Tabrz Iran. W9609 A B S T R A C T Accurate end condtons
More informationDIOPHANTINE EQUATIONS WITH BINOMIAL COEFFICIENTS AND PERTURBATIONS OF SYMMETRIC BOOLEAN FUNCTIONS
DIOPHANTINE EQUATIONS WITH BINOMIAL COEFFICIENTS AND PERTURBATIONS OF SYMMETRIC BOOLEAN FUNCTIONS FRANCIS N CASTRO, OSCAR E GONZÁLEZ, AND LUIS A MEDINA Abstract Ths work presents a study of perturbatons
More informationn-step cycle inequalities: facets for continuous n-mixing set and strong cuts for multi-module capacitated lot-sizing problem
n-step cyce nequates: facets for contnuous n-mxng set and strong cuts for mut-modue capactated ot-szng probem Mansh Bansa and Kavash Kanfar Department of Industra and Systems Engneerng, Texas A&M Unversty,
More informationQuantum Runge-Lenz Vector and the Hydrogen Atom, the hidden SO(4) symmetry
Quantum Runge-Lenz ector and the Hydrogen Atom, the hdden SO(4) symmetry Pasca Szrftgser and Edgardo S. Cheb-Terrab () Laboratore PhLAM, UMR CNRS 85, Unversté Le, F-59655, France () Mapesoft Let's consder
More informationNice plotting of proteins II
Nce plottng of protens II Fnal remark regardng effcency: It s possble to wrte the Newton representaton n a way that can be computed effcently, usng smlar bracketng that we made for the frst representaton
More informationRoot Structure of a Special Generalized Kac- Moody Algebra
Mathematcal Computaton September 04, Volume, Issue, PP8-88 Root Structu of a Specal Generalzed Kac- Moody Algebra Xnfang Song, #, Xaox Wang Bass Department, Bejng Informaton Technology College, Bejng,
More information3. Stress-strain relationships of a composite layer
OM PO I O U P U N I V I Y O F W N ompostes ourse 8-9 Unversty of wente ng. &ech... tress-stran reatonshps of a composte ayer - Laurent Warnet & emo Aerman.. tress-stran reatonshps of a composte ayer Introducton
More informationApplied Mathematics Letters
Appled Matheatcs Letters 2 (2) 46 5 Contents lsts avalable at ScenceDrect Appled Matheatcs Letters journal hoepage: wwwelseverco/locate/al Calculaton of coeffcents of a cardnal B-splne Gradr V Mlovanovć
More informationMEM 255 Introduction to Control Systems Review: Basics of Linear Algebra
MEM 255 Introducton to Control Systems Revew: Bascs of Lnear Algebra Harry G. Kwatny Department of Mechancal Engneerng & Mechancs Drexel Unversty Outlne Vectors Matrces MATLAB Advanced Topcs Vectors A
More informationSupport Vector Machines. Vibhav Gogate The University of Texas at dallas
Support Vector Machnes Vbhav Gogate he Unversty of exas at dallas What We have Learned So Far? 1. Decson rees. Naïve Bayes 3. Lnear Regresson 4. Logstc Regresson 5. Perceptron 6. Neural networks 7. K-Nearest
More informationModeling curves. Graphs: y = ax+b, y = sin(x) Implicit ax + by + c = 0, x 2 +y 2 =r 2 Parametric:
Modelng curves Types of Curves Graphs: y = ax+b, y = sn(x) Implct ax + by + c = 0, x 2 +y 2 =r 2 Parametrc: x = ax + bxt x = cos t y = ay + byt y = snt Parametrc are the most common mplct are also used,
More information