Convergence of a Generalized Midpoint Iteration

Similar documents
Course 212: Academic Year Section 1: Metric Spaces

M17 MAT25-21 HOMEWORK 6

7 Complete metric spaces and function spaces

1 Topology Definition of a topology Basis (Base) of a topology The subspace topology & the product topology on X Y 3

METRIC SPACES KEITH CONRAD

Problem Set 2: Solutions Math 201A: Fall 2016

U e = E (U\E) e E e + U\E e. (1.6)

Some Background Material

MATH 51H Section 4. October 16, Recall what it means for a function between metric spaces to be continuous:

Mathematics for Economists

Lecture Notes in Advanced Calculus 1 (80315) Raz Kupferman Institute of Mathematics The Hebrew University

The small ball property in Banach spaces (quantitative results)

Introduction to Real Analysis Alternative Chapter 1

Chapter 2 Metric Spaces

Set, functions and Euclidean space. Seungjin Han

R N Completeness and Compactness 1

Maths 212: Homework Solutions

Definition 2.1. A metric (or distance function) defined on a non-empty set X is a function d: X X R that satisfies: For all x, y, and z in X :

CLASS NOTES FOR APRIL 14, 2000

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi

Metric Spaces and Topology

Mid Term-1 : Practice problems

LECTURE 15: COMPLETENESS AND CONVEXITY

MATH 4200 HW: PROBLEM SET FOUR: METRIC SPACES

Math LM (24543) Lectures 01

A VERY BRIEF REVIEW OF MEASURE THEORY

MATH 202B - Problem Set 5

Chapter II. Metric Spaces and the Topology of C

Continuity. Matt Rosenzweig

Part III. 10 Topological Space Basics. Topological Spaces

REVIEW OF ESSENTIAL MATH 346 TOPICS

The Heine-Borel and Arzela-Ascoli Theorems

Econ Lecture 3. Outline. 1. Metric Spaces and Normed Spaces 2. Convergence of Sequences in Metric Spaces 3. Sequences in R and R n

Real Analysis. Joe Patten August 12, 2018

01. Review of metric spaces and point-set topology. 1. Euclidean spaces

Economics 204 Fall 2012 Problem Set 3 Suggested Solutions

Analysis Comprehensive Exam Questions Fall F(x) = 1 x. f(t)dt. t 1 2. tf 2 (t)dt. and g(t, x) = 2 t. 2 t

Analysis Finite and Infinite Sets The Real Numbers The Cantor Set

Math 320-2: Final Exam Practice Solutions Northwestern University, Winter 2015

EC 521 MATHEMATICAL METHODS FOR ECONOMICS. Lecture 1: Preliminaries

What to remember about metric spaces

Math 320-2: Midterm 2 Practice Solutions Northwestern University, Winter 2015

Math 421, Homework #7 Solutions. We can then us the triangle inequality to find for k N that (x k + y k ) (L + M) = (x k L) + (y k M) x k L + y k M

Characterisation of Accumulation Points. Convergence in Metric Spaces. Characterisation of Closed Sets. Characterisation of Closed Sets

THEOREMS, ETC., FOR MATH 515

SOLUTIONS TO SOME PROBLEMS

MAS331: Metric Spaces Problems on Chapter 1

2. Metric Spaces. 2.1 Definitions etc.

CHAPTER 1. Metric Spaces. 1. Definition and examples

2 Topology of a Metric Space

MAT 570 REAL ANALYSIS LECTURE NOTES. Contents. 1. Sets Functions Countability Axiom of choice Equivalence relations 9

Existence and Uniqueness

Contents Ordered Fields... 2 Ordered sets and fields... 2 Construction of the Reals 1: Dedekind Cuts... 2 Metric Spaces... 3

Math 201 Topology I. Lecture notes of Prof. Hicham Gebran

THE CONTRACTION MAPPING THEOREM, II

4 Countability axioms

Math 4317 : Real Analysis I Mid-Term Exam 1 25 September 2012

Chapter 2. Metric Spaces. 2.1 Metric Spaces

II - REAL ANALYSIS. This property gives us a way to extend the notion of content to finite unions of rectangles: we define

