Moving frame and integrable system of the discrete centroaffine curves in R 3

Size: px
Start display at page:

Download "Moving frame and integrable system of the discrete centroaffine curves in R 3"

Transcription

1 Movig frame ad itegrable system of the discrete cetroaffie curves i R 3 Yu Yag, Yahua Yu Departmet of Mathematics, Northeaster Uiversity, Sheyag 0004, P R Chia arxiv: v2 [mathdg] 27 Nov 206 Abstract Ay two equivalet discrete curves must have the same ivariats at the correspodig poits uder a affie trasformatio I this paper, we costruct the movig frame ad ivariats for the discrete cetroaffie curves, which could be used to discrimiate the same discrete curves from differet graphics, ad estimate whether a polygo flow is stable or periodically stable I fact, usig the similar method as the Freet-Serret frame, a discrete curve ca be uiquely idetified by its cetroaffie curvatures ad torsios I 878, Darboux studied the problem of midpoit iteratio of polygos[2] Berlekamp et al studied this problem i detail[2] Now, through the cetroaffie curvatures ad torsios, the iteratio process ca be clearly quatified Exactly, we describe the whole iteratio process by usig cetroaffie curvatures ad torsios, ad its periodicity could be directly exhibited As a applicatio, we would obtai some stable discrete space curves with chageless curvatures ad torsios after multistep iteratio For the petagram map of a polygo, the affiely regular polygos are stable Furthermore, we fid the covex hexagos with parallel ad equi-legth opposite sides are periodically stable, ad some covex parallel ad equi-legth opposite sides octagos are also periodically stable The proofs of these results are obtaied usig the structure equatios of the discrete cetroaffie curves ad the itegrable coditios of its flows MSC 200: 52C07, 53A5 Key Words: Discrete differetial geometry, affie trasformatio, discrete curvature flow Itroductio Discrete differetial geometry has attracted much attetio recetly, maily due to the growth of computer graphics Oe of the mai issues i discrete differetial geometry is to defie suitable discrete aalogous of the cocepts of smooth differetial geometry[4, 6] More recetly, the expasio of computer graphics ad applicatios i mathematical physics have give a great impulse to the issue of givig discrete equivalets of affie differetial geometric objects[3, 9, 0] I [5] a cosistet defiitio of discrete affie spheres is proposed, both for defiite ad idefiite metrics ad i [2] a similar costructio is doe i the cotext of improper affie spheres This work was supported by NSFC(Nos ad 37080) Correspodig author addresses: yagyu@maileueduc (Y Yag), yuyahua@maileueduc(y Yu)

2 Group based movig frames have a wide rage of applicatios, from the classical equivalece problems i differetial geometry to more moder applicatios such as computer visio[20, 25] The first results for the computatio of discrete ivariats usig group based movig frames were give by Olver[22] who calls them joit ivariats; moder applicatios to date iclude computer visio[23] ad umerical schemes for systems with a Lie symmetry[8] Movig frames for discrete applicatios as formulated by Olver do give geeratig sets of discrete ivariats, ad the recursio formulas for differetial ivariats are so successful for the applicatio of movig frames to calculus-based applicatios Recet developmet of a theory of discrete equivariat movig frames has bee applied to itegrable differetial-differece systems[20, 25] Followig the ideas of Klei, preseted i his famous lecture at Erlage, several geometers i the early 20th cetury proposed the study of curves ad surfaces with respect to differet trasformatio groups I geometry, a affie trasformatio, affie map or a affiity is a fuctio betwee affie spaces which preserves poits, straight lies ad plaes Also, sets of parallel lies remai parallel after a affie trasformatio A affie trasformatio does ot ecessarily preserve agles betwee lies or distaces betwee poits, though it does preserve ratios of distaces betwee poits lyig o a straight lie Examples of affie trasformatios iclude traslatio, scalig, homothety, similarity trasformatio, reflectio, rotatio, shear mappig, ad compositios of them i ay combiatio ad sequece A cetroaffie trasformatio is othig but a geeral liear trasformatio R x Ax R, where A GL(,R) I 907 Tzitzéica foud that for a surface i Euclidea 3-space the property that the ratio of the Gauss curvature to the fourth power of the distace of the taget plae from the origi is costat is ivariat uder a cetroaffie trasformatio The surfaces with this property tur out to be what are ow called Tzitzéica surfaces, or proper affie spheres with ceter at the origi I cetroaffie differetial geometry, the theory of hypersurfaces has a log history The otio of cetroaffie miimal hypersurfaces was itroduced by Wag [3] as extremals for the area itegral of the cetroaffie metric See also [33, 34] for the classificatio results about cetroaffie traslatiosurfaces adcetroaffie ruled surfaces ir 3 Several authors studied the curves uder cetroaffie trasformatio group usig some methods([3, 5], etc), ad Liu defie cetroaffie ivariat arc legth ad cetroaffie curvature fuctios of a curve i affie -space directly by the parameter trasformatios ad the cetroaffie trasformatios[8] Usig the equivariat method of movig frames, Olver costructed the explicit formulas for the geeratig differetial ivariats ad ivariat differetial operators for curves i 2-dimesioal ad 3-dimesioal cetro-equi-affie ad cetroaffie geometry ad surfaces i 3-dimesioal cetro-equi-affie geometry[24] The study of discrete itegrable systems is rather ew It bega with discretizig cotiuous itegrable systems i 970s The best kow discretizatio of the Korteweg-de Vries equatio (KdV) is the Toda lattice[30] Aother famous itegrable discretizatio of the KdV equatio is the Volterra lattice[6, 9] I this paper, the defiitios ad costructios of discrete itegrable systems are ature ad useful It arises as aalogues of curvature flows for polygo evolutios I mathematics, curvature refers to ay of a umber of loosely related cocepts i differet areas of geometry Ituitively, curvature is the amout by which a geometric object deviates from beig flat, but this is defied i differet ways depedig o the cotext The arragemet of the paper is as follows: I Sect 2 we recall the basic theory ad otios for cetroaffie differetial geometry, the basic otatios for discrete curves ad cetroaffie curves There are some results with cetroaffie curvatures ad torsios for cetroaffie plaar curves ad space curves i Sect 3 ad Sect 4 I Sect 5, we exted the curve shorteig flow to the discrete cetroaffie curve I Sects 6 ad 7, we study the trasversal flow ad taget flow for a discrete cetroaffie curve respectively Some iterest examples for polygo iteratio 2

3 are show Fially, Sect 8 describes some applicatios of cetroaffie curvatures ad torsios We coclude with idicatios of future work 2 Affie mappigs ad trasformatio groups, basic otatios If X ad Y are affie spaces, the every affie trasformatio f : X Y is of the form x M x+ b, where M is a liear trasformatio o X ad b is a vector i Y Ulike a purely liear trasformatio, a affie map eed ot preserve the zero poit i a liear space Thus, every liear trasformatio is affie, but ot every affie trasformatio is liear For may purposes a affie space ca be thought of as Euclidea space, though the cocept of affie space is far more geeral (ie, all Euclidea spaces are affie, but there are affie spaces that are o-euclidea) I affie coordiates, which iclude Cartesia coordiates i Euclidea spaces, each output coordiate of a affie map is a liear fuctio (i the sese of calculus) of all iput coordiates Aother way to deal with affie trasformatios systematically is to select a poit as the origi; the, ay affie trasformatio is equivalet to a liear trasformatio (of positio vectors) followed by a traslatio It is well kow that the set of all automorphisms of a vector space V of dimesio m forms a group We use the followig stadard otatios for this group ad its subgroups([7]): GL(m,R) := {L : V V L isomorphism}; SL(m,R) := {L GL(m,R) detl = } Correspodigly, for a affie space A,dimA = m, we have the followig affie trasformatio groups A(m) := {α : A A L α regular} is the regular affie group S(m) := {α A detα = } is the uimodular(equiaffie) group Z p (m) := {α A α(p) = p} is the cetroaffie group with ceter p A τ(m) := {α : A A there exists b(α) V,st pα(p) = b(α), p A} is the group of trasformatios o A Let G be oe of the groups above ad S,S 2 A subsets The S ad S 2 are called equivalet modulo G if there exists a α G such that S 2 = αs The stadard properties of affie mappigs come from the properties of the associated liear mappig Recall i particular: (i) Parallelism is ivariat uder affie mappigs (ii) The partitio ratio of 3 poits is ivariat uder affie mappigs (iii) The ratio of the volumes of two parallelepipeds is affiely ivariat Moreover, covexity is a affie property 3

4 Theorem 2 α : A A is a regular affie trasformatio if ad oly if α is bijective, cotiuous ad preserves covexity (For a more geeral result see [32] ) I cetroaffie geometry we fix a poit i A(the origi O A without loss of geerality) ad cosider the geometric properties i variat uder the cetroaffie group Z p Thus the mappig π 0 : A V idetifies A with the vector space V ad Z O with GL(m,R) We will start very simple, by discretizig the otio of a smooth curve That is, we wat to defie a discrete aalog to a smooth map from a iterval I R to R By discrete we mea here that the map should ot be defied o a iterval i R but o a discrete (ordered) set of poits therei It turs out that this is basically all we eed to demad i this case: Defiitio 22 Let I Z be a iterval (the itersectio of a iterval i R with Z, possibly ifiite) A map r : I R is called a discrete curve, whe we put the startig poit of the vector r to the origi O R Obviously, a discrete curve is a polygo A discrete curve r is said to be periodic (or closed) if I = Z ad if there is a p Z such that r(k) = r(k +p) for all k I The smallest possible value of p is called the period I fact, we ca defie (2) r(t) = (t k) r(k)+(k + t) r(k +), t (k,k +),k Z The simplicity of a smooth curve ca be geeralized to the discrete case Defiitio 23 A closed curve r is simple if it has o further self-itersectios i oe period; that is, if t,t 2 [k,k + p),t t 2,k Z, the r(t ) r(t 2 ) We call a o closed curve is simple if t,t 2 [a,b),t t 2, the r(t ) r(t 2 ) For coveiece, sometimes we will write r k = r(k) ad eve r = r k, r = r k+, ad r = r k Defiitio 24 The edge taget vector of a discrete curve r : I R is defied as the forward differece t k := r k+ r k We could have writte t := r r as well The lies passig through the termial poits of r k ad r k+, k I Z are called taget lies of the discrete curve r With these preparatios, it is ature to give the defiitio of covexity for a discrete plaar curve as follows Defiitio 25 A covex discrete curve is the curve i the plae which lies completely o oe side of each ad every oe of its taget lies I classical differetial geometry of curve, the followig results are well-kow, ad we will have some similar results for a discrete curve with cetroaffie curvatures ad torsios i the ext two sectios Theorem 26 A closed regular plaar simple curve C is covex if ad oly if its curvature is either always o-egative or always o-positive, ie, if ad oly if the turig agle (the agle of the taget to the curve) is a weakly mootoe fuctio of the parametrizatio of the curve(for details see [4]) 4

