The Bernstein polynomial basis: a centennial retrospective
|
|
- Muriel Foster
- 6 years ago
- Views:
Transcription
1 The Bernstein polynomial basis: a centennial retrospective a sociological study in the evolution of mathematical ideas Rida T. Farouki Department of Mechanical & Aerospace Engineering, University of California, Davis
2 synopsis 1912: Sergei Natanovich Bernstein constructive proof of Weierstrass approximation theorem 1960s: Paul de Faget de Casteljau, Pierre Étienne Bézier and the origins of computer-aided geometric design elucidation of Bernstein basis properties and algorithms 1980s: intrinsic numerical stability of the Bernstein form algorithms & representations for computer-aided design diversification of applications in scientific computing
3 Weierstrass approximation theorem Given any continuous function f(x) on an interval [ a, b ] and a tolerance ɛ > 0, a polynomial p n (x) of sufficiently high degree n exists, such that f(x) p n (x) ɛ for x [ a, b ]. Polynomials can uniformly approximate any continuous f(x), x [ a, b ]. Original (1885) proof by Weierstrass is existential in nature begins by expressing f(x) as a convolution 1 + [ f(x) = lim f(t) exp k 0 πk (t x)2 k 2 ] dt with a Dirac delta function, and relies heavily on analytic limit arguments.
4 Sergei Natanovich Bernstein ( ) (photo: Russian Academy of Sciences)
5 academic career of S. N. Bernstein 1904: Sorbonne PhD thesis, on analytic nature of PDE solutions (worked with Hilbert at Göttingen during academic year) 1913: Kharkov PhD thesis (polynomial approximation of functions) 1912: Comm. Math. Soc. Kharkov paper (2 pages): constructive proof of Weierstrass theorem introduction of Bernstein basis : Professor in Kharkov & Director of Mathematical Institute political purge: moved to USSR Academy of Sciences (Leningrad) : escapes to Kazakhstan during the siege of Leningrad : Steklov Math. Institute, Russian Academy of Sciences, Moscow edited complete works of Chebyshev (died 1968) collected works of Bernstein published in 4 volumes,
6 Bernstein s proof of Weierstrass theorem Russian school of approximation theory, founded by Chebyshev, favors constructive approximation methods over existential proofs given f(t) continuous on t [ 0, 1 ] define p n (t) = f(k/n) b n k(t), b n k(t) = k=0 ( ) n (1 t) n k t k k p n (t) = convex combination of sampled values f(0), f( 1 n ),..., f(1) f(t) p n (t) = O( 1 n ) for t [ 0, 1 ] = p n (t) converges uniformly to f(t) as n derivatives of p n (t) also converge to those of f(t) as n
7 S. N. Bernstein, Comm. Kharkov Math. Soc. (1912)
8 connection with probability theory basis function b n k (t) = ( n k ) (1 t) n k t k probability of k successes in n trials of random process with individual probability of success t in each trial binomial probability distribution non-negativity & partition-of-unity properties of b n k (t)
9 slow convergence of Bernstein approximations Bernstein polynomial approximations of degree n = 10, 30, 100, 300, 1000 to a triangular wave This fact seems to have precluded any numerical application of Bernstein polynomials from having been made. Perhaps they will find application when the properties of the approximant in the large are of more importance than the closeness of the approximation. Philip J. Davis, Interpolation and Approximation (1963)
10 Paul de Casteljau & Pierre Bézier (1960s) an emerging application computer-aided design Paul de Faget de Casteljau theory of courbes et surfaces à pôles developed at André Citroën, SA in the early 1960s de Casteljau s work unpublised (regarded as proprietary by Citroën) revealed to outside world by Wolfgang Böhm in mid-1980s Pierre Étienne Bézier implemented methods for computer-aided design and manufacturing at Renault during 1960s and 1970s Bézier published numerous articles and books describing his ideas basic problem: provide intuitive & interactive means for design and manipulation of free form curves and surfaces by computer, in the automotive, aerospace, and related industries identification of de Casteljau s and Bézier s ideas with Bernstein form of polynomials came later, through work of Forrest, Riesenfeld, et al.
11 de Casteljau Courbes et surfaces à pôles (Société Anonyme André Citroën, 1963) pôles = pilot points (interpolation of polynomials with polar forms)
12 de Casteljau barycentric coordinates de Casteljau s λ and µ = interval barycentric coordinates, with λ + µ = 1 example for t [ 0, 1 ] take λ = 1 t and µ = t, and expand (λ + µ) n 1 = [ (1 t) + t ] n = k=0 ( ) n (1 t) n k t k = k b n k(t) k=0 Bernstein basis {b n k (t)} is non negative and forms a partition of unity de Casteljau also considers extension to barycentric coordinates and multivariate polynomial bases on triangular and simplex domains
13 computer-aided design in the early 60s... the designers were astonished and scandalized. Was it some kind of joke? It was considered nonsense to represent a car body mathematically. It was enough to please the eye, the word accuracy had no meaning... reaction at Citroën to de Casteljau s ideas Citroën s first attempts at digital shape representation used a Burroughs E101 computer featuring 128 program steps, a 220-word memory, and a 5 kw power consumption! De Casteljau s insane persistence led to an increased adoption of computer-aided design methods in Citroën from 1963 onwards. My stay at Citroën was not only an adventure for me, but also an adventure for Citroën! P. de Casteljau
14 correspondence with de Casteljau (1991)
15 Bézier s point-vector form of a polynomial curve specify degree-n curve by initial point p 0 and n vectors a 1,..., a n r(t) = p 0 + k=1 a k f n k (t), f n k (t) = ( 1)k (k 1)! dk 1 (1 t) n 1 tk dt k 1 t a 2 a r(t) f 3 1(t) f 3 2(t) a f 3 3(t) 0.2 p Left: Bézier point-vector specification of a cubic curve. Right: cubic basis functions f 3 1 (t), f 3 2 (t), f 3 3 (t) associated with the vectors a 1 a 2, a 3. mischievous Bézier f n 1 (t),..., f n n (t) = basis of Onésime Durand!
16 control-point form of a Bézier curve n Forrest (1972) : fi n(t) = b n k(t), b n k(t) = k=i ( ) n (1 t) n k t k k p 1 p 2 p3 p 0 re-write as r(t) = p k b n k(t), k=0 p k = p k 1 + a k manipulate cure shape by moving control points p 0,..., p n
17 convex-hull, variation-diminishing, degree-elevation properties of the Bézier form r(t) = p k b n k(t), control points p 0,..., p n k=0
18 de Casteljau algorithm evaluates & subdivides r(t) initialize set t = τ and p 0 k = p k for k = 0,..., n for r = 1,..., n for j = r,..., n { p r j = (1 τ) pr 1 j 1 + τ pr 1 j } generates a triangular array of points { p r j } p 0 0 p 0 1 p 0 2 p 0 n p 1 1 p 1 2 p 1 n p 2 2 p 2 n p n n
19 p n n = evaluated point r(τ) on curve p 0 0, p 1 1,..., p n 1 n 1, pn n = control points for subsegment t [ 0, τ ] of r(t) p n n, p n 1 n,..., p 1 n, p 0 n = control points for subsegment t [ τ, 1 ] of r(t) p 5 p 0 p 4 p 1 r(τ) p 3 p 2
20 interlude... lost in translation warning sign on bathroom door in Beijing hotel
21 English on vacation in a Bucharest hotel lobby The elevator is being fixed for the next day. During that time we regret that you will be unbearable. in a Paris hotel elevator Please leave your values at the front desk. in a Zurich hotel Because of the impropriety of entertaining guests of the opposite sex in your bedroom, it is suggested that the lobby be used for this purpose. in an Acapulco restaurant The manager has personally passed all the water served here. in Germany s Schwarzwald It is strictly forbidden on our Black Forest camping site that people of different sex for instance, men and women live together in one tent unless they are married with each other for that purpose.
22 in an Athens hotel Guests are expected to complain at the office between 9 and 11 am daily. instructions for AC in Japanese hotel If you want just condition of warm in your room, please control yourself. in a Yugoslav hotel The flattening of underwear with pleasure is the job of the chambermaid. in a Japanese hotel You are invited to take advantage of the chambermaid. on the menu of a Swiss restaurant Our wines leave you nothing to hope for. in a Bangkok dry cleaners Drop your trousers here for best results. Japanese rental car instructions When passenger of foot heave in sight, tootle the horn. Trumpet him melodiously at first, but if he still obstacles your passage, then tootle him with vigor.
23 Bernstein basis functions b n k(t) t b n k (t) = ( n k ) (1 t) n k t k roots of multiplicity k and n k at t = 0 and t = 1
24 properties of the Bernstein basis b n k (t) = ( n k ) (1 t) n k t k, k = 0,..., n partition of unity : k=0 b n k(t) 1 non-negativity : b n k (t) 0 for t [ 0, 1 ] symmetry : b n k (t) = bn n k (1 t) recursion : b n+1 k (t) = t b n k 1 (t) + (1 t) bn k (t) unimodality : b n k (t) has maximum at t = k/n
25 properties of Bernstein form, p(t) = k=0 c k b n k(t) end-point values : p(0) = c 0 and p(1) = c n lower & upper bounds : min k c k p(t) max k c k variation diminishing : # roots = signvar(c 0,..., c n ) 2m derivatives & integrals : coefficients of p (t) & p(t) dt = differences & partial sums of c 0,..., c n recursive algorithms for subdivision, degree elevation, arithmetic operations, composition, resultants, etc. root isolation (subdivision & variation-diminishing property)
26 the plague of numerical instability... or, the temptation to kick the computer Do you ever want to kick the computer? Does it iterate endlessly on your newest algorithm that should have converged in three iterations? And does it finally come to a crashing halt with the insulting message that you divided by zero? These minor trauma are, in fact, the ways the computer manages to kick you and, unfortunately, you almost always deserve it! For it is a sad fact that most of us can more readily compute than think...
27 numerical stability of polynomials p(t) has coefficients c 0,..., c n in basis Φ = {φ 0 (t),..., φ n (t)} p(t) = c k φ k (t) k=0 how sensitive is a value or root of p(t) to perturbations of maximum relative magnitude ɛ in the coefficients c 0,..., c n? condition number for value of p(t) : δp(t) C Φ (p(t)) ɛ, C Φ (p(t)) = c k φ k (t) k=0 condition number for root τ of p(t) : δτ C Φ (τ) ɛ, C Φ (τ) = 1 p (τ) c k φ k (t) k=0
28 condition numbers for power and Bernstein forms c j = p(t) = a k t k = c k b n k(t) k=0 k=0 ( j j ( j ( k) n a k, t k) k = ( k) n b k) n j (t) k=0 j=k Theorem. C B (p(t)) C P (p(t)) for any polynomial p(t) and all t [ 0, 1 ]. Proof (triangle inequality). C B (p(t)) = c j b n j j (t) = j=0 j=0 k=0 ( j a k ( k) n b k) n j (t) = k=0 j=k k=0 ( j ( k) n a k b k) n j (t) a k t k = C P (p(t)).
29 Wilkinson s perfidious polynomial problem: compute the roots of the degree 20 polynomial p(t) = (t 1)(t 2) (t 20) = 20 k=0 a k t k using (software) floating-point arithmetic J. H. Wilkinson (1959), The evaluation of the zeros of ill-conditioned polynomials, Parts I & II, Numerische Mathematik 1, & The cosy relationship that mathematicians enjoyed with polynomials suffered a severe setback in the early fifties when electronic computers came into general use. Speaking for myself, I regard it as the most traumatic experience in my career as a numerical analyst. J. H. Wilkinson, The Perfidious Polynomial, in Studies in Numerical Analysis (1984)
30 root condition numbers for Wilkinson polynomial root power basis Bernstein basis
31 perturbed roots of Wilkinson polynomial ɛ = root power basis Bernstein basis ± i ± i ± i ± i ± i
32 evaluating Wilkinson s t = a 0 = a 1 t = a 2 t 2 = a 3 t 3 = a 4 t 4 = a 5 t 5 = a 6 t 6 = a 7 t 7 = a 8 t 8 = a 9 t 9 = a 10 t 10 = a 11 t 11 = a 12 t 12 = a 13 t 13 = a 14 t 14 = a 15 t 15 = a 16 t 16 = a 17 t 17 = a 18 t 18 = a 19 t 19 = a 20 t 20 = p(t) =
33 perturbation regions for p(t) = (t 1 6 ) (t 1) ε = ε = perturbed Bernstein form perturbed power form
34 optimal stability of Bernstein basis Ψ = {ψ 0 (t),..., ψ n (t)} and Φ = {φ 0 (t),..., φ n (t)} non negative on [ a, b ] Theorem. If ψ j (t) = M jk φ k (t) with M jk 0, k=0 then the condition numbers for the value of any degree n polynomial p(t) at any point t [ a, b ] in the bases Φ and Ψ satisfy C Φ (p(t)) C Ψ (p(t)). We say that the Φ basis is systematically more stable than the Ψ basis. Example: Φ = {b n 0(t),..., b n n(t)} and Ψ = {1, t,..., t n } in fact, the Bernstein basis is optimally stable (it is impossible to construct a basis on [ 0, 1 ] that is systematically more stable).
35 optimal stability (sketch) P n = set of all non negative bases for degree-n polynomials on [ a, b ]. For Φ, Ψ P n we write Φ Ψ if Ψ = M Φ for a non negative matrix M. The relation is a partial ordering of the set of non-negative bases P n. Theorem. Φ Ψ C Φ (p(t)) C Ψ (p(t)) for all p(t) P n and t [ a, b ]. Definition. Φ is a minimal basis in P n if no Ψ exists, such that Ψ Φ. A minimal basis in P n is optimally stable it is impossible to construct a non negative basis on [ a, b ] that is systematically more stable. Theorem. The Bernstein basis is minimal in P n, and is optimally stable. It is the only minimal basis whose basis functions have no roots in (a, b).
36 condition numbers can be very large!
37 least-squares polynomial approximation minimize 1 0 orthogonal basis [ f(t) p n (t) ] 2 dt, p n (t) = 1 0 φ j (t) φ k (t) dt = { = a k = 1 1 f(t) φ k (t) dt β k 0 β k a k φ k (t) k=0 j = k 0 j k permanence of coefficients: a 0,..., a n unchanged when n n + 1 orthogonality impossible for non negative bases, but Bernstein basis is intimately related to Legendre basis
38 Legendre and Bernstein bases on t [ 0, 1 ] recurrence relation L 0 (t) = 1, L 1 (t) = 2t 1 (k + 1)L k+1 (t) = (2k + 1)(2t 1) L k (t) k L k 1 (t) Rodrigues formula Bernstein form L k (t) = L k (t) = ( 1)k k! d k [ (1 t)t ]k dtk k ( ) k ( 1) k+i b k i (t) i i=0 Pascal s triangle with alternating signs!!
39 Bernstein form of the Legendre polynomials L 0 (t) = 1 b 0 0(t), L 1 (t) = 1 b 1 0(t) + 1 b 1 1(t), L 2 (t) = 1 b 2 0(t) 2 b 2 1(t) + 1 b 2 2(t), L 3 (t) = 1 b 3 0(t) + 3 b 3 1(t) 3 b 3 2(t) + 1 b 3 3(t), L 4 (t) = 1 b 4 0(t) 4 b 4 1(t) + 6 b 4 2(t) 4 b 4 3(t) + 1 b 4 4(t), L 5 (t) = 1 b 5 0(t) + 5 b 5 1(t) 10 b 5 2(t) + 10 b 5 3(t) 5 b 5 4(t) + 1 b 5 5(t), Bernstein form of Legendre polynomial derivatives e.g., L 4 (t) L 4 (t) = 1 b 4 0(t) 4 b 4 1(t) + 6 b 4 2(t) 4 b 4 3(t) + 1 b 4 4(t), 1 2 L 4(t) = 5 b 3 0(t) 10 b 3 1(t) + 10 b 3 2(t) 5 b 3 3(t), 1 3 L 4(t) = 15 b 2 0(t) 20 b 2 1(t) + 15 b 2 2(t), 1 3 L 4 (t) = 35 b 1 0(t) 35 b 1 1(t), 1 5 L 4 (t) = 70 b 0 0(t),
40 Legendre Bernstein basis transformations p(t) = c j = M jk = ( 1 ) n k a k L k (t) = k=0 M jk a k, a j = k=0 min(j,k) i=max(0,j+k n) M 1 jk = 2j + 1 n + j + 1 c k b n k(t) k=0 k=0 ( 1) k+i ( j i ( ) j n ( 1) j+i k i=0 M 1 jk c k )( k i )( ) n k j i ( )( ) j j i i ( ) n + j k + i condition number C p (M) = M p M 1 p, C 1 (M) = 2 n > C (M)
41 condition numbers for basis transformations matrix condition number Bernstein Legendre Bernstein power Bernstein Hermite Hermite power polynomial degree
42 extension to rational forms rational Bézier curve r(t) = w k p k b n k(t) k=0 w k b n k(t) k=0 defined by control points p 0,..., p n and scalar weights w 0,..., w n set of rational curves is closed under projective transformations conic segments as rational quadratic Bézier curves (w 0 = w 2 = 1) w 1 < 1 (ellipse) w 1 = 1 (parabola) w 1 > 1 (hyperbola)
43 bivariate & multivariate generalizations p p T 2 T 3 T 1 p p 1 p 2 + p barycentric coordinates: (u, v, w) = (area(t 1), area(t 2 ), area(t 3 )) area(t ) 1 = (u + v + w) n = b n ijk(u, v, w), b n ijk(u, v, w) = n! i!j!k! ui v j w k i+j+k=n triangular surface patch: r(u, v, w) = i+j+k=n p ijk b n ijk(u, v, w)
44 bivariate de Casteljau algorithm generates tetrahedral array evaluates and subdivides r(u, v, w) r = r = r = r =
45 generalization to B-spline basis continuous domain [ 0, 1 ] partitioned domain [ t 0, t 1,..., t N+n+1 ] B 3 0(t) B 3 1(t) B 3 2(t) B 3 3(t) B 3 N 2(t) B 3 N 1(t) B 3 N(t) t 0 =t 1 =t 2 =t 3 t 4 t 5 t N t N+1 =t N+2 =t N+3 =t N+4 retain partition-of-unity, non-negativity, recursion properties + compact support & control over continuity (multiple knots) local shape modification C 2 embedded linear segment
46 scientific computing applications real solutions of systems of algebraic equations; identifying extrema or bounds on constrained or unconstrained polynomial functions in one or several variables (optimization) using Bernstein basis properties robust stability of dynamic systems with uncertain physical parameters (Kharitonov generalization of Routh-Hurwitz criterion) definition of barycentric coordinates and partition-of-unity polynomial basis functions over general polygon or polytope domains for use in the finite-element and meshless analysis methods modelling of inter-molecular potential energy surfaces; design of filters for signal processing applications; inputs to neurofuzzy networks modelling non-linear dynamical systems; reconstruction of 3D models and calibration of optical range sensors
47 closure 100 years have elapsed since introduction of Bernstein basis Bernstein form was limited to theory, rather than practice, of polynomial approximation for 50 years after its introduction applications in design, rather than approximation, pioneered 50 years ago by de Casteljau and Bézier now universally adopted as a fundamental representation for computer-aided geometric design applications optimally stable basis for polynomials defined over finite domains Bernstein basis intimately related to Legendre orthogonal basis increasing adoption in diverse scientific computing applications In theory, there is no difference between theory and practice. In practice, there is.... Yogi Berra
Interpolation and polynomial approximation Interpolation
Outline Interpolation and polynomial approximation Interpolation Lagrange Cubic Splines Approximation B-Splines 1 Outline Approximation B-Splines We still focus on curves for the moment. 2 3 Pierre Bézier
More informationCurves, Surfaces and Segments, Patches
Curves, Surfaces and Segments, atches The University of Texas at Austin Conics: Curves and Quadrics: Surfaces Implicit form arametric form Rational Bézier Forms and Join Continuity Recursive Subdivision
More informationM2R IVR, October 12th Mathematical tools 1 - Session 2
Mathematical tools 1 Session 2 Franck HÉTROY M2R IVR, October 12th 2006 First session reminder Basic definitions Motivation: interpolate or approximate an ordered list of 2D points P i n Definition: spline
More informationThe Essentials of CAGD
The Essentials of CAGD Chapter 4: Bézier Curves: Cubic and Beyond Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/essentials-cagd c 2000
More information1.1. The analytical denition. Denition. The Bernstein polynomials of degree n are dened analytically:
DEGREE REDUCTION OF BÉZIER CURVES DAVE MORGAN Abstract. This paper opens with a description of Bézier curves. Then, techniques for the degree reduction of Bézier curves, along with a discussion of error
More informationBernstein polynomials of degree N are defined by
SEC. 5.5 BÉZIER CURVES 309 5.5 Bézier Curves Pierre Bézier at Renault and Paul de Casteljau at Citroën independently developed the Bézier curve for CAD/CAM operations, in the 1970s. These parametrically
More informationApproximation of Circular Arcs by Parametric Polynomials
Approximation of Circular Arcs by Parametric Polynomials Emil Žagar Lecture on Geometric Modelling at Charles University in Prague December 6th 2017 1 / 44 Outline Introduction Standard Reprezentations
More informationSubdivision Matrices and Iterated Function Systems for Parametric Interval Bezier Curves
International Journal of Video&Image Processing Network Security IJVIPNS-IJENS Vol:7 No:0 Subdivision Matrices Iterated Function Systems for Parametric Interval Bezier Curves O. Ismail, Senior Member,
More informationComputer Aided Design. B-Splines
1 Three useful references : R. Bartels, J.C. Beatty, B. A. Barsky, An introduction to Splines for use in Computer Graphics and Geometric Modeling, Morgan Kaufmann Publications,1987 JC.Léon, Modélisation
More informationVisualizing Bezier s curves: some applications of Dynamic System Geogebra
Visualizing Bezier s curves: some applications of Dynamic System Geogebra Francisco Regis Vieira Alves Instituto Federal de Educação, Ciência e Tecnologia do Estado do Ceará IFCE. Brazil fregis@ifce.edu.br
More informationComputergrafik. Matthias Zwicker Universität Bern Herbst 2016
Computergrafik Matthias Zwicker Universität Bern Herbst 2016 2 Today Curves Introduction Polynomial curves Bézier curves Drawing Bézier curves Piecewise curves Modeling Creating 3D objects How to construct
More informationCMSC427 Parametric curves: Hermite, Catmull-Rom, Bezier
CMSC427 Parametric curves: Hermite, Catmull-Rom, Bezier Modeling Creating 3D objects How to construct complicated surfaces? Goal Specify objects with few control points Resulting object should be visually
More informationAn O(h 2n ) Hermite approximation for conic sections
An O(h 2n ) Hermite approximation for conic sections Michael Floater SINTEF P.O. Box 124, Blindern 0314 Oslo, NORWAY November 1994, Revised March 1996 Abstract. Given a segment of a conic section in the
More informationMA 323 Geometric Modelling Course Notes: Day 12 de Casteljau s Algorithm and Subdivision
MA 323 Geometric Modelling Course Notes: Day 12 de Casteljau s Algorithm and Subdivision David L. Finn Yesterday, we introduced barycentric coordinates and de Casteljau s algorithm. Today, we want to go
More informationOn-Line Geometric Modeling Notes
On-Line Geometric Modeling Notes CUBIC BÉZIER CURVES Kenneth I. Joy Visualization and Graphics Research Group Department of Computer Science University of California, Davis Overview The Bézier curve representation
More informationLecture 23: Hermite and Bezier Curves
Lecture 23: Hermite and Bezier Curves November 16, 2017 11/16/17 CSU CS410 Fall 2017, Ross Beveridge & Bruce Draper 1 Representing Curved Objects So far we ve seen Polygonal objects (triangles) and Spheres
More informationEXPLICIT ERROR BOUND FOR QUADRATIC SPLINE APPROXIMATION OF CUBIC SPLINE
J. KSIAM Vol.13, No.4, 257 265, 2009 EXPLICIT ERROR BOUND FOR QUADRATIC SPLINE APPROXIMATION OF CUBIC SPLINE YEON SOO KIM 1 AND YOUNG JOON AHN 2 1 DEPT OF MATHEMATICS, AJOU UNIVERSITY, SUWON, 442 749,
More informationBézier Curves and Splines
CS-C3100 Computer Graphics Bézier Curves and Splines Majority of slides from Frédo Durand vectorportal.com CS-C3100 Fall 2017 Lehtinen Before We Begin Anything on your mind concerning Assignment 1? CS-C3100
More informationTHE PYTHAGOREAN THEOREM
THE STORY SO FAR THE PYTHAGOREAN THEOREM USES OF THE PYTHAGOREAN THEOREM USES OF THE PYTHAGOREAN THEOREM SOLVE RIGHT TRIANGLE APPLICATIONS USES OF THE PYTHAGOREAN THEOREM SOLVE RIGHT TRIANGLE APPLICATIONS
More informationCSE 167: Lecture 11: Bézier Curves. Jürgen P. Schulze, Ph.D. University of California, San Diego Fall Quarter 2012
CSE 167: Introduction to Computer Graphics Lecture 11: Bézier Curves Jürgen P. Schulze, Ph.D. University of California, San Diego Fall Quarter 2012 Announcements Homework project #5 due Nov. 9 th at 1:30pm
More informationO1 History of Mathematics Lecture VIII Establishing rigorous thinking in analysis. Monday 30th October 2017 (Week 4)
O1 History of Mathematics Lecture VIII Establishing rigorous thinking in analysis Monday 30th October 2017 (Week 4) Summary French institutions Fourier series Early-19th-century rigour Limits, continuity,
More informationCGT 511. Curves. Curves. Curves. What is a curve? 2) A continuous map of a 1D space to an nd space
Curves CGT 511 Curves Bedřich Beneš, Ph.D. Purdue University Department of Computer Graphics Technology What is a curve? Mathematical ldefinition i i is a bit complex 1) The continuous o image of an interval
More informationNovel polynomial Bernstein bases and Bézier curves based on a general notion of polynomial blossoming
See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/283943635 Novel polynomial Bernstein bases and Bézier curves based on a general notion of polynomial
More informationIntroduction to Curves. Modelling. 3D Models. Points. Lines. Polygons Defined by a sequence of lines Defined by a list of ordered points
Introduction to Curves Modelling Points Defined by 2D or 3D coordinates Lines Defined by a set of 2 points Polygons Defined by a sequence of lines Defined by a list of ordered points 3D Models Triangular
More informationSmooth Path Generation Based on Bézier Curves for Autonomous Vehicles
Smooth Path Generation Based on Bézier Curves for Autonomous Vehicles Ji-wung Choi, Renwick E. Curry, Gabriel Hugh Elkaim Abstract In this paper we present two path planning algorithms based on Bézier
More informationIntroduction to Computer Graphics. Modeling (1) April 13, 2017 Kenshi Takayama
Introduction to Computer Graphics Modeling (1) April 13, 2017 Kenshi Takayama Parametric curves X & Y coordinates defined by parameter t ( time) Example: Cycloid x t = t sin t y t = 1 cos t Tangent (aka.
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 7 Interpolation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationMAT300/500 Programming Project Spring 2019
MAT300/500 Programming Project Spring 2019 Please submit all project parts on the Moodle page for MAT300 or MAT500. Due dates are listed on the syllabus and the Moodle site. You should include all neccessary
More informationLecture 20: Bezier Curves & Splines
Lecture 20: Bezier Curves & Splines December 6, 2016 12/6/16 CSU CS410 Bruce Draper & J. Ross Beveridge 1 Review: The Pen Metaphore Think of putting a pen to paper Pen position described by time t Seeing
More informationRational Bézier Patch Differentiation using the Rational Forward Difference Operator
Rational Bézier Patch Differentiation using the Rational Forward Difference Operator Xianming Chen, Richard F. Riesenfeld, Elaine Cohen School of Computing, University of Utah Abstract This paper introduces
More informationCheck boxes of Edited Copy of Sp Topics (was 217-pilot)
Check boxes of Edited Copy of 10024 Sp 11 213 Topics (was 217-pilot) College Algebra, 9th Ed. [open all close all] R-Basic Algebra Operations Section R.1 Integers and rational numbers Rational and irrational
More informationNumerische Mathematik
Numer. Math. (997) 76: 479 488 Numerische Mathematik c Springer-Verlag 997 Electronic Edition Exponential decay of C cubic splines vanishing at two symmetric points in each knot interval Sang Dong Kim,,
More informationInterpolation and Deformations A short cookbook
Interpolation and Deformations A short cookbook 600.445 Fall 2000; Updated: 29 September 205 Linear Interpolation p ρ 2 2 = [ 40 30 20] = 20 T p ρ = [ 0 5 20] = 5 p 3 = ρ3 =? 0?? T [ 20 20 20] T 2 600.445
More informationGeometric meanings of the parameters on rational conic segments
Science in China Ser. A Mathematics 005 Vol.48 No.9 09 09 Geometric meanings of the parameters on rational conic segments HU Qianqian & WANG Guojin Department of Mathematics, Zhejiang University, Hangzhou
More informationKatholieke Universiteit Leuven Department of Computer Science
Interpolation with quintic Powell-Sabin splines Hendrik Speleers Report TW 583, January 2011 Katholieke Universiteit Leuven Department of Computer Science Celestijnenlaan 200A B-3001 Heverlee (Belgium)
More informationPrentice Hall Mathematics, Algebra Correlated to: Achieve American Diploma Project Algebra II End-of-Course Exam Content Standards
Core: Operations on Numbers and Expressions Priority: 15% Successful students will be able to perform operations with rational, real, and complex numbers, using both numeric and algebraic expressions,
More informationLecture 20: Lagrange Interpolation and Neville s Algorithm. for I will pass through thee, saith the LORD. Amos 5:17
Lecture 20: Lagrange Interpolation and Neville s Algorithm for I will pass through thee, saith the LORD. Amos 5:17 1. Introduction Perhaps the easiest way to describe a shape is to select some points on
More informationInterpolation. Chapter Interpolation. 7.2 Existence, Uniqueness and conditioning
76 Chapter 7 Interpolation 7.1 Interpolation Definition 7.1.1. Interpolation of a given function f defined on an interval [a,b] by a polynomial p: Given a set of specified points {(t i,y i } n with {t
More informationUtah Core State Standards for Mathematics - Precalculus
A Correlation of A Graphical Approach to Precalculus with Limits A Unit Circle Approach 6 th Edition, 2015 to the Resource Title: with Limits 6th Edition Publisher: Pearson Education publishing as Prentice
More informationArsène Pérard-Gayot (Slides by Piotr Danilewski)
Computer Graphics - Splines - Arsène Pérard-Gayot (Slides by Piotr Danilewski) CURVES Curves Explicit y = f x f: R R γ = x, f x y = 1 x 2 Implicit F x, y = 0 F: R 2 R γ = x, y : F x, y = 0 x 2 + y 2 =
More informationBÉZIER FITTING TO ALMOST OVAL GEAR DEVICES
BÉZIER FITTING TO ALMOST OVAL GEAR DEVICES Karol Gajda (a), Piotr Krawiec (b), Adam Marlewsi (c) (a) Institute of Mathematics, Poznan University of Technology, Poland (b) Chair of Basics of Machine Design,
More informationMCPS Algebra 2 and Precalculus Standards, Categories, and Indicators*
Content Standard 1.0 (HS) Patterns, Algebra and Functions Students will algebraically represent, model, analyze, and solve mathematical and real-world problems involving functional patterns and relationships.
More informationPre-Calculus and Trigonometry Capacity Matrix
Review Polynomials A1.1.4 A1.2.5 Add, subtract, multiply and simplify polynomials and rational expressions Solve polynomial equations and equations involving rational expressions Review Chapter 1 and their
More informationCheck boxes of Edited Copy of Sp Topics (was 261-pilot)
Check boxes of Edited Copy of 10023 Sp 11 253 Topics (was 261-pilot) Intermediate Algebra (2011), 3rd Ed. [open all close all] R-Review of Basic Algebraic Concepts Section R.2 Ordering integers Plotting
More informationO1 History of Mathematics Lecture VIII Establishing rigorous thinking in analysis. Monday 31st October 2016 (Week 4)
O1 History of Mathematics Lecture VIII Establishing rigorous thinking in analysis Monday 31st October 2016 (Week 4) Summary French institutions Fourier series Early-19th-century rigour Limits, continuity,
More informationSMSU Mathematics Course Content
Southwest Minnesota State University Department of Mathematics SMSU Mathematics Course Content 2012-2013 Thefollowing is a list of possibletopics and techniques to cover in your SMSU College Now course.
More informationSupplement. Applications of the Maximum Modulus Theorem to Polynomials
Supplement: Applications of the Maximum Modulus Theorem 1 Supplement. Applications of the Maximum Modulus Theorem to Polynomials Note. These notes are a supplement to Section 4.54 ( The Maximum Principle
More informationCHINO VALLEY UNIFIED SCHOOL DISTRICT INSTRUCTIONAL GUIDE ALGEBRA II
CHINO VALLEY UNIFIED SCHOOL DISTRICT INSTRUCTIONAL GUIDE ALGEBRA II Course Number 5116 Department Mathematics Qualification Guidelines Successful completion of both semesters of Algebra 1 or Algebra 1
More informationDRAFT. Pre-calculus Curriculum Map Quarter 1 Chapters P 2. Extraneous Critical numbers Test intervals
Quarter 1 Chapters P 2 Plot points in the coordinate plane and use distance and midpoint formulas. Sketch graphs of equations. Find and use slope of a line to write and graph linear equations. Solve equations:
More informationax 2 + bx + c = 0 where
Chapter P Prerequisites Section P.1 Real Numbers Real numbers The set of numbers formed by joining the set of rational numbers and the set of irrational numbers. Real number line A line used to graphically
More information2. Polynomial interpolation
2. Polynomial interpolation Contents 2. POLYNOMIAL INTERPOLATION... 1 2.1 TYPES OF INTERPOLATION... 1 2.2 LAGRANGE ONE-DIMENSIONAL INTERPOLATION... 2 2.3 NATURAL COORDINATES*... 15 2.4 HERMITE ONE-DIMENSIONAL
More information1 ** The performance objectives highlighted in italics have been identified as core to an Algebra II course.
Strand One: Number Sense and Operations Every student should understand and use all concepts and skills from the pervious grade levels. The standards are designed so that new learning builds on preceding
More informationMATH-1420 Review Concepts (Haugen)
MATH-40 Review Concepts (Haugen) Unit : Equations, Inequalities, Functions, and Graphs Rational Expressions Determine the domain of a rational expression Simplify rational expressions -factor and then
More information1 f(u)b n i (u)du, (1.1)
AMRX Applied Mathematics Research express 25, No. 4 On a Modified Durrmeyer-Bernstein Operator and Applications Germain E. Randriambelosoa 1 Introduction Durrmeyer [6] has introduced a Bernstein-type operator
More informationOutline. 1 Interpolation. 2 Polynomial Interpolation. 3 Piecewise Polynomial Interpolation
Outline Interpolation 1 Interpolation 2 3 Michael T. Heath Scientific Computing 2 / 56 Interpolation Motivation Choosing Interpolant Existence and Uniqueness Basic interpolation problem: for given data
More informationCurves - Foundation of Free-form Surfaces
Crves - Fondation of Free-form Srfaces Why Not Simply Use a Point Matrix to Represent a Crve? Storage isse and limited resoltion Comptation and transformation Difficlties in calclating the intersections
More informationHome Page. Title Page. Contents. Bezier Curves. Milind Sohoni sohoni. Page 1 of 27. Go Back. Full Screen. Close.
Bezier Curves Page 1 of 27 Milind Sohoni http://www.cse.iitb.ac.in/ sohoni Recall Lets recall a few things: 1. f : [0, 1] R is a function. 2. f 0,..., f i,..., f n are observations of f with f i = f( i
More informationIntermediate Algebra
Intermediate Algebra 978-1-63545-084-2 To learn more about all our offerings Visit Knewton.com Source Author(s) (Text or Video) Title(s) Link (where applicable) Openstax Lyn Marecek, MaryAnne Anthony-Smith
More informationPART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435
PART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435 Professor Biswa Nath Datta Department of Mathematical Sciences Northern Illinois University DeKalb, IL. 60115 USA E mail: dattab@math.niu.edu
More informationII.3. POLYNOMIAL ANALYSIS. Didier HENRION
II.3. POLYNOMIAL ANALYSIS Didier HENRION henrion@laas.fr EECI Graduate School on Control Supélec - Spring 2010 Polynomial methods Based on the algebra of polynomials and polynomial matrices, typically
More informationCollege Algebra with Corequisite Support: Targeted Review
College Algebra with Corequisite Support: Targeted Review 978-1-63545-056-9 To learn more about all our offerings Visit Knewtonalta.com Source Author(s) (Text or Video) Title(s) Link (where applicable)
More informationL(x; 0) = f(x) This family is equivalent to that obtained by the heat or diffusion equation: L t = 1 2 L
Figure 1: Fine To Coarse Scales. The parameter t is often associated with the generation of a scale-space family, where t = 0 denotes the original image and t > 0 increasing (coarser) scales. Scale-Space
More informationarxiv: v4 [cs.gr] 23 Nov 2015
A de Casteljau Algorithm for Bernstein type Polynomials based on (p, q)-integers arxiv:507.040v4 [cs.gr 3 Nov 05 Khalid Khan, ) D.K. Lobiyal ) and Adem ilicman ) ) School of Computer and System Sciences,
More informationTennessee s State Mathematics Standards Precalculus
Tennessee s State Mathematics Standards Precalculus Domain Cluster Standard Number Expressions (N-NE) Represent, interpret, compare, and simplify number expressions 1. Use the laws of exponents and logarithms
More informationMATH Spring 2010 Topics per Section
MATH 101 - Spring 2010 Topics per Section Chapter 1 : These are the topics in ALEKS covered by each Section of the book. Section 1.1 : Section 1.2 : Ordering integers Plotting integers on a number line
More informationLagrange Interpolation and Neville s Algorithm. Ron Goldman Department of Computer Science Rice University
Lagrange Interpolation and Neville s Algorithm Ron Goldman Department of Computer Science Rice University Tension between Mathematics and Engineering 1. How do Mathematicians actually represent curves
More informationSample Exam 1 KEY NAME: 1. CS 557 Sample Exam 1 KEY. These are some sample problems taken from exams in previous years. roughly ten questions.
Sample Exam 1 KEY NAME: 1 CS 557 Sample Exam 1 KEY These are some sample problems taken from exams in previous years. roughly ten questions. Your exam will have 1. (0 points) Circle T or T T Any curve
More informationNonstationary Subdivision Schemes and Totally Positive Refinable Functions
Nonstationary Subdivision Schemes and Totally Positive Refinable Functions Laura Gori and Francesca Pitolli December, 2007 Abstract In this paper we construct a class of totally positive refinable functions,
More informationPractical Linear Algebra: A Geometry Toolbox
Practical Linear Algebra: A Geometry Toolbox Third edition Chapter 2: Here and There: Points and Vectors in 2D Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/pla
More informationAlgebra 2. Chapter 4 Exponential and Logarithmic Functions. Chapter 1 Foundations for Functions. Chapter 3 Polynomial Functions
Algebra 2 Chapter 1 Foundations for Chapter 2 Quadratic Chapter 3 Polynomial Chapter 4 Exponential and Logarithmic Chapter 5 Rational and Radical Chapter 6 Properties and Attributes of Chapter 7 Probability
More informationGEOMETRIC MODELLING WITH BETA-FUNCTION B-SPLINES, I: PARAMETRIC CURVES
International Journal of Pure and Applied Mathematics Volume 65 No. 3 2010, 339-360 GEOMETRIC MODELLING WITH BETA-FUNCTION B-SPLINES, I: PARAMETRIC CURVES Arne Lakså 1, Børre Bang 2, Lubomir T. Dechevsky
More informationAlgebra & Trigonometry for College Readiness Media Update, 2016
A Correlation of Algebra & Trigonometry for To the Utah Core Standards for Mathematics to the Resource Title: Media Update Publisher: Pearson publishing as Prentice Hall ISBN: SE: 9780134007762 TE: 9780133994032
More informationExtrapolation Methods for Approximating Arc Length and Surface Area
Extrapolation Methods for Approximating Arc Length and Surface Area Michael S. Floater, Atgeirr F. Rasmussen and Ulrich Reif March 2, 27 Abstract A well-known method of estimating the length of a parametric
More informationPRECALCULUS BISHOP KELLY HIGH SCHOOL BOISE, IDAHO. Prepared by Kristina L. Gazdik. March 2005
PRECALCULUS BISHOP KELLY HIGH SCHOOL BOISE, IDAHO Prepared by Kristina L. Gazdik March 2005 1 TABLE OF CONTENTS Course Description.3 Scope and Sequence 4 Content Outlines UNIT I: FUNCTIONS AND THEIR GRAPHS
More informationPre-Calculus and Trigonometry Capacity Matrix
Pre-Calculus and Capacity Matri Review Polynomials A1.1.4 A1.2.5 Add, subtract, multiply and simplify polynomials and rational epressions Solve polynomial equations and equations involving rational epressions
More informationNonlinear Stationary Subdivision
Nonlinear Stationary Subdivision Michael S. Floater SINTEF P. O. Box 4 Blindern, 034 Oslo, Norway E-mail: michael.floater@math.sintef.no Charles A. Micchelli IBM Corporation T.J. Watson Research Center
More informationQuantile Textbook Report
Quantile Textbook Report Algebra 2 Author Charles, Randall I., et al StateEdition West Virginia Grade Algebra 2 1 Expressions, Equations, and Inequalities 1.1 Patterns and Expressions 930Q 1.2 Properties
More informationHigh-speed cornering by CNC machines under prescribed bounds on axis accelerations and toolpath contour error
High-speed cornering by CNC machines under prescribed bounds on axis accelerations and toolpath contour error Charlie A. Ernesto and Rida T. Farouki Department of Mechanical and Aerospace Engineering,
More informationPennsylvania Academic Standards for Mathematics Grade 11
correlated to the Pennsylvania Academic Standards for Mathematics Grade 11 1/2003 2001 Algebra 1: Concepts and Skills 2001 correlated to the Pennsylvania Academic Standards for Mathematics Grade 11 2.1
More informationMath 160 Final Exam Info and Review Exercises
Math 160 Final Exam Info and Review Exercises Fall 2018, Prof. Beydler Test Info Will cover almost all sections in this class. This will be a 2-part test. Part 1 will be no calculator. Part 2 will be scientific
More informationNonlinear Means in Geometric Modeling
Nonlinear Means in Geometric Modeling Michael S. Floater SINTEF P. O. Box 124 Blindern, 0314 Oslo, Norway E-mail: Michael.Floater@math.sintef.no Charles A. Micchelli IBM Corporation T.J. Watson Research
More informationSTUDY GUIDE Math 20. To accompany Intermediate Algebra for College Students By Robert Blitzer, Third Edition
STUDY GUIDE Math 0 To the students: To accompany Intermediate Algebra for College Students By Robert Blitzer, Third Edition When you study Algebra, the material is presented to you in a logical sequence.
More informationHigh School Mathematics Math IV
High School Mathematics Math IV The fundamental purpose of Mathematics IV is to generalize and abstract learning accumulated through previous courses and to provide the final springboard to calculus. Students
More informationVolusia County Mathematics Curriculum Map. Pre-Calculus. Course Number /IOD
Volusia County Mathematics Curriculum Map Pre-Calculus Course Number 1202340/IOD Mathematics Department Volusia County Schools Revised June 9, 2012 Pre- Calculus Curriculum Map 120234/IOD COMPONENTS OF
More informationGauss-Lobatto to Bernstein polynomials transformation
Gauss-Lobatto to Bernstein polynomials transformation Loredana Coluccio, Alfredo Eisinberg, Giuseppe Fedele Dip. Elettronica Informatica e Sistemistica, Università degli Studi della Calabria, 87036, Rende
More informationApplications. More Counting Problems. Complexity of Algorithms
Recurrences Applications More Counting Problems Complexity of Algorithms Part I Recurrences and Binomial Coefficients Paths in a Triangle P(0, 0) P(1, 0) P(1,1) P(2, 0) P(2,1) P(2, 2) P(3, 0) P(3,1) P(3,
More informationElementary and Intermediate Algebra
Elementary and Intermediate Algebra 978-1-63545-106-1 To learn more about all our offerings Visit Knewton.com Source Author(s) (Text or Video) Title(s) Link (where applicable) OpenStax Lynn Marecek, Santa
More informationMiller Objectives Alignment Math
Miller Objectives Alignment Math 1050 1 College Algebra Course Objectives Spring Semester 2016 1. Use algebraic methods to solve a variety of problems involving exponential, logarithmic, polynomial, and
More informationMath II. Number and Quantity The Real Number System
MATHEMATICS Math II The high school mathematics curriculum is designed to develop deep understanding of foundational math ideas. In order to allow time for such understanding, each level focuses on concepts
More informationMathematics High School Advanced Mathematics Plus Course
Mathematics High School Advanced Mathematics Plus Course, a one credit course, specifies the mathematics that students should study in order to be college and career ready. The Advanced Mathematics Plus
More informationCS 536 Computer Graphics NURBS Drawing Week 4, Lecture 6
CS 536 Computer Graphics NURBS Drawing Week 4, Lecture 6 David Breen, William Regli and Maxim Peysakhov Department of Computer Science Drexel University 1 Outline Conic Sections via NURBS Knot insertion
More informationCollege Algebra with Corequisite Support: A Blended Approach
College Algebra with Corequisite Support: A Blended Approach 978-1-63545-058-3 To learn more about all our offerings Visit Knewtonalta.com Source Author(s) (Text or Video) Title(s) Link (where applicable)
More informationMyMathLab for School Precalculus Graphical, Numerical, Algebraic Common Core Edition 2016
A Correlation of MyMathLab for School Precalculus Common Core Edition 2016 to the Tennessee Mathematics Standards Approved July 30, 2010 Bid Category 13-090-10 , Standard 1 Mathematical Processes Course
More information9th and 10th Grade Math Proficiency Objectives Strand One: Number Sense and Operations
Strand One: Number Sense and Operations Concept 1: Number Sense Understand and apply numbers, ways of representing numbers, the relationships among numbers, and different number systems. Justify with examples
More informationNotes on uniform convergence
Notes on uniform convergence Erik Wahlén erik.wahlen@math.lu.se January 17, 2012 1 Numerical sequences We begin by recalling some properties of numerical sequences. By a numerical sequence we simply mean
More informationApproximation of Continuous Functions
Approximation of Continuous Functions Francis J. Narcowich October 2014 1 Modulus of Continuit Recall that ever function continuous on a closed interval < a x b < is uniforml continuous: For ever ɛ > 0,
More informationAnalysis-3 lecture schemes
Analysis-3 lecture schemes (with Homeworks) 1 Csörgő István November, 2015 1 A jegyzet az ELTE Informatikai Kar 2015. évi Jegyzetpályázatának támogatásával készült Contents 1. Lesson 1 4 1.1. The Space
More informationCollege Algebra with Corequisite Support: A Compressed Approach
College Algebra with Corequisite Support: A Compressed Approach 978-1-63545-059-0 To learn more about all our offerings Visit Knewton.com Source Author(s) (Text or Video) Title(s) Link (where applicable)
More informationORF 363/COS 323 Final Exam, Fall 2018
Name: Princeton University ORF 363/COS 323 Final Exam, Fall 2018 January 16, 2018 Instructor: A.A. Ahmadi AIs: Dibek, Duan, Gong, Khadir, Mirabelli, Pumir, Tang, Yu, Zhang 1. Please write out and sign
More information