NiraDynandElzaFarkhi. The application of spline subdivision schemes to data consisting of convex compact sets,

Size: px
Start display at page:

Download "NiraDynandElzaFarkhi. The application of spline subdivision schemes to data consisting of convex compact sets,"

Transcription

1 Splne Subdvson Schemes for Convex Compact Sets NraDynandElzaFarkh The applcaton of splne subdvson schemes to data consstng of convex compact sets, wth addton replaced by Mnkowsk sums of sets, s nvestgated. These methods generate n the lmt set-valued functons, whch can be expressed explctly n terms of lnear combnatons of nteger shfts of B-splnes wth the ntal data as coecents. The subdvson technques are used to conclude that these lmt set-valued splne functons have shape preservng propertes smlar to those of the usual splne functons. Ths extenson of subdvson methods from the scalar settng to the set-valued case has applcaton n the approxmate reconstructon of 3-D bodes from nte collectons of ther parallel cross-sectons. Keywords: convex sets, support functons, Mnkowsk addton, set-valued functons, splne subdvson, shape preservaton, approxmaton. 1 Introducton Subdvson schemes are recursve methods for the generaton of smooth functons from dscrete data. By these methods at each recurson step, new dscrete values on a ner grd are computed by weghted sums of the already exstng dscrete values. In the lmt of the recursve process, data s dened on a dense set of ponts. Consderng ths data as functon values, under certan condtons, a lmt contnuous functon s dened by ths process. The theory of subdvson processes s presented n e.g. [4], [6]. In ths work we apply a class of subdvson methods wth postve weghts to data consstng of convex compact sets, replacng the addton by Mnkowsk sums of sets, and obtan n the lmt set-valued functons. Ths extenson of subdvson methods from the scalar settng to the set-valued case has applcaton n the approxmate reconstructon of 3-D bodes from nte collectons of ther parallel cross-sectons. The class of subdvson methods consdered here conssts of methods whch generate splne functons n the scalar settng. These methods have shape preservng propertes and approxmaton propertes n the set-valued case, smlar to those n the scalar case. The proof of the shape preservng propertes of the lmt set-valued functons reles on the subdvson technque. Yet the lmt multfuncton (set-valued functon) has a smple explct form n terms of the ntal sets and nteger shfts of B-splnes, and there s no need to compute t recursvely. Ths explct form also yelds the smoothness and the approxmaton propertes of the lmt set-valued functon. Among the splne subdvson schemes, there s a scheme whch generates n the lmt pecewse lnear nterpolants to the gven data of sets. In [11] pecewse lnear approxmaton to convex set-valued functons s studed. Ths approxmaton conssts of a pecewswe lnear nterpolant to samples of a multfuncton, wth addton replaced by Mnkowsk sums. Thus, one can regard the present paper as an extenson of [11]. 1

2 The mathematcal tools used for analysng set-valued functons nclude the support functon technque for descrbng convex compact sets (see e.g. [10]) and methods of embeddng the cone of convex compact subsets of R n n a lnear normed space, wth addton dened as the Mnkowsk sum of sets ([2], [7], [8], [9]). In any such lnear normed space, we ntroduce a partal order generated by the set ncluson order n the cone of convex compact sets. Thus amultfuncton wth convex compact mages from R to R n s consdered as an abstract functon wth values n a partally ordered normed lnear space. Monotoncty and convexty of such abstract functons are easly expressed by the postvty of ther rst and second nte derences. The paper s organsed as follows. In secton 2 basc facts about convex compact sets, ther support functons and ther embeddng n a normed lnear space wth a partal order, are presented. Secton 3 presents a smple example of a shape-preservng set-valued subdvson scheme. The man results about splne subdvson methods appled to convex compact sets are derved n Secton 4. The explct form of the lmt multvalued splne functon s obtaned together wth ts approxmaton and shape-preservng propertes. 2 Prelmnares Denote by C(R n ) the cone of nonempty convex compact subsets of R n : Recall the dentons of Mnkowsk sum and multplcaton by scalars of sets A B 2C(R n ): A + B = f a + b j a 2 A b 2 B g A = f a j a 2 A g: Snce the technque of support functons s central to ths text, we recall the denton and the basc propertes of these functons (see e.g. [7], [10]). For a set A 2C(R n ) ts support functon (A ) :R n! R s dened as follows: (A l) = maxhl a l a2a 2 Rn where h s the Eucldean nner product. Note that for any xedl 2 R n (A l) s nte. The followng propertes of are well known ([10]): 1 o : (A l) = (A l) 0: 2 o : (A l 1 + l 2 ) (A l 1 )+ (A l 2 ): 3 o : (A ) s Lpschtz contnuous: wth a constant kak = max kak where kk s the a2a Eucldean norm, namely j (A l 1 ) ; (A l 2 )jkakkl 1 ; l 2 k,forl 1 l 2 2 R n. 4 o : A scalar functon : R n! R s a support functon of a convex compact set t satses 1 o 2 o (see e.g. [10], Theorem 13.2 and ts corollares). 5 o : (A + B ) = (A )+ (B ): 6 o : (A ) = (A ) 0 7 o : A B () (A l) (B l) for each l 2 R n : 2