5 I classical differetial geometry, the Freet-Serret formulas describe the kiematic properties of a particle movig alog a cotiuous, differetiable curve i three-dimesioal Euclidea space R 3, or the geometric properties of the curve itself irrespective of ay motio More specifically, the formulas describe the derivatives of the so-called taget, ormal, ad biormal uit vectors i terms of each other The taget, ormal, ad biormal uit vectors, ofte called T, N, ad B, or collectively the Freet-Serret frame, together form a orthoormal basis spaig R 3 ad are defied as follows: T is the uit vector taget to the curve, poitig i the directio of motio N is the ormal uit vector, the derivative of T with respect to the arclegth parameter of the curve, divided by its legth B is the biormal uit vector, the cross product of T ad N The Freet-Serret formulas are: ( dt ds, dn ds, db ds ) = (T,N,B) 0 κ 0 κ 0 τ 0 τ 0 where d is the derivative with respect to arc legth, κ is the curvature, ad τ is the torsio of ds the curve The two scalars κ ad τ effectively defie the curvature ad torsio of a space curve The associated collectio, T, N, B, κ ad τ, is called the Freet-Serret apparatus Ituitively, curvature measures the failure of a curve to be a straight lie, while torsio measures the failure of a curve to be plaar At the same time we have Theorem 27 (Fudametal theorem of space curves) Let r (s) ad r 2 (s) be two vectorvalued fuctios that represet the space curves C ad C 2 respectively, ad suppose that these curves have the same o-vaishig curvature κ(s) ad the same torsio τ(s) The C ad C 2 are cogruet such that each ca be rigidly shifted/rotated so that every poit o C coicides with every poit o C 2 (for details see [4]) Before startig ext sectio, we shall eed a defiitio of the discrete cetroaffie curves i R 2 ad R 3 Defiitio 28 A discrete plaar curve r : I R 2 is called a cetroaffie plaar curve if the edge taget vector t k is ot parallel to positio vectors r(k) ad r(k +), ad A discrete curve r : I R 3 is called a cetroaffie curve if the edge taget vectors t k, t k ad the positio vector r(k) are ot coplaar 3 Discrete plaar curves uder the affie trasformatio I this sectio we wat to cosider two ivariats ad their geometrical properties uder the affie trasformatio, although we call them the first ad secod cetroaffie curvatures just for uity, that because i the ext sectio we will use them together with cetroaffie torsios for a space discrete cetroaffie curve, oly uder the cetroaffie trasformatio Now let vectorvalued fuctio r : I Z R 2 represet a discrete curve C Below we give the defiitio of the cetroaffie curvatures by usig the otatios of the previous sectio Defiitio 3 For the discrete plaar curve C, if [ t k, t k ] = 0, we call its first cetroaffie curvature κ k = 0 at poit r(k), which implies the curve is a straight lie locally to r(k), where [ ] deotes the stadard determiat i R 2 If [ t k, t k ] 0, the first ad secod cetroaffie curvatures at the poit r(k) are defied by (3) κ k = [ t k, t k+ ] [ t k, t k ], κ k = [ t k, t k+ ] [ t k, t k ] 5,

6 Uder a affie trasformatio r = A r + b, we have t k = A t k So [ t k, t k+ ] [ t k, t k ] = (deta)[ t k, t k+] (deta)[ t k, t k ] = [ t k, t k+ ] [ t k, t k ], [ t k, t k+ ] [ t k, t k ] = (deta)[ t k, t k+ ] (deta)[ t k, t k ] = [ t k, t k+ ], [ t k, t k ] which implies κ k ad κ k are affie ivariats, of course, ad also cetroaffie ivariats Figure : The sig of the first curvature Now let us explai the geometrical meaig of the cetroaffie curvatures I Figure, P k deotes the ed poit of vector r(k) Therefore, the determiat [ t k, t k ] i R 2 exactly represets twice the oriet area of the triagle Pk P k P k+, ad [ t k, t k+ ] represets twice the oriet area of the triagle Pk P k+ P k+2 As a result, if the poits P k+2 ad P k lie o the same side of the straight lie P k P k+, the first cetroaffie curvature κ k takes positive value, ad if they lie o differet sides of the straight lie P k P k+, the first cetroaffie curvature κ k takes egative value I details, we ca see the left oe i Figure 2 Sice the triagles Pk P k P k+ ad Pk P k+ P k+2 have a edge P k P k+ i commo, accordig to the relatio of height ad area of a triagle, if poit P k+2 lies o differet lies, the cetroaffie curvature κ k is differet, ad if it lies o the same lie, the cetroaffie curvature κ k is same I the right of Figure 2, let P k+ Q = t k, P k+ R = t k Thus, [t k,t k ] represets the oriet area of the triagle Pk+ QR, ad [ t k, t k+ ] represets twice the oriet area of the triagle Pk+ QP k+2 These two triagles have a commo edge P k+ Q Similarly, by Eq (3), we ca also obtai the differet secod cetroaffie curvatures o the differet lies Figure 2: Geometric meaig of the curvatures Hece, we ca make use of above coclusios to decide the positio of the poit P k+2 I Figure 3, whe the poit P k+2 lies at differet itersectios of two group parallel lies, the cetroaffie curvature pair {κ k, κ k } is differet O the cotrary, give a pair {κ k, κ k }, the poit 6

7 P k+2 ca be uiquely decided, which lies at the itersectio of two straight lies κ k, κ k For example, i Figure 3, at the itersectio of two straight lies κ k = 2, κ k =, the poit P k+2 is uiquely determied Figure 3: Decidig a curve by the curvatures κ k ad κ k I fact, from Eq (3), we ca obtai the chai structure (32) r k+2 r k+ = κ k ( r k r k )+ κ k ( r k+ r k ) This shows that (33) r k+2 = κ k r k +( κ k κ k ) r k +(+ κ k ) r k+ Usig a simple matrix multiplicatio, it is coveiet to express Eqs (32) ad (33) by ( ) κk (34) ( r k+2 r k+, r k+ r k ) = ( r k+ r k, r k r k ) κ k 0 ad (35) ( r k+2, r k+, r k ) = ( r k+, r k, r k ) + κ k 0 κ k κ k 0 κ k 0 0 Exactly, this expressio ca be cosidered as a state trasitio process similar to a Markov chai, ad we otice that the sum of every colum of the trasitio matrix is If κ k 0, the matrices ( ) κk κ k 0 ad + κ k 0 κ k κ k 0 κ k 0 0 are reversible The from Eqs (34) ad (35) it follows that 0 (36) ( r k+ r k, r k r k ) = ( r k+2 r k+, r k+ r k ) κ k κ k 7 κ k

8 ad 0 0 (37) ( r k+, r k, r k ) = ( r k+2, r k+, r k ) 0 κ k 0 + κ k κ k κ k + κ k κ k, which are the iverse chais of Eqs (34) ad (35) I fact, the sum of every colum of these two trasitio matrices also is We shall ow start a discussio about two discrete plaar curves with same cetroaffie curvatures uder a affie trasformatio or a cetroaffie trasformatio, ad fid how the cetroaffie curvatures affect a discrete plaar curve by the chai structure (35) The followig two propositios tell us that two discrete plaar curves with same cetroaffie curvatures are affie equivalet, ad they are cetroaffie equivalet up to a traslatio trasformatio Propositio 32 Give two sequeces of umber {κ,κ 2, } ad { κ, κ 2, }, where κ k 0, κ k 0, k Z, up to a cetroaffie trasformatio, there exists a discrete plaar curve r : I R 2 such that κ k ad κ k is the first ad secod cetroaffie curvature of the curve r Proof I a plae, give ay two liearly idepedet vector groups { t 0, t } ad { e 0, e }, there must exists a ivertible matrix A of size 2 such that ( e 0, e ) = A( t 0, t ) Sice the cetroaffie curvatures are ivariat uder a affie trasformatio, of course, they also are ivariat uder a cetroaffie trasformatio, we ca choose two fixed liearly idepedet vectors { e 0, e } as the first two taget vectors { t 0 = r() r(0), t = r(2) r()} From Eq (32) we ca obtai e 2, e 3, i tur The we make e 0, e, head ad tail dockig, which is a discrete curve with cetroaffie curvatures κ k ad κ k,(k = 0,,2, ) The above propositio tells us a discrete plaar curve ca be determied by cetroaffie curvatures with respect to a cetroaffie trasformatio I fact, with differet startig poit of e 0, the curve is differet uder the cetroaffie trasformatio I the followig propositio, we kow they are affie equivalet Propositio 33 Assume two discrete plaar curves r(k), r(k) have same cetroaffie curvatures o the correspodig poits, the, there exist a o-degeerate matrix A of size 2 ad a costat vector C such that r(k) = A r(k)+ C, for all k Z, that is, these two curves are affie equivalet Proof Clearly, there is a o-degeerate matrix A of size 2 satisfyig that From Eq (32), we get Oe after aother, it follows that { r() r(0), r(2) r()} = A { r() r(0), r(2) r() } r(3) r(2) = A( r(3) r(2)) r(k +) r(k) = A( r(k +) r(k)),k Z 8

9 A costat vector C is give by C = r(0) A r(0), It is easily see that r(0) = A r(0)+ C, r() = r(0) + r() r(0) = A r(0)+ C +A( r() r(0)) = A r()+ C, r(2) = r() + r(2) r() = A r()+ C +A( r(2) r()) = A r(2)+ C, Clearly, we obtai r(k) = A r(k)+ C, k Z, which completes the proof of the propositio Remark 34 Sice the cetroaffie trasformatio does ot iclude the traslatio trasformatio, the differet choice of the startig poit will geerate differet discrete cetroaffie curves However, with give cetroaffie curvatures, there exists oe ad oly oe curve uder affie trasformatio, that is, two discrete plaar curves with same cetroaffie curvatures are affie equivalet If two discrete plaar curves r(k), r(k) are cetroaffie equivalet, that is, there exists a cetroaffie trasformatio A such that r(k) = A r(k), k Z, these two curves have same cetroaffie curvatures o the correspodig poits O the other had, if we fix a iitial vector r 0, with give cetroaffie curvatures, there exists oly oe discrete plaar curve I the followig, we wat to describe some results that belog to the global differetial geometry of a discrete plaar curve Observe Eq (35), we obtai (38) ( r p+2, r p+, r p ) = ( r 2, r, r 0 ) + κ 0 κ κ p 0 κ κ p 0 κ p κ p 0 κ p 0 0, p Z For a discrete closed curve with period p Z, we have ( r p+2, r p+, r p ) = ( r 2, r, r 0 ) Thus the followig lemma is obvious Lemma 35 A discrete cetroaffie curve is a closed curve with period p Z if ad oly if + κ 0 κ κ 0 κ 0 0 where E is the idetity matrix of size 3 + κ 2 0 κ 2 κ 2 0 κ κ p 0 κ p κ p 0 κ p 0 0 = E,

10 Notice that immediately, we have det + κ k 0 κ k κ k 0 κ k 0 0 = κ k, Corollary 36 If a discrete plaar curve is closed with period p, its first cetroaffie curvature satisfies that κ κ 2 κ p = As we kow, i classical differetial geometry, a closed plaar curve with costat curvature is a circle Naturally, it is iterestig to cosider the similar problems for a discrete plaar curve Firstly, let us give the followig defiitio Defiitio 37 A discrete cetroaffie curve is called costat curvature cetroaffie curve if its first ad secod cetroaffie curvatures are costat I the particular case, accordig to the defiitio of the cetroaffie curvature, a discrete plaar curve withκ = 0 isalie Hece, fromowowe assume κ 0 The followig propositioshows how to get a discrete closed plaar curve depedig o the costat cetroaffie curvatures Propositio 38 A discrete plaar curve with costat cetroaffie curvatures is closed if ad oly if the curvatures κ =, κ 2, ad there exist θ R,p,l Z, such that cosθ = κ, pθ = 2lπ, 2 where p is the period of the discrete closed curve, p ad l are coprime Proof It is easy to see from Lemma 35 ad Corollary 36, that a discrete plaar curve with costat cetroaffie curvatures κ ad κ is closed if ad oly if there exists a iteger p Z such that (39) To take the fact det ito accout, easily we get κ = or κ = + κ 0 κ κ 0 κ κ 0 κ κ 0 κ 0 0 p = E = κ Now if κ =, the trasitio matrix ca be rewritte as + κ 0 κ By a simple calculatio, we obtai its eigevalues λ =,λ 2,3 = κ κ2 2 ± +4 Sice λ 2, λ 3, it is 2 p + κ 0 impossible to fid a iteger p satisfyig that κ 0 = E κ 0 O the other had, if κ =, the trasitio matrix is κ 0, ad its eigevalues 0 0 are λ =,λ 2,3 = κ κ2 2 ± 4 Accordig to Eq (39) we obtai λ = λ 2 = λ 3 =, which 2 0

