A Special Case Solution to the Perspective 3-Point Problem William J. Wolfe California State University Channel Islands
|
|
- Elinor Walters
- 5 years ago
- Views:
Transcription
1 A Secial Case Solution to the Persective -Point Problem William J. Wolfe California State University Channel Islands Abstract In this aer we address a secial case of the ersective -oint roblem. We assume that the oints form an equilateral triangle, and that the center of ersective lies on a lane erendicular to the triangle and containing one of its altitudes. We then resent a geometric roof that accounts for all ossible solution sets, from to 4 solutions. We emhasize the shared-oint and shared-side solution airs, as first develoed in reference [5]. Introduction Given the ersective view of non-colinear control oints (i.e.: known ΔABC we wish to determine the ossible camera-triangle configurations. This roblem, called the ersective -oint roblem (PP, has been addressed by several authors. Reference [] shows that there can be at most 4 solutions, and a secific 4-solution case is constructed. Reference [] shows that there are exactly 4 solutions whenever the center of ersective (CP is directly over the orthocenter of the triangle. Reference [] delineates the history of the roblem and discusses the relative effectiveness of the various algebraic solutions that have been roosed in the literature. Reference [4] rovides a comlete algebraic solution along with some geometric techniques to hel visualize configurations with,, or 4 solutions. Reference [5] focuses on the secial case where there is at least one solution with CP on a lane erendicular to the triangle and containing one of the its altitudes, and shows that there are 4 solutions, rovided that CP is within the "danger cylinder". Reference [6] identifies the "danger cylinder" as the region of instability caused by the degeneracy of the Jacobian of the nonlinear equations that define the camera-triangle constraints. Reference [7] suggests that the instability is related to the existence of exactly solutions on the boundary of the danger cylinder. Here we address the same secial case as develoed in [5] but we add an additional assumtion: the known triangle is an equilateral triangle. Then, as should be exected, we can ush the interesting results from [5] even further. In articular, a "Key Theorem" facilitates the derivation of a relatively straightforward geometric solution. To our knowledge, this is the most comlex case for which a comlete geometric characterization of all solutions is unambiguously rovided. The Secial Case We assume a in-hole camera model and that the corresondence between the image oints and the corners of the triangle (A, B, C is known. Figure a shows the notation
2 for the assumed equilateral triangle. As in [5], we assume that there is at least one solution with CP on a lane orthogonal to the triangle and containing one of its altitudes, as shown in figure b. Without any further loss of generality let side AB be the base of that altitude, and label the orthogonal lane π AB. We also know that CPm is a erendicular bisector of AB, and that CPm bisects = ACPB. Figure. The triangle is equilateral with side s and altitude h. O is the orthocenter of the triangle, and m is the midoint of AB. CP is on the lane π AB, which is erendicular to the lane of the triangle and contains the altitudecm. Aroach We start with the triangle in any of the configurations that satisfy our assumtions, such as the configuration shown in figure b. To get multile solutions we maneuver the triangle into other configurations that generate the same ersective measurements. From our assumtions it is clear that the camera's ersective measurements can be used to calculate ( ACPB and φ ( CCPm, as labeled in figure b (and in figure. We think of φ as the "angle of elevation" of the line CPC abovecpm. Because of our assumtions the roblem reduces to rather than measured angles. Generally, the ersective measurements consist of the angles subtended by each of the three sides of the triangle, but in our secial case the roblem reduces to the measurement of, the angle subtended by side AB, and φ, the angle elevation of C. We will use the geometry of the initial configuration to generate all ossible cameratriangle configurations that are consistent with the measured values of and φ. In each
3 case we will rovide geometric constructs and exlicit formulas for the locations of the triangle with resect to the camera. Shared-Side Solution Pair Given the initial configuration (figure b, there is an easy way to generate a second solution. Rotate about AB until C hits the same line of ersective again. This creates a "shared-side solution air" similar to those that are defined in [5] (see figure. This air of solutions shares side AB. We sometimes refer to this air as the "front-back air" because one solution is tyically leaning frontward (toward CP while the other is leaning backward (away from CP. There are a coule of excetions. If the starting solution is such that the triangle's altitude (Cm is erendicular to the line of ersective running through C (CPC, then rotating about AB will not roduce another solution (i.e.: because of tangency. In this case the shared-side "air" reduces to a single solution. This oint of tangency also identifies the largest value of φ that is feasible for the given : φ = arcsin( tan max The other excetion is when the new intersection is behind the camera (i.e.: when the distance from CP to m is less than h. These situations are easy to diagnose and deal with. Figure : The "shared-side" solution air. By rotating about AB the oint C can meet the line of ersective in laces. Both configurations have the same ersective measurements. For our secial case, the ersective measurements reduce to the measurement of angles and φ.
4 Shared-Point Solution Pair Starting with the initial configuration (figure b, while keeing fixed, slide AB toward CP. abel the new location of the triangle A'B'C' (see figure. If A' is closer to CP (as it is in figure then there is a symmetric location with B' closer to CP. This symmetric air will form the "shared-oint solution air" as also described in [5]. We sometimes refer to this air as the "left-right" air. After moving AB to one of these locations we can rotate about A ' B' until C' intersects π AB (see figure 4. Now, we rove the Key Theorem. Figure : While keeing fixed, Δ ABC has been slid to a new location, labeled Δ A'B'C'. Then, the triangle is rotated about A ' B' to move C' into π AB. Figure 4: The triangle can be slid, keeing fixed, to many ossible locations, but in each case we can rotate about A ' B' until C' intersects π AB. On the left side of the figure the triangle has been slid as close to CP as ossible (A', or B', coincident with CP. On the right side of the figure the triangle has been slid as far away from CP as ossible, in which case the left and right handed symmetric air reduces to a single triangle that leans back at the same angle as any of the other configurations shown in the figure. The Key Theorem roves that the height of the intersection (z, as well as the leaning angle, is the same for each of these configurations. 4
5 Key Theorem The height at which C' intersects the orthogonal lane in a shared-oint air, as shown in figures and 4, is given by: tan = h z Proof: Refer to the construction in figure 5. Using a sequence of similar right triangles we have: R B R A Δ R B Δ ( cos( ( sin( sin( ( cos( ( t t + = + = ( ( sin( sin( 5 s x x x B B s s R + = = + Δ Δ Using one more air of similar right triangles, and substituting the results for t and x, we get: tan( sin( 4 s R y t x y B m = = Δ Δ Now, y and z are the legs of a right triangle with hyotenuse h, as shown in figure. Therefore: tan = h z QED. 5
6 Figure 5: ooking down on the lane formed by A, B, and CP. The side AB has been slid, keeing fixed, to an arbitrary osition. For any such osition there is a symmetric osition. The figure labels one of them A B and the other A R BR ( for "left" and R for "right". The Key Theorem shows that the variable y deends only on, indeendent of and. Note that the only roerty of the triangle that is used in the roof of the theorem is that CPm is a erendicular bisector of AB. Therefore the theorem is also true if the triangle is isosceles (with base side AB. 6
7 It is surrising that z deends only on, and not on and. Since z deends only on, we can lot the locus of oints in π AB at height z abovecpm. Each of these oints is the otential location of a share-oint solution air. Consider the line λ, in π AB, at height z above CPm (see figure 4. Whenever the line of ersective (from CP to C' intersects λ we have the otential location of a shared-oint solution air. Figure 6 looks erendicularly at π AB. The circle in the figure has radius h, and it is centered at m, the midoint of AB. The side AB is shown only as a single oint, labeled "m" in the figure. This circle defines the locations of share-side solutions (as was already shown in figure. The shared-side solutions corresond to the intersections of the line of ersective and the circle. The share-oint solutions corresond to the intersections of the line of ersective with λ. Notice that the only intersections on λ that roduce feasible shared-oint candidates are those between c and c. If the intersection is to the left of c then the configuration violates our assumtions (φ > φ max. And, as was reviously exlained, c is the furthest osition that the lean-back solution can reach for our given configuration. Therefore, the only candidates on λ for feasible shared-oint solutions occur between c and c. Figure 6: An orthogonal view of the π AB. ines of ersective emanating from CP that intersect the circle at two laces secify the locations of shared-side solution airs. A line 7
8 of ersective that intersects λ, anywhere between c and c, is the otential location of a shared-oint solution air. Combining Shared-Side and Shared-Point Solutions Position in figure 6 shows a side view of the triangle as it leans back at the angle: arctan tan This corresonds to the far right art of figure 4, where the triangle is in the osition with AB erendicular to CPm and leaning back at the same angle as each of the sharedoint configurations. That is, this is where the "left" and "right" members of the sharedoint air become a single solution. In articular, lines of ersective that intersect λ beyond Position do not corresond to alternate solutions since the triangle cannot be that far away from CP and still have AB erendicular to CPm and subtending the angle. Position corresonds to the smallest value of φ for which shared-oint solutions are ossible. ines of ersective emanating from CP and intersecting the circle below Position corresond to configurations with exactly solutions, coming from shared-side airs. These are somewhat extreme front-back airings, with the extreme case being φ = 0 (CP in the lan of the triangle. As the line of ersective rises (i.e.: φ increases from Position to Position we have both shared-side and shared-oint solution airs, as indicated by the intersection twice with the circle and once with λ. This results in configurations with exactly 4 solutions. Position is where CP is directly over the orthocenter of the triangle. Notice its symmetric ositioning with resect to Position (a consequence of the equilateral triangle assumtion. At Position, each of the shared-oint solutions share a common side with one of the shared-side solutions (the "lean-back" solution, which is consistent with the solutions derived in []. Between Position and Position we still have 4 solutions coming from distinct sharedside and shared-oint airs, but the shared-side solutions get closer and closer as we aroach Position. At Position the line of ersective is tangent to the circle (φ = φ max, in which case the shared-side configuration becomes a single solution, but the shared-oint solutions are still resent as is indicated by the intersection of the line of ersective with λ at c. Therefore, there are solutions at Position. Given the angles φ and, the exact ositions of each triangle in solutions described above can be derived in a straightforward manner from the equations in the Key Theorem and the labels in figures 5 and 6. 8
9 Varying That accounts for all the feasible solutions when CP is at least a distance h from m, as it is in figure 6. If CP gets closer to the triangle, the same tye of sketch shows that there are situations with exactly solution because the "front" member of the shared-side solution is now behind the camera. However, for large enough value of φ, it is ossible to have shared-oint solutions, which in combination with the single shared-side solution roduces exactly solutions. For a systematic look at such ossibilities, we let vary. et d = CPm be the distance from CP to m. From figure 5 we see that d is inversely related to : d = h tan et φ, φ and φ be the elevation angles of Position, Position and Position resectively, in figure 6. Given any value of d, we can comute the values of φ, φ and φ (from figure 6 as functions of : T T φ = arctan( + T T T φ = arctan( T φ = arctan( T = tan T T From these equations we can lot φ, φ and φ as a function of the distance d, as shown in figure 7. If φ, the angle of elevation of the line of ersective, is less than or equal to φ then we will have exactly solutions, the shared-side air. If φ < φ < φ there are exactly 4 solutions, the combined shared-oint and shared-side airs. When φ = φ, CP is directly over the orthocenter. φ cannot be greater than φ (= φ max. As d increases, decreases gradually to 0 while Positions, and in figure 6 move toward the to of the circle (converging to a single oint at the to of the circle when d =, = 0. When d = h, the oint of tangency (Position is coincident with CP and we have φ = π/ and = π /. 9
10 Figures 7 and 8 show the ossible solution sets for d h. As d decreases, Positions,, and move down the circle. When d = h, CP is coincident with the oint of tangency (Position and φ = π/, = π /. When d = h, CP can no longer be directly over the orthocenter because φ = π/. At d = h/ the shared-oint solutions are no longer geometrically ossible (φ = 0, z = 0, = π/. See Chart for a summary of combinations of d, and φ. Figure 8 shows the same grah a figure 7, but this time with labels indicating the exact number of solutions associated with each configuration. Figure 8 omits φ since it simly tells us when CP is directly over the orthocenter, and does not cause a change in the number of solutions. Figure 7: When CP is at a distance h from m, Position (from figure 6 is coincident with CP. At this distance we lose one of the shared-side solutions because the leanforward solution is behind the camera. At d = h the CP can no longer be over the orthocenter (φ = π/. At d = h/ the shared-oint solutions are lying flat (φ = 0, z = 0, = π/, and if d < h/ ( > π/ the shared-oint solutions are no longer ossible. d φ φ φ h π tan ( tan ( π π tan ( π * π 0 * * 0 π * * * h h 0
11 Chart : Values of, φ, φ and φ corresonding to samled values of d between h and 0. The *'s indicate infeasible configurations. For examle, when d = h, φ = π/, but when d < h Position goes behind the camera, so φ is ignored. Figure 8: Same as figure 7, but with φ omitted, and labeled with the number of solutions associated with each region of values for φ and d. When d is between 0 to h/ we get a single solution: the lean-back member of the shared-side air. No other solution is geometrically ossible. When the distance is between h/ and there are two ossibilities. If φ φ, then we again have only one solution, the lean-back member of the shared-side air. If φ > φ then the shared-oint solutions are ossible and together with the lean-back solution make distinct solutions. For d > h, we have two ossibilities. If φ φ we get exactly solutions, the shared side air. If φ > φ we get exactly 4 solutions u until φ = φ, which is where the tangency occurs, so we lose one of the shared-side solutions, resulting in exactly solutions. Fixed Triangle Figure 9 shows the same configurations as figure 6 but this time with the triangle fixed s and CP moving along a circle in the orthogonal lane. This circle has radius r =, tan and it is centered at m. Included in the icture is the location of the "danger cylinder", the right cylinder whose base is the circle that circumscribes the triangle. Positions and corresond to being on the boundary of the danger cylinder, and Position corresonds to being directly over the triangle's orthocenter.
12 Figure 9: The same configurations as figure 6, but with a fixed triangle, and including the location of the danger cylinder. Positions and are on the boundary of the cylinder and Position is directly over the orthocenter. Position is clearly the lace where the triangle's altitude is tangent to the line of ersective. Shared-oint solutions are ossible when CP is between Position and Position, as are the shared-side solutions, for a total of 4 solutions. It is easy to see the symmetry between Positions and since they are located directly above oints at distance h/ from m, which is how we know that Position is on the boundary of the cylinder. If CP is on the same circular trajectory but outside the cylinder, then there are exactly solutions, the front-back air. Conclusions For the secial case of an equilateral triangle, and the center of ersective in one of the orthogonal lanes, we have rovided a thorough geometric solution. Some of the analysis can be directly extended to the case of an isosceles triangle, but then the results would be restricted to the secific orthogonal lane that is orthogonal to the base of the triangle, whereas the results resented here aly to any of the three orthogonal lanes. Another extension of these results would be to consider a small erturbation off the orthogonal lane. If CP were close to, but not on, the orthogonal lane then the sharedside solutions would be eliminated. That is, the shared-side solutions deend heavily on the symmetry created by the orthogonal lane. Once that symmetry is gone, it aears that the only ossible solutions are those that are small adjustments to the share-oint solutions, which would no longer share a oint but would be close to their configurations as in the orthogonal case. This leads us to the conjecture that -solution case is the norm
13 and that the 4-solution case results from geometrically rare (set of measure 0 symmetries. References [] Martin A. Fischler, Robert C. Bolles, "Random samle consensus: a aradigm for model fitting with alications to image analysis and automated cartograhy", Communications of the ACM, Volume 4, Issue 6 (June 98, Pages: 8-95 [] William J. Wolfe, Donald Mathis, Cheryl Weber Sklair, Michael Magee, "The Persective View of Three Points", IEEE Transactions on Pattern Analysis and Machine Intelligence, Volume, Issue (January 99, Pages: [] Robert M. Haralick, Chung-Nan ee, Karsten Ottenberg, Michael Nölle, "Review and analysis of solutions of the three oint ersective ose estimation roblem", International Journal of Comuter Vision, Volume, Issue (December 994, Pages: [4] Xiao-Shan Gao, Xiao-Rong Hou, Jianliang Tang, Hang-Fei Cheng, "Comlete Solution Classification for the Persective-Three-Point Problem", IEEE Transactions on Pattern Analysis and Machine Intelligence, Volume 5, Issue 8 (August 00, Pages: [5] Cai-Xia Zhang, Zhan-Yi Hu, "A general sufficient condition of four ositive solutions of the PP roblem", Journal of Comuter Science and Technology, Volume 0, Issue 6 (November 005, Pages: [6] E. H. Thomson "SPACE RESECTION : FAIURE CASES", The Photogrammetric Record 5 (7, 0 07, 966. [7] ZHANG Cai-Xia, HU Zhan-Yi, "Why is the Danger Cylinder Dangerous in the PP Problem?", Acta Automatica Sinica, 006, (4:
CMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationTHE ERDÖS - MORDELL THEOREM IN THE EXTERIOR DOMAIN
INTERNATIONAL JOURNAL OF GEOMETRY Vol. 5 (2016), No. 1, 31-38 THE ERDÖS - MORDELL THEOREM IN THE EXTERIOR DOMAIN PETER WALKER Abstract. We show that in the Erd½os-Mordell theorem, the art of the region
More information19th Bay Area Mathematical Olympiad. Problems and Solutions. February 28, 2017
th Bay Area Mathematical Olymiad February, 07 Problems and Solutions BAMO- and BAMO- are each 5-question essay-roof exams, for middle- and high-school students, resectively. The roblems in each exam are
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationarxiv:cond-mat/ v2 25 Sep 2002
Energy fluctuations at the multicritical oint in two-dimensional sin glasses arxiv:cond-mat/0207694 v2 25 Se 2002 1. Introduction Hidetoshi Nishimori, Cyril Falvo and Yukiyasu Ozeki Deartment of Physics,
More informationSystem Reliability Estimation and Confidence Regions from Subsystem and Full System Tests
009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract
More informationAn Improved Calibration Method for a Chopped Pyrgeometer
96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum
More informationApplications to stochastic PDE
15 Alications to stochastic PE In this final lecture we resent some alications of the theory develoed in this course to stochastic artial differential equations. We concentrate on two secific examles:
More informationSYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY
SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY FEDERICA PASQUOTTO 1. Descrition of the roosed research 1.1. Introduction. Symlectic structures made their first aearance in
More informationNUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS
NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS Tariq D. Aslam and John B. Bdzil Los Alamos National Laboratory Los Alamos, NM 87545 hone: 1-55-667-1367, fax: 1-55-667-6372
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationOXFORD UNIVERSITY. MATHEMATICS, JOINT SCHOOLS AND COMPUTER SCIENCE WEDNESDAY 4 NOVEMBER 2009 Time allowed: hours
OXFORD UNIVERSITY MATHEMATICS, JOINT SCHOOLS AND COMPUTER SCIENCE WEDNESDAY 4 NOVEMBER 2009 Time allowed: 2 1 2 hours For candidates alying for Mathematics, Mathematics & Statistics, Comuter Science, Mathematics
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More informationMath 121: Calculus 1 - Fall 2012/2013 Review of Precalculus Concepts
Introduction Math : Calculus - Fall 0/0 Review of Precalculus Concets Welcome to Math - Calculus, Fall 0/0! This roblems in this acket are designed to hel you review the toics from Algebra and Precalculus
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More information1. INTRODUCTION. Fn 2 = F j F j+1 (1.1)
CERTAIN CLASSES OF FINITE SUMS THAT INVOLVE GENERALIZED FIBONACCI AND LUCAS NUMBERS The beautiful identity R.S. Melham Deartment of Mathematical Sciences, University of Technology, Sydney PO Box 23, Broadway,
More informationComplex Analysis Homework 1
Comlex Analysis Homework 1 Steve Clanton Sarah Crimi January 27, 2009 Problem Claim. If two integers can be exressed as the sum of two squares, then so can their roduct. Proof. Call the two squares that
More informationA M,ETHOD OF MEASURING THE RESISTIVITY AND HALL' COEFFICIENT ON LAMELLAE OF ARBITRARY SHAPE
'. ' 220 HILlS TECHNICAL REVIEW VOLUME 20 A M,ETHOD OF MEASURING THE RESISTIVITY AND HALL' COEFFICIENT ON LAMELLAE OF ARBITRARY SHAE 621.317.331:538.632.083 Resistivity and Hall-coefficient measurements
More informationStatics and dynamics: some elementary concepts
1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and
More informationCentral Force Motion Challenge Problems
Central Force Motion Challenge Problems Problem 1: Ellitic Orbit A satellite of mass m s is in an ellitical orbit around a lanet of mass m which is located at one focus of the ellise. The satellite has
More informationChapter 7 Rational and Irrational Numbers
Chater 7 Rational and Irrational Numbers In this chater we first review the real line model for numbers, as discussed in Chater 2 of seventh grade, by recalling how the integers and then the rational numbers
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationElliptic Curves and Cryptography
Ellitic Curves and Crytograhy Background in Ellitic Curves We'll now turn to the fascinating theory of ellitic curves. For simlicity, we'll restrict our discussion to ellitic curves over Z, where is a
More informationBrownian Motion and Random Prime Factorization
Brownian Motion and Random Prime Factorization Kendrick Tang June 4, 202 Contents Introduction 2 2 Brownian Motion 2 2. Develoing Brownian Motion.................... 2 2.. Measure Saces and Borel Sigma-Algebras.........
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More informationUse of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek
Use of Transformations and the Reeated Statement in PROC GLM in SAS Ed Stanek Introduction We describe how the Reeated Statement in PROC GLM in SAS transforms the data to rovide tests of hyotheses of interest.
More informationTowards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK
Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)
More informationPrinciples of Computed Tomography (CT)
Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel
More informationThe Fekete Szegő theorem with splitting conditions: Part I
ACTA ARITHMETICA XCIII.2 (2000) The Fekete Szegő theorem with slitting conditions: Part I by Robert Rumely (Athens, GA) A classical theorem of Fekete and Szegő [4] says that if E is a comact set in the
More informationParticipation Factors. However, it does not give the influence of each state on the mode.
Particiation Factors he mode shae, as indicated by the right eigenvector, gives the relative hase of each state in a articular mode. However, it does not give the influence of each state on the mode. We
More informationarxiv:math/ v4 [math.gn] 25 Nov 2006
arxiv:math/0607751v4 [math.gn] 25 Nov 2006 On the uniqueness of the coincidence index on orientable differentiable manifolds P. Christoher Staecker October 12, 2006 Abstract The fixed oint index of toological
More informationMulti-Operation Multi-Machine Scheduling
Multi-Oeration Multi-Machine Scheduling Weizhen Mao he College of William and Mary, Williamsburg VA 3185, USA Abstract. In the multi-oeration scheduling that arises in industrial engineering, each job
More informationA Closed-Form Solution to the Minimum V 2
Celestial Mechanics and Dynamical Astronomy manuscrit No. (will be inserted by the editor) Martín Avendaño Daniele Mortari A Closed-Form Solution to the Minimum V tot Lambert s Problem Received: Month
More informationRobustness of classifiers to uniform l p and Gaussian noise Supplementary material
Robustness of classifiers to uniform l and Gaussian noise Sulementary material Jean-Yves Franceschi Ecole Normale Suérieure de Lyon LIP UMR 5668 Omar Fawzi Ecole Normale Suérieure de Lyon LIP UMR 5668
More informationResearch Article Circle Numbers for Star Discs
International Scholarly Research Network ISRN Geometry Volume 211, Article ID 479262, 16 ages doi:1.542/211/479262 Research Article Circle Numbers for Star Discs W.-D. Richter Institute of Mathematics,
More informationA Note on Guaranteed Sparse Recovery via l 1 -Minimization
A Note on Guaranteed Sarse Recovery via l -Minimization Simon Foucart, Université Pierre et Marie Curie Abstract It is roved that every s-sarse vector x C N can be recovered from the measurement vector
More informationAn Analysis of Reliable Classifiers through ROC Isometrics
An Analysis of Reliable Classifiers through ROC Isometrics Stijn Vanderlooy s.vanderlooy@cs.unimaas.nl Ida G. Srinkhuizen-Kuyer kuyer@cs.unimaas.nl Evgueni N. Smirnov smirnov@cs.unimaas.nl MICC-IKAT, Universiteit
More informationFinding Shortest Hamiltonian Path is in P. Abstract
Finding Shortest Hamiltonian Path is in P Dhananay P. Mehendale Sir Parashurambhau College, Tilak Road, Pune, India bstract The roblem of finding shortest Hamiltonian ath in a eighted comlete grah belongs
More informationStrong Matching of Points with Geometric Shapes
Strong Matching of Points with Geometric Shaes Ahmad Biniaz Anil Maheshwari Michiel Smid School of Comuter Science, Carleton University, Ottawa, Canada December 9, 05 In memory of Ferran Hurtado. Abstract
More informationImplementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis
CST0 191 October, 011, Krabi Imlementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis Chakrit Suvanjumrat and Ekachai Chaichanasiri* Deartment of Mechanical Engineering, Faculty
More information#A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS
#A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS Norbert Hegyvári ELTE TTK, Eötvös University, Institute of Mathematics, Budaest, Hungary hegyvari@elte.hu François Hennecart Université
More informationLower bound solutions for bearing capacity of jointed rock
Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying
More informationPHYS 301 HOMEWORK #9-- SOLUTIONS
PHYS 0 HOMEWORK #9-- SOLUTIONS. We are asked to use Dirichlet' s theorem to determine the value of f (x) as defined below at x = 0, ± /, ± f(x) = 0, - < x
More informationNotes on Instrumental Variables Methods
Notes on Instrumental Variables Methods Michele Pellizzari IGIER-Bocconi, IZA and frdb 1 The Instrumental Variable Estimator Instrumental variable estimation is the classical solution to the roblem of
More information7. The olynomial P (x) =x, 9x 4 +9x, 4x +4x, 4 can be written in the form P (x) =(x, ) n Q(x), where n is a ositive integer and Q(x) is not divisible
. Which of the following intervals contains all of the real zeros of the function g(t) =t, t + t? (a) (,; ) (b) ( ; ) (c) (,; ) (d) (; ) (e) (,; ) 4 4. The exression ((a), + b), is equivalent to which
More informationUniform Law on the Unit Sphere of a Banach Space
Uniform Law on the Unit Shere of a Banach Sace by Bernard Beauzamy Société de Calcul Mathématique SA Faubourg Saint Honoré 75008 Paris France Setember 008 Abstract We investigate the construction of a
More informationANALYTIC NUMBER THEORY AND DIRICHLET S THEOREM
ANALYTIC NUMBER THEORY AND DIRICHLET S THEOREM JOHN BINDER Abstract. In this aer, we rove Dirichlet s theorem that, given any air h, k with h, k) =, there are infinitely many rime numbers congruent to
More informationwhere x i is the ith coordinate of x R N. 1. Show that the following upper bound holds for the growth function of H:
Mehryar Mohri Foundations of Machine Learning Courant Institute of Mathematical Sciences Homework assignment 2 October 25, 2017 Due: November 08, 2017 A. Growth function Growth function of stum functions.
More informationε i (E j )=δj i = 0, if i j, form a basis for V, called the dual basis to (E i ). Therefore, dim V =dim V.
Covectors Definition. Let V be a finite-dimensional vector sace. A covector on V is real-valued linear functional on V, that is, a linear ma ω : V R. The sace of all covectors on V is itself a real vector
More informationRANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES
RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES AARON ZWIEBACH Abstract. In this aer we will analyze research that has been recently done in the field of discrete
More informationMultiplicative group law on the folium of Descartes
Multilicative grou law on the folium of Descartes Steluţa Pricoie and Constantin Udrişte Abstract. The folium of Descartes is still studied and understood today. Not only did it rovide for the roof of
More informationResearch Article An iterative Algorithm for Hemicontractive Mappings in Banach Spaces
Abstract and Alied Analysis Volume 2012, Article ID 264103, 11 ages doi:10.1155/2012/264103 Research Article An iterative Algorithm for Hemicontractive Maings in Banach Saces Youli Yu, 1 Zhitao Wu, 2 and
More informationINJECTIVITY RADIUS AND DIAMETER OF THE MANIFOLDS OF FLAGS IN THE PROJECTIVE PLANES
INJECTIVITY RADIUS AND DIAMETER OF THE MANIFOLDS OF FLAGS IN THE PROJECTIVE PLANES THOMAS PÜTTMANN Abstract. The manifolds of flags in the rojective lanes RP, CP, HP, and OP are among the very few comact
More informationarxiv: v2 [math.ac] 5 Jan 2018
Random Monomial Ideals Jesús A. De Loera, Sonja Petrović, Lily Silverstein, Desina Stasi, Dane Wilburne arxiv:70.070v [math.ac] Jan 8 Abstract: Insired by the study of random grahs and simlicial comlexes,
More informationTopic 7: Using identity types
Toic 7: Using identity tyes June 10, 2014 Now we would like to learn how to use identity tyes and how to do some actual mathematics with them. By now we have essentially introduced all inference rules
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationarxiv: v1 [physics.data-an] 26 Oct 2012
Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch
More informationSolved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.
Solved Problems Solved Problems P Solve the three simle classification roblems shown in Figure P by drawing a decision boundary Find weight and bias values that result in single-neuron ercetrons with the
More informationLilian Markenzon 1, Nair Maria Maia de Abreu 2* and Luciana Lee 3
Pesquisa Oeracional (2013) 33(1): 123-132 2013 Brazilian Oerations Research Society Printed version ISSN 0101-7438 / Online version ISSN 1678-5142 www.scielo.br/oe SOME RESULTS ABOUT THE CONNECTIVITY OF
More informationSingular Frégier Conics in Non-Euclidean Geometry
Singular Frégier onics in on-euclidean Geometry Hans-Peter Schröcker University o Innsbruck, Austria arxiv:1605.07437v1 [math.mg] 24 May 2016 May 25, 2016 The hyotenuses o all right triangles inscribed
More informationThe Graph Accessibility Problem and the Universality of the Collision CRCW Conflict Resolution Rule
The Grah Accessibility Problem and the Universality of the Collision CRCW Conflict Resolution Rule STEFAN D. BRUDA Deartment of Comuter Science Bisho s University Lennoxville, Quebec J1M 1Z7 CANADA bruda@cs.ubishos.ca
More informationOn the stability and integration of Hamilton-Poisson systems on so(3)
Rendiconti di Matematica, Serie VII Volume 37, Roma (016), 1 4 On the stability and integration of Hamilton-Poisson systems on so(3) R. M. ADAMS R. BIGGS W. HOLDERBAUM C. C. REMSING Abstract: We consider
More informationA Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition
A Qualitative Event-based Aroach to Multile Fault Diagnosis in Continuous Systems using Structural Model Decomosition Matthew J. Daigle a,,, Anibal Bregon b,, Xenofon Koutsoukos c, Gautam Biswas c, Belarmino
More informationOn the Toppling of a Sand Pile
Discrete Mathematics and Theoretical Comuter Science Proceedings AA (DM-CCG), 2001, 275 286 On the Toling of a Sand Pile Jean-Christohe Novelli 1 and Dominique Rossin 2 1 CNRS, LIFL, Bâtiment M3, Université
More informationNetwork Configuration Control Via Connectivity Graph Processes
Network Configuration Control Via Connectivity Grah Processes Abubakr Muhammad Deartment of Electrical and Systems Engineering University of Pennsylvania Philadelhia, PA 90 abubakr@seas.uenn.edu Magnus
More informationAP Calculus Testbank (Chapter 10) (Mr. Surowski)
AP Calculus Testbank (Chater 1) (Mr. Surowski) Part I. Multile-Choice Questions 1. The grah in the xy-lane reresented by x = 3 sin t and y = cost is (A) a circle (B) an ellise (C) a hyerbola (D) a arabola
More informationB8.1 Martingales Through Measure Theory. Concept of independence
B8.1 Martingales Through Measure Theory Concet of indeendence Motivated by the notion of indeendent events in relims robability, we have generalized the concet of indeendence to families of σ-algebras.
More informationOn the relationship between sound intensity and wave impedance
Buenos Aires 5 to 9 Setember, 16 Acoustics for the 1 st Century PROCEEDINGS of the nd International Congress on Acoustics Sound Intensity and Inverse Methods in Acoustics: Paer ICA16-198 On the relationshi
More informationPrincipal Components Analysis and Unsupervised Hebbian Learning
Princial Comonents Analysis and Unsuervised Hebbian Learning Robert Jacobs Deartment of Brain & Cognitive Sciences University of Rochester Rochester, NY 1467, USA August 8, 008 Reference: Much of the material
More informationdn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential
Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system
More informationLINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL
LINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL Mohammad Bozorg Deatment of Mechanical Engineering University of Yazd P. O. Box 89195-741 Yazd Iran Fax: +98-351-750110
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationAn Inverse Problem for Two Spectra of Complex Finite Jacobi Matrices
Coyright 202 Tech Science Press CMES, vol.86, no.4,.30-39, 202 An Inverse Problem for Two Sectra of Comlex Finite Jacobi Matrices Gusein Sh. Guseinov Abstract: This aer deals with the inverse sectral roblem
More informationSession 5: Review of Classical Astrodynamics
Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve
More informationReal Analysis 1 Fall Homework 3. a n.
eal Analysis Fall 06 Homework 3. Let and consider the measure sace N, P, µ, where µ is counting measure. That is, if N, then µ equals the number of elements in if is finite; µ = otherwise. One usually
More informationPreliminary Round Question Booklet
First Annual Pi Day Mathematics Cometition Question Booklet 016 goes on and on, and e is just as cursed. Iwonder,howdoes begin When its digits are reversed? -MartinGardner Pi Day Mathematics Cometition
More informationQuantitative estimates of propagation of chaos for stochastic systems with W 1, kernels
oname manuscrit o. will be inserted by the editor) Quantitative estimates of roagation of chaos for stochastic systems with W, kernels Pierre-Emmanuel Jabin Zhenfu Wang Received: date / Acceted: date Abstract
More informationImprovement on the Decay of Crossing Numbers
Grahs and Combinatorics 2013) 29:365 371 DOI 10.1007/s00373-012-1137-3 ORIGINAL PAPER Imrovement on the Decay of Crossing Numbers Jakub Černý Jan Kynčl Géza Tóth Received: 24 Aril 2007 / Revised: 1 November
More informationNumerical Linear Algebra
Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and
More informationWolfgang POESSNECKER and Ulrich GROSS*
Proceedings of the Asian Thermohysical Proerties onference -4 August, 007, Fukuoka, Jaan Paer No. 0 A QUASI-STEADY YLINDER METHOD FOR THE SIMULTANEOUS DETERMINATION OF HEAT APAITY, THERMAL ONDUTIVITY AND
More informationarxiv: v2 [math.na] 6 Apr 2016
Existence and otimality of strong stability reserving linear multiste methods: a duality-based aroach arxiv:504.03930v [math.na] 6 Ar 06 Adrián Németh January 9, 08 Abstract David I. Ketcheson We rove
More informationMTH234 Chapter 15 - Multiple Integrals Michigan State University
MTH24 Chater 15 - Multile Integrals Michigan State University 6 Surface Area Just as arc length is an alication of a single integral, surface area is an alication of double integrals. In 15.6 we comute
More informationCERIAS Tech Report The period of the Bell numbers modulo a prime by Peter Montgomery, Sangil Nahm, Samuel Wagstaff Jr Center for Education
CERIAS Tech Reort 2010-01 The eriod of the Bell numbers modulo a rime by Peter Montgomery, Sangil Nahm, Samuel Wagstaff Jr Center for Education and Research Information Assurance and Security Purdue University,
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More information16.2. Infinite Series. Introduction. Prerequisites. Learning Outcomes
Infinite Series 6.2 Introduction We extend the concet of a finite series, met in Section 6., to the situation in which the number of terms increase without bound. We define what is meant by an infinite
More informationSolutions to Problem Set 5
Solutions to Problem Set Problem 4.6. f () ( )( 4) For this simle rational function, we see immediately that the function has simle oles at the two indicated oints, and so we use equation (6.) to nd the
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell February 10, 2010 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationHENSEL S LEMMA KEITH CONRAD
HENSEL S LEMMA KEITH CONRAD 1. Introduction In the -adic integers, congruences are aroximations: for a and b in Z, a b mod n is the same as a b 1/ n. Turning information modulo one ower of into similar
More informationOn Wald-Type Optimal Stopping for Brownian Motion
J Al Probab Vol 34, No 1, 1997, (66-73) Prerint Ser No 1, 1994, Math Inst Aarhus On Wald-Tye Otimal Stoing for Brownian Motion S RAVRSN and PSKIR The solution is resented to all otimal stoing roblems of
More informationCONVOLVED SUBSAMPLING ESTIMATION WITH APPLICATIONS TO BLOCK BOOTSTRAP
Submitted to the Annals of Statistics arxiv: arxiv:1706.07237 CONVOLVED SUBSAMPLING ESTIMATION WITH APPLICATIONS TO BLOCK BOOTSTRAP By Johannes Tewes, Dimitris N. Politis and Daniel J. Nordman Ruhr-Universität
More informationMinimal Surfaces in R 3
Minimal Suaces in R 3 AJ Vargas May 10, 018 1 Introduction The goal of this note is to give a brief introduction to the theory of minimal suaces in R 3, and to show how one would go about generalizing
More informationLecture 1.2 Pose in 2D and 3D. Thomas Opsahl
Lecture 1.2 Pose in 2D and 3D Thomas Osahl Motivation For the inhole camera, the corresondence between observed 3D oints in the world and 2D oints in the catured image is given by straight lines through
More informationQUIZ ON CHAPTER 4 - SOLUTIONS APPLICATIONS OF DERIVATIVES; MATH 150 FALL 2016 KUNIYUKI 105 POINTS TOTAL, BUT 100 POINTS = 100%
QUIZ ON CHAPTER - SOLUTIONS APPLICATIONS OF DERIVATIVES; MATH 150 FALL 016 KUNIYUKI 105 POINTS TOTAL, BUT 100 POINTS = 100% = x + 5 1) Consider f x and the grah of y = f x in the usual xy-lane in 16 x
More informationSpin as Dynamic Variable or Why Parity is Broken
Sin as Dynamic Variable or Why Parity is Broken G. N. Golub golubgn@meta.ua There suggested a modification of the Dirac electron theory, eliminating its mathematical incomleteness. The modified Dirac electron,
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationNode-voltage method using virtual current sources technique for special cases
Node-oltage method using irtual current sources technique for secial cases George E. Chatzarakis and Marina D. Tortoreli Electrical and Electronics Engineering Deartments, School of Pedagogical and Technological
More informationThe inverse Goldbach problem
1 The inverse Goldbach roblem by Christian Elsholtz Submission Setember 7, 2000 (this version includes galley corrections). Aeared in Mathematika 2001. Abstract We imrove the uer and lower bounds of the
More informationTHE 2D CASE OF THE BOURGAIN-DEMETER-GUTH ARGUMENT
THE 2D CASE OF THE BOURGAIN-DEMETER-GUTH ARGUMENT ZANE LI Let e(z) := e 2πiz and for g : [0, ] C and J [0, ], define the extension oerator E J g(x) := g(t)e(tx + t 2 x 2 ) dt. J For a ositive weight ν
More informationElementary Analysis in Q p
Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some
More information8 STOCHASTIC PROCESSES
8 STOCHASTIC PROCESSES The word stochastic is derived from the Greek στoχαστικoς, meaning to aim at a target. Stochastic rocesses involve state which changes in a random way. A Markov rocess is a articular
More information