3 8 o : The Hausdor dstance between two sets A B 2C(R n ) s gven, n terms of the support functons of these two sets, by where S n;1 s the unt sphere n R n : haus(a B) = max j (A l) ; (B l)j l2s n;1 Note that by 1 o and 2 o, (A ) s postvely homogeneous and convex. Also, functons satsfyng 1 o,2 o are called sublnear. As we noted n 3 o, (see e.g. [10]), sublnear functons dened on all R n are Lpschtz contnuous. The next proposton wll be used n what follows. Proposton 2.1 Let A 1 A 2 B 1 B 2 2C(R n ) and A 1 B 1 : Then the equalty A 1 + B 2 = A 2 + B 1 mples A 2 B 2 : Proof: The proof follows from propertes 5 o 7 o : A 1 + B 2 = A 2 + B 1 () (A 1 )+ (B 2 ) = (A 2 )+ (B 1 ): From A 1 B 1 follows (A 1 ) (B 1 ): Ths, combned wth the above equalty mples (B 2 ) (A 2 ) =) B 2 A 2 : There are varous ways to construct a lnear normed vector space D(R n ) relatve to the Mnkowsk sum, n whch the cone C(R n )sembedded by anembeddng J : C(R n )! D(R n ) wth the followng propertes(see e.g. [7], [9], [8]): () J(A + B) =J(A)+J(B): () J(A) =J(A) 0: () J(A) =J(B) () A = B: (v) kj(a) ; J(B)k =haus(a B). These propertes mply J(f0g) =, where s the zero element ofd(r n ). A smple embeddng s J(A) = (A ). It s easy to check, usng the above stated propertes of support functons, that ths embeddng has the requred four propertes. Havng an embeddng of the cone C(R n )nto a lnear normed vector space D(R n ), we can ntroduce a partal order n D(R n ) and therefore n C(R n ). For that, we dene the followng cone n D(R n ): K = n C 2D(R n )j C = J(A) ; J(B) A B2C(R n ) A B The cone K determnes the followng partal order n D(R n ): For A B 2D(R n ) A B () B ; A 2K: (2) Remark 2.1 Here are some observatons regardng the partal order ntroduced above: 1. By () and () K s a convex cone. o (1) 3

4 2. Note that (1) does not depend on the choce of the sets A and B n the followng sense: f C = J(A 1 ) ; J(B 1 )=J(A 2 ) ; J(B 2 ) where A 1 B 1 then by () and () we get A 1 + B 2 = A 2 + B 1 and by Proposton 2.1 A 2 B 2 : 3. By the prevous observaton, the order dened nd(r n ) by (2) nduces the regular ncluson order n C(R n ) namely for A B 2C(R n ) A B () J(A) J(B): Thus we denote for A B 2C(R n ) A B A B. 4. Our denton of the postve cone K n D(R n ) concdes wth the postve cone n some concrete lnear spaces n whch C(R n ) s embedded, lke the space ofthepars of convex compact sets and the space of derences of support functons (see e.g. [8]). Remark 2.2 It follows from the facts that K s a cone and D(R n ) s a lnear space that for A B (a) ;A ;B (b) A + C B + C for every C 2D(R n ) (c) If C D, then A + C B + D for 0: The prevous remark justes the notons of postve and negatve elements of D(R n ). The element A 2D(R n )scallednonnegatve when A.e.A 2K: If A A 6= then A s called postve. The element B 2D(R n )scallednonpostve f ;B s nonnegatve and B s negatve when ;B s postve. We wll call a convex compact set A nonnegatve when J(A),.e. 0 2 A: A s postve 0 2 A and A 6= f0g. We are nterested n ths work n set-valued mappngs from R to C(R n ) (called also multmaps or multfunctons), and n partcular n multmaps of the form F (t) = N =1 A f (t) (3) where A 2C(R n ) f : R 7! R f (t) 0forallt 2 R. LetS be the cone of maps of the form (3). We saythatf 2Ss C k f n (3) f 2 C k for =1 ::: N. Denton 2.1 A mappng F : R!D(R n ) s called (A) monotone ncreasng f t 1 t 2 mples F (t 1 ) F (t 2 ) (.e. F (t 2 ) ; F (t 1 ) 2K). (B) monotone decreasng f F (;t) s monotone ncreasng. (C) convex f F (t 1 +(1; )t 2 ) F (t 1 )+(1; )F (t 2 ) for each 2 [0 1] t 1 t 2 2 R: (4) 4

5 (D) concave f F (;t) s convex. Denton 2.2 Dene for a gven functon F : R!D(R n ) the k-th forward nte derence at the pont t wth a step h>0 Forasequence ff g we dene k h F (t) = k j=0 ( k F ) = (;1) k;j k j k j=0! (;1) k;j k j F (t + jh):! F +j : In case 1 F s nonnegatve for all t h > 0 (), the functon (sequence) s monotone ncreasng. If 2 F s nonpostve for all t h > 0 (), the functon (sequence) s called convex. Remark 2.3 For the sake of smplcty we sometmes dentfy the set A 2 C(R n ) wth ts embedded mage J(A) and for sets A B 2C(R n ) we denote by A ; B the derence J(A) ; J(B). Geometrcally, a monotone ncreasng set-valued map F : R!C(R n ) has a growng mage as the argument t ncreases. Formally, F s monotone when for a gven h>0 the rst derence 1 h (t) =F (t + h) ; F (t) s of a constant sgn. Smlarly, the convexty means that the second derence 2 hf (t) =F (t)+f(t +2h) ; 2F (t + h) s nonpostve (.e. F (t)+f(t +2h) 2F (t + h)). The nequalty (4) s opposte to the common denton of convex scalar functons. We choose t ths way n order to ensure the convexty of the graph of F: Clearly, n the case of a convex map F : R!C(R n ) F (t 1 +(1; )t 2 ) F (t 1 )+(1; )F (t 2 ) for each 2 [0 1] t 1 t 2 2 R whch means that the graph of F s a convex set n R n+1 : The last ncluson and property 7 o of support functons mply that for each gven drecton l the support functon (F () l) of the convex multmap F s a concave scalar functon. 3 Chakn Subdvson Scheme for Convex Compact Sets Let F 0 =0 ::: N be convex compact sets n Rn : We seek a set-valued functon n S whch has a smlar structure to the the pecewse lnear multfuncton F : R! R n dened by F ( + ) =(1; )F 0 + F =0 1 ::: N ; 1 but s smoother n the sense of (3). Note that every pecewse lnear multfuncton s n S and s C 0. 5