11 shows that κ 2, at the same time, there is a θ satisfyig that κ κ2 2 ± 4 = cosθ± siθ 2 AgaifromEq (39), wekow(cosθ± siθ) p = exp(± pθ) =, whichimplies pθ = 2lπ, where l Z, p ad l are coprime The the trasitio equatio (35) tells us this is a closed discrete plaar curve with period p O the cotrary, if these coditios hold, Eq (39) is satisfied Together with Eq (35), we ca obtai a discrete closed plaar curve I the classical differetial geometry, it is well kow that a closed regular plaar simple curve is covex if ad oly if its curvature is either always o-egative or always o-positive, ie, if ad oly if the turig agle (the agle of the taget to the curve) is a weakly mootoe fuctio of the parametrizatio of the curve Similarly, for a discrete plaar curve we obtai Propositio 39 A discrete closed plaar simple curve C is covex if ad oly if its first cetroaffie curvature κ k > 0 Proof I fact, a discrete closed plaar simple curve is a plaar polygo Exactly, covexity ad cetroaffie curvatures are ivariat uder the affie trasformatio i the plae Therefore, we ca cosider the problem by affiely trasformig the polygo to a fixed polygo o the Euclidea plae As we kow, a polygo is covex if ad oly if each of its iterior agles has a measure that is strictly less tha π From Eq (3), we obtai κ k = t k+ siα k+ t k siα k, where α k is the k th iterior agle of the polygo Sice 0 < α i < π(i =,2,,p), we get every κ k > 0 O the other had, if κ i > 0(i =,2,,p), we obtai siα,siα 2,,siα p have the same sig Hece, 0 < α i < π(i =,2,,p), which implies the polygo is covex By a simple calculatio, the followig result is obvious Remark 30 If a discrete closed curve is a triagle, the its cetroaffie curvatures satisfy that κ = κ 2 = κ 3 =, κ = κ 2 = κ 3 = If a discrete closed curve is a parallelogram, we have κ = κ 2 = κ 3 = κ 4 =, κ = κ 2 = κ 3 = κ 4 = 0 If a covex polygo has p sides, the its iterior agle sum is give by the equatio (p 2)π Obviously, usig the same method as above ad the graph show i the right of Figure 2, it is easy to prove the followig result Corollary 3 If a p polygo except parallelogram is covex, where p > 3, its secod cetroaffie curvature κ k >, ad there are o more tha 2 o positive secod cetroaffie curvatures κ Furthermore, if κ k 0 ad κ l 0, we must have k l The followig two propositios illustrate how to get a covex closed curve depedig o costat cetroaffie curvatures ad estimate whether it has self-itersectios Propositio 32 A discrete plaar curve with costat cetroaffie curvatures is the simple covex closed curve if there exist θ ad p, where θ R,p Z, satisfyig that cosθ = κ, pθ = 2π 2 ad κ =

12 Proof From Propositio 38, the curve is a closed curve with period p ad p 3 Obviously, we have κ < 2 Sice cosθ = κ, pθ = 2π, p 3, we have κ = or 0 κ < 2 If κ =, it 2 follows that θ = 2π,p = 3, which is a triagle Certaily this is a simple covex closed curve 3 If κ = 0, we obtai θ = π,p = 4 By a simple calculatio, we get κ =, κ = 0 for the 2 square show i Figure 4, which is a simple covex closed curve Accordig to Propositio 33 ad Theorem 2, it is clearly that a discrete plaar curve with κ =, κ = 0 is a simple covex closed curve with period 4 Figure 4: Simple covex closed curves If 0 < κ < 2 ad there exist θ ad p, such that cosθ = κ, pθ = 2π, we make a equilateral 2 polygo of size p as i Figure 4, where poit A is the ceter of the equilateral polygo, P i are the ed poits of vector r(i),i = 0,,, It is easy to check that the first cetroaffie curvature κ = Now let us calculate its secod cetroaffie curvature κ c By usig the otatio T = ( cosθ siθ siθ cosθ ), we have AP k = AP k T, AP k+ = AP k T 2, 2 AP k+2 = AP k T 3

13 The κ c = [ t k, t k+ ] [ t k, t k ] = [ AP k AP k, AP k+2 AP k+ ] [ AP k AP k, AP k+ AP k ] = [(T E) AP k,t 2 (T E) AP k ] [(T E) AP k,t(t E) AP k ] = [ AP k,t 2 AP k ] [ AP k,tap k ] = [ AP k, AP k+ ] [ AP k, AP k ] = Area( AP k P k+ ) Area( APk P k ) = 2cosθ = κ So this equilateral polygo is a simple covex closed curve with costat cetroaffie curvature κ = ad κ Agai from Propositio 33 ad Theorem 2, a discrete plaar curve with costat cetroaffie curvatures κ = ad above κ is a simple covex closed curve Figure 5: Closed curves with self-itersectios Propositio 33 A discrete plaar curve with costat cetroaffie curvatures is closed curve with self-itersectios if κ =, κ 2, θ R,p,l Z, such that cosθ = κ ad pθ = 2lπ,l >, 2 where p ad l are relatively prime Proof Usig the similar methods as Propositio 32, we ca obtai the cetroaffie curvatures of the equilateral polygo of size p show i Figure 5 By a direct calculatio, we obtai κ =, κ = 2 cos θ Clearly, this closed curve has self-itersectios From Propositio 33 ad Theorem 2, a discrete plaar curve with costat cetroaffie curvature κ = ad above κ is a closed 3

14 curve with self-itersectios Accordig to the proofs of Propositio 32 ad Propositio 33, we kow if the period p is eve, A is the symmetric ceter of the curve Sice the partitio ratio of 3 poits is ivariat uder affie mappigs, it is immediate to get Corollary 34 A discrete plae closed curve with costat cetroaffie curvatures is cetrosymmetric if ad oly if its period is eve Fially, we ca defie the affiely regular polygo usig the affie curvatures Defiitio 35 A plaar polygo is a affiely regular polygo with period p if ad oly if it have costat affie curvatures κ =, κ = 2cos 2lπ 2l, where p ad l are relatively prime ad < p p Especially, l =, it is a affiely regular simple polygos (a simple polygo is oe that does ot itersect itself aywhere) 4 Discrete cetroaffie space curve i R 3 Let curve r : I Z R 3 be a cetroaffie discrete curve deoted by C, ad the by the defiitio 28, we have [ r k, t k, t k ] 0, where [ ] deotes the stadard determiat i R 3 I the followig the cetroaffie curvatures ad cetroaffie torsios of a cetroaffie discrete space curve i R 3 will be defied Defiitio 4 The first, secod cetroaffie curvatures ad cetroaffie torsios of the discrete cetroaffie curve r at poit r(k) are defied by (4) κ k := [ r k+, t k, t k+ ] [ r k, t k, t k ], κ k := [ r k+, t k, t k+ ], τ k := [ t k, t k, t k+ ] [ r k, t k, t k ] [ r k, t k, t k ] ByDefiitio4,uderacetroaffietrasformatioR 3 x A x R 3,whereA GL(3,R),it is easy to see that the first, secod cetroaffie curvatures ad cetroaffie torsios are ivariat However, uder a affie trasformatio x A x + C, where C R 3 is a costat vector, the first, secod cetroaffie curvatures ad cetroaffie torsios may chage Hece, we have Propositio 42 The first, secod cetroaffie curvatures ad cetroaffie torsios are cetroaffie ivariats ad ot affie ivariats I fact, [ r k, t k, t k ] = [ r k, r k, r k+ ] The the cetroaffie curvatures ad torsios ca be rewritte as (42) κ k = [ r k, r k+, r k+2 ] [ r k, r k, r k+ ], κ k = [ r k+, t k, r k+2 ] [ r k, r k, r k+ ], τ k = [ t k, t k, t k+ ] [ r k, r k, r k+ ] By a direct calculatio, it follows that (43) r k+2 = κ k r k +( κ k κ k ) r k +(τ k + κ k +) r k+, k Z, ad (44) ( r k+2, r k+, r k ) = ( r k+, r k, r k ) τ k ++ κ k 0 κ k κ k 0 κ k 0 0, 4

15 which are the three dimesioal curve chai structures This formula is called the Freet-Serret formula of a discrete cetroaffie curve O the other had, whe τ = 0, from Eq (43) we get the chai of edge taget vector (45) t k+ = κ k t k + κ k t k, k Z, which is coicidet with Eq (32) If κ k 0, we otice that the iverse chai ca be represeted as (46) ( r k+, r k, r k ) = ( r k+2, r k+, r k ) 0 0 κ k 0 τ k++ κ k κ k 0 + κ k κ k I this sectio, we oly cosider the discrete cetroaffie space curve C uder the cetroaffie trasformatio Firstly, the followig propositio states that with the give first, secod cetroaffie curvatures ad cetroaffie torsios, the curve is oly determied up to a cetroaffie trasformatio Propositio 43 Two curve C ad C are cetroaffie equivalet if ad oly if they have same cetroaffie curvatures κ k, κ k ad torsios τ k, for all k I Z Proof It is easy to see from Eq (4) that if the curves C ad C satisfy that r(k) = A r(k), where A is a 3 3 matrix, they have same curvatures κ k, κ k ad torsios τ k O the other had, if curves C ad C have same cetroaffie curvatures ad torsios at correspodig poits, we eed to show they are cetroaffie equivalet Obviously there exist a matrix A of size 3 such that ( r(0), r(), r(2)) = A( r(0), r(), r(2)) From Eq (43), by the same cetroaffie curvatures κ k, κ k ad torsios τ k, it is simple to prove that r(k) = A r(k), k I Z This meas the curves C ad C are cetroaffie equivalet Next, we will cosider the geometric iterpretatio for the cetroaffie curvatures ad cetroaffie torsios by figures We deote the ed poit of vector r(k) by P k ad the plaar icludig the poits P k,p k ad P k+ by π 0 I Figure 6, accordig to Eq (4) we kow if the poit P k+2 lies differet plaes which parallel to the plae π 0, the torsios τ k are differet If the poit P k+2 lies differet place i a same plae which parallels to the plae π 0, the torsios τ k are same Similarly, i the left of Figure 7, let π represet the plae cotaiig the poits P k,p k+ ad the origi O By Eq (4) we obtai if P k+2 lies differet plaes which parallel to the plae π, the curvatures κ k are differet If P k+2 lies differet place i a same plae which parallels to the plae π, the cetroaffie curvatures κ k are same I the right of Figure 7, we give a vector P k+ P k = P k P k The [ OP k, OP k, OP k+ ] = [ OP k, OP k+, OP k ] Let π 2 represet the plae cotaiig the poits P k,p k+ ad the origi O Agai from Eq (4), we coclude that if P k+2 lies differet plaes which parallel to the plae π 2, the curvatures κ k are differet If P k+2 lies differet place i a same plae which parallels to the plae π 2, the cetroaffie curvatures κ k are same 5

16 Figure 6: Cetroaffie torsios i differet plaes Figure 7: Cetroaffie curvatures i differet plaes 6

