B-spline curves. 1. Properties of the B-spline curve. control of the curve shape as opposed to global control by using a special set of blending

Similar documents
L5 Polynomial / Spline Curves

Geometric Modeling

Curves - Foundation of Free-form Surfaces

CS5620 Intro to Computer Graphics

2.160 System Identification, Estimation, and Learning Lecture Notes No. 17 April 24, 2006

An Expansion of the Derivation of the Spline Smoothing Theory Alan Kaylor Cline

Application of Legendre Bernstein basis transformations to degree elevation and degree reduction

Alternating Direction Implicit Method

Computer Graphics. Geometric Modeling. Geometric Modeling. Page. Copyright Gotsman, Elber, Barequet, Karni, Sheffer Computer Science - Technion

DISTURBANCE TERMS. is a scalar and x i

Five-axis Spline Interpolation Algorithm for Digital Manufacturing System

( ) ( ) ( ( )) ( ) ( ) ( ) ( ) ( ) = ( ) ( ) + ( ) ( ) = ( ( )) ( ) + ( ( )) ( ) Review. Second Derivatives for f : y R. Let A be an m n matrix.

Taylor s Series and Interpolation. Interpolation & Curve-fitting. CIS Interpolation. Basic Scenario. Taylor Series interpolates at a specific

Computer Graphics. Shi-Min Hu. Tsinghua University

Third handout: On the Gini Index

Processing of Information with Uncertain Boundaries Fuzzy Sets and Vague Sets

u(x, t) = u 0 (x ct). This Riemann invariant u is constant along characteristics λ with x = x 0 +ct (u(x, t) = u 0 (x 0 )):

Algorithms Theory, Solution for Assignment 2

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS. Aysegul Akyuz Dascioglu and Nese Isler

Lecture 07: Poles and Zeros

Fundamentals of Regression Analysis

A unified matrix representation for degree reduction of Bézier curves

Numerical Analysis Formulae Booklet

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

A scalar t is an eigenvalue of A if and only if t satisfies the characteristic equation of A: det (A ti) =0

On the convergence of derivatives of Bernstein approximation

Assignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class)

Bezier curve and its application

ρ < 1 be five real numbers. The

Ideal multigrades with trigonometric coefficients

The Control Vector Scheme for Design of Planar Primitive PH curves

Functions of Random Variables

Open and Closed Networks of M/M/m Type Queues (Jackson s Theorem for Open and Closed Networks) Copyright 2015, Sanjay K. Bose 1

Log1 Contest Round 2 Theta Complex Numbers. 4 points each. 5 points each

Chapter 4 Multiple Random Variables

Lecture 5: Interpolation. Polynomial interpolation Rational approximation

Lecture 3 Probability review (cont d)

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II

Consumer theory. A. The preference ordering B. The feasible set C. The consumption decision. A. The preference ordering. Consumption bundle

2. Independence and Bernoulli Trials

Meromorphic Solutions of Nonlinear Difference Equations

Chapter 14 Logistic Regression Models

Discrete Adomian Decomposition Method for. Solving Burger s-huxley Equation

Econometric Methods. Review of Estimation

TESTS BASED ON MAXIMUM LIKELIHOOD

A NEW FINITE ELEMENT CONSIDERING SHEAR LAG

Motion Estimation Based on Unit Quaternion Decomposition of the Rotation Matrix

Summary of the lecture in Biostatistics

Transforms that are commonly used are separable

Journal of Mathematical Analysis and Applications

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

Part 4b Asymptotic Results for MRR2 using PRESS. Recall that the PRESS statistic is a special type of cross validation procedure (see Allen (1971))

arxiv:math/ v1 [math.gm] 8 Dec 2005

Overview of the weighting constants and the points where we evaluate the function for The Gaussian quadrature Project two

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

Special Instructions / Useful Data

Lecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES

Multiple Choice Test. Chapter Adequacy of Models for Regression

International Journal of Scientific & Engineering Research, Volume 5, Issue 9, September ISSN

Multivariate Transformation of Variables and Maximum Likelihood Estimation

u(x, t) = u 0 (x ct). This Riemann invariant u is constant along characteristics λ with x = x 0 +ct (u(x, t) = u 0 (x 0 )):

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

HYPERBOLIC PARABOLOID SPLINE FOR SPATIAL DATA FITTING AND APPROXIMATION

Lecture 2: The Simple Regression Model

A tighter lower bound on the circuit size of the hardest Boolean functions

For combinatorial problems we might need to generate all permutations, combinations, or subsets of a set.

18.413: Error Correcting Codes Lab March 2, Lecture 8

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

Lecture 3. Sampling, sampling distributions, and parameter estimation

PROJECTION PROBLEM FOR REGULAR POLYGONS

EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM

Simple Linear Regression