Overview of normed linear spaces

Math 341: Convex Geometry. Xi Chen

Introduction to Dynamical Systems

Convex Analysis and Economic Theory Winter 2018

Simple Abelian Topological Groups. Luke Dominic Bush Hipwood. Mathematics Institute

Topological properties of Z p and Q p and Euclidean models

Real Analysis Notes. Thomas Goller

*Room 3.13 in Herschel Building

Part V. 17 Introduction: What are measures and why measurable sets. Lebesgue Integration Theory

Introduction to Topology

Functional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...

Topology. Xiaolong Han. Department of Mathematics, California State University, Northridge, CA 91330, USA address:

Lecture 3. Econ August 12

Fourth Week: Lectures 10-12

Appendix B Convex analysis

Sequence of maximal distance codes in graphs or other metric spaces

The Arzelà-Ascoli Theorem

15 Real Analysis II Sequences and Limits , Math for Economists Fall 2004 Lecture Notes, 10/21/2004

Topology, Math 581, Fall 2017 last updated: November 24, Topology 1, Math 581, Fall 2017: Notes and homework Krzysztof Chris Ciesielski

Metric Spaces Math 413 Honors Project

1 Lecture 4: Set topology on metric spaces, 8/17/2012

Metric Spaces Math 413 Honors Project

MA677 Assignment #3 Morgan Schreffler Due 09/19/12 Exercise 1 Using Hölder s inequality, prove Minkowski s inequality for f, g L p (R d ), p 1:

Completeness and quasi-completeness. 1. Products, limits, coproducts, colimits

Discrete Geometry. Problem 1. Austin Mohr. April 26, 2012

Near convexity, metric convexity, and convexity

MATH 54 - TOPOLOGY SUMMER 2015 FINAL EXAMINATION. Problem 1

BALL-POLYHEDRA. 1. Introduction

Math 118B Solutions. Charles Martin. March 6, d i (x i, y i ) + d i (y i, z i ) = d(x, y) + d(y, z). i=1

Chapter 3: Baire category and open mapping theorems

Locally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem

Functional Analysis Exercise Class

Indeed, if we want m to be compatible with taking limits, it should be countably additive, meaning that ( )

Optimization and Optimal Control in Banach Spaces

A FIXED POINT THEOREM FOR GENERALIZED NONEXPANSIVE MULTIVALUED MAPPINGS

converges as well if x < 1. 1 x n x n 1 1 = 2 a nx n

2 Metric Spaces Definitions Exotic Examples... 3

Euclidean Models of the p-adic Integers

Math Tune-Up Louisiana State University August, Lectures on Partial Differential Equations and Hilbert Space

Measure Theory on Topological Spaces. Course: Prof. Tony Dorlas 2010 Typset: Cathal Ormond

MATH 31BH Homework 1 Solutions

Transcription:

J. Able, D. Bradley, A.S. Moon under the supervision of Dr. Xingping Sun REU Final Presentation July 31st, 2014

Preliminary Words O Rourke s conjecture We begin with a motivating question concerning the convergence of sequences of centroid polyhedra. Conjecture Given any convex polyhedron P (0) in R 3, the iteration process which takes the convex hull c(p (n) ) of the centroids of faces of P (n) converges to an ellipsoid as n.

Preliminary Words

Preliminary Words O Rourke s conjecture To work towards a solution of O Rourke s conjecture, we examine a simpler iteration of midpoints in R 2 similar to the midpoint iteration of convex polygons.

Preliminary Words Notation Let d(, ) denote the Euclidean distance between two points. Let denote the Euclidean norm. Let c( ) denote the convex hull of a set. Let cl( ) denote set closure (closure includes limit points). Let B(x, ε) denote the open ball of radius ε about x. Let diam( ) denote the diameter of a set.

Preliminary Words Hausdorff distance For two nonempty compact subsets X, Y of (R 2, d), we denote the Hausdorff distance by d H (X, Y ) = max { sup x X inf d(x, y), sup y Y y Y inf d(y, x) x X }.