17 From the defiitio ad above geometric iterpretatio, we kow if τ k = 0, the curve is a plaar curve show i Figure 8 I this case, by compariso, we ca fid the defiitios of cetroaffie curvatures are coicidet with the plae situatio i the above sectio I Figure 8, let P k+2 P = P k P k Because the poits P k,p k,p k+,p k+2 are i the same plae π, so is the poit P Hece, the first cetroaffie curvature κ k = [ r k+, t k, t k+ ] is the ratio of the oriet area [ r k, t k, t k ] of the triagle Pk P k+ P k+2 ad the oriet area of the triagle Pk P k P k+ At the same time the secod cetroaffie curvature κ k = [ r k+, t k, t k+ ] is the ratio of the oriet area of the triagle [ r k, t k, t k ] Pk+ P k+2 P ad the oriet area of the triagle P k P k P k+ Exactly, their geometry meaig is same as a plaar curve defied i the above sectio Therefore we have Figure 8: Cetroaffie plaar curve i R 3 Remark 44 If τ k = 0, k Z, the curve C is a plaar curve ad the cetroaffie curvatures i Eq (4) are as same as defied i Eq (3) for the plaar curve A discrete space curve with period p is closed if ad oly if r(k +p) = r(k), k Z Through the three dimesioal curve chai, that is, Eq (44), it is easy to see Lemma 45 A discrete cetroaffie space curve is a closed curve with period p if ad oly if the cetroaffie curvatures ad torsios satisfy that + κ +τ 0 + κ 2 +τ κ p +τ p 0 κ κ 0 κ 2 κ 2 0 κ p κ p 0 = E, κ 0 0 κ κ p 0 0 where E is the idetity matrix of size 3 Observe that It is immediate to obtai det + κ k +τ k 0 κ k κ k 0 κ k 0 0 = κ k Corollary 46 If a discrete cetroaffie curve is closed with period p, the κ κ 2 κ p = I the above sectio, we have cosidered a discrete plaar curve with costat cetroaffie curvatures Similarly, we ca obtai 7

18 Propositio 47 Let the curve C be a discrete cetroaffie closed p polygo with costat cetroaffie curvatures κ, κ ad torsios τ, the we have κ = or κ = Moreover, if κ =, the curve is plaar curve If κ =, we have 0 < κ < 4,τ < 0, τ = 2 κ, ad there exist a real umber θ ad a iteger l which is relatively prime to p satisfyig that cosθ = κ,pθ = 2lπ, 2 where p is eve umber Proof Sice C is a discrete cetroaffie closed curve with costat cetroaffie curvatures κ, κ ad torsios τ, we have from Corollary 46 that κ = or κ = Obviously, if κ =, p is eve umber + κ+τ 0 Assume λ,λ 2,λ 3 are the eigevalues of the matrix κ κ 0 The eigevalue κ 0 0 equatio of this matrix is (47) λ 3 (τ + κ+)λ 2 +(κ+ κ)λ k = 0 Hece we obtai (48) λ +λ 2 +λ 3 = τ + κ+, λ λ 2 +λ 2 λ 3 +λ 3 λ = κ+ κ, λ λ 2 λ 3 = κ p + κ+τ 0 That the curve C is closed, that is, κ κ 0 = E, implies λ,λ 2,λ 3 are ot equal κ 0 0 to each other ad (49) λ p = λ p 2 = λ p 3 = If κ =, from Eqs (48) ad (49) we ca assume λ =, λ 2 = cosθ+ siθ, λ 3 = cosθ siθ, where pθ = 2lπ, p, l are relatively prime Immediately Eq (48) geerates λ 2 +λ 3 = τ + κ, λ 2 +λ 3 = κ, which implies τ = 0 ad the curve C is a plaar curve If κ =, from Eqs (48) ad (49) we ca assume λ =, λ 2 = cosθ+ siθ, λ 3 = cosθ siθ, where pθ = 2lπ ad p,l are relatively prime Agai from Eq (48) we get Fially we obtai λ 2 +λ 3 = τ + κ+2, λ 2 +λ 3 = 2 κ τ = 2 κ, cosθ = κ 2 The we complete the proof From the above proof, it is clear that if κ =, r(k + p ) = r(k), which implies the curve is 2 cetrosymmetric ad the ceter of symmetry is the origi O Hece, we have Corollary 48 If a discrete cetroaffie closed curve with costat cetroaffie curvatures ad torsios is ot a plaar curve, it is symmetric aroud the origi O Corollary 49 If a discrete cetroaffie curve with costat cetroaffie curvatures ad torsios is closed with period p, where p is a eve umber, the it is cetrosymmetric 8

19 5 Flows o curves Curve-shorteig flowisthesimplest exampleofacurvatureflow Itmoves eachpoitoaplaar curve γ i the iwards ormal directio ν with speed proportioal to the siged curvature κ at that poit, as described by the equatio γ = κ ν The ame curve-shorteig comes from t the fact that the curve is always movig so as to decrease its legth as efficietly as possible I this sectio, the geeral form of the flow o curves is cosidered Furthermore, it ca be exteded to the discrete curves See [7] for more details 5 Flows o smooth curves The motio of a curve r : I R N i space could be described by applyig some vector field v I geeral v might deped o the whole curve Exactly, if v depeds oly o a small eighborhood at each poit of the curve, we call the geerated flow a local flow With these coditios, the evolutio process of r uder the flow geerated by v ca be described by a differetial equatio (5) t r = v( r, r, r, ) A oe-parameter family of curves (52) r : I J R N which is a solutio of Eq (5) i the sese that (53) t r(s,t) = v( r(s,t), r (s,t), r (s,t), ) =: v(s,t) for all (s,t) I J is called the evolutio of the curve r 0 (s) = r(s,0) uder the flow give by v For this particular iitial curve r 0, the vector field v becomes a o-parameter family of vector fields alog the parametrizatio (54) v : I J R N, ad Eq (5) becomes (55) t r(s,t) = v(s,t) The map (56) Φ : (t, r) Φ t r, where Φ t r(s) = r(s,t) is the evolutio of r is called the curve flow give by v Note that i geeral Φ might ot be well defied due to lack of existece ad uiqueess of solutios of Eq (5) for arbitrary r 0 Additioally if oe might wat the flow to be geometric, ie oly deped o the shape of the curve It should be ivariat with respect to Euclidea motios, reparametrizatio of the curve The flow is the well defied o the correspodig equivalece classes of parametrized curves Example(plaar geometric flow) For plaar curves these two coditios ca be realized by the asatz: v = v(κ,κ,κ, ) = α(κ,κ,κ, ) T +β(κ,κ,κ, ) N Exactly, the flow ca be geeralized i affie space if the flow is ivariat with respect to affie trasformatio 9

20 52 Flows o discrete curves The previous sectio was devoted to the cotiuous case We ow study the discrete case For I := [0,,] Z fiite iterval, I = Z := Z/Z ad I = Z, we defie the space (57) C I := { r : I R N } of fiite,fiite closed ad ifiite curves respectively A flow of discrete curves is give by a vector field (58) v : C I TC I, r v [ r ] T r C I, or o some submaifold U C I I the fiite case we have T r C I = (R N ) So v gives a directio i R N at every vertex k which possibly depeds o the whole curve r We state this relatio as (59) v k [ r ] R N For a give iitial curve r : I R N the vector field v o C I becomes a oe-parameter family of vector fields alog the parametrizatio of the curve (50) v : I J R N, where J R is a ope iterval The actio of the flow leads to a cotiuous deformatio of the curve r k = r k (0), (5) r : I J R N satisfyig that (52) t r(t) = v k (t) By a local flow we mea a flow which at every vertex k oly depeds o the curve at the adjacet vertices, ie r k, r k, r k+ : (53) v k [ r ] = v( r k, r k, r k+ ) I affie geometry, the result of a flow is very differet A tiy curve segmet ca be affiely equivalet to a huge oe So we defie its stability depedig o its affie ivariats, such as cetraffie curvatures ad torsios Defiitio 5 A discrete curve is stable if it remais its cetroaffie curvatures ad the torsios uchaged i the ext descedat A discrete curve is periodically stable if oe of its descedats have the same cetroaffie curvatures ad torsios as that 6 Trasversal flow o discrete cetroaffie space curves Now we exted the flow to a cetroaffie space curve, which implies the flow is ivariat with respect to the affie trasformatio or cetroaffie trasformatio As metioed previously, the cetroaffie curvatures are affie ivariat whe the discrete curve lies i 2-dimesioal plae Whe the curve is plaar, the results are suitable for the affie trasformatio Sice the cetroaffie curvatures are coicidet for a discrete plaar curve o matter is i 2-dimesioal plae or i 3-dimesioal space, here we cosider flows o the discrete curve i 3-dimesioal 20

21 space We draw a aalogy betwee the two cases by replacig a curve with a discrete curve, ad the time with a discrete time For a cetroaffie space curve, its positio vector field is always trasversal to the osculatio plae The discrete trasversal flow ca be defied by (6) t r k = v k := α k r k, where α k = α k [ r ] should deped o the curve i ay way ad be cetroaffiely ivariat Let r m R3 be a discrete cetroaffie curve, where is the idex of the vertices ad m is the discrete deformatio parameter For a local discrete trasversal flow r m, we have (62) r m+ r m = αm r m, where α m is cetroaffiely ivariat This equatio ca also be writte as (63) r m+ = (+α m ) r m By takig the otatio (64) β m = +α m, the trasfer equatio is easily show (65) ( r m+ m+ +, r, r m+ ) = ( r m +, r m, r m )Mm, β m where M m = β m β m I fact, we have obtaied the structure equatio (35) for a discrete cetroaffie space curve i the previous sectio By a simplificatio, the followig otatio is used (66) L m = + κ m +τm 0 κ m κm 0 κ m 0 0, (L m ) = 0 0 κ m 0 + κm +τm 0 κ m κ m + κ m κ m Now, by usig the compatibility coditio of the liear system (35) ad (65) it is o difficult to get L m M m + = M m L m+, (67) (68) (69) β+2 m (+ κm +τm ) = βm + (+ κm+ +τ m+ ), β m +2(κ m + κ m ) = β m (κ m+ + κ m+ ), β m +2 κm = β m κm+ 2

22 Obviously, if β m m+ = 0, by Eq (63), we have r = 0, which is cotrary to the defiitio of the discrete cetroaffie curve So i the followig we should assume β 0 By a direct computatio, from Eqs (67)-(69), it shows the cetroaffie curvatures ad torsios of ext geeratio vertex r m+ are (60) (6) τ m+ = βm +2 β m κ m+ = βm +2 β m κ m + βm +2 β m + κ m, (+ κ m +τ m ) βm +2 (κ m β m + κ m ), (62) κ m+ = βm +2 β+ m Therefore, we may coclude that (κ m + κ m ) βm +2 κ m β m Propositio 6 Uder a discrete trasversal motio, the curve r m ad r m+ have the relatio (60)-(62) I order to be more clear, we could chage Eqs (60)-(62) to a matrix form (63) τ m+ κ m+ κ m+ = β m +2 β m + 0 β+2 m βm +2 β+ m β m β m +2 β m β+2 m βm +2 β m β m β+2 m βm β+ m +2 β m β m +2 β m β m +2 β m τ m κ m κ m Now, the followig propositio shows that if β is costat, the trasversal flow of a discrete curve would keep stable Propositio 62 If β is costat, the the discrete trasversal motio r m+ of a discrete curve r m is cetroaffiely equivalet to the curve r m Whe we assume κ m +κ m 0, β is costat if ad oly if the discrete trasversal motio r m+ of a discrete curve r m is cetroaffiely equivalet to the curve r m Proof If β is costat, from Eq (63) it is easy to get τ m+ = τ m, κ m+ = κ m,κ m+ = κ m, which implies the curve r m+ is cetroaffie equivalet to the curve r m O the other had, if the curve r m+ is cetroaffie equivalet to the curve r m, the τ m+ = τ m, κ m+ = κ m,κ m+ = κ m From Eq (6), we get β m = β+2 m The by Eq (62) ad κ m +κ m 0, we have β+ m = βm +2 Agai from Eq (60) we get βm = βm + Hece, we obtai β m is costat The by observig Eq (60), we ca state the followig propositio 22