6 Consder the followng teratve procedure of reconstructng F known as Chakn algorthm when appled to scalar functons ([5], [6]). At level k +1 k 0 of the procedure we calculate for =0 ::: 2 k (N ; 1) the sets 2 = 1 4 F k F k (5) 2 = 3 4 F k F k : (6) Thanks to the postvty of the coecents, we obtan convex compact sets at each stage of the process. Moreover, snce the coecents form a convex combnaton, the scheme has two notceable propertes: 1. It preserves monotoncty,.e. If for all F k ;1 F k then 2;1 2 2 for all. Ths follows drectly from the equaltes 2 ; 2;1 = 1 4 (F k ; F;1) k 2 ; 2 = 1 2 (F k ; F k ) (7) whch means that the rst nte derences reman n K at each stage of the process f they are n K for k =0: 2. It preserves convexty, namely the second derences reman n K at each teraton f they are n t for k = 0, snce ; 2 2 = 1 4 (F k + F k +2 ; 2F) k 2;1 + 2 ; 2 2 = 1 4 (F k ;1 + F k ; 2F k ): Atthek-th teraton (k 1) we construct a pecewse lnear multmap F k :[0 N;1+2 ;k ]! C(R n ) satsfyng F k (t) = tk ; t t k ; t k ;1 F k ;1 + t ; tk ;1 t k ; t k F k tk ;1 t t k : ;1 Lemma 3.1 The sequence of set-valued functons ff k (t)g 1 k=1 converges (unformly n the nterval [0 N ; 1]) to a Lpschtz contnuous multfuncton F 1 (t) wth convex values. Proof: Denote the support functons k (l) = (F k l)andk (l t) = (F k (t) l): Then by (5), (6) and propertes 5 o 6 o of support functons t follows that for each xed drecton l 2 R n k+1 2 (l) = 1 4 k (l)+ 3 4 k (l) k+1 2 (l) =3 4 k (l)+1 4 k (l): 6

7 It means that for each gven l the Chakn subdvson procedure s realzed on the scalar values k (l). Hence by the well known convergence of the Chakn algorthm for scalar subdvson (see e.g. [6]) for each xedl2r n there exsts a lmt functon 1 (l t) = lm k!1 k (l t): In the followng we prove that 1 (l t) s Lpschtz contnuous n l and t. Let us x a pont t 2 [0 N ; 1]: By propertes 1 o 2 o of the support functons k ( t), t follows that the lmt functon 1 ( t) satses 1 o 2 o. Moreover, snce the ntal sets F 0 k =0 :::N are unformly bounded, then by (5),(6)F =0 :::2k (N ; 1) + 1 and F k (t) are unformly bounded by the same constant. Hence by property 3 o the functons k ( t) are Lpschtz contnuous wth an absolute constant, ndependent of t and k, and therefore the lmt functon 1 ( t)slpschtz contnuous wth the same constant and the convergence s unform wth respect to l. Moreover, propertes 1 o 2 o characterzng each support functon hold for the lmt functon 1 ( t) for xed t: Therefore, by 4 o 1 ( t) s the support functon of a convex compact set whch sdenotedby F 1 (t): That F 1 (t) = lm F k (t) can be concluded from k!1 haus(f 1 (t) F k (t)) = max l2s 1 j 1 (l t) ; k (l t)j!0 as k!1: To see that ths convergence s unform for t 2 [0 N ; 1], we observeby (7) and by property (v) of J that haus(f k Fk ) ;1 2;k max haus(f 0 j+1 F0 j ): 0jN;1 Ths means that the Lpschtz constants of the pecewse lnear maps F k (t), and therefore of F 1 (t), ( or of the functons k (l ) 1 (l ) ) do not exceed the constant max haus(f 0 1N F0 ;1): Hence k (l ) andf k () are unformly Lpschtz contnuous wth the same constant nthe nterval [0 N;1] and converge unformly to 1 (l )andf 1 ()) respectvely, onthsnterval. The shape preservng propertes of F 1 (t) and ts smoothness follow from the analyss of more general schemes, done n the next secton. 4 A Class of Shape Preservng Subdvson Schemes Let the ntal sequence F 0 =0 ::: N of convex compact sets n Rn be gven. For convenence we dene F 0 = f0g for 2 Zn f0 1 ::: Ng. Consder a ntely supported subdvson scheme gven by = a [m] ;2j F k j 2 Z k =0 1 2 ::: (8) wth the splne weghts a [m] m+1 = =2 m = 0 1 ::: m+ 1 and a [m] f0 1 ::: m +1g: Note that Chakn algorthm s the specal case m =2. = 0 for 2 Zn 7