Propositions About Compact Sets Propositions about compact sets We identify a polygon P with the convex hull of a finite set of vertices in R 2. Recall the Heine-Borel characterization of compact subsets of R 2 (compact iff bounded and closed). Then P is a bounded and closed set of the plane. There are many useful results about nested sequences of compact subsets which can be applied to polygons on the plane.

Propositions About Compact Sets Propositions about compact sets Proposition Suppose K n+1 K n for all n N, where all K n are compact, nonempty sets. Then the intersection I = n N K n. is nonempty and compact. Proof. I is closed and bounded, since arbitrary intersections of closed sets are closed and I K 1. Define a sequence {a p } such that a p K p \ K p 1. Then there exists by compactness a subsequence {a nk } with limit in I.

Propositions About Compact Sets Propositions about compact sets Proposition Let I = n N K n be the intersection of K n. If lim n diam(k n ) = 0, then I = {x 0 } for some x 0 K 1. Proof. If I is not a singleton, then diam(i ) > 0. But I K n for all n, implying lim diam(k n ) > 0.

Propositions About Compact Sets Propositions about compact sets Proposition Let {K n } n N be a sequence of compact subsets such that K n+1 K n for all n N, and I their infinite intersection n N K n. Then lim d H (K n, I ) = 0.

Propositions About Compact Sets Proof. By compactness, choose sequences {x n }, {y n } s.t. x n K n \ I and y n I for all n N, such that {x n } α, {y n } β and d(x n, y n ) = sup inf d(x, y). x K n y I By a contradiction, α I, so that And we know: d(x n, y n ) = inf y I d(x n, y) d(x n, α) d(x n, α) 0 as n

Propositions About Compact Sets Propositions about compact sets Proposition Suppose {A n } is a sequence of compact sets such that A n A n+1 and A n N K for all n N and some compact K. Let U = n A n. Then lim A n = cl(u). Proof. Similar to previous proof concerning intersections. Identify convergent subsequences of relevant sequences and bound the Hausdorff distance between the sequences so that this distance approaches 0.

Midpoint Polygon Sequences Background First solved by the French mathematician J.G. Darboux, the midpoint polygon problem has been examined using diverse techniques, including finite Fourier analysis and matrix products [1] [2] [3] [4]. We present an elementary analysis proof of the midpoint polygon problem.

Midpoint Polygon Sequences Proposition The original problem statement from the American Mathematical Monthly reads Proposition Let Π = P 0 P 1 P k 1 be a closed polygon in the plane (k 2). Denote by P 0, P 1,..., P k 1 the midpoints of the sides P 0 P 1, P 1 P 2,..., P k 2 P k 1, respectively, obtaining the fist derived polygon Π = P 0 P 1 P k 1. Repeat the same construction on Π obtaining the second derived polygon Π, and finally, after n constructions, obtain the nth derived polygon Π (n) = P (n) 0 P(n) 1 P (n) k 1. Show that the vertices of Π(n) converge, as n, to the centroid of the original points P 0, P 1,..., P k 1.

Midpoint Polygon Sequences Midpoint problem In the language of elementary analysis and Hausdorff convergence, Proposition Given an initial set of vertices { } P (0) = v (0) 1,..., v (0) Q in general linear position, the sequence of convex hulls c(p (n) ) produced by midpoint iteration converge in the Hausdorff metric to the centroid of c(p (0) ). Proof. Without loss of generality, suppose the centroid of P (0) is the origin. The Hausdorff limit lim P (n) is equal to the intersection n N c(p(n) ), which we know is nonempty and compact.

Midpoint Polygon Sequences Midpoint problem Generally, the edge lengths are bounded: sup v (n) i v (n) i+1 (1 i 2 )n 1 diam(c(p (0) )) And the perimeter bounds the diameter of a closed convex polygon, so diam(c(p (k) )) Q i=1 v (n) i v (n) i+1 Q(1 2 )n 1 diam(c(p (0) )). Conclude lim diam(c(p (n) )) = 0. Then lim c(p (n) ) = {x 0 } for some x 0. But by calculation, we see that the centroid of P (n+1) is identically the centroid of P (n), the origin. So 0 c(p (n) ) is the limit of { c(p (n) ) }