23 Propositio 63 If the discrete trasversal motio r m+ of a discrete plaar curve r m is still plaar, that is, τ m+ = τ m = 0, the coefficiets of motio satisfy that (64) or (65) β m +2 β m +2 = = ( β m β m + β m + (κ m + κ m )+, β m, = κ m ( β m β m + β m κ m + (+ κ m β ) + m )(+ κ m, κm κm,κm )Tra, β m ) κ m ( β m ) β m Especially, they ca be represeted as the trasitio trasformatio + κ m 0 (66) (,, ) = (,, ) κ m β+2 m β+ m β m β+ m β m β m κm 0 κ m 0 0 ad (67) ( κ m+ κ m+ ) = β m +2 β m + 0 β+2 m βm +2 β+ m β m β+2 m β m ( κ m κ m ) Furthermore, accordig to Eq (64), the followig propositio is obvious Propositio 64 For a plaar p polygo, if its ext geeratio is still plaar uder the trasversal motio, the coefficiets of motio β m,β m 2,,β m p should satisfy that (68) S β m β m 2 β m 3 β m p = 0, where S = κ m κm + κ m 0 0 κ m κ m 2 κ m 2 κm 2 + κ m κ m 3 κ m 3 κ m 3 + κ m κ m 4 κ m 4 κm κ m p κ m p + κ m p + κ m p 0 0 κ m p κ m p κ m p I order to obtai the solutios of the above equatio, firstly, we eed to study the rak of coefficiet matrix S Propositio 65 The rak of matrix S is p 3 23

24 Proof Whe τ = 0, from Eq (44), we have S ( r m, r m 2,, r m p ) Tra = 0, which implies the colum vectors of ( r m, r m 2,, r m p )Tra are three solutios of liear equatio Sx = 0 Sice r m k, r m k+, r m k+,(k =,2,,p 2) are liearly idepedet, it is easy to see that the rak of ( r m, r m 2,, r ) m p is 3 If there exists a vector group ( r, r 2,, r p ) satisfyig that S ( ) Tra r, r 2,, r p = 0 From Propositio 43, we kow there must be a matrix A of size 3 satisfyig that r k = A r m k (k =,2,,p) By the traspositio, it follows that ( r, r 2,, r p ) Tra = ( r m, r m 2,, r m p ) TraA Tra, which implies arbitrary solutio of liear equatio Sx = 0 ca be liearly represeted by colum vectors of ( r m, r m 2,, r m p ) Tra If we assume the colum vectors of ( r m, r m 2,, r m p ) Tra are V,V 2,V 3, the vector group {V,V 2,V 3 } is a system of fudametal solutio for liear equatio Sx = 0 Hece, it is clear that the rak of matrix S is p 3 Now from Eq (68) ad the proof of Propositio 65, we have β m β m 2 β m 3 β m p = a V +a 2 V 2 +a 3 V 3, where a,a 2,a 3 are arbitrary costat Ifact, β m should becetroaffie ivariat, but V,V 2,V 3 areot So we eed to fidasystem of fudametal solutios for liear equatio Sx = 0, which are cetroaffie ivariat As we kow, there exists a matrix A, such that A r m = (,0,0)Tra,A r m 2 = (0,,0)Tra,A r m 3 = (0,0,)Tra Hece, we obtai a stadard vector group r m i = A r m i,(i =,2,,p) which are ivariat uder cetroaffie trasformatio Now assume the colum vectors of ( r m, r m 2,, r m p ) Tra are V, V 2, V 3, ad the vector group { V, V 2, V 3 } is a system of fudametal solutios for liear equatio Sx = 0 ad cetraffie ivariat Thus we obtai Corollary 66 If a trasversal motio of a plaar p polygo r m remais plaar, the coefficiets β m satisfy that β m β m 2 β m 3 β m p = a m V +a m V 2 2 +a m V 3 3, where a m,a m 2,a m 3 are some cetroaffie ivariats that ca esure β m i 0(i =,2,,p) 24

25 Hece, if a m = am 2 = am 3, we have βm = β2 m = = βp m = a m Accordig to Propositio 62, this implies the discrete curve r m is stable uder the above trasversal motio Computer experimets show that there are some plaar polygos whose trasversal flows will reach stable I the followig, we give a example, which shows after multistep iteratio the coefficiets a m,am 2,am 3 will be same Example of stable trasversal flow Let us take a m = p i= κm i p, a m 2 = p i= κm i p, a m 3 = am +am 2 2 It is obvious that a m,a m 2,a m 3 are cetroaffie ivariat Usig computer experimets, we choose a covex plaar petago, which will keep plaar durig the iteratio process I Table, we ca see that, after the thirty-eighth step, this trasversal cetroaffie flow will reach stable That is, from the o, the mea value of the first cetroaffie curvatures is equal to that of the secod cetroaffie curvatures Table : A stable trasversal flow Iitial poits (0,22,) (8,2,) (2,0,) (37,2,) (48,28,) a a 2 κ κ κ κ κ κ κ κ Taget flow o discrete cetroaffie curves Sice taget vectors of a discrete cetroaffie curve r live o edges, it is ot istatly clear what the taget directio at a vertex k should be If we wat it to deped oly o the eighborig taget vectors, a obvious symmetric choice would be t k + t k However, if the motios of the discrete curve always lie o the taget plae, we ca defie the motios as the discrete taget flow of discrete cetroaffie curves, that is, (7) t r k = v k := α k t k +β k t k, where α k = α k [ r ] ad β k = β k [ r ] should deped o the curve i ay way ad be cetroaffie ivariat Now let r m R 3 be a discrete cetroaffie space curve, where is the idex of the vertices ad m is the discrete deformatio parameter For a local discrete taget flow r m, it is ature to see (72) r m+ r m = αm t m +β m t m = αm r m + +(βm αm +) r m βm r m, where α m ad β m are cetroaffie ivariat This equatio ca also be writte as (73) r m+ = ( r m +, r m, r m )(αm, αm +βm, βm )Tra 25

26 The by Eqs (35) ad (66) we get ad r m+ + = ( r m +2, r m +, r m )(α m +, α m + +β m +, β m +) Tra = ( r m +, r m, r m )Lm (αm +, αm + +βm +, βm + )Tra, r m+ = ( r m, r m, r m 2 )(αm, αm +βm, βm )Tra = ( r m +, r m, r m )(Lm ) (α m, αm +βm, βm )Tra Obviously, the three equatios above ca be regarded as state trasitio equatios (74) ( r m+ m+ +, r, r m+ ) = ( r m +, r m, r m )Mm, where M m = +β+ m +αm + (τm + κm ) αm βm κ m α+ m (κm + κm ) βm + β m αm + αm + βm κ m (+τ m + κm ) α m + κm β m α m βm κm κ m Similarly, the compatibility coditio of the liear system (35) ad (74) L m Mm + = Mm Lm+ yields (75) M m + κ m+ +τ m+ κ m+ κ m+ κ m+ = L m +α m +2(τ m + + κ m +)+β m +2 β m +2 αm +2 (κm + + κm + ) α m +2 κm + By Eq (74) it is easy to see M m is o-degeerate So it follows that + κ m+ +τ m+ (76) = (M m ) L m κ m+ κ m+ κ m+ More clearly, it ca be writte as (77) τ m+ κ m+ κ m+ Ideed, we have =α+2 m (M m ) L m (M m ) L m +α m +2 (τm + + κm + )+βm +2 β m +2 α m +2(κ m + + κ m +) α m +2 κm β m +2 β m +2 0 τ m + κ m + Remark 7 By the defiitio of discrete taget flow, it is easy to see the taget flows of a plaar discrete curve are still plaar O the other had, this result ca also be show by usig that the sum of elemets of every colum vector is i Eq (76) Hece, if τ m = 0, it is clear that τ m+ = κ m +

27 7 Iteratio of defiite proportioal divisio poit The midpoit map is perhaps the simplest polygo iteratio Startig with a N polygo P, we create a ew N polygo P 2 whose vertices are midpoits of the edges of P For almost every choice of P, if we iterate this process, the obtaied sequece of polygos P k will coverge to its cetroid Furthermore, Berlekamp et al cocluded that the descedats of a plaar polygo approach a affiely regular polygo The results ca be exteded to the more geeral trasformatio with the formula r m+ = a 0 r m +a r m + + +a N r m +N, where a 0,a,a N are ay costats ad where the N subscripts o the r m +k are to be computed modulo N[2] For o-plaar polygos, they obtaied almost all o-plaar polygos lack plaar descedats If the first descedat is o-plaar the so are all the rest O the other had, all o-plaar polygos have descedats which differ arbitrarily little (relative to their size) from plaar polygos Ad almost all o-plaar polygos have descedats which differ arbitrarily little (relative to their size) from plaar covex polygos[2] Now we ca visually show these results by its cetroaffie curvatures ad torsios By Defiitio35adRemark44, τ = 0represets theplaardiscrete curve, adκ=, κ = 2cos 2lπ N,τ = 0 represet the plaar affiely regular polygos To make the results above more clearly, we ca use the taget flow o a polygo with the related cetroaffie curvatures ad torsios Exactly, the iteratio of polygos with defiite proportioal divisio poit ca be described by r m+ = ( α) r m +α r m +, where α is costat ad 0 < α < Usig Eq (7), we have α m = α, βm = 0 By usig Eq (77), we get the iteratio of torsios ad curvatures for the taget flow r m, which may be writte as (78) (79) (70) τ m+ κ m+ κ m+ = (α α2 ) κ m +α2 κ m +( α)2 α(α ) α( α)2 τ m + +(α α 2 ) κ m + +α 2 κ m + +( α) 2 α( α) 2 τ m +(α α 2 ) κ m +α 2 κ m +( α) 2 α( α)τm + +(α α 2 ) κ m + +α 2 κ m + +( α) 2, α(α ) = ( κ m + α τ m α κm )α( α)2 + +(α α2 ) κ m + +α2 κ m + +( α)2 α( α) 2 τ m +(α α2 ) κ m +α2 κ m +( α)2 α = κ m α κm +, α( α) 2 τ m + +(α α2 ) κ m + +α2 κ m + +( α)2 α( α) 2 τ m +(α α2 ) κ m +α2 κ m +( α)2 It ca be cocluded from Eqs (78), (79) ad (70) that Propositio 72 Uder iteratios of defiite proportioal divisio poit, a space discrete cetroaffie curve with costat cetroaffie curvatures ad torsios is stable Especially, a affiely regular polygo is stable By previous coclusios[2], we have 27

CAUCHY'S FORMULA AND EIGENVAULES (PRINCIPAL STRESSES) IN 3-D

CAUCHY'S FORMULA AND EIGENVAULES (PRINCIPAL STRESSES) IN 3-D GG303 Lecture 19 11/5/0 1 CAUCHY'S FRMULA AN EIGENVAULES (PRINCIPAL STRESSES) IN 3- I II Mai Topics A Cauchy s formula Pricipal stresses (eigevectors ad eigevalues) Cauchy's formula A Relates tractio vector

More information

5. Best Unbiased Estimators

5. Best Unbiased Estimators Best Ubiased Estimators http://www.math.uah.edu/stat/poit/ubiased.xhtml 1 of 7 7/16/2009 6:13 AM Virtual Laboratories > 7. Poit Estimatio > 1 2 3 4 5 6 5. Best Ubiased Estimators Basic Theory Cosider agai

More information

EXERCISE - BINOMIAL THEOREM

EXERCISE - BINOMIAL THEOREM BINOMIAL THOEREM / EXERCISE - BINOMIAL THEOREM LEVEL I SUBJECTIVE QUESTIONS. Expad the followig expressios ad fid the umber of term i the expasio of the expressios. (a) (x + y) 99 (b) ( + a) 9 + ( a) 9

More information

A New Constructive Proof of Graham's Theorem and More New Classes of Functionally Complete Functions

A New Constructive Proof of Graham's Theorem and More New Classes of Functionally Complete Functions A New Costructive Proof of Graham's Theorem ad More New Classes of Fuctioally Complete Fuctios Azhou Yag, Ph.D. Zhu-qi Lu, Ph.D. Abstract A -valued two-variable truth fuctio is called fuctioally complete,

More information

The Limit of a Sequence (Brief Summary) 1

The Limit of a Sequence (Brief Summary) 1 The Limit of a Sequece (Brief Summary). Defiitio. A real umber L is a it of a sequece of real umbers if every ope iterval cotaiig L cotais all but a fiite umber of terms of the sequece. 2. Claim. A sequece