8 The scheme (8) when appled to scalar values ff 0 g, s unformly convergent and ts lmt functon f 1 () s of the form (see e.g. [6]) f 1 (t) = 2Z f 0 B m(t ; ) (9) where the functon B m () s a B-splne of degree m, wth nteger knots and support [0 m+1]: In the followng we obtan a set-valued analog of (9). As n the prevous secton, at the k-th teraton (k 1) we construct a pecewse lnear multmap F k : R!C(R n ) satsfyng F k (t k )=F k for t k =2 ;k : Let the generatng functon of the sequence fa [m] g be a [m] (z) = m+1 =0 a [m] z : Then a [m] () hastheform Clearly the coecents a [m] a [m] (z) = (1 + z)m+1 2 m : are nonnegatve numbers and satsfy 2Z a [m] 2 =1 2Z a [m] 2 =1: (10) Denote the derences G k = J(F k ) ; J(F k ;1 ) 1 k 0: Note that Gk 2D(Rn ). It s easy to show as n the case of scalar functons (see e.g.[6]) that Proposton 4.1 The derences G k satsfy G k+1 = j b [m] ;2j Gk j 2 Z k =0 1 ::: (11) where the generatng functon b [m] () of the sequence fb [m] g s b [m] (z) = 2Z b [m] z = a[m] (z) 1+z = (1 + z)m 2 m : (12) It s clear from the last Proposton that the coecents b [m] are nonnegatve and [ m 2 ] =0 [ m b [m] 2 ] 2 = =0 b [m] 2 = 1 2 : (13) Ths means that f the ntal sets ff 0 g N =0 form a monotone ncreasng sequence (.e. the derences G 0 =1 ::: N are n the cone K, then the derences Gk =1 ::: N reman n ths cone at each stage of the subdvson process (11),.e. for each k the sequence ff k g s monotone ncreasng. 8

9 By the same reasonng the second derences H k = G k ; G k ;1 can be obtaned by a subdvson scheme wth a generatng functon c [m] (z) = b[m] (z) 1+z = (1 + z)m;1 2 m : (14) Snce c [m] (z) has nonnegatve coecents, the second derences reman n K at each teraton k, provded they are n K for k = 0. Therefore the scheme (8) s shape preservng. By the above argumentaton t s easy to conclude that the subdvson schemes (8) preserve the sgn of the derences F k of order 1 m.e. ( F k ) belongs to K, provded ( F 0 ) 2K for all. Thus we have proved the followng Proposton 4.2 The subdvson scheme (8) s monotoncty and convexty preservng,.e. F k j F k j+1 for all j =) F k+1 for all (15) F k + F k 2F k j;1 j+1 j for all j =) + ;1 k+1 2F + for all : (16) Theorem 4.1 The set-valued mappngs ff k ()g 1 k=1 multmap F 1 () wth convex mages of the form converge unformly on R to a splne F 1 (t) = 2Z F 0 B m(t ; ) for each t 2 R: (17) Proof: Fx a drecton l 2 R n and denote k (l) = (F k l), k (l t) = (F k (t) l). Then (8) mples k+1 (l) = a [m] ;2j k j (l) hence by (9) for each l 1 (l t) = 2Z 0 (l)b m (t ; ): Snce the last sum s nte for any t, and snce the coecents B m (t ; ) are nonnegatve, t follows that for xed t the scalar functon 1 ( t) s a support functon of the set F 1 (t) = 2Z F 0 B m(t ; ): (18) The fact that F 1 (t) = lm k!1 F k (t) can be proved by the same argumentaton as n the proof of Lemma 3.1 wth (7) replaced by (11) and (13). The next corollary follows from the fact that B m 2 C m;1. Corollary 4.1 F 1 2S and F 1 2 C m;1. For m =1, F 1 s pecewse lnear multfuncton nterpolatng ff 0 g.e. satsfyng F 1 () = F 0 F1 (t) =(t ; )F 0 0 +( +1; t)f for t +1 2 Z. For m = 2 (the Chakn scheme), F 1 2 C 1 and t s pecewse quadratc n the sense that for each xedl the support functon (l ) s a lnear combnaton of quadratc B-splnes. The followng shape preservng propertes of F 1 follow from Proposton 4.2 and the dscusson above t. 9

10 Corollary 4.2 Let ff 0 g N =1 be an ntal sequence ofcompact convex sets. (a) If the ntal sequence ff 0 g N =1 s monotone ncreasng,.e. F 0 F 0 =1 ::: N; 1, then the map F 1 (t) s monotone ncreasng n the sense of Remark 2.3. (b) If the ntal sequence sconvex,.e. F 0 + F 0 F 0 +2 =0 :::N; 2, then the map F 1 (t) has a convex graph. Proof: (a) Wth the notatons of the prevous theorem, t follows from (15) that for every l 2 R n k (l) k (l): Snce a pecewse lnear nterpolatng scalar functon of monotone data s tself monotone, t follows that for each l 2 R n k (l t + h) k (l t) for every t and h>0. The last nequalty mples F k (t + h) F k (t) for all t and h>0. Therefore the set-valued map F k () s monotone ncreasng. The last concluson holds also for the lmt F 1 (t) = lm k!1 F k (t). The proof of (b) s smlar, based on (16) and the fact that a scalar pecewse lnear nterpolant of convex (concave) data s convex (concave). Note that n ths case the setvalued functons are convex, but ther support functons are concave for each xed drecton l (see Remark 2.3). The lmt multmaps generated by the above subdvson processes approxmate smooth multmaps n S wth the same approxmaton order as the lmt functons generated by the correspondng scalar schemes approxmate smooth scalar functons. Proposton 4.3 Let F 2S\C r wth r>0. Then haus(f (t) F (jh)b m ( t h ; j)) chs where s = mnfr 2g and the constant c depends only on F. Proof: Let F (t)= P A f (t), where I s a nte subset of Z, andf 2 C r. Denote 2I ~f (t) = P f (jh)b m ( t ; j). Then by awell-known approxmaton result for scalar functons h (see e.g. [3]), we obtan kf ; f ~ k 1 c h s : P wth the constant c dependng only on f.let~ F(t) = A f ~ (t), then 2I ~F(t) = B m ( t ; j)f (jh): h 10

11 Hence for each t haus(f (t) F(t)) ~ = kj( A f (t)) ; J( A ~ f (t))k ka kkf ; f ~ k 1 ch s (19) 2I 2I 2I where we used property (v) of the embeddng J. For general multmaps (not from S) whch are Hausdor contnuous, we get a weaker approxmaton result. Proposton 4.4 Let F be a general multmap wth values n C(R n ) whch s Hausdor contnuous, then haus(f (t) F (jh)b m ( t ; j)) = o(1): h Proof: Foraxedt there s only a nte number of terms n P F (jh)b m ( t ; j) : Now h snce P P B m ( t ; j) =1 we can wrte F (t) = F (t)b h m ( t ; j). Usng an argument smlar h to (19) and by the fact that F s contnuous, we obtan haus(f (t) F (jh)b m ( t h ; j)) = k j2( h t ;m h t )\Z J(F (t)) ; J(F (jh)) B m ( t ; j)k = o(1): h If we assume enough smoothness on the support functons of the sets F (t) for all t,we get a smlar result to that of Proposton 4.3. Proposton 4.5 Let F be a multmap dened onr wth values n C(R n ) such that the t-dependent support functon (F () l) hasasecond dervatve n t unformly bounded nt and n l 2 S n;1. Then haus(f (t) F (jh)b m ( t h ; j)) = O(h2 ): Then Proof: Let us denote the splne multfuncton S h m F (t) = Denote for any l 2 R n F (jh)b m ( t h ; j) = j2( t h ;m t h )\Z F (jh)b m ( t h ; j): haus(f (t) (S h m F )(t)) = k (F (t) ) ; ((S h m F )(t) )k 1 S n;1 : l (t) =(t l)= (F (t) l). For scalar functons t s known that jf(t) ; (S h f)(t)j 1 m 2 sup j d2 f(t) jh 2 : (20) t dt 2 Set f(t) = l (t). Snce Sm h (F () l)= (Sm h F () l), we get (F (t) l) ; ((S h m F )(t) l)j 1 2 k d2 11 dt 2 l()k 1 h 2 :