Midpoint Polygon Sequences Regularity in the limit Proposition Let e (n) j denote the jth edge length v (n) j v (n) j+1 of the nth iteration of the initial polygon of Q vertices. Then lim n e (n) i e (n) j = 0 for any i, j indices modulo Q. Proof. e (n) i e (n) j e (n) i + e (n) j 2 diam(p (n) ) since every edge is bounded by the diameter. Now lim diam(c(p (n) )) = 0, so lim n e (n) i e (n) j = 0.

Midpoint Polygon Sequences Regularity in the limit By the previous proposition, given ε > 0, there is a large enough N such that n > N implies e (n) i for all i, j indices taken modulo Q. e (n) j < ε

Generalized Midpoint Iteration Generalized Midpoint Iteration Instead of restricting to midpoints between adjacent vertices, we now consider midpoints of every pair of vertices. This resembles the O Rourke iteration of centroids of faces of polyhedra because each iteration contains more points than the previous. We will find that this process can be characterized by fixed points.

Generalized Midpoint Iteration Generalized Midpoint Iteration Informally, the Generalized Midpoint Iteration (GMI) on P (0) produces vertex sets P (k) containing every possible midpoint combination of the previous vertex set P (k 1), and records how many times p P (k) occurs as a midpoint by the multiplicity function µ k : P (k) N.

Generalized Midpoint Iteration

Generalized Midpoint Iteration

Generalized Midpoint Iteration

Generalized Midpoint Iteration

Generalized Midpoint Iteration

Generalized Midpoint Iteration Subsystems Within the GMI, we can consider individual sets of three points and their direct descendants for all iterations. The centroid of these subsystems will be in the limit by the classical midpoint iteration on a polygon. This tells us that that GMI does not converge to a single point like the classical midpoint iteration.

Generalized Midpoint Iteration Fixed points A point is fixed if it remains in all succeeding iterations. Under the GMI, a point becomes fixed if it has multiplicity greater than 2 (since there are three ways of pairing a multiplicity 3 point with itself). Interestingly, we can guarantee, under the GMI on polygons with more than 3 vertices, the existence of explicit fixed points past the second iteration. Also, since the midpoint of a fixed point with any other point is fixed (has multiplicity greater than 2), we can eventually fill up the object of convergence through one initial fixed point.

Generalized Midpoint Iteration Fixed points Given any four points, we can find a fixed point of multiplicity 3 after just two iterations using these combinations: 1 ( ) 3p (0) 1 + 2p (0) 2 + 2p (0) 3 + p (0) 4 4 = 1 ( 1 2 2 (p(0) 1 + p (0) 2 ) + 1 2 (p(0) 1 + p (0) 3 )) + 1 2 (1 2 (p(0) 1 + p (0) 2 ) + 1 ) 2 (p(0) 3 + p (0) 4 ) = 1 ( 1 2 2 (p(0) 1 + p (0) 2 ) + 1 2 (p(0) 1 + p (0) 4 )) + 1 2 (1 2 (p(0) 1 + p (0) 3 ) + 1 ) 2 (p(0) 2 + p (0) 3 ) = 1 ( 1 2 2 (p(0) 1 + p (0) 2 ) + 1 2 (p(0) 2 + p (0) 3 )) + 1 2 (1 2 (p(0) 1 + p (0) 3 ) + 1 ) 2 (p(0) 1 + p (0) 4 )

Generalized Midpoint Iteration GMI Examples

Hausdorff Convergence Results for the GMI Goal We now prove several convergence results involving the sets of fixed points F (n), iteration midpoints P (n), and their convex hulls c(f (n) ) and c(p (n) ) respectively. lim d H (F (n), c(f (n) )) = 0 lim d H (P (n), F (n) ) = 0 lim d H (P (n), c(f (n) )) = 0

Hausdorff Convergence Results for the GMI Intuitition

Hausdorff Convergence Results for the GMI Existence of limits We observe that the sequence of convex hulls { c(p (n) ) } n N form a nested sequence of compact sets. Proposition lim c(p (n) ) = c(p (n) ) Proof. This is a direct consequence of the earlier proposition about the intersection of nested compact sets, satisfied because the c(p (n) ) are certainly nested and compact.