More information

Sequences and Series

Sequences and Series Sequeces ad Series Matt Rosezweig Cotets Sequeces ad Series. Sequeces.................................................. Series....................................................3 Rudi Chapter 3 Exercises........................................

More information

point estimator a random variable (like P or X) whose values are used to estimate a population parameter

point estimator a random variable (like P or X) whose values are used to estimate a population parameter Estimatio We have oted that the pollig problem which attempts to estimate the proportio p of Successes i some populatio ad the measuremet problem which attempts to estimate the mea value µ of some quatity

More information

A random variable is a variable whose value is a numerical outcome of a random phenomenon.

A random variable is a variable whose value is a numerical outcome of a random phenomenon. The Practice of Statistics, d ed ates, Moore, ad Stares Itroductio We are ofte more iterested i the umber of times a give outcome ca occur tha i the possible outcomes themselves For example, if we toss

More information

Chapter 5: Sequences and Series

Chapter 5: Sequences and Series Chapter 5: Sequeces ad Series 1. Sequeces 2. Arithmetic ad Geometric Sequeces 3. Summatio Notatio 4. Arithmetic Series 5. Geometric Series 6. Mortgage Paymets LESSON 1 SEQUENCES I Commo Core Algebra I,

More information

Rafa l Kulik and Marc Raimondo. University of Ottawa and University of Sydney. Supplementary material

Rafa l Kulik and Marc Raimondo. University of Ottawa and University of Sydney. Supplementary material Statistica Siica 009: Supplemet 1 L p -WAVELET REGRESSION WITH CORRELATED ERRORS AND INVERSE PROBLEMS Rafa l Kulik ad Marc Raimodo Uiversity of Ottawa ad Uiversity of Sydey Supplemetary material This ote

More information

The Valuation of the Catastrophe Equity Puts with Jump Risks

The Valuation of the Catastrophe Equity Puts with Jump Risks The Valuatio of the Catastrophe Equity Puts with Jump Risks Shih-Kuei Li Natioal Uiversity of Kaohsiug Joit work with Chia-Chie Chag Outlie Catastrophe Isurace Products Literatures ad Motivatios Jump Risk

More information

Estimating Proportions with Confidence