12 Now, by our assumptons sup l2s n;1 j d2 dt 2 l ()j L, andwe obtan for every t k (F (t) ) ; ((S h m F )(t) )k 1 S n;1 1 2 Lh2 : Note that f F s dened n the nte nterval [0 N], then the estmate near the boundary of the nterval s O(h). Ths follows from the correspondng result n the scalar case. Before concludng the paper, we state a conjecture that was nspred by the famous theorem of R. Aumann [1] statng that the ntegral of a compact-valued multmap s a convex set. Snce the Remann ntegral s the lmt of Remann sums, whch are n essence averages wth postve weghts, we expect that the repeated applcaton of (8) wth the splne weghts generates n the lmt a convex- valued multfuncton, even when the ntal sets are compact but not convex. Conjecture: Any of the subdvson methods (8) wth m>1 appled to ntal compact sets, ffj 0 g generates as a lmt a convex-valued multfuncton, gven by F 1 (t) = cof 0 j B m(t ; j) where coa denotes the convex hull of A. References [1] R. J. Aumann,Integrals of set-valued functons, J. of Math. Analyss and Applcatons, 12 (1965) 1{12. [2] R. Baer and E. Farkh,Drected Sets and Derences of Convex Compact Sets, n: M.P. Pols, A.L. Dontchev, P. Kall, I. Lasecka and A.W. Olbrot, eds., Systems Modellng and Optmzaton, Proc. of the 18th IFIP TC7 Conference, CRC Research Notes n Mathematcs, (Chapman and Hall,1999) 135{143. [3] C. deboor,a Practcal Gude to Splnes (Sprnger Verlag, New York,1978). [4] A. S. Cavaretta, W. Dahmen and C.A. Mcchell,Statonary Subdvson, Memors of AMS, No. 453 (1991). [5] G. M. Chakn,An algorthm for hgh speed curve generaton, Computer Graphcs and Image Processng, 3 (1974) 346{349. [6] N. Dyn,Subdvson schemes n Computer-Aded Geometrc Desgn, n: W. Lght, ed., Advances n Numercal Analyss, Vol. II, Wavelets, Subdvson Algorthms and Radal Bass Functons (Clarendon Press, Oxford, 1992) 36{104. [7] L. Hormander,Sur la foncton d'appu des ensembles convexes dans un espace localement convexe, Arkv for Matematk, 3 (1954)181{

13 [8] B. Margols,Compact, Convex Sets n R n and a new Banach Lattce, I.- Theory, Numer. Funct. Anal. Optmz. 11 (1990) 555{576. [9] H. Radstrom,An embeddng theorem for spaces of convex sets, Proceedngs of the Amercan Mathematcal Socety 3 (1952) 165{169. [10] R. T. Rockafellar,Convex Analyss (Prnceton Unversty Press, Prnceton, 1970). [11] R. A. Vtale,Approxmatons of convex set-valued functons, J. Approx. Theory 26 (1979) 301{316. Nra Dyn nradyn@math.tau.ac.l Elza Farkh elza@math.tau.ac.l School of Mathematcal Scences Sackler Faculty of Exact Scences Tel Avv Unversty, Tel Avv, Israel 13

Convexity preserving interpolation by splines of arbitrary degree

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

More metrics on cartesian products

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

Approximations of Set-Valued Functions Based on the Metric Average

Approximations of Set-Valued Functions Based on the Metric Average Approxmatons of Set-Valued Functons Based on the Metrc Average Nra Dyn, Alona Mokhov School of Mathematcal Scences Tel-Avv Unversty, Israel Abstract. Ths paper nvestgates the approxmaton of set-valued

More information

Appendix B. Criterion of Riemann-Stieltjes Integrability

Appendix B. Criterion of Riemann-Stieltjes Integrability Appendx B. Crteron of Remann-Steltes Integrablty Ths note s complementary to [R, Ch. 6] and [T, Sec. 3.5]. The man result of ths note s Theorem B.3, whch provdes the necessary and suffcent condtons for

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Approximation of Univariate Set-Valued Functionsan Overview

Approximation of Univariate Set-Valued Functionsan Overview Approxmaton of Unvarate Set-Valued Functonsan Overvew ra Dyn, Elza Farkh, Alona Mokhov School of Mathematcal Scences Tel-Avv Unversty, Israel Abstract The paper s an updated survey of our work on the approxmaton

More information

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS BOUNDEDNESS OF THE IESZ TANSFOM WITH MATIX A WEIGHTS Introducton Let L = L ( n, be the functon space wth norm (ˆ f L = f(x C dx d < For a d d matrx valued functon W : wth W (x postve sem-defnte for all

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

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

APPENDIX A Some Linear Algebra

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

PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 125, Number 7, July 1997, Pages 2119{2125 S (97) THE STRONG OPEN SET CONDITION

PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 125, Number 7, July 1997, Pages 2119{2125 S (97) THE STRONG OPEN SET CONDITION PROCDINGS OF TH AMRICAN MATHMATICAL SOCITY Volume 125, Number 7, July 1997, Pages 2119{2125 S 0002-9939(97)03816-1 TH STRONG OPN ST CONDITION IN TH RANDOM CAS NORBRT PATZSCHK (Communcated by Palle. T.

More information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

DIEGO AVERNA. A point x 2 X is said to be weakly Pareto-optimal for the function f provided

DIEGO AVERNA. A point x 2 X is said to be weakly Pareto-optimal for the function f provided WEAKLY PARETO-OPTIMAL ALTERNATIVES FOR A VECTOR MAXIMIZATION PROBLEM: EXISTENCE AND CONNECTEDNESS DIEGO AVERNA Let X be a non-empty set and f =f 1 ::: f k ):X! R k a functon. A pont x 2 X s sad to be weakly

More information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter

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

DEGREE REDUCTION OF BÉZIER CURVES USING CONSTRAINED CHEBYSHEV POLYNOMIALS OF THE SECOND KIND

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 information

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

General viscosity iterative method for a sequence of quasi-nonexpansive mappings Avalable onlne at www.tjnsa.com J. Nonlnear Sc. Appl. 9 (2016), 5672 5682 Research Artcle General vscosty teratve method for a sequence of quas-nonexpansve mappngs Cuje Zhang, Ynan Wang College of Scence,

More information

Exercise Solutions to Real Analysis

Exercise Solutions to Real Analysis xercse Solutons to Real Analyss Note: References refer to H. L. Royden, Real Analyss xersze 1. Gven any set A any ɛ > 0, there s an open set O such that A O m O m A + ɛ. Soluton 1. If m A =, then there

More information

arxiv: v1 [math.na] 8 Aug 2018

arxiv: v1 [math.na] 8 Aug 2018 Convergence analyss of corner cuttng algorthms refnng ponts and refnng nets of functons Costanza Cont a,, Nra Dyn b, Luca Roman c a Dpartmento d Ingegnera Industrale, Unverstà d Frenze, Italy b School

More information

Some basic inequalities. Definition. Let V be a vector space over the complex numbers. An inner product is given by a function, V V C

Some basic inequalities. Definition. Let V be a vector space over the complex numbers. An inner product is given by a function, V V C Some basc nequaltes Defnton. Let V be a vector space over the complex numbers. An nner product s gven by a functon, V V C (x, y) x, y satsfyng the followng propertes (for all x V, y V and c C) (1) x +

More information

Maximizing the number of nonnegative subsets

Maximizing the number of nonnegative subsets Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum

More information

A new Approach for Solving Linear Ordinary Differential Equations

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

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

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

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

Lecture 12: Discrete Laplacian

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

MINIMAL MULTI-CONVEX PROJECTIONS ONTO SUBSPACES OF POLYNOMIALS JOANNA MEISSNER 1

MINIMAL MULTI-CONVEX PROJECTIONS ONTO SUBSPACES OF POLYNOMIALS JOANNA MEISSNER 1 MINIMAL MULTI-CONVEX PROJECTIONS ONTO SUBSPACES OF POLYNOMIALS JOANNA MEISSNER IMUJ Preprnt 9/9 Abstract. Let X C N [, ], ), where N 3 and let V be a lnear subspace of Π N, where Π N denotes the space

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

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

Bezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0

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

Report on Image warping

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

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

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

1 Convex Optimization

1 Convex Optimization Convex Optmzaton We wll consder convex optmzaton problems. Namely, mnmzaton problems where the objectve s convex (we assume no constrants for now). Such problems often arse n machne learnng. For example,

More information

SUCCESSIVE MINIMA AND LATTICE POINTS (AFTER HENK, GILLET AND SOULÉ) M(B) := # ( B Z N)

SUCCESSIVE MINIMA AND LATTICE POINTS (AFTER HENK, GILLET AND SOULÉ) M(B) := # ( B Z N) SUCCESSIVE MINIMA AND LATTICE POINTS (AFTER HENK, GILLET AND SOULÉ) S.BOUCKSOM Abstract. The goal of ths note s to present a remarably smple proof, due to Hen, of a result prevously obtaned by Gllet-Soulé,

More information

2 S. S. DRAGOMIR, N. S. BARNETT, AND I. S. GOMM Theorem. Let V :(d d)! R be a twce derentable varogram havng the second dervatve V :(d d)! R whch s bo

2 S. S. DRAGOMIR, N. S. BARNETT, AND I. S. GOMM Theorem. Let V :(d d)! R be a twce derentable varogram havng the second dervatve V :(d d)! R whch s bo J. KSIAM Vol.4, No., -7, 2 FURTHER BOUNDS FOR THE ESTIMATION ERROR VARIANCE OF A CONTINUOUS STREAM WITH STATIONARY VARIOGRAM S. S. DRAGOMIR, N. S. BARNETT, AND I. S. GOMM Abstract. In ths paper we establsh

More information

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

Bé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 information

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

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

20. Mon, Oct. 13 What we have done so far corresponds roughly to Chapters 2 & 3 of Lee. Now we turn to Chapter 4. The first idea is connectedness.

20. Mon, Oct. 13 What we have done so far corresponds roughly to Chapters 2 & 3 of Lee. Now we turn to Chapter 4. The first idea is connectedness. 20. Mon, Oct. 13 What we have done so far corresponds roughly to Chapters 2 & 3 of Lee. Now we turn to Chapter 4. The frst dea s connectedness. Essentally, we want to say that a space cannot be decomposed

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

On a nonlinear compactness lemma in L p (0, T ; B).

On a nonlinear compactness lemma in L p (0, T ; B). On a nonlnear compactness lemma n L p (, T ; B). Emmanuel Matre Laboratore de Matématques et Applcatons Unversté de Haute-Alsace 4, rue des Frères Lumère 6893 Mulouse E.Matre@ua.fr 3t February 22 Abstract

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES

TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES SVANTE JANSON Abstract. We gve explct bounds for the tal probabltes for sums of ndependent geometrc or exponental varables, possbly wth dfferent

More information

Perron Vectors of an Irreducible Nonnegative Interval Matrix

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

2 MADALINA ROXANA BUNECI subset G (2) G G (called the set of composable pars), and two maps: h (x y)! xy : G (2)! G (product map) x! x ;1 [: G! G] (nv

2 MADALINA ROXANA BUNECI subset G (2) G G (called the set of composable pars), and two maps: h (x y)! xy : G (2)! G (product map) x! x ;1 [: G! G] (nv An applcaton of Mackey's selecton lemma Madalna Roxana Bunec Abstract. Let G be a locally compact second countable groupod. Let F be a subset of G (0) meetng each orbt exactly once. Let us denote by df

More information

Journal of Universal Computer Science, vol. 1, no. 7 (1995), submitted: 15/12/94, accepted: 26/6/95, appeared: 28/7/95 Springer Pub. Co.

Journal of Universal Computer Science, vol. 1, no. 7 (1995), submitted: 15/12/94, accepted: 26/6/95, appeared: 28/7/95 Springer Pub. Co. Journal of Unversal Computer Scence, vol. 1, no. 7 (1995), 469-483 submtted: 15/12/94, accepted: 26/6/95, appeared: 28/7/95 Sprnger Pub. Co. Round-o error propagaton n the soluton of the heat equaton by

More information

Geometry of Müntz Spaces

Geometry of Müntz Spaces WDS'12 Proceedngs of Contrbuted Papers, Part I, 31 35, 212. ISBN 978-8-7378-224-5 MATFYZPRESS Geometry of Müntz Spaces P. Petráček Charles Unversty, Faculty of Mathematcs and Physcs, Prague, Czech Republc.

More information

Approximate D-optimal designs of experiments on the convex hull of a finite set of information matrices

Approximate D-optimal designs of experiments on the convex hull of a finite set of information matrices Approxmate D-optmal desgns of experments on the convex hull of a fnte set of nformaton matrces Radoslav Harman, Mára Trnovská Department of Appled Mathematcs and Statstcs Faculty of Mathematcs, Physcs

More information

On the Multicriteria Integer Network Flow Problem

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

Self-complementing permutations of k-uniform hypergraphs

Self-complementing permutations of k-uniform hypergraphs Dscrete Mathematcs Theoretcal Computer Scence DMTCS vol. 11:1, 2009, 117 124 Self-complementng permutatons of k-unform hypergraphs Artur Szymańsk A. Paweł Wojda Faculty of Appled Mathematcs, AGH Unversty

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Anti-van der Waerden numbers of 3-term arithmetic progressions.

Anti-van der Waerden numbers of 3-term arithmetic progressions. Ant-van der Waerden numbers of 3-term arthmetc progressons. Zhanar Berkkyzy, Alex Schulte, and Mchael Young Aprl 24, 2016 Abstract The ant-van der Waerden number, denoted by aw([n], k), s the smallest

More information

MATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2)

MATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2) 1/16 MATH 829: Introducton to Data Mnng and Analyss The EM algorthm (part 2) Domnque Gullot Departments of Mathematcal Scences Unversty of Delaware Aprl 20, 2016 Recall 2/16 We are gven ndependent observatons

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

n α 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

n α 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 information

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS Journal of Mathematcal Scences: Advances and Applcatons Volume 25, 2014, Pages 1-12 A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS JIA JI, WEN ZHANG and XIAOFEI QI Department of Mathematcs

More information

Modeling curves. Graphs: y = ax+b, y = sin(x) Implicit ax + by + c = 0, x 2 +y 2 =r 2 Parametric:

Modeling 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

On some variants of Jensen s inequality

On some variants of Jensen s inequality On some varants of Jensen s nequalty S S DRAGOMIR School of Communcatons & Informatcs, Vctora Unversty, Vc 800, Australa EMMA HUNT Department of Mathematcs, Unversty of Adelade, SA 5005, Adelade, Australa

More information

Research Article Relative Smooth Topological Spaces

Research Article Relative Smooth Topological Spaces Advances n Fuzzy Systems Volume 2009, Artcle ID 172917, 5 pages do:10.1155/2009/172917 Research Artcle Relatve Smooth Topologcal Spaces B. Ghazanfar Department of Mathematcs, Faculty of Scence, Lorestan

More information

Counterexamples to the Connectivity Conjecture of the Mixed Cells

Counterexamples to the Connectivity Conjecture of the Mixed Cells Dscrete Comput Geom 2:55 52 998 Dscrete & Computatonal Geometry 998 Sprnger-Verlag New York Inc. Counterexamples to the Connectvty Conjecture of the Mxed Cells T. Y. L and X. Wang 2 Department of Mathematcs

More information

Lecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem.

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

ON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES. 1. Introduction

ON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES. 1. Introduction ON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES TAKASHI ITOH AND MASARU NAGISA Abstract We descrbe the Haagerup tensor product l h l and the extended Haagerup tensor product l eh l n terms of

More information

Polynomial Reproduction by Symmetric Subdivision Schemes

Polynomial Reproduction by Symmetric Subdivision Schemes Polynomal Reproducton by Symmetrc Subdvson Schemes Nra Dyn School of Mathematcal Scences Tel Avv Unversty Malcolm A. Sabn Computer Laboratory Unversty of Cambrdge Ka Hormann Department of Informatcs Clausthal

More information

REGULAR POSITIVE TERNARY QUADRATIC FORMS. 1. Introduction

REGULAR POSITIVE TERNARY QUADRATIC FORMS. 1. Introduction REGULAR POSITIVE TERNARY QUADRATIC FORMS BYEONG-KWEON OH Abstract. A postve defnte quadratc form f s sad to be regular f t globally represents all ntegers that are represented by the genus of f. In 997

More information

arxiv: v1 [math.co] 12 Sep 2014

arxiv: v1 [math.co] 12 Sep 2014 arxv:1409.3707v1 [math.co] 12 Sep 2014 On the bnomal sums of Horadam sequence Nazmye Ylmaz and Necat Taskara Department of Mathematcs, Scence Faculty, Selcuk Unversty, 42075, Campus, Konya, Turkey March

More information

Ballot Paths Avoiding Depth Zero Patterns

Ballot Paths Avoiding Depth Zero Patterns Ballot Paths Avodng Depth Zero Patterns Henrch Nederhausen and Shaun Sullvan Florda Atlantc Unversty, Boca Raton, Florda nederha@fauedu, ssull21@fauedu 1 Introducton In a paper by Sapounaks, Tasoulas,

More information

arxiv: v1 [math.co] 1 Mar 2014

arxiv: v1 [math.co] 1 Mar 2014 Unon-ntersectng set systems Gyula O.H. Katona and Dánel T. Nagy March 4, 014 arxv:1403.0088v1 [math.co] 1 Mar 014 Abstract Three ntersecton theorems are proved. Frst, we determne the sze of the largest

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

CHAPTER III Neural Networks as Associative Memory

CHAPTER III Neural Networks as Associative Memory CHAPTER III Neural Networs as Assocatve Memory Introducton One of the prmary functons of the bran s assocatve memory. We assocate the faces wth names, letters wth sounds, or we can recognze the people

More information

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton

More information

Memoirs on Dierential Equations and Mathematical Physics Volume 11, 1997, 67{88 Guram Kharatishvili and Tamaz Tadumadze THE PROBLEM OF OPTIMAL CONTROL

Memoirs on Dierential Equations and Mathematical Physics Volume 11, 1997, 67{88 Guram Kharatishvili and Tamaz Tadumadze THE PROBLEM OF OPTIMAL CONTROL Memors on Derental Equatons and Mathematcal Physcs Volume 11, 1997, 67{88 Guram Kharatshvl and Tamaz Tadumadze THE PROBLEM OF OPTIMAL CONTROL FOR NONLINEAR SYSTEMS WITH VARIABLE STRUCTURE, DELAYS AND PIECEWISE

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

More information

Remarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence

Remarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence Remarks on the Propertes of a Quas-Fbonacc-lke Polynomal Sequence Brce Merwne LIU Brooklyn Ilan Wenschelbaum Wesleyan Unversty Abstract Consder the Quas-Fbonacc-lke Polynomal Sequence gven by F 0 = 1,

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.070J Fall 013 Lecture 1 10/1/013 Martngale Concentraton Inequaltes and Applcatons Content. 1. Exponental concentraton for martngales wth bounded ncrements.

More information

A Hybrid Variational Iteration Method for Blasius Equation

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

The lower and upper bounds on Perron root of nonnegative irreducible matrices

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

REAL ANALYSIS I HOMEWORK 1

REAL ANALYSIS I HOMEWORK 1 REAL ANALYSIS I HOMEWORK CİHAN BAHRAN The questons are from Tao s text. Exercse 0.0.. If (x α ) α A s a collecton of numbers x α [0, + ] such that x α

More information

Lecture 4: Universal Hash Functions/Streaming Cont d

Lecture 4: Universal Hash Functions/Streaming Cont d CSE 5: Desgn and Analyss of Algorthms I Sprng 06 Lecture 4: Unversal Hash Functons/Streamng Cont d Lecturer: Shayan Oves Gharan Aprl 6th Scrbe: Jacob Schreber Dsclamer: These notes have not been subjected

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space.

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space. Lnear, affne, and convex sets and hulls In the sequel, unless otherwse specfed, X wll denote a real vector space. Lnes and segments. Gven two ponts x, y X, we defne xy = {x + t(y x) : t R} = {(1 t)x +

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

NUMERICAL DIFFERENTIATION

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

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP C O L L O Q U I U M M A T H E M A T I C U M VOL. 80 1999 NO. 1 FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP BY FLORIAN K A I N R A T H (GRAZ) Abstract. Let H be a Krull monod wth nfnte class

More information

First day August 1, Problems and Solutions

First day August 1, Problems and Solutions FOURTH INTERNATIONAL COMPETITION FOR UNIVERSITY STUDENTS IN MATHEMATICS July 30 August 4, 997, Plovdv, BULGARIA Frst day August, 997 Problems and Solutons Problem. Let {ε n } n= be a sequence of postve

More information

Another converse of Jensen s inequality

Another converse of Jensen s inequality Another converse of Jensen s nequalty Slavko Smc Abstract. We gve the best possble global bounds for a form of dscrete Jensen s nequalty. By some examples ts frutfulness s shown. 1. Introducton Throughout

More information

Monotonic Interpolating Curves by Using Rational. Cubic Ball Interpolation

Monotonic Interpolating Curves by Using Rational. Cubic Ball Interpolation Appled Mathematcal Scences, vol. 8, 204, no. 46, 7259 7276 HIKARI Ltd, www.m-hkar.com http://dx.do.org/0.2988/ams.204.47554 Monotonc Interpolatng Curves by Usng Ratonal Cubc Ball Interpolaton Samsul Arffn

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

A note on almost sure behavior of randomly weighted sums of φ-mixing random variables with φ-mixing weights

A note on almost sure behavior of randomly weighted sums of φ-mixing random variables with φ-mixing weights ACTA ET COMMENTATIONES UNIVERSITATIS TARTUENSIS DE MATHEMATICA Volume 7, Number 2, December 203 Avalable onlne at http://acutm.math.ut.ee A note on almost sure behavor of randomly weghted sums of φ-mxng

More information

Online Classification: Perceptron and Winnow

Online Classification: Perceptron and Winnow E0 370 Statstcal Learnng Theory Lecture 18 Nov 8, 011 Onlne Classfcaton: Perceptron and Wnnow Lecturer: Shvan Agarwal Scrbe: Shvan Agarwal 1 Introducton In ths lecture we wll start to study the onlne learnng

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

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

Restricted divisor sums

Restricted divisor sums ACTA ARITHMETICA 02 2002) Restrcted dvsor sums by Kevn A Broughan Hamlton) Introducton There s a body of work n the lterature on varous restrcted sums of the number of dvsors of an nteger functon ncludng

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

NP-Completeness : Proofs

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

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v10157-01-00-y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION

More information

SL n (F ) Equals its Own Derived Group

SL n (F ) Equals its Own Derived Group Internatonal Journal of Algebra, Vol. 2, 2008, no. 12, 585-594 SL n (F ) Equals ts Own Derved Group Jorge Macel BMCC-The Cty Unversty of New York, CUNY 199 Chambers street, New York, NY 10007, USA macel@cms.nyu.edu

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

On Finite Rank Perturbation of Diagonalizable Operators

On Finite Rank Perturbation of Diagonalizable Operators Functonal Analyss, Approxmaton and Computaton 6 (1) (2014), 49 53 Publshed by Faculty of Scences and Mathematcs, Unversty of Nš, Serba Avalable at: http://wwwpmfnacrs/faac On Fnte Rank Perturbaton of Dagonalzable

More information

Homogenization of reaction-diffusion processes in a two-component porous medium with a non-linear flux-condition on the interface

Homogenization of reaction-diffusion processes in a two-component porous medium with a non-linear flux-condition on the interface Homogenzaton of reacton-dffuson processes n a two-component porous medum wth a non-lnear flux-condton on the nterface Internatonal Conference on Numercal and Mathematcal Modelng of Flow and Transport n

More information