Hausdorff Convergence Results for the GMI Existence of limits Proposition lim c(f (n) ) exists, with respect to the Hausdorff metric. Proof. By the definition of fixed points, F (n) F (n+1). This implies the convex hulls c(f (n) ) c(f (n+1) ) form a reverse-nested sequence of subsets. Lastly, c(f (n) ) c(p (0) ) for all n. So { c(f (n) ) } satisfies the conditions of Proposition 3.9. As a remark, this implies { c(f (n) ) } are a Cauchy sequence.

Hausdorff Convergence Results for the GMI Hausdorff distance Recall the Hausdorff distance between two sets is given by d H (X, Y ) = max { sup x X inf d(x, y), sup y Y y Y inf d(y, x) x X }.

Hausdorff Convergence Results for the GMI Midpoint Triangle Property

Hausdorff Convergence Results for the GMI Propositions Proposition lim d H (F (n), c(f (n) )) = 0 Proof. Note F (n) c(f (n) ), so we consider the distance from c(f (n) ) to F (n). Let x c(f (n) ). There are three noncollinear points in F (n) such that x is in their convex hull. Taking n iterations, we have 4 n congruent triangles, since these points are all fixed, and x is in one of these 4 n triangles.

Hausdorff Convergence Results for the GMI Propositions

Hausdorff Convergence Results for the GMI Propositions Let T n represent a triangle of the nth subdivision containing x. By the Schoenberg midpoint iteration, for arbitrary ε > 0, there exists an N such that for n > N, the vertices {v 1, v 2, v 3 } of T n are within ε distance of the centroid c of T n. Then This implies inf i d H (F (n), c(f (n) )) = d(x, v i ) sup d(c, v j ) < ε j sup x c(f (n) ) inf d(x, y) < ε. y F (n)

Hausdorff Convergence Results for the GMI Propositions Proposition lim d H (P (n), F (n) ) = 0 Proof. Note F (n) P (n), so the Hausdorff distance is equal to d H (P (n), F (n) ) = sup x P (n) inf d(x, y). y F (n)

Hausdorff Convergence Results for the GMI Propositions Let p P (n+1) arbitrarily. By definition of the GMI, there exists at least one pair p, q P (n) such that p+q 2 = p. Then by the Midpoint Triangle Property, d(p, F (n+1) ) 1 2 d(p, F (n) ) 1 2 sup x P (n) inf d(x, y). y F (n)

Hausdorff Convergence Results for the GMI Propositions Since p was arbitrary and P (n) and F (n) are finite sets, we may choose p such that sup d(z, F (n+1) ) = d(p, F (n+1) ). z P (n+1) This is precisely the inequality d H (P (n+1), F (n+1) ) 1 2 d H(P (n), F (n) ).

Hausdorff Convergence Results for the GMI Propositions Proposition lim d H (P (n), c(f (n) )) = 0 Proof. Triangle inequality: d H (P (n), c(f (n) )) d H (P (n), F (n) ) + d H (F (n), c(f (n) )) < 1 2 ε + 1 2 ε = ε for appropriate N.

Hausdorff Convergence Results for the GMI Main result Theorem The sequence of convex hulls of fixed points F (n) and the sequence of convex hulls of P (n) have the same limit. That is, lim c(f (n) ) = lim c(p (n) ) = c(p (n) ) And consequently lim d H (P (n), c(p (n) )) = 0.

Hausdorff Convergence Results for the GMI Main result Proof. For notation, let F = lim c(f (n) ) and P = lim c(p (n) ). Seeking contradiction, suppose F P. Then d H (F, P) 0. By proposition 5.2, the convergent sequence { c(f (n) ) } is also n N a Cauchy sequence. Therefore, for all ε there exists N such that for n, m > N, d H (c(f (n) ), c(f (m) )) < ε.

Hausdorff Convergence Results for the GMI Main result By compactness, choose p c(p (N) ) and f c(f (N) ) such that d(p, f ) = d H (c(p (N) ), c(f (N) )). Since F P, d(p, f ) is necessarily positive. Now, the midpoint of f and p is in F (N+1) by proposition 4.3. Denote g = p+f 2. Then d H (c(f (N) ), c(f (N+1) )) Passing to the limit, inf d(g, y) 1 d(p, f ) > 0. y c(f (N) ) 2 lim N d H (c(f (N) ), c(f (N+1) )) > 0. contradicting the fact that { c(f (k) ) } is a Cauchy sequence. Conclude F = P necessarily.