Estimating Proportions with Confidence Aoucemets: Discussio today is review for midterm, o credit. You may atted more tha oe discussio sectio. Brig sheets of otes ad calculator to midterm. We will provide Scatro form. Homework: (Due Wed Chapter

More information

Lecture 4: Parameter Estimation and Confidence Intervals. GENOME 560 Doug Fowler, GS

Lecture 4: Parameter Estimation and Confidence Intervals. GENOME 560 Doug Fowler, GS Lecture 4: Parameter Estimatio ad Cofidece Itervals GENOME 560 Doug Fowler, GS (dfowler@uw.edu) 1 Review: Probability Distributios Discrete: Biomial distributio Hypergeometric distributio Poisso distributio

More information

1 Estimating sensitivities

1 Estimating sensitivities Copyright c 27 by Karl Sigma 1 Estimatig sesitivities Whe estimatig the Greeks, such as the, the geeral problem ivolves a radom variable Y = Y (α) (such as a discouted payoff) that depeds o a parameter

More information

Statistics for Economics & Business

Statistics for Economics & Business Statistics for Ecoomics & Busiess Cofidece Iterval Estimatio Learig Objectives I this chapter, you lear: To costruct ad iterpret cofidece iterval estimates for the mea ad the proportio How to determie

More information

Limits of sequences. Contents 1. Introduction 2 2. Some notation for sequences The behaviour of infinite sequences 3

Limits of sequences. Contents 1. Introduction 2 2. Some notation for sequences The behaviour of infinite sequences 3 Limits of sequeces I this uit, we recall what is meat by a simple sequece, ad itroduce ifiite sequeces. We explai what it meas for two sequeces to be the same, ad what is meat by the -th term of a sequece.

More information

We learned: $100 cash today is preferred over $100 a year from now

We learned: $100 cash today is preferred over $100 a year from now Recap from Last Week Time Value of Moey We leared: $ cash today is preferred over $ a year from ow there is time value of moey i the form of willigess of baks, busiesses, ad people to pay iterest for its

More information

Combining imperfect data, and an introduction to data assimilation Ross Bannister, NCEO, September 2010

Combining imperfect data, and an introduction to data assimilation Ross Bannister, NCEO, September 2010 Combiig imperfect data, ad a itroductio to data assimilatio Ross Baister, NCEO, September 00 rbaister@readigacuk The probability desity fuctio (PDF prob that x lies betwee x ad x + dx p (x restrictio o

More information

Models of Asset Pricing

Models of Asset Pricing APPENDIX 1 TO CHAPTER 4 Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see

More information

5 Statistical Inference

5 Statistical Inference 5 Statistical Iferece 5.1 Trasitio from Probability Theory to Statistical Iferece 1. We have ow more or less fiished the probability sectio of the course - we ow tur attetio to statistical iferece. I statistical

More information

Today: Finish Chapter 9 (Sections 9.6 to 9.8 and 9.9 Lesson 3)

Today: Finish Chapter 9 (Sections 9.6 to 9.8 and 9.9 Lesson 3) Today: Fiish Chapter 9 (Sectios 9.6 to 9.8 ad 9.9 Lesso 3) ANNOUNCEMENTS: Quiz #7 begis after class today, eds Moday at 3pm. Quiz #8 will begi ext Friday ad ed at 10am Moday (day of fial). There will be

More information

Models of Asset Pricing

Models of Asset Pricing APPENDIX 1 TO CHAPTER4 Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see

More information

Random Sequences Using the Divisor Pairs Function

Random Sequences Using the Divisor Pairs Function Radom Sequeces Usig the Divisor Pairs Fuctio Subhash Kak Abstract. This paper ivestigates the radomess properties of a fuctio of the divisor pairs of a atural umber. This fuctio, the atecedets of which

More information

APPLICATION OF GEOMETRIC SEQUENCES AND SERIES: COMPOUND INTEREST AND ANNUITIES

APPLICATION OF GEOMETRIC SEQUENCES AND SERIES: COMPOUND INTEREST AND ANNUITIES APPLICATION OF GEOMETRIC SEQUENCES AND SERIES: COMPOUND INTEREST AND ANNUITIES Example: Brado s Problem Brado, who is ow sixtee, would like to be a poker champio some day. At the age of twety-oe, he would

More information

Maximum Empirical Likelihood Estimation (MELE)

Maximum Empirical Likelihood Estimation (MELE) Maximum Empirical Likelihood Estimatio (MELE Natha Smooha Abstract Estimatio of Stadard Liear Model - Maximum Empirical Likelihood Estimator: Combiatio of the idea of imum likelihood method of momets,

More information

Lecture 4: Probability (continued)

Lecture 4: Probability (continued) Lecture 4: Probability (cotiued) Desity Curves We ve defied probabilities for discrete variables (such as coi tossig). Probabilities for cotiuous or measuremet variables also are evaluated usig relative

More information

INTERVAL GAMES. and player 2 selects 1, then player 2 would give player 1 a payoff of, 1) = 0.

INTERVAL GAMES. and player 2 selects 1, then player 2 would give player 1 a payoff of, 1) = 0. INTERVAL GAMES ANTHONY MENDES Let I ad I 2 be itervals of real umbers. A iterval game is played i this way: player secretly selects x I ad player 2 secretly ad idepedetly selects y I 2. After x ad y are

More information

Models of Asset Pricing

Models of Asset Pricing 4 Appedix 1 to Chapter Models of Asset Pricig I this appedix, we first examie why diversificatio, the holdig of may risky assets i a portfolio, reduces the overall risk a ivestor faces. The we will see

More information

1 Random Variables and Key Statistics

1 Random Variables and Key Statistics Review of Statistics 1 Radom Variables ad Key Statistics Radom Variable: A radom variable is a variable that takes o differet umerical values from a sample space determied by chace (probability distributio,

More information

Bayes Estimator for Coefficient of Variation and Inverse Coefficient of Variation for the Normal Distribution

Bayes Estimator for Coefficient of Variation and Inverse Coefficient of Variation for the Normal Distribution Iteratioal Joural of Statistics ad Systems ISSN 0973-675 Volume, Number 4 (07, pp. 7-73 Research Idia Publicatios http://www.ripublicatio.com Bayes Estimator for Coefficiet of Variatio ad Iverse Coefficiet

More information

An Empirical Study of the Behaviour of the Sample Kurtosis in Samples from Symmetric Stable Distributions

An Empirical Study of the Behaviour of the Sample Kurtosis in Samples from Symmetric Stable Distributions A Empirical Study of the Behaviour of the Sample Kurtosis i Samples from Symmetric Stable Distributios J. Marti va Zyl Departmet of Actuarial Sciece ad Mathematical Statistics, Uiversity of the Free State,

More information

Outline. Plotting discrete-time signals. Sampling Process. Discrete-Time Signal Representations Important D-T Signals Digital Signals

Outline. Plotting discrete-time signals. Sampling Process. Discrete-Time Signal Representations Important D-T Signals Digital Signals Outlie Discrete-Time Sigals-Itroductio. Plottig discrete-time sigals. Samplig Process. Discrete-Time Sigal Represetatios Importat D-T Sigals Digital Sigals Discrete-Time Sigals-Itroductio The time variable

More information

x satisfying all regularity conditions. Then

x satisfying all regularity conditions. Then AMS570.01 Practice Midterm Exam Sprig, 018 Name: ID: Sigature: Istructio: This is a close book exam. You are allowed oe-page 8x11 formula sheet (-sided). No cellphoe or calculator or computer is allowed.

More information

Sampling Distributions and Estimation

Sampling Distributions and Estimation Cotets 40 Samplig Distributios ad Estimatio 40.1 Samplig Distributios 40. Iterval Estimatio for the Variace 13 Learig outcomes You will lear about the distributios which are created whe a populatio is

More information

Hopscotch and Explicit difference method for solving Black-Scholes PDE

Hopscotch and Explicit difference method for solving Black-Scholes PDE Mälardale iversity Fiacial Egieerig Program Aalytical Fiace Semiar Report Hopscotch ad Explicit differece method for solvig Blac-Scholes PDE Istructor: Ja Röma Team members: A Gog HaiLog Zhao Hog Cui 0

More information

CHAPTER 8: CONFIDENCE INTERVAL ESTIMATES for Means and Proportions

CHAPTER 8: CONFIDENCE INTERVAL ESTIMATES for Means and Proportions CHAPTER 8: CONFIDENCE INTERVAL ESTIMATES for Meas ad Proportios Itroductio: We wat to kow the value of a parameter for a populatio. We do t kow the value of this parameter for the etire populatio because

More information

43. A 000 par value 5-year bod with 8.0% semiaual coupos was bought to yield 7.5% covertible semiaually. Determie the amout of premium amortized i the 6 th coupo paymet. (A).00 (B).08 (C).5 (D).5 (E).34

More information

Math 312, Intro. to Real Analysis: Homework #4 Solutions

Math 312, Intro. to Real Analysis: Homework #4 Solutions Math 3, Itro. to Real Aalysis: Homework #4 Solutios Stephe G. Simpso Moday, March, 009 The assigmet cosists of Exercises 0.6, 0.8, 0.0,.,.3,.6,.0,.,. i the Ross textbook. Each problem couts 0 poits. 0.6.

More information

EVEN NUMBERED EXERCISES IN CHAPTER 4

EVEN NUMBERED EXERCISES IN CHAPTER 4 Joh Riley 7 July EVEN NUMBERED EXERCISES IN CHAPTER 4 SECTION 4 Exercise 4-: Cost Fuctio of a Cobb-Douglas firm What is the cost fuctio of a firm with a Cobb-Douglas productio fuctio? Rather tha miimie

More information

Calculation of the Annual Equivalent Rate (AER)

Calculation of the Annual Equivalent Rate (AER) Appedix to Code of Coduct for the Advertisig of Iterest Bearig Accouts. (31/1/0) Calculatio of the Aual Equivalet Rate (AER) a) The most geeral case of the calculatio is the rate of iterest which, if applied

More information

1 The Power of Compounding

1 The Power of Compounding 1 The Power of Compoudig 1.1 Simple vs Compoud Iterest You deposit $1,000 i a bak that pays 5% iterest each year. At the ed of the year you will have eared $50. The bak seds you a check for $50 dollars.

More information

Section Mathematical Induction and Section Strong Induction and Well-Ordering

Section Mathematical Induction and Section Strong Induction and Well-Ordering Sectio 4.1 - Mathematical Iductio ad Sectio 4. - Strog Iductio ad Well-Orderig A very special rule of iferece! Defiitio: A set S is well ordered if every subset has a least elemet. Note: [0, 1] is ot well

More information

Notes on Expected Revenue from Auctions

Notes on Expected Revenue from Auctions Notes o Epected Reveue from Auctios Professor Bergstrom These otes spell out some of the mathematical details about first ad secod price sealed bid auctios that were discussed i Thursday s lecture You

More information

FINM6900 Finance Theory How Is Asymmetric Information Reflected in Asset Prices?

FINM6900 Finance Theory How Is Asymmetric Information Reflected in Asset Prices? FINM6900 Fiace Theory How Is Asymmetric Iformatio Reflected i Asset Prices? February 3, 2012 Referece S. Grossma, O the Efficiecy of Competitive Stock Markets where Traders Have Diverse iformatio, Joural

More information

Chapter 4 - Consumer. Household Demand and Supply. Solving the max-utility problem. Working out consumer responses. The response function

Chapter 4 - Consumer. Household Demand and Supply. Solving the max-utility problem. Working out consumer responses. The response function Almost essetial Cosumer: Optimisatio Chapter 4 - Cosumer Osa 2: Household ad supply Cosumer: Welfare Useful, but optioal Firm: Optimisatio Household Demad ad Supply MICROECONOMICS Priciples ad Aalysis

More information

MA Lesson 11 Section 1.3. Solving Applied Problems with Linear Equations of one Variable

MA Lesson 11 Section 1.3. Solving Applied Problems with Linear Equations of one Variable MA 15200 Lesso 11 Sectio 1. I Solvig Applied Problems with Liear Equatios of oe Variable 1. After readig the problem, let a variable represet the ukow (or oe of the ukows). Represet ay other ukow usig

More information

Binomial Model. Stock Price Dynamics. The Key Idea Riskless Hedge

Binomial Model. Stock Price Dynamics. The Key Idea Riskless Hedge Biomial Model Stock Price Dyamics The value of a optio at maturity depeds o the price of the uderlyig stock at maturity. The value of the optio today depeds o the expected value of the optio at maturity

More information

Overlapping Generations

Overlapping Generations Eco. 53a all 996 C. Sims. troductio Overlappig Geeratios We wat to study how asset markets allow idividuals, motivated by the eed to provide icome for their retiremet years, to fiace capital accumulatio

More information

CHAPTER 8: CONFIDENCE INTERVAL ESTIMATES for Means and Proportions

CHAPTER 8: CONFIDENCE INTERVAL ESTIMATES for Means and Proportions CHAPTER 8: CONFIDENCE INTERVAL ESTIMATES for Meas ad Proportios Itroductio: I this chapter we wat to fid out the value of a parameter for a populatio. We do t kow the value of this parameter for the etire

More information

Subject CT5 Contingencies Core Technical. Syllabus. for the 2011 Examinations. The Faculty of Actuaries and Institute of Actuaries.

Subject CT5 Contingencies Core Technical. Syllabus. for the 2011 Examinations. The Faculty of Actuaries and Institute of Actuaries. Subject CT5 Cotigecies Core Techical Syllabus for the 2011 Examiatios 1 Jue 2010 The Faculty of Actuaries ad Istitute of Actuaries Aim The aim of the Cotigecies subject is to provide a groudig i the mathematical

More information

Monetary Economics: Problem Set #5 Solutions

Monetary Economics: Problem Set #5 Solutions Moetary Ecoomics oblem Set #5 Moetary Ecoomics: oblem Set #5 Solutios This problem set is marked out of 1 poits. The weight give to each part is idicated below. Please cotact me asap if you have ay questios.

More information

Unbiased estimators Estimators

Unbiased estimators Estimators 19 Ubiased estimators I Chapter 17 we saw that a dataset ca be modeled as a realizatio of a radom sample from a probability distributio ad that quatities of iterest correspod to features of the model distributio.

More information

18.S096 Problem Set 5 Fall 2013 Volatility Modeling Due Date: 10/29/2013

18.S096 Problem Set 5 Fall 2013 Volatility Modeling Due Date: 10/29/2013 18.S096 Problem Set 5 Fall 2013 Volatility Modelig Due Date: 10/29/2013 1. Sample Estimators of Diffusio Process Volatility ad Drift Let {X t } be the price of a fiacial security that follows a geometric

More information

Lecture 9: The law of large numbers and central limit theorem

Lecture 9: The law of large numbers and central limit theorem Lecture 9: The law of large umbers ad cetral limit theorem Theorem.4 Let X,X 2,... be idepedet radom variables with fiite expectatios. (i) (The SLLN). If there is a costat p [,2] such that E X i p i i=

More information

LESSON #66 - SEQUENCES COMMON CORE ALGEBRA II

LESSON #66 - SEQUENCES COMMON CORE ALGEBRA II LESSON #66 - SEQUENCES COMMON CORE ALGEBRA II I Commo Core Algebra I, you studied sequeces, which are ordered lists of umbers. Sequeces are extremely importat i mathematics, both theoretical ad applied.

More information

Monopoly vs. Competition in Light of Extraction Norms. Abstract

Monopoly vs. Competition in Light of Extraction Norms. Abstract Moopoly vs. Competitio i Light of Extractio Norms By Arkadi Koziashvili, Shmuel Nitza ad Yossef Tobol Abstract This ote demostrates that whether the market is competitive or moopolistic eed ot be the result

More information

Online appendices from The xva Challenge by Jon Gregory. APPENDIX 10A: Exposure and swaption analogy.

Online appendices from The xva Challenge by Jon Gregory. APPENDIX 10A: Exposure and swaption analogy. APPENDIX 10A: Exposure ad swaptio aalogy. Sorese ad Bollier (1994), effectively calculate the CVA of a swap positio ad show this ca be writte as: CVA swap = LGD V swaptio (t; t i, T) PD(t i 1, t i ). i=1

More information

DESCRIPTION OF MATHEMATICAL MODELS USED IN RATING ACTIVITIES

DESCRIPTION OF MATHEMATICAL MODELS USED IN RATING ACTIVITIES July 2014, Frakfurt am Mai. DESCRIPTION OF MATHEMATICAL MODELS USED IN RATING ACTIVITIES This documet outlies priciples ad key assumptios uderlyig the ratig models ad methodologies of Ratig-Agetur Expert

More information

Chapter 8: Estimation of Mean & Proportion. Introduction

Chapter 8: Estimation of Mean & Proportion. Introduction Chapter 8: Estimatio of Mea & Proportio 8.1 Estimatio, Poit Estimate, ad Iterval Estimate 8.2 Estimatio of a Populatio Mea: σ Kow 8.3 Estimatio of a Populatio Mea: σ Not Kow 8.4 Estimatio of a Populatio

More information

A New Approach to Obtain an Optimal Solution for the Assignment Problem

A New Approach to Obtain an Optimal Solution for the Assignment Problem Iteratioal Joural of Sciece ad Research (IJSR) ISSN (Olie): 231-7064 Idex Copericus Value (2013): 6.14 Impact Factor (2015): 6.31 A New Approach to Obtai a Optimal Solutio for the Assigmet Problem A. Seethalakshmy

More information

Online appendices from Counterparty Risk and Credit Value Adjustment a continuing challenge for global financial markets by Jon Gregory

Online appendices from Counterparty Risk and Credit Value Adjustment a continuing challenge for global financial markets by Jon Gregory Olie appedices from Couterparty Risk ad Credit Value Adjustmet a APPENDIX 8A: Formulas for EE, PFE ad EPE for a ormal distributio Cosider a ormal distributio with mea (expected future value) ad stadard

More information

Introduction to Probability and Statistics Chapter 7

Introduction to Probability and Statistics Chapter 7 Itroductio to Probability ad Statistics Chapter 7 Ammar M. Sarha, asarha@mathstat.dal.ca Departmet of Mathematics ad Statistics, Dalhousie Uiversity Fall Semester 008 Chapter 7 Statistical Itervals Based

More information

Subject CT1 Financial Mathematics Core Technical Syllabus

Subject CT1 Financial Mathematics Core Technical Syllabus Subject CT1 Fiacial Mathematics Core Techical Syllabus for the 2018 exams 1 Jue 2017 Subject CT1 Fiacial Mathematics Core Techical Aim The aim of the Fiacial Mathematics subject is to provide a groudig

More information

ACTUARIAL RESEARCH CLEARING HOUSE 1990 VOL. 2 INTEREST, AMORTIZATION AND SIMPLICITY. by Thomas M. Zavist, A.S.A.

ACTUARIAL RESEARCH CLEARING HOUSE 1990 VOL. 2 INTEREST, AMORTIZATION AND SIMPLICITY. by Thomas M. Zavist, A.S.A. ACTUARIAL RESEARCH CLEARING HOUSE 1990 VOL. INTEREST, AMORTIZATION AND SIMPLICITY by Thomas M. Zavist, A.S.A. 37 Iterest m Amortizatio ad Simplicity Cosider simple iterest for a momet. Suppose you have

More information

ON DIFFERENTIATION AND HARMONIC NUMBERS

ON DIFFERENTIATION AND HARMONIC NUMBERS ON DIFFERENTIATION AND HARMONIC NUMBERS ERIC MORTENSON Abstract. I a paper of Adrews ad Uchimura [AU, it is show how differetiatio applied to hypergeometric idetities produces formulas for harmoic ad q-harmoic

More information

Chapter 8. Confidence Interval Estimation. Copyright 2015, 2012, 2009 Pearson Education, Inc. Chapter 8, Slide 1

Chapter 8. Confidence Interval Estimation. Copyright 2015, 2012, 2009 Pearson Education, Inc. Chapter 8, Slide 1 Chapter 8 Cofidece Iterval Estimatio Copyright 2015, 2012, 2009 Pearso Educatio, Ic. Chapter 8, Slide 1 Learig Objectives I this chapter, you lear: To costruct ad iterpret cofidece iterval estimates for

More information

Research Article The Probability That a Measurement Falls within a Range of n Standard Deviations from an Estimate of the Mean

Research Article The Probability That a Measurement Falls within a Range of n Standard Deviations from an Estimate of the Mean Iteratioal Scholarly Research Network ISRN Applied Mathematics Volume 0, Article ID 70806, 8 pages doi:0.540/0/70806 Research Article The Probability That a Measuremet Falls withi a Rage of Stadard Deviatios

More information

ECON 5350 Class Notes Maximum Likelihood Estimation

ECON 5350 Class Notes Maximum Likelihood Estimation ECON 5350 Class Notes Maximum Likelihood Estimatio 1 Maximum Likelihood Estimatio Example #1. Cosider the radom sample {X 1 = 0.5, X 2 = 2.0, X 3 = 10.0, X 4 = 1.5, X 5 = 7.0} geerated from a expoetial

More information

Appendix 1 to Chapter 5

Appendix 1 to Chapter 5 Appedix 1 to Chapter 5 Models of Asset Pricig I Chapter 4, we saw that the retur o a asset (such as a bod) measures how much we gai from holdig that asset. Whe we make a decisio to buy a asset, we are

More information

FOUNDATION ACTED COURSE (FAC)

FOUNDATION ACTED COURSE (FAC) FOUNDATION ACTED COURSE (FAC) What is the Foudatio ActEd Course (FAC)? FAC is desiged to help studets improve their mathematical skills i preparatio for the Core Techical subjects. It is a referece documet

More information

Twitter: @Owe134866 www.mathsfreeresourcelibrary.com Prior Kowledge Check 1) State whether each variable is qualitative or quatitative: a) Car colour Qualitative b) Miles travelled by a cyclist c) Favourite