Suggested Answers, Problem Set 4 ECON The R 2 for the unrestricted model is by definition u u u u

Analysis of Lagrange Interpolation Formula

N-dimensional Auto-Bäcklund Transformation and Exact Solutions to n-dimensional Burgers System

PTAS for Bin-Packing

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin

A Study on Generalized Generalized Quasi hyperbolic Kac Moody algebra QHGGH of rank 10

Polyphase Filters. Section 12.4 Porat

CS 1675 Introduction to Machine Learning Lecture 12 Support vector machines

Mechanical Design Technology (Free-form Surface) April 28, /12

The internal structure of natural numbers, one method for the definition of large prime numbers, and a factorization test

Analysis of Variance with Weibull Data

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions.

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS

Subdivision of Uniform ωb-spline Curves and Two Proofs of Its C k 2 -Continuity

Chapter 4 Multiple Random Variables

ECONOMETRIC THEORY. MODULE VIII Lecture - 26 Heteroskedasticity

Introduction to local (nonparametric) density estimation. methods

Computational Geometry

A conic cutting surface method for linear-quadraticsemidefinite

Given a table of data poins of an unknown or complicated function f : we want to find a (simpler) function p s.t. px (

Chapter 16 Measurement Error Models

LINEAR RECURRENT SEQUENCES AND POWERS OF A SQUARE MATRIX

h-analogue of Fibonacci Numbers

EXPECTATION IDENTITIES OF GENERALIZED RECORD VALUES FROM NEW WEIBULL-PARETO DISTRIBUTION AND ITS CHARACTERIZATION

Probability and. Lecture 13: and Correlation

Research Article Gauss-Lobatto Formulae and Extremal Problems

Transcription:

B-sple crve Copyrght@, YZU Optmal Desg Laboratory. All rghts reserved. Last pdated: Yeh-Lag Hs (--9). ote: Ths s the corse materal for ME Geometrc modelg ad compter graphcs, Ya Ze Uversty. art of ths materal s adapted from CAD/CAM Theory ad ractce, by Ibrahm Zed, McGraw-Hll, 99. Ths materal s be sed strctly for teachg ad learg of ths corse. B-sple crves. ropertes of the B-sple crve B-sple crves are powerfl geeralzato of Bezer crves. They provde local cotrol of the crve shape as opposed to global cotrol by sg a specal set of bledg fctos that provde local flece. They also provde the ablty to add cotrol pots wthot creasg the degree of the crve. B-sple crves have the ablty to terpolate or approxmate a set of gve data pots. Iterpolato s sefl dsplayg desg or egeerg reslts sch as stress or dsplacemet dstrbto a part whle approxmato s good to desg free-form crves. Iterpolato s also sefl f the desger has measred data pots had that mst le o the resltg crve. Assgmet Use yor CAD software to draw several B-sple crves. Use the same cotrol pots to draw both approxmato crves ad terpolato crves. Report yor observatos. The theory of B-sple crves separates the degree of the resltg crve from the mber of the gve cotrol pots. The B-sple crve defed by + cotrol pots s gve by http://desger.mech.yz.ed.tw/

B-sple crve where,, max (), are the B-sple fctos sggested by Cox ad de Boor 97. Frst, the parameter cotrols the degree (-) of the resltg B-sple crve ad s sally depedet of the mber of cotrol pots. The B-sple fctos have the followg propertes: artto of ty:, ostvty:, Local spport:, f, Cotty:, s (-) tmes cotosly dfferetable The frst property esres that the relatoshp betwee the crve ad ts defg cotrol pots s varat der affe trasformatos. The secod property garatees that the crve segmet les completely wth the covex hll of. The thrd property dcates that each segmet of a B-sple crve s fleced by oly cotrol pots or each cotrol pot affects oly crve segmets, as show Fgre. It s sefl to otce that the Berste polyomal, B,, has the same frst two propertes metoed above. http://desger.mech.yz.ed.tw/

B-sple crve * ** Fgre. Local cotrol of B-sple crves. Assgmet Assme a set of two-dmesoal cotrol pots the x-y plae, the B-sple crve defed by ths set of cotrol pots s,. Assme a set of coeffcets the followg affe trasformato dscssed Chapter to trasform () to get ˆ : x a y a xx yx x a x a xy yy y b x y b y Use the same affe trasformato to trasform to get crve defed by ths set of cotrol pots s ˆ =, ad the B-sple,. Show that sg partto of ty dscssed above. Ths also esres that the relatoshp betwee the crve ad ts defg cotrol pots s varat der affe trasformatos. Ths s a hard qesto. Yo do ot have to prove t for geeral cases. Yo jst have to assme a set of cotrol pots ad coeffcets for the affe trasformato, ad show that ths s tre for the mbers yo assmed. http://desger.mech.yz.ed.tw/

B-sple crve Assgmet Use yor CAD software to draw crves smlar to those Fgre. Idcate the local cotrol of the crves.. The B-sple fcto The B-sple fcto also has the property of recrso, whch s defed as,,, () where,, (), otherwse The are called parametrc ots or ot vales. For a ope crve,, j j,, j j j () where ad the rage of s j () () Relato () shows that (++) ots are eeded to create a (-) degree crve defed by (+) cotrol pots. (7) Ths relato shows that a mmm of two, three, ad for cotrol pots are reqred to defe a lear, qadratc, ad cbc B-sple crve respectvely. A cbc B-sple s sffcet for a large mber of applcatos. Fgre shows the shapes of the B-sple fctos. http://desger.mech.yz.ed.tw/

B-sple crve. (a) Lear fcto ( = ).7.. (b) Qadratc fcto ( = ). (c) Cbc fcto ( = ) Fgre. B-sple fctos.. A Example for Dervg the B-sple Fcto To derstad the recrsve atre of the B-sple fcto, the B-sple fcto for a cbc sple crve defed by cotrol pots s derved. Ths cbc sple has ad. Eght ots are eeded to calclate the B-sple fctos. Eqato () gves the ot vector as 7 The rage of (Eqato ()) s. Eqato () gves, (8),,,, To calclate the above B-sple fctos, se Eqato () ad () together wth the ot vector as follows: http://desger.mech.yz.ed.tw/

B-sple crve,,,,, elsewhere,,, elsewhere,,,,, elsewhere,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 7 7,,,,,,,,,,,,,,,,,,,,, http://desger.mech.yz.ed.tw/

B-sple crve,,,,, 7 7,,,,,,,,,,,,,, 7,, 7 Sbstttg, to Eqato (8) gves,, Sbstttg, to ths eqato gves the crve eqato as, Ths eqato s the same as the oe for the Bezer crve wth the same cotrol pots. Ths the cbc B-sple crve defed by for cotrol pots s detcal to the cbc Bezer crve defed by the same pots. Ths fact ca be geeralzed for a (-)-degree crve defed by cotrol pots. Assgmet Derve the B-sple fctos for =, =, ad =, =. Do t oce yor lfetme. Derve the B-sple fctos for =, =, ad shows ts local cotrol capablty. The gap betwee the eghborg ot vales defe Eqato () s always form, wth the vale of. These ots are called form ots ad a B-sple crve based o form ots s called a form B-sple crve. As we modfy the shape of a crve, we ofte add or delete ot vales ad so prodce o-form gaps betwee the ots, whch reslt o-form B-sple crves. o-form B-sples are 7 http://desger.mech.yz.ed.tw/

B-sple crve cosdered the geeral form of B-sple crves. Ths most CAD systems provde the capablty of creatg ad modfyg o-form B-sple.. Ratoal Crves A ratoal crve s defed by the algebrac rato of two polyomals whle a oratoal crve s defed by oe polyomal. I a ratoal crve, the cotrol pots, y, z are gve homogeeos coordates x h, y h, z h, h x. Ths the coordates of a pot o a ratoal crve the homogeeos space, x h, y h, z h, h s obtaed from x h y h z h h h x h y h z, (9), (), () h, () Ths a ratoal B-sple crve ca be expressed as R,, max () where R, are the ratoal B-sple bass fctos ad are gve by h, R, () h, The above eqato shows that R fctos,. If we sbsttte bass fctos R oratoal B-sple coterparts., are a geeralzato of the oratoal bass h the eqato, R. The ratoal,,, have early all the aalytc ad geometrc characterstcs of ther 8 http://desger.mech.yz.ed.tw/

B-sple crve Assgmet Show that the ratoal B-sple bass fctos propertes: partto of ty ad postvty. R, have the followg to se ote that The ma dfferece betwee ratoal ad oratoal B-sple crves s the ablty h at each cotrol pot to cotrol the behavor of the ratoal B-sples. h s do ot have to be tegers. The greater the vale of a partclar closer the crve s plled toward the cotrol pot. h, the Ratoal sples provde a exact represetato for qadrc crves (cocs), sch as crcles ad ellpses. oratoal sples, whch are polyomals, ca oly approxmate cocs. Ths allows graphcs pacages to model all crve shapes wth oe represetato, ratoal sples, wthot eedg a lbrary of crve fctos to hadle dfferet desg shapes. For example, sg three cotrol pots,,,,, homogeeos coordates h =h =, h cos ad ther, a crclar arc ca be defed as: () Ad x y () The we ca easly show that x y, that s, Assgmet s a exact crclar arc. Wrte a program Matlab to draw a fll crcle sg Eqatos () ad (). Draw a crcle the tradtoal way ad compare the two crcles. Show yor Matlab program too. 9 http://desger.mech.yz.ed.tw/