Hausdorff Convergence Results for the GMI Main result

Hausdorff Convergence Results for the GMI Main result Consequently, d H (P (n), c(p (n) )) d H (P (n), c(f (n) )) + d H (c(f (n) ), c(p (n) )) lim d H (P (n), c(p (n) )) = 0.

Hausdorff Convergence Results for the GMI Main result The takeaway from this result is that the P (n) become dense in their convex hulls c(p (n) ), and, since the convex hulls c(p (n) ) converge to their intersection n N c(p(n) ), the P (n) become dense in the limit. The c(f (n) ) are approaching the limit from the inside, the c(p (n) ) are approaching the limit from the outside, and the P (n) are dense in the limit where they meet.

Unit Square Example Unit square example We still don t know anything about the shape of the limit, whether it becomes smooth or elliptical. So we considered the simplest case of the GMI on a unit square.

Unit Square Example Unit square example It turns out that the unit square limits on the boundary to a 44-gon! (The inside will be filled up because of our major theorem on density.) This is unexpected behavior because we expected a circle, or at least an ellipse, or at the very least something smooth. While this behavior makes the GMI unlikely to directly generalize to O Rourke s conjecture, it is still very interesting.

Unit Square Example Unit square example Surely if we ignore multiplicity and fixed points, the unit square would limit to a circle, right? It looks like it!

Unit Square Example Unit square example But actually our calculations in Mathematica suggest that this limit is not a circle. After five iterations, the area of the limit is less than the area of the smallest circle containing all the known limit points. This suggests that O Rourke s conjecture may be false. However, there are many differences between the GMI and the O Rourke iteration, including using centroids instead of midpoints, using R 3 instead of R 2, and taking only centroids on the boundary instead of all centroids. In general, however, this sort of iteration probably does not produce strictly elliptical bodies.

Unit Square Example

Conclusion What we showed In the GMI, the points become dense in the limit, which is equivalently the union of the convex hulls of fixed points or the intersection of the convex hulls of all points. This process turns a discrete set into a subset of R 2.

Conclusion What we didn t show We weren t able to prove O Rourke s conjecture about the centroid-of-faces iteration converging to an ellipsoid. In fact, our results do not include a characterization of the shape of the limit set. But the theory we developed still may lead in some way to a solution of O Rourke s conjecture.

Conclusion Multiplicity and fixed points We developed the concept of multiplicity in order to find some structure in this ever-increasing set of points under the GMI. However, we saw that fixed points from multiplicity lead to a strange convergence result on the unit square, which does not become smooth but converges to a 44-gon. Since fixed points do not lead to the convergence results we expected, they may not be the right tool to work towards a proof of O Rourke s conjecture (although they are interesting in their own right!). However, the proof techniques may still apply.

Conclusion Future work Prove a similar density result (lim d H (P (n), c(p (n) ) = 0) without counting multiplicity and without the fixed points that come from multiplicity. Prove results on the shape of the limit does the GMI converge to an ellipse, some sort of generalized ellipse, or a union of ellipses? (Or just to some blob) Consider a generalized centroid iteration, taking centroids of three points instead of midpoints. More generally, we could iterate by any convex combination of points and hope for some sort of result. Extend to three dimensions or to general R n. Some of our proofs rely on the structure of the plane. (Or to general normed linear spaces eg consider functions instead of points.)

Conclusion Darboux, G., Sur un probleme de geometrie elementaire, Bulletin des sciences mathematiques et astronomiques 2e serie, 2 (1878), 298-304. Ouyang, Charles, A Problem Concerning the Dynamic Geometry of Polygons (2013). Schoenberg, I. J., The finite Fourier series and elementary geometry, Amer. Math. Monthly, 57 (1950), 390404. Sun, Xingping and Eric Hintikka, Convergence of Sequences of Polygons, (2014).