More information

of Asset Pricing R e = expected return

of Asset Pricing R e = expected return Appedix 1 to Chapter 5 Models of Asset Pricig EXPECTED RETURN I Chapter 4, we saw that the retur o a asset (such as a bod) measures how much we gai from holdig that asset. Whe we make a decisio to buy

More information

Proceedings of the 5th WSEAS Int. Conf. on SIMULATION, MODELING AND OPTIMIZATION, Corfu, Greece, August 17-19, 2005 (pp )

Proceedings of the 5th WSEAS Int. Conf. on SIMULATION, MODELING AND OPTIMIZATION, Corfu, Greece, August 17-19, 2005 (pp ) Proceedigs of the 5th WSEAS It. Cof. o SIMULATION, MODELING AND OPTIMIZATION, Corfu, Greece, August 7-9, 005 (pp488-49 Realized volatility estimatio: ew simulatio approach ad empirical study results JULIA

More information

AMS Portfolio Theory and Capital Markets

AMS Portfolio Theory and Capital Markets AMS 69.0 - Portfolio Theory ad Capital Markets I Class 6 - Asset yamics Robert J. Frey Research Professor Stoy Brook iversity, Applied Mathematics ad Statistics frey@ams.suysb.edu http://www.ams.suysb.edu/~frey/

More information

0.1 Valuation Formula:

0.1 Valuation Formula: 0. Valuatio Formula: 0.. Case of Geeral Trees: q = er S S S 3 S q = er S S 4 S 5 S 4 q 3 = er S 3 S 6 S 7 S 6 Therefore, f (3) = e r [q 3 f (7) + ( q 3 ) f (6)] f () = e r [q f (5) + ( q ) f (4)] = f ()

More information

of Asset Pricing APPENDIX 1 TO CHAPTER EXPECTED RETURN APPLICATION Expected Return

of Asset Pricing APPENDIX 1 TO CHAPTER EXPECTED RETURN APPLICATION Expected Return APPENDIX 1 TO CHAPTER 5 Models of Asset Pricig I Chapter 4, we saw that the retur o a asset (such as a bod) measures how much we gai from holdig that asset. Whe we make a decisio to buy a asset, we are

More information

Math 124: Lecture for Week 10 of 17

Math 124: Lecture for Week 10 of 17 What we will do toight 1 Lecture for of 17 David Meredith Departmet of Mathematics Sa Fracisco State Uiversity 2 3 4 April 8, 2008 5 6 II Take the midterm. At the ed aswer the followig questio: To be revealed

More information

Anomaly Correction by Optimal Trading Frequency

Anomaly Correction by Optimal Trading Frequency Aomaly Correctio by Optimal Tradig Frequecy Yiqiao Yi Columbia Uiversity September 9, 206 Abstract Uder the assumptio that security prices follow radom walk, we look at price versus differet movig averages.

More information

Neighboring Optimal Solution for Fuzzy Travelling Salesman Problem

Neighboring Optimal Solution for Fuzzy Travelling Salesman Problem Iteratioal Joural of Egieerig Research ad Geeral Sciece Volume 2, Issue 4, Jue-July, 2014 Neighborig Optimal Solutio for Fuzzy Travellig Salesma Problem D. Stephe Digar 1, K. Thiripura Sudari 2 1 Research

More information

CHAPTER 2 PRICING OF BONDS

CHAPTER 2 PRICING OF BONDS CHAPTER 2 PRICING OF BONDS CHAPTER SUARY This chapter will focus o the time value of moey ad how to calculate the price of a bod. Whe pricig a bod it is ecessary to estimate the expected cash flows ad

More information

Forecasting bad debt losses using clustering algorithms and Markov chains

Forecasting bad debt losses using clustering algorithms and Markov chains Forecastig bad debt losses usig clusterig algorithms ad Markov chais Robert J. Till Experia Ltd Lambert House Talbot Street Nottigham NG1 5HF {Robert.Till@uk.experia.com} Abstract Beig able to make accurate

More information

ii. Interval estimation:

ii. Interval estimation: 1 Types of estimatio: i. Poit estimatio: Example (1) Cosider the sample observatios 17,3,5,1,18,6,16,10 X 8 X i i1 8 17 3 5 118 6 16 10 8 116 8 14.5 14.5 is a poit estimate for usig the estimator X ad

More information

Problem Set 1a - Oligopoly

Problem Set 1a - Oligopoly Advaced Idustrial Ecoomics Sprig 2014 Joha Steek 6 may 2014 Problem Set 1a - Oligopoly 1 Table of Cotets 2 Price Competitio... 3 2.1 Courot Oligopoly with Homogeous Goods ad Differet Costs... 3 2.2 Bertrad

More information

Solutions to Problem Sheet 1

Solutions to Problem Sheet 1 Solutios to Problem Sheet ) Use Theorem.4 to prove that p log for all real x 3. This is a versio of Theorem.4 with the iteger N replaced by the real x. Hit Give x 3 let N = [x], the largest iteger x. The,

More information

Institute of Actuaries of India Subject CT5 General Insurance, Life and Health Contingencies

Institute of Actuaries of India Subject CT5 General Insurance, Life and Health Contingencies Istitute of Actuaries of Idia Subject CT5 Geeral Isurace, Life ad Health Cotigecies For 2017 Examiatios Aim The aim of the Cotigecies subject is to provide a groudig i the mathematical techiques which

More information

CAPITAL PROJECT SCREENING AND SELECTION

CAPITAL PROJECT SCREENING AND SELECTION CAPITAL PROJECT SCREEIG AD SELECTIO Before studyig the three measures of ivestmet attractiveess, we will review a simple method that is commoly used to scree capital ivestmets. Oe of the primary cocers

More information

A New Second-Order Corrector Interior-Point Algorithm for P (κ)-lcp

A New Second-Order Corrector Interior-Point Algorithm for P (κ)-lcp Filomat 3:0 07, 6379 639 https://doi.org/0.9/fil70379k Published by Faculty of Scieces ad Mathematics, Uiversity of Niš, Serbia Available at: http://www.pmf.i.ac.rs/filomat A New Secod-Order Corrector

More information

. (The calculated sample mean is symbolized by x.)

. (The calculated sample mean is symbolized by x.) Stat 40, sectio 5.4 The Cetral Limit Theorem otes by Tim Pilachowski If you have t doe it yet, go to the Stat 40 page ad dowload the hadout 5.4 supplemet Cetral Limit Theorem. The homework (both practice

More information

Mixed and Implicit Schemes Implicit Schemes. Exercise: Verify that ρ is unimodular: ρ = 1.

Mixed and Implicit Schemes Implicit Schemes. Exercise: Verify that ρ is unimodular: ρ = 1. Mixed ad Implicit Schemes 3..4 The leapfrog scheme is stable for the oscillatio equatio ad ustable for the frictio equatio. The Euler forward scheme is stable for the frictio equatio but ustable for the

More information

A Bayesian perspective on estimating mean, variance, and standard-deviation from data

A Bayesian perspective on estimating mean, variance, and standard-deviation from data Brigham Youg Uiversity BYU ScholarsArchive All Faculty Publicatios 006--05 A Bayesia perspective o estimatig mea, variace, ad stadard-deviatio from data Travis E. Oliphat Follow this ad additioal works

More information

Topic-7. Large Sample Estimation

Topic-7. Large Sample Estimation Topic-7 Large Sample Estimatio TYPES OF INFERENCE Ò Estimatio: É Estimatig or predictig the value of the parameter É What is (are) the most likely values of m or p? Ò Hypothesis Testig: É Decidig about

More information

NORMALIZATION OF BEURLING GENERALIZED PRIMES WITH RIEMANN HYPOTHESIS

NORMALIZATION OF BEURLING GENERALIZED PRIMES WITH RIEMANN HYPOTHESIS Aales Uiv. Sci. Budapest., Sect. Comp. 39 2013) 459 469 NORMALIZATION OF BEURLING GENERALIZED PRIMES WITH RIEMANN HYPOTHESIS We-Bi Zhag Chug Ma Pig) Guagzhou, People s Republic of Chia) Dedicated to Professor

More information

r i = a i + b i f b i = Cov[r i, f] The only parameters to be estimated for this model are a i 's, b i 's, σe 2 i

r i = a i + b i f b i = Cov[r i, f] The only parameters to be estimated for this model are a i 's, b i 's, σe 2 i The iformatio required by the mea-variace approach is substatial whe the umber of assets is large; there are mea values, variaces, ad )/2 covariaces - a total of 2 + )/2 parameters. Sigle-factor model:

More information

SUPPLEMENTAL MATERIAL

SUPPLEMENTAL MATERIAL A SULEMENTAL MATERIAL Theorem (Expert pseudo-regret upper boud. Let us cosider a istace of the I-SG problem ad apply the FL algorithm, where each possible profile A is a expert ad receives, at roud, a

More information

Inferential Statistics and Probability a Holistic Approach. Inference Process. Inference Process. Chapter 8 Slides. Maurice Geraghty,

Inferential Statistics and Probability a Holistic Approach. Inference Process. Inference Process. Chapter 8 Slides. Maurice Geraghty, Iferetial Statistics ad Probability a Holistic Approach Chapter 8 Poit Estimatio ad Cofidece Itervals This Course Material by Maurice Geraghty is licesed uder a Creative Commos Attributio-ShareAlike 4.0

More information

MODIFICATION OF HOLT S MODEL EXEMPLIFIED BY THE TRANSPORT OF GOODS BY INLAND WATERWAYS TRANSPORT

MODIFICATION OF HOLT S MODEL EXEMPLIFIED BY THE TRANSPORT OF GOODS BY INLAND WATERWAYS TRANSPORT The publicatio appeared i Szoste R.: Modificatio of Holt s model exemplified by the trasport of goods by ilad waterways trasport, Publishig House of Rzeszow Uiversity of Techology No. 85, Maagemet ad Maretig

More information

AUTOMATIC GENERATION OF FUZZY PAYOFF MATRIX IN GAME THEORY

AUTOMATIC GENERATION OF FUZZY PAYOFF MATRIX IN GAME THEORY AUTOMATIC GENERATION OF FUZZY PAYOFF MATRIX IN GAME THEORY Dr. Farha I. D. Al Ai * ad Dr. Muhaed Alfarras ** * College of Egieerig ** College of Coputer Egieerig ad scieces Gulf Uiversity * Dr.farha@gulfuiversity.et;

More information

14.30 Introduction to Statistical Methods in Economics Spring 2009

14.30 Introduction to Statistical Methods in Economics Spring 2009 MIT OpeCourseWare http://ocwmitedu 430 Itroductio to Statistical Methods i Ecoomics Sprig 009 For iformatio about citig these materials or our Terms of Use, visit: http://ocwmitedu/terms 430 Itroductio

More information

B = A x z

B = A x z 114 Block 3 Erdeky == Begi 6.3 ============================================================== 1 / 8 / 2008 1 Correspodig Areas uder a ormal curve ad the stadard ormal curve are equal. Below: Area B = Area

More information

Parametric Density Estimation: Maximum Likelihood Estimation

Parametric Density Estimation: Maximum Likelihood Estimation Parametric Desity stimatio: Maimum Likelihood stimatio C6 Today Itroductio to desity estimatio Maimum Likelihood stimatio Itroducto Bayesia Decisio Theory i previous lectures tells us how to desig a optimal

More information