Surface Curvatures. (Com S 477/577 Notes) Yan-Bin Jia. Oct 23, Curve on a Surface: Normal and Geodesic Curvatures
|
|
- Valerie Horn
- 6 years ago
- Views:
Transcription
1 Surface Curvatures (Com S 477/577 Notes) Yan-Bin Jia Oct 23, Curve on a Surface: Normal and Geodesic Curvatures One way to examine how much a surface bends is to look at the curvature of curves on the surface. Let γ(t) = σ(u(t),v(t)) be a unit-speed curve in a surface patch σ. Thus, γ is a unit tangent vector to σ, and it is perpendicular to the surface normal ˆn at the same point. The three vectors γ, ˆn γ, and ˆn form a local coordinate frame by the right-hand rule. Differentiating γ γ yields that γ is orthogonal to γ. Hence γ is a linear combination of ˆn and ˆn γ: γ = κ nˆn+κ gˆn γ (1) Here κ n is called the normal curvature and κ g is the geodesic curvature of γ. ˆn ˆn γ γ φ γ κ n γ κ g σ Since ˆn and ˆn γ are orthogonal to each other, (1) implies that Since γ is unit-speed, its curvature is κ n = γ ˆn and κ g = γ (ˆn γ). κ = γ = κ nˆn+κ gˆn γ = κ 2 n +κ2 g. (2) The material is adapted from the book Elementary Differential Geometry by Andrew Pressley, Springer-Verlag, 1
2 Let ψ be the angle between the principal normal ˆn γ of the curve, in the direction of γ, and the surface normal ˆn. We have κ n = κˆn γ ˆn = κcosψ. (3) Therefore, equation (2) implies κ g = ±κsinψ. If γ is regular but arbitrary-speed, the normal and geodesic curvatures of γ are defined to be those of a unit-speed reparametrization of the same curve. A unit-speed curve γ is a geodesic if κ g = 0. By (1), its acceleration γ is always normal to the surface. Geodesics have many applications that we will devote one lecture to the topic later on. 2 Darboux Frame on a Curve On the unit-speed surface curve γ, the frame formed by the unit vectors γ, ˆn γ, and ˆn is called the Darboux frame on the curve. This frame is different from the Frenet frame on the curve defined by γ, the principal normal γ/ γ, and the binormal γ γ/ γ. Let us rename the three unit vectors γ, ˆn γ, ˆn as T, V, U, respectively. Their derivatives must be respectively orthogonal to themselves. Differentiation of U T = 0 yields U T = U T = U (κ n U +κ g V) (by (1)) = κ n. Thevector U has a second component along V. Thegeodesic torsion, definedto beτ g = U V, describes the negative rate of change of U in the direction of V. With τ g we characterize this change rate completely as below: U = κ g T τ g V. Similarly, we express V in terms of T,U,V, and combine it with the above equation and (1) into the following compact form (where the vectors are viewed as scalars ): T V U 0 κ g κ n = κ g 0 τ g κ n τ g 0 T V U (4) The formulas (4) describe the geometry of a curve at a point in a local frame suiting the curve as well as a surface on which it lies, whereas the Frenet formulas describe its geometry at the point in a local frame best suiting the curve alone. The curve γ is asymptotic provided its tangent γ always points in a direction in which κ n = 0. In some sense, the surface bends less along γ than it does along a general curve. In (4), T = κ g V. Thus, κ = κ g and V is aligned with the principal normal T / T (assuming κ g 0). Consequently, the Darboux frame T-V-U coincides with the Frenet frame everywhere. 2
3 3 The Second Fundamental Form Let σ be a surface patch in R 3 with standard unit normal ˆn = σ u σ v σ u σ v = σ u σ v EG F 2. (5) With an increase ( u, v) in the parameter values, the movement of the point is described by Taylor s series below: σ(u+ u,v + v) σ(u,v) = σ u u+σ v v + 1 ( ( σ uu ( u) 2 +2σ uv u v +σ vv ( v) 2 )+O ( u+ v) 3). 2 The first order terms are tangent to the surface, hence perpendicular to ˆn. The terms of order higher than two tend to zero as ( u) 2 +( v) 2 does so. The deviation of σ from the tangent plane is determined by the dot product of the second order term with the surface normal ˆn, namely, it is where 1 ( L( u) 2 +2M u v +N( v) 2), (6) 2 L = σ uu ˆn = σ uu (σ u σ v ) σ u σ v = det(σ uu σ u σ v ), EG F 2 by (5) M = σ uv ˆn = det(σ uv σ u σ v ), (7) EG F 2 N = σ vv ˆn = det(σ vv σ u σ v ). EG F 2 Recall that a unit-speed space curve α(s) can be approximated around s = 0 up to the second order as α(0) + sˆt(0) κ(0)s2ˆn(0), where ˆt(0) and ˆn(0) are the tangent and principal normal, respectively, and κ(0) the curvature. The expression(6) for the surface is analogous to the curvature term 1 2 κ(0)s2ˆn(0) for a curve. In particular, the expression Ldu 2 +2Mdudv +Ndv 2 is the second fundamental form of σ. While the first fundamental form permits the calculation of metric properties such as length and area on a surface patch, the second fundamental form captures how curved a surface patch is. The roles of the two fundamental forms on describing the local geometry are analogous to those of speed and acceleration for a parametric curve. Just as a unit-speed space curve is determined up to a rigid motion by its curvature and torsion, a surface patch is determined up to a rigid motion by its first and second fundamental forms. 3
4 Example 1. Consider a surface of revolution σ(u,v) = (f(u)cosv,f(u)sinv,g(u)), where f(u) > 0 always holds. The figure below plots a catenoid where ( u f(u) = 2cosh and g(u) = u 2) with (u,v) [ 2,2] [0,2π]. Assume the profile curve u (f(u), 0, g(u)) is unit-speed, i.e., f 2 + ġ 2 = 1, where a dot denotes differentiation with respect to u. We perform the following calculations: σ u = ( f cosv, f sinv,ġ), σ v = ( f sinv,f cosv,0), E = σ u σ u = f 2 +ġ 2 = 1, F = σ u σ v = 0, G = σ v σ v = f 2, σ u σ v = ( fġcosv, fġsinv,ff), σ u σ v = f, σ u σ v ˆn = σ u σ v = ( ġcosv, ġsinv, f), σ uu = ( f cosv, f sinv, g), = ( f sinv, f cosv,0), σ uv σ vv = ( f cosv, f sinv,0), L = σ uu ˆn = f g fġ, M = σ uv ˆn = 0, N = σ vv ˆn = fġ. Thus, the second fundamental form is ( f g fġ)du 2 +fġdv 2. 4
5 4 Principal Curvatures Letγ(t) = σ(u(t),v(t)) beaunit-speedcurveonasurfacepatchσ withthestandardunitnormal ˆn. Below we obtain the normal curvature of the curve: κ n = γ ˆn = d γ ˆn dt = ˆn d dt (σ u u+σ v v) ( ) = ˆn σ u ü+σ v v +(σ uu u+σ uv v) u+(σ uv u+σ vv v) v = L u 2 +2M u v +N v 2, (8) where L,M,N are coefficients of the second fundamental form defined (7). In the last step above, we used the fact that ˆn is normal to both σ u and σ v, i.e., ˆn σ u = ˆn σ v = 0. Equation (8) states that the normal curvature of γ(t) depends on u, v, u, and v. The first two quantities specify the location of the point on the surface σ, and are thus curve independent. Let û be the unit tangent vector so that γ = û. Take the dot products of the equation û = σ u u+σ v v with σ u and σ v separately, yielding E u+f v = û σ u, F u+g v = û σ v, where E,F,G are the coefficients of the first fundamental form of the surface patch. Since σ u and σ v are linearly independent, the coefficient matrix in the above linear system in u and v is non-singular. We solve the system to obtain ) ( ( u E F = v F G ) 1 ) (û σu û σ v This implies that u and v are independent of the parametrization of γ. By (8), we conclude that any two unit-speed curves passing through the same point in the same direction û must have the same normal curvature at this point. Subsequently, we refer to κ n as the normal curvature of the surfaceσ atthepointp = σ(u,v)inthetangentdirectionofû. It measures the curving of the surface in that direction. Generally, the surface bends at different rates in different tangent directions. ˆn The tangent û and the normal ˆn at p defines a plane that cuts a û curve α out of the patch. This curve is called the normal section of σ in the û direction. Since the principal normal ˆn γ of the α normal section is related to the surface normal ˆn by ˆn γ = ±ˆn. Equation (3) implies that the curvature of the normal section is σ the normal curvature κ n at the point or its opposite, depending on the choice of the surface normal ˆn at the point. 5
6 To analyze the normal curvature κ n further, we make use of the first and second fundamental forms: Edu 2 +2Fdudv +Gdv 2 and Ldu 2 +2Mdudv +Ndv 2. For convenience, we introduce two symmetric matrices ( ) E F F 1 = and F F G 2 = ( L M M N ) The tangent vector of the unit-speed curve γ(t) = σ(u(t),v(t)) is γ = uσ u + vσ v. Introducing T = ( u, v) t, where t denotes the transpose operator, equation (8) can be rewritten as κ n = T t F 2 T. (9) Since σ u and σ v span the tangent plane, consider two tangent vectors: We obtain their inner product: where ˆt 1 = ξ 1 σ u +η 1 σ v and ˆt 2 = ξ 2 σ u +η 2 σ v. ˆt 1 ˆt 2 = (ξ 1 σ u +η 1 σ v ) (ξ 2 σ u +η 2 σ v ) T 1 = = Eξ 1 ξ 2 +F(ξ 1 η 2 +ξ 2 η 1 )+Gη 1 η 2 ( ) = T1 t E F T F G 2, (10) ( ξ1 η 1 ) and T 2 = The principal curvatures of the surface patch σ are the roots of the equation det(f 2 κf 1 ) = L κe M κf M κf N κg = 0. (11) From the linear independence of σ u and σ v, it is easy to show that matrix F 1 is always invertible. Equation (11) essentially states that the principal curvatures are the eigenvalues of F1 1 F 2. Let κ be a principal curvature of F1 1 F 2, and T = (ξ,η) T the corresponding eigenvector. That (F1 1 F 2)T = κt implies (F 2 κf 1 )T = 0. (12) The unit tangent vector ˆt in the direction of ξσ u +ησ v is called the principal vector corresponding to the principal curvature κ. Theorem 1 Let κ 1 and κ 2 be the principal curvatures at a point p of a surface patch σ. Then (i) κ 1,κ 2 R; ( ξ2 (ii) if κ 1 = κ 2 = κ, then F 2 = κf 1 and every tangent vector at p is a principal vector. (iii) if κ 1 κ 2, then the two corresponding principal vectors are perpendicular to each other η 2 ). 6
7 For a proof of the theorem, we refer the reader to [2, ]. Intuitively, the principal vectors give the directions of maximum and minimum bending of the surface at the point, and the principal curvatures measure the bending rates. In case (ii), the point is umbilic. In this case, the surface bends the same amount in all directions at p (thus all directions are principal). Example 2. A sphere bends the same amount in every direction. Take the unit sphere in Example 9 in the notes Surfaces, for instance, with the parametrization We found previously that σ(θ,φ) = (cosθcosφ,cosθsinφ,sinθ), E = 1, F = 0, G = cos 2 θ. Since a sphere is a surface of revolution, we can plug in the result from Example 1, with f(θ) = cosθ and g(u) = sinθ: L = f g fġ = sinθ( sinθ) ( cosθ)cosθ = 1, M = 0, N = fġ = cos 2 θ. Hence the principal curvatures are the roots of det(f 2 kf 1 ) = 1 κ 0 0 cos 2 θ κcos 2 θ = 0. Hence κ = 1. And every tangent direction is a principal vector. Example 3. Consider a cylinder with the z-axis as its axis and circular cross sections of unit radius. The parametrization is given as σ(u,v) = (cosv,sinv,u). ˆn ˆt 1 ˆt 2 p The coefficients of the first and second fundamental forms can be computed as E = 1, F = 0, G = 1, L = 0, M = 0, N = 1. The principal curvatures are roots of 0 κ κ 7
8 So we obtain κ 1 = 0 and κ 2 = 1. The eigenvectors T i = (ξ i,η i ), i = 1,2 of F 1 1 F 2 are found from solving the equation (F 2 κ i F 1 )T i = 0. The results are T 1 = λ 1 (1,0) t and T 2 = λ 2 (0,1) t for any non-zero λ 1,λ 2 R. Hence the principal vector ˆt 1 is along the direction of 1σ u +0σ v, i.e., ˆt 1 = (0,0,1). The principal vector ˆt 2 is along the direction of 0σ u +1σ v = ( sinv,cosv,0), i.e., ˆt 2 = ( sinv,cosv,0). A curve γ on the surface σ is a principal curve if its velocity γ always points in a principal direction, that is, the direction of a principal vector. At every point on a principal curve, the normal curvature is a maximum or minimum. The next figure shows some principal curves on the ellipsoid x y2 5 +z2 = 1. 5 Euler s Formula Suppose the two principal curvatures κ 1 κ 2 at p on the surface σ. Then by Theorem 1(iii), the two corresponding principal vectors ˆt 1 = ξ 1 σ u + η 1 σ v and ˆt 2 = ξ 2 σ u + η 2 σ v must be orthogonal to each other. Denote by T 1 = (ξ 1,η 1 ) t and T 2 = (ξ 2,η 2 ) t. Replace the T in (12) with T j, j = 1,2, multiply both sides of the equation by Ti t to the left, and move the second resulting term to the right hand side of the equation. This yields T t i F 2T j = κ j T t i F 1T j, i,j = 1,2. Meanwhile, the orthogonality of the two principal vectors implies that ˆt 1 ˆt 2 = T t 1F 1 T 2 = 0, from (10). 8
9 To summarize, we have { TiF t 2 T j = κ i TiF t κi, if i = j, 1 T j = 0, otherwise (13) With the principal curvatures and vectors at p, we can evaluate the normal curvature in any direction. Theorem 2 Let κ 1,κ 2 be the principal curvatures, and ˆt 1,ˆt 2 the two corresponding principal vectors of a patch σ at p. The normal curvature of σ in the direction û = cosθˆt 1 +sinθˆt 2 is κ n = κ 1 cos 2 θ+κ 2 sin 2 θ. Proof Let û = ξσ u + ησ v and T = (ξ,η) t. We first look at the special case κ 1 = κ 2 = κ. By Theorem 1(ii), û = ξσ u +ησ v is a principal vector. The normal curvature in the direction û is κ n = T t F 2 T (by (9)) = κt t F 1 T (by (12)) = κû û = κ. (14) Meanwhile, we have κ 1 cos 2 θ+κ 2 sin 2 θ = κ(cos 2 θ+sin 2 θ) = κ. So the theorem holds when the point is umbilic. Assume κ 1 κ 2. Therefore by Theorem 1(iii), ˆt 1 and ˆt 2 are perpendicular to each other. Let Thus, So we have û = ξσ u +ησ v, where ˆt i = ξ i σ u +η i σ v, and T i = (ξ i,η i ) t. û = cosθ(ξ 1 σ u +η 1 σ v )+sinθ(ξ 2 σ u +η 2 σ v ) = (ξ 1 cosθ +ξ 2 sinθ)σ u +(η 1 cosθ+η 2 sinθ)σ v. ξ = ξ 1 cosθ+ξ 2 sinθ, η = η 1 cosθ+η 2 sinθ, The above is written succinctly as T = cosθt 1 +sinθt 2. By equation (9) the normal curvature in the û direction is κ n = (cosθt t 1 +sinθtt 2 )F 2(cosθT 1 +sinθt 2 ) = cos 2 θt t 1F 2 T 1 +cosθsinθ(t t 1F 2 T 2 +T t 2F 2 T 1 )+sin 2 θt t 2F 2 T 2 = κ 1 cos 2 θ+κ 2 sin 2 θ. The last step above followed from the equation (13). Theorem 2 implies that κ 1 and κ 2 are the maximum and minimum of any normal curvatures at the point. Equivalently, among all tangent directions at the point, the geometry varies the most in one principal direction while the least in the other. 9
10 References [1] B. O Neill. Elementary Differential Geometry. Academic Press, Inc., [2] A. Pressley. Elementary Differential Geometry. Springer-Verlag London,
Curves on a Surface. (Com S 477/577 Notes) Yan-Bin Jia. Oct 24, 2017
Curves on a Surface (Com S 477/577 Notes) Yan-Bin Jia Oct 24, 2017 1 Normal and Geodesic Curvatures One way to examine how much a surface bends is to look at the curvature of curves on the surface. Let
More informationPrincipal Curvatures
Principal Curvatures Com S 477/577 Notes Yan-Bin Jia Oct 26, 2017 1 Definition To furtheranalyze thenormal curvatureκ n, we make useof the firstandsecond fundamental forms: Edu 2 +2Fdudv +Gdv 2 and Ldu
More informationRiemannian Geometry, Key to Homework #1
Riemannian Geometry Key to Homework # Let σu v sin u cos v sin u sin v cos u < u < π < v < π be a parametrization of the unit sphere S {x y z R 3 x + y + z } Fix an angle < θ < π and consider the parallel
More information5.1 Gauss Remarkable Theorem
5.1 Gauss Remarkable Theorem Recall that, for a surface M, its Gauss curvature is defined by K = κ 1 κ 2 where κ 1, κ 2 are the principal curvatures of M. The principal curvatures are the eigenvalues of
More information1.17 The Frenet-Serret Frame and Torsion. N(t) := T (t) κ(t).
Math 497C Oct 1, 2004 1 Curves and Surfaces Fall 2004, PSU Lecture Notes 7 1.17 The Frenet-Serret Frame and Torsion Recall that if α: I R n is a unit speed curve, then the unit tangent vector is defined
More informationSmarandache Curves on S 2. Key Words: Smarandache Curves, Sabban Frame, Geodesic Curvature. Contents. 1 Introduction Preliminaries 52
Bol. Soc. Paran. Mat. (s.) v. (04): 5 59. c SPM ISSN-75-88 on line ISSN-00787 in press SPM: www.spm.uem.br/bspm doi:0.569/bspm.vi.94 Smarandache Curves on S Kemal Taşköprü and Murat Tosun abstract: In
More informationk-type null slant helices in Minkowski space-time
MATHEMATICAL COMMUNICATIONS 8 Math. Commun. 20(2015), 8 95 k-type null slant helices in Minkowski space-time Emilija Nešović 1, Esra Betül Koç Öztürk2, and Ufuk Öztürk2 1 Department of Mathematics and
More informationHomework JWR. Feb 6, 2014
Homework JWR Feb 6, 2014 1. Exercise 1.5-12. Let the position of a particle at time t be given by α(t) = β(σ(t)) where β is parameterized by arclength and σ(t) = α(t) is the speed of the particle. Denote
More informationThe Smarandache Curves on H 0
Gazi University Journal of Science GU J Sci 9():69-77 (6) The Smarandache Curves on H Murat SAVAS, Atakan Tugkan YAKUT,, Tugba TAMIRCI Gazi University, Faculty of Sciences, Department of Mathematics, 65
More informationA Characterization for Bishop Equations of Parallel Curves according to Bishop Frame in E 3
Bol. Soc. Paran. Mat. (3s.) v. 33 1 (2015): 33 39. c SPM ISSN-2175-1188 on line ISSN-00378712 in press SPM: www.spm.uem.br/bspm doi:10.5269/bspm.v33i1.21712 A Characterization for Bishop Equations of Parallel
More informationEuler Savary s Formula On Complex Plane C
Applied Mathematics E-Notes, 606, 65-7 c ISSN 607-50 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Euler Savary s Formula On Complex Plane C Mücahit Akbıyık, Salim Yüce Received
More informationDENSITY OF PERIODIC GEODESICS IN THE UNIT TANGENT BUNDLE OF A COMPACT HYPERBOLIC SURFACE
DENSITY OF PERIODIC GEODESICS IN THE UNIT TANGENT BUNDLE OF A COMPACT HYPERBOLIC SURFACE Marcos Salvai FaMAF, Ciudad Universitaria, 5000 Córdoba, Argentina. e-mail: salvai@mate.uncor.edu Abstract Let S
More informationCURVATURE AND TORSION FORMULAS FOR CONFLICT SETS
GEOMETRY AND TOPOLOGY OF CAUSTICS CAUSTICS 0 BANACH CENTER PUBLICATIONS, VOLUME 6 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 004 CURVATURE AND TORSION FORMULAS FOR CONFLICT SETS MARTIJN
More informationON GENERALIZED NULL MANNHEIM CURVES IN MINKOWSKI SPACE-TIME. Milica Grbović, Kazim Ilarslan, and Emilija Nešović
PUBLICATIONS DE L INSTITUT MATHÉMATIQUE Nouvelle série, tome 99(113 (016, 77 98 DOI: 10.98/PIM1613077G ON GENERALIZED NULL MANNHEIM CURVES IN MINKOWSKI SPACE-TIME Milica Grbović, Kazim Ilarslan, and Emilija
More information1-TYPE CURVES AND BIHARMONIC CURVES IN EUCLIDEAN 3-SPACE
International Electronic Journal of Geometry Volume 4 No. 1 pp. 97-101 (2011) c IEJG 1-TYPE CURVES AND BIHARMONIC CURVES IN EUCLIDEAN 3-SPACE (Communicated by Shyuichi Izumiya) Abstract. We study 1-type
More informationDifferential Geometry: Curvature, Maps, and Pizza
Differential Geometry: Curvature, Maps, and Pizza Madelyne Ventura University of Maryland December 8th, 2015 Madelyne Ventura (University of Maryland) Curvature, Maps, and Pizza December 8th, 2015 1 /
More informationResearch Article The Smarandache Curves on S 2 1 and Its Duality on H2 0
Hindawi Publishing Corporation Journal of Applied Mathematics Volume 04, Article ID 93586, pages http://dx.doi.org/0.55/04/93586 Research Article The Smarandache Curves on S and Its Duality on H 0 Atakan
More informationTangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.
Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey
More informationMath F412: Homework 3 Solutions February 14, Solution: By the Fundamental Theorem of Calculus applied to the coordinates of α we have
1. Let k() be a mooth function on R. Let α() = ( θ() = co(θ(u)) du, k(u) du in(θ(u)) du). Show that α i a mooth unit peed curve with igned curvature κ p () = k(). By the Fundamental Theorem of Calculu
More informationarxiv: v1 [math.dg] 31 Mar 2014 Generalized Similar Frenet Curves
arxiv:14037908v1 [mathdg] 31 Mar 2014 Generalize Similar Frenet Curves Fatma GÖKÇELİK, Seher KAYA, Yusuf YAYLI, an F Nejat EKMEKCİ Abstract The paper is evote to ifferential geometric invariants etermining
More information4: SINGLE-PERIOD MARKET MODELS
4: SINGLE-PERIOD MARKET MODELS Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 4: Single-Period Market Models 1 / 87 General Single-Period
More informationPart 3: Trust-region methods for unconstrained optimization. Nick Gould (RAL)
Part 3: Trust-region methods for unconstrained optimization Nick Gould (RAL) minimize x IR n f(x) MSc course on nonlinear optimization UNCONSTRAINED MINIMIZATION minimize x IR n f(x) where the objective
More informationIMPA Commodities Course : Forward Price Models
IMPA Commodities Course : Forward Price Models Sebastian Jaimungal sebastian.jaimungal@utoronto.ca Department of Statistics and Mathematical Finance Program, University of Toronto, Toronto, Canada http://www.utstat.utoronto.ca/sjaimung
More informationSaddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance
Saddlepoint Approximation Methods for Pricing Financial Options on Discrete Realized Variance Yue Kuen KWOK Department of Mathematics Hong Kong University of Science and Technology Hong Kong * This is
More informationSquare-Root Measurement for Ternary Coherent State Signal
ISSN 86-657 Square-Root Measurement for Ternary Coherent State Signal Kentaro Kato Quantum ICT Research Institute, Tamagawa University 6-- Tamagawa-gakuen, Machida, Tokyo 9-86, Japan Tamagawa University
More informationSNELL S LAW AND UNIFORM REFRACTION. Contents
SNELL S LAW AND UNIFORM REFRACTION CRISTIAN E. GUTIÉRREZ Contents 1. Snell s law of refraction 1 1.1. In vector form 1 1.2. κ < 1 2 1.3. κ > 1 3 1.4. κ = 1 4 2. Uniform refraction 4 2.1. Surfaces with
More informationsymmys.com 3.2 Projection of the invariants to the investment horizon
122 3 Modeling the market In the swaption world the underlying rate (3.57) has a bounded range and thus it does not display the explosive pattern typical of a stock price. Therefore the swaption prices
More information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationSystems of Ordinary Differential Equations. Lectures INF2320 p. 1/48
Systems of Ordinary Differential Equations Lectures INF2320 p. 1/48 Lectures INF2320 p. 2/48 ystems of ordinary differential equations Last two lectures we have studied models of the form y (t) = F(y),
More informationChapter 2 Portfolio Management and the Capital Asset Pricing Model
Chapter 2 Portfolio Management and the Capital Asset Pricing Model In this chapter, we explore the issue of risk management in a portfolio of assets. The main issue is how to balance a portfolio, that
More informationInvariant Variational Problems & Integrable Curve Flows. Peter J. Olver University of Minnesota olver
Invariant Variational Problems & Integrable Curve Flows Peter J. Olver University of Minnesota http://www.math.umn.edu/ olver Cocoyoc, November, 2005 1 Variational Problems x = (x 1,..., x p ) u = (u 1,...,
More informationELEMENTS OF MATRIX MATHEMATICS
QRMC07 9/7/0 4:45 PM Page 5 CHAPTER SEVEN ELEMENTS OF MATRIX MATHEMATICS 7. AN INTRODUCTION TO MATRICES Investors frequently encounter situations involving numerous potential outcomes, many discrete periods
More informationM5MF6. Advanced Methods in Derivatives Pricing
Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................
More informationProject 1: Double Pendulum
Final Projects Introduction to Numerical Analysis II http://www.math.ucsb.edu/ atzberg/winter2009numericalanalysis/index.html Professor: Paul J. Atzberger Due: Friday, March 20th Turn in to TA s Mailbox:
More informationApplication of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem
Application of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem Malgorzata A. Jankowska 1, Andrzej Marciniak 2 and Tomasz Hoffmann 2 1 Poznan University
More informationOn Toponogov s Theorem
On Toponogov s Theorem Viktor Schroeder 1 Trigonometry of constant curvature spaces Let κ R be given. Let M κ be the twodimensional simply connected standard space of constant curvature κ. Thus M κ is
More informationPortfolio Choice. := δi j, the basis is orthonormal. Expressed in terms of the natural basis, x = j. x j x j,
Portfolio Choice Let us model portfolio choice formally in Euclidean space. There are n assets, and the portfolio space X = R n. A vector x X is a portfolio. Even though we like to see a vector as coordinate-free,
More informationMECHANICS OF MATERIALS
CHAPTER 7 Transformations MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf Lecture Notes: J. Walt Oler Texas Tech Universit of Stress and Strain 006 The McGraw-Hill Companies,
More informationThe Margins of Global Sourcing: Theory and Evidence from U.S. Firms by Pol Antràs, Teresa C. Fort and Felix Tintelnot
The Margins of Global Sourcing: Theory and Evidence from U.S. Firms by Pol Antràs, Teresa C. Fort and Felix Tintelnot Online Theory Appendix Not for Publication) Equilibrium in the Complements-Pareto Case
More informationChapter 5 Finite Difference Methods. Math6911 W07, HM Zhu
Chapter 5 Finite Difference Methods Math69 W07, HM Zhu References. Chapters 5 and 9, Brandimarte. Section 7.8, Hull 3. Chapter 7, Numerical analysis, Burden and Faires Outline Finite difference (FD) approximation
More informationFinancial Market Models. Lecture 1. One-period model of financial markets & hedging problems. Imperial College Business School
Financial Market Models Lecture One-period model of financial markets & hedging problems One-period model of financial markets a 4 2a 3 3a 3 a 3 -a 4 2 Aims of section Introduce one-period model with finite
More informationUniform Refraction in Negative Refractive Index Materials
Haverford College Haverford Scholarship Faculty Publications Mathematics 2015 Uniform Refraction in Negative Refractive Index Materials Eric Stachura Haverford College, estachura@haverford.edu Cristian
More informationLeader or Follower? A Payoff Analysis in Quadratic Utility Harsanyi Economy
Leader or Follower? A Payoff Analysis in Quadratic Utility Harsanyi Economy Sai Ma New York University Oct. 0, 015 Model Agents and Belief There are two players, called agent i {1, }. Each agent i chooses
More information9.1 Principal Component Analysis for Portfolios
Chapter 9 Alpha Trading By the name of the strategies, an alpha trading strategy is to select and trade portfolios so the alpha is maximized. Two important mathematical objects are factor analysis and
More informationLecture Quantitative Finance Spring Term 2015
implied Lecture Quantitative Finance Spring Term 2015 : May 7, 2015 1 / 28 implied 1 implied 2 / 28 Motivation and setup implied the goal of this chapter is to treat the implied which requires an algorithm
More informationEE 521 Instrumentation and Measurements Fall 2007 Solutions for homework assignment #2
Problem 1 (1) EE 51 Instrumentation and Measurements Fall 007 Solutions for homework assignment # f(x) = 1 σ π e () If the height of the peaks in the distribution as drawn are assumed to be 1, then the
More informationChapter 10: Mixed strategies Nash equilibria, reaction curves and the equality of payoffs theorem
Chapter 10: Mixed strategies Nash equilibria reaction curves and the equality of payoffs theorem Nash equilibrium: The concept of Nash equilibrium can be extended in a natural manner to the mixed strategies
More informationLecture 7: Bayesian approach to MAB - Gittins index
Advanced Topics in Machine Learning and Algorithmic Game Theory Lecture 7: Bayesian approach to MAB - Gittins index Lecturer: Yishay Mansour Scribe: Mariano Schain 7.1 Introduction In the Bayesian approach
More informationNon replication of options
Non replication of options Christos Kountzakis, Ioannis A Polyrakis and Foivos Xanthos June 30, 2008 Abstract In this paper we study the scarcity of replication of options in the two period model of financial
More informationA THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES
Proceedings of ALGORITMY 01 pp. 95 104 A THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES BEÁTA STEHLÍKOVÁ AND ZUZANA ZÍKOVÁ Abstract. A convergence model of interest rates explains the evolution of the
More informationA No-Arbitrage Theorem for Uncertain Stock Model
Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe
More informationSupplementary online material to Information tradeoffs in dynamic financial markets
Supplementary online material to Information tradeoffs in dynamic financial markets Efstathios Avdis University of Alberta, Canada 1. The value of information in continuous time In this document I address
More informationMethod of Characteristics
The Ryan C. Trinity University Partial Differential Equations January 22, 2015 Linear and Quasi-Linear (first order) PDEs A PDE of the form A(x,y) u x +B(x,y) u y +C 1(x,y)u = C 0 (x,y) is called a (first
More informationLaurence Boxer and Ismet KARACA
SOME PROPERTIES OF DIGITAL COVERING SPACES Laurence Boxer and Ismet KARACA Abstract. In this paper we study digital versions of some properties of covering spaces from algebraic topology. We correct and
More informationMath F412: Homework 4 Solutions February 20, κ I = s α κ α
All prts of this homework to be completed in Mple should be done in single worksheet. You cn submit either the worksheet by emil or printout of it with your homework. 1. Opre 1.4.1 Let α be not-necessrily
More informationSmooth estimation of yield curves by Laguerre functions
Smooth estimation of yield curves by Laguerre functions A.S. Hurn 1, K.A. Lindsay 2 and V. Pavlov 1 1 School of Economics and Finance, Queensland University of Technology 2 Department of Mathematics, University
More informationA Numerical Approach to the Estimation of Search Effort in a Search for a Moving Object
Proceedings of the 1. Conference on Applied Mathematics and Computation Dubrovnik, Croatia, September 13 18, 1999 pp. 129 136 A Numerical Approach to the Estimation of Search Effort in a Search for a Moving
More informationInformation, Interest Rates and Geometry
Information, Interest Rates and Geometry Dorje C. Brody Department of Mathematics, Imperial College London, London SW7 2AZ www.imperial.ac.uk/people/d.brody (Based on work in collaboration with Lane Hughston
More informationCharacterization of the Optimum
ECO 317 Economics of Uncertainty Fall Term 2009 Notes for lectures 5. Portfolio Allocation with One Riskless, One Risky Asset Characterization of the Optimum Consider a risk-averse, expected-utility-maximizing
More informationQuantitative Risk Management
Quantitative Risk Management Asset Allocation and Risk Management Martin B. Haugh Department of Industrial Engineering and Operations Research Columbia University Outline Review of Mean-Variance Analysis
More informationMartingale Pricing Theory in Discrete-Time and Discrete-Space Models
IEOR E4707: Foundations of Financial Engineering c 206 by Martin Haugh Martingale Pricing Theory in Discrete-Time and Discrete-Space Models These notes develop the theory of martingale pricing in a discrete-time,
More informationClassifying Solvency Capital Requirement Contribution of Collective Investments under Solvency II
Classifying Solvency Capital Requirement Contribution of Collective Investments under Solvency II Working Paper Series 2016-03 (01) SolvencyAnalytics.com March 2016 Classifying Solvency Capital Requirement
More informationSimulating more interesting stochastic processes
Chapter 7 Simulating more interesting stochastic processes 7. Generating correlated random variables The lectures contained a lot of motivation and pictures. We'll boil everything down to pure algebra
More informationLikelihood Methods of Inference. Toss coin 6 times and get Heads twice.
Methods of Inference Toss coin 6 times and get Heads twice. p is probability of getting H. Probability of getting exactly 2 heads is 15p 2 (1 p) 4 This function of p, is likelihood function. Definition:
More informationMATH 5510 Mathematical Models of Financial Derivatives. Topic 1 Risk neutral pricing principles under single-period securities models
MATH 5510 Mathematical Models of Financial Derivatives Topic 1 Risk neutral pricing principles under single-period securities models 1.1 Law of one price and Arrow securities 1.2 No-arbitrage theory and
More informationHomework # 8 - [Due on Wednesday November 1st, 2017]
Homework # 8 - [Due on Wednesday November 1st, 2017] 1. A tax is to be levied on a commodity bought and sold in a competitive market. Two possible forms of tax may be used: In one case, a per unit tax
More informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationModelling strategies for bivariate circular data
Modelling strategies for bivariate circular data John T. Kent*, Kanti V. Mardia, & Charles C. Taylor Department of Statistics, University of Leeds 1 Introduction On the torus there are two common approaches
More informationAsymmetric Information: Walrasian Equilibria, and Rational Expectations Equilibria
Asymmetric Information: Walrasian Equilibria and Rational Expectations Equilibria 1 Basic Setup Two periods: 0 and 1 One riskless asset with interest rate r One risky asset which pays a normally distributed
More informationCEO Attributes, Compensation, and Firm Value: Evidence from a Structural Estimation. Internet Appendix
CEO Attributes, Compensation, and Firm Value: Evidence from a Structural Estimation Internet Appendix A. Participation constraint In evaluating when the participation constraint binds, we consider three
More information3.2 No-arbitrage theory and risk neutral probability measure
Mathematical Models in Economics and Finance Topic 3 Fundamental theorem of asset pricing 3.1 Law of one price and Arrow securities 3.2 No-arbitrage theory and risk neutral probability measure 3.3 Valuation
More information3. The Discount Factor
3. he Discount Factor Objectives Eplanation of - Eistence of Discount Factors: Necessary and Sufficient Conditions - Positive Discount Factors: Necessary and Sufficient Conditions Contents 3. he Discount
More informationarxiv: v1 [q-fin.mf] 16 Jan 2019
arxiv:1901.05113v1 [q-fin.mf] 16 Jan 2019 Instantaneous Arbitrage and the CAPM Lars Tyge Nielsen Department of Mathematics Columbia University January 2019 Abstract This paper studies the concept of instantaneous
More informationLabor Economics Field Exam Spring 2011
Labor Economics Field Exam Spring 2011 Instructions You have 4 hours to complete this exam. This is a closed book examination. No written materials are allowed. You can use a calculator. THE EXAM IS COMPOSED
More informationOnline Appendices to Financing Asset Sales and Business Cycles
Online Appendices to Financing Asset Sales usiness Cycles Marc Arnold Dirk Hackbarth Tatjana Xenia Puhan August 22, 2017 University of St. allen, Unterer raben 21, 9000 St. allen, Switzerl. Telephone:
More informationCHOICE THEORY, UTILITY FUNCTIONS AND RISK AVERSION
CHOICE THEORY, UTILITY FUNCTIONS AND RISK AVERSION Szabolcs Sebestyén szabolcs.sebestyen@iscte.pt Master in Finance INVESTMENTS Sebestyén (ISCTE-IUL) Choice Theory Investments 1 / 65 Outline 1 An Introduction
More informationThe mean-variance portfolio choice framework and its generalizations
The mean-variance portfolio choice framework and its generalizations Prof. Massimo Guidolin 20135 Theory of Finance, Part I (Sept. October) Fall 2014 Outline and objectives The backward, three-step solution
More informationROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices
ROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices Bachelier Finance Society Meeting Toronto 2010 Henley Business School at Reading Contact Author : d.ledermann@icmacentre.ac.uk Alexander
More informationLinearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing
Linearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing Liuren Wu, Baruch College Joint work with Peter Carr and Xavier Gabaix at New York University Board of
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationTechniques for Calculating the Efficient Frontier
Techniques for Calculating the Efficient Frontier Weerachart Kilenthong RIPED, UTCC c Kilenthong 2017 Tee (Riped) Introduction 1 / 43 Two Fund Theorem The Two-Fund Theorem states that we can reach any
More information3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors
3.4 Copula approach for modeling default dependency Two aspects of modeling the default times of several obligors 1. Default dynamics of a single obligor. 2. Model the dependence structure of defaults
More informationWhat Can Rational Investors Do About Excessive Volatility and Sentiment Fluctuations?
What Can Rational Investors Do About Excessive Volatility and Sentiment Fluctuations? Bernard Dumas INSEAD, Wharton, CEPR, NBER Alexander Kurshev London Business School Raman Uppal London Business School,
More informationMATH362 Fundamentals of Mathematical Finance. Topic 1 Mean variance portfolio theory. 1.1 Mean and variance of portfolio return
MATH362 Fundamentals of Mathematical Finance Topic 1 Mean variance portfolio theory 1.1 Mean and variance of portfolio return 1.2 Markowitz mean-variance formulation 1.3 Two-fund Theorem 1.4 Inclusion
More informationFluctuations. Shocks, Uncertainty, and the Consumption/Saving Choice
Fluctuations. Shocks, Uncertainty, and the Consumption/Saving Choice Olivier Blanchard April 2005 14.452. Spring 2005. Topic2. 1 Want to start with a model with two ingredients: Shocks, so uncertainty.
More informationGirsanov s Theorem. Bernardo D Auria web: July 5, 2017 ICMAT / UC3M
Girsanov s Theorem Bernardo D Auria email: bernardo.dauria@uc3m.es web: www.est.uc3m.es/bdauria July 5, 2017 ICMAT / UC3M Girsanov s Theorem Decomposition of P-Martingales as Q-semi-martingales Theorem
More informationQuantitative Techniques (Finance) 203. Derivatives for Functions with Multiple Variables
Quantitative Techniques (Finance) 203 Derivatives for Functions with Multiple Variables Felix Chan October 2006 1 Introduction In the previous lecture, we discussed the concept of derivative as approximation
More informationPermutation Factorizations and Prime Parking Functions
Permutation Factorizations and Prime Parking Functions Amarpreet Rattan Department of Combinatorics and Optimization University of Waterloo Waterloo, ON, Canada N2L 3G1 arattan@math.uwaterloo.ca June 10,
More informationStatistics 13 Elementary Statistics
Statistics 13 Elementary Statistics Summer Session I 2012 Lecture Notes 5: Estimation with Confidence intervals 1 Our goal is to estimate the value of an unknown population parameter, such as a population
More informationMATH4512 Fundamentals of Mathematical Finance. Topic Two Mean variance portfolio theory. 2.1 Mean and variance of portfolio return
MATH4512 Fundamentals of Mathematical Finance Topic Two Mean variance portfolio theory 2.1 Mean and variance of portfolio return 2.2 Markowitz mean-variance formulation 2.3 Two-fund Theorem 2.4 Inclusion
More informationSHORT-TERM RELATIVE ARBITRAGE IN VOLATILITY-STABILIZED MARKETS
SHORT-TERM RELATIVE ARBITRAGE IN VOLATILITY-STABILIZED MARKETS ADRIAN D. BANNER INTECH One Palmer Square Princeton, NJ 8542, USA adrian@enhanced.com DANIEL FERNHOLZ Department of Computer Sciences University
More informationPDE Methods for the Maximum Drawdown
PDE Methods for the Maximum Drawdown Libor Pospisil, Jan Vecer Columbia University, Department of Statistics, New York, NY 127, USA April 1, 28 Abstract Maximum drawdown is a risk measure that plays an
More informationMFE Macroeconomics Week 3 Exercise
MFE Macroeconomics Week 3 Exercise The first row in the figure below shows monthly data for the Federal Funds Rate and CPI inflation for the period 199m1-18m8. 1 FFR CPI inflation 8 1 6 4 1 199 1995 5
More informationParametric Inference and Dynamic State Recovery from Option Panels. Nicola Fusari
Parametric Inference and Dynamic State Recovery from Option Panels Nicola Fusari Joint work with Torben G. Andersen and Viktor Todorov July 2012 Motivation Under realistic assumptions derivatives are nonredundant
More informationA Confluent Hypergeometric System Associated with Ф 3 and a Confluent Jordan-Pochhammer Equation. Shun Shimomura
Research Report KSTS/RR-98/002 Mar. 11, 1998 A Confluent Hypergeometric System Associated with Ф 3 and a Confluent Jordan-Pochhammer Equation by Department of Mathematics Keio University Department of
More informationElements of Economic Analysis II Lecture II: Production Function and Profit Maximization
Elements of Economic Analysis II Lecture II: Production Function and Profit Maximization Kai Hao Yang 09/26/2017 1 Production Function Just as consumer theory uses utility function a function that assign
More informationFinancial Modelling Using Discrete Stochastic Calculus
Preprint typeset in JHEP style - HYPER VERSION Financial Modelling Using Discrete Stochastic Calculus Eric A. Forgy, Ph.D. E-mail: eforgy@yahoo.com Abstract: In the present report, a review of discrete
More informationConsumption- Savings, Portfolio Choice, and Asset Pricing
Finance 400 A. Penati - G. Pennacchi Consumption- Savings, Portfolio Choice, and Asset Pricing I. The Consumption - Portfolio Choice Problem We have studied the portfolio choice problem of an individual
More informationStochastic Volatility (Working Draft I)
Stochastic Volatility (Working Draft I) Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu 1 Introduction When using the Black-Scholes-Merton model to price derivative
More informationExercises in Growth Theory and Empirics
Exercises in Growth Theory and Empirics Carl-Johan Dalgaard University of Copenhagen and EPRU May 22, 2003 Exercise 6: Productive government investments and exogenous growth Consider the following growth
More information