Fitting a Natural Spline to Samples of the Form (t, f(t))
|
|
- Diana Todd
- 5 years ago
- Views:
Transcription
1 Fitting a Natural Spline to Samples of the Form (t, f(t)) David Eberly, Geometric Tools, Redmond WA This work is licensed under the Creative Commons Attribution 4 International License To view a copy of this license, visit or send a letter to Creative Commons, PO Box 1866, Mountain View, CA 9442, USA Created: January 5, 28 Last Modified: January 5, 215 Contents 1 Introduction 2 11 Constraints for Free Splines 2 12 Constraints for Clamped Splines 3 13 Constraints for Periodic Splines 3 2 Constructing Free Splines 3 3 Constructing Clamped Splines 5 4 Constructing Periodic Splines 6 41 Three Samples 6 42 Four Samples 7 43 The General Case 9 1
2 1 Introduction Consider a function f : D IR m, whose domain is D = [t min, t max ] and whose range consists of m-tuples with m 1 The only information we have about the function is a set of samples, {(t i, f i )} n i=, f i = f(t i ) and t min t < t 1 < < t n < t n t max We want to fit the function with a piecewise cubic polynomial spline S (t), t [t, t 1 ] S 1 (t), t [t 1, t 2 ] S(t) = S n 2 (t), t [t n 2, t n ] S n (t), t [t n, t n ] S i (t) = a i + b i (t t i ) + c i (t t i ) 2 + d i (t t i ) 3 (2) for i n 1 Moreover, we want S(t) to be a C 2 function; that is, S(t), S (t), and S (t) must be continuous functions of t on the interval [t, t n ] Each of S i (t) is a C 2 function on the open interval (t i, t i+1 ), so the only points of concern to be C 2 every are the sample points We have n polynomials to determine, each having 4 unknown coefficients Thus, we have 4n unknowns Imposing the C 2 condition leads to linear constraints on the unknowns We need the polynomials, their first derivatives, and their second derivatives to match at the sample points t 1 through t n Specifically, S i (t i+1 ) = S i+1 (t i ), S i(t i+1 ) = S i+1(t i ), S i (t i+1 ) = S i+1(t i ), i n 2 (3) This is a set of 3(n 1) constraints on the unknowns We also need the polynomials to pass through the sample points, so S i (t i ) = f i, i n 1, and S n (t n ) = f n (4) This is a set of n + 1 constraints on the unknowns Equations (3) and (4) represent 4n 2 constraints on the 4n unknowns We have two degrees of freedom Three standard choices are the following (1) 11 Constraints for Free Splines We may allow the endpoints to have whatever slope necessary by specifying the second derivatives at those points to be zero: S (t ) =, S n(t n ) = (5) The resulting curve is said to be a free spline These two constraints and the previously mentioned ones give us 4n linear equations in 4n unknowns; that is, the coefficients of the cubic polynomials are obtained by solving a linear system of equations It turns out that the system is tridiagonal, so a specialized linear solver may be used that is O(n) General linear solvers are O(n 3 ), so the specialized solver saves computational time for a large number of sample points 2
3 12 Constraints for Clamped Splines We may clamp the endpoints in the sense that we specify the first derivatives at those points: S (t ) = σ, S n(t n ) = σ 1 (6) for user-specified constants σ and σ 1 The resulting curve is said to be a clamped spline These two constraints and the previously mentioned ones give us 4n linear equations in 4n unknowns As for the case of free splines, the linear system is tridiagonal and may be solved by a specialized linear solver 13 Constraints for Periodic Splines We may require that the first and last samples have the same function values, namely, that f n = f Moreover, the C 2 condition requires us to impose the derivative equality constraints S (t ) = S n(t n ), S (t ) = S n(t n ) (7) The resulting curve is said to be a periodic spline These two constraints and the previously mentioned ones produce a system of 4n linear equations in 4n unknowns, but this system is not tridiagonal However, the system may be solved by a specialized linear solver that is O(n) 2 Constructing Free Splines This section describes how to set up and solve the tridiagonal linear system of equations for the polynomial coefficients when the endpoints are required to have zero second-order derivatives The ratio of differences of function values and time values occurs often in the derivation, so let us first define q i = f i+1 f i i (8) Define i = t i+1 t i for all relevant i The conditions S i (t i+1 ) = S i+1 (t i ) and S i (t i ) = f i are equivalent to The conditions S i (t i+1) = S i+1 (t i) are equivalent to The conditions S i (t i+1) = S i+1 (t i) are equivalent to i b i + 2 i c i + 3 i d i = f i+1 f i (9) 2 i c i i d i = b i+1 b i (1) 3 i d i = c i+1 c i (11) Substitute Equation (11) into Equations (9) and (1) to obtain ( ) i b i + 2 ci+1 + 2c i i = f i+1 f i (12) 3 and i (c i+1 + c i ) = b i+1 b i (13) 3
4 Increment the index in Equation (12) to obtain ( ) i+1 b i ci+2 + 2c i+1 i+1 = f i+2 f i+1 (14) 3 Compute i times Equation (14) and subtract from it i+1 times Equation (12) to obtain ( ) i i+1 (b i+1 b i ) + i 2 ci+2+2c i+1 i+1 3 i+1 2 i = i (f i+2 f i+1 ) i+1 (f i+1 f i ) ( ) ci+1+2c i 3 (15) Substitute Equation (13) into Equation (15), collect terms, and divide by i i+1 to obtain i+1 c i+2 + 2( i+1 + i )c i+1 + i c i = 3(q i+1 q i ) (16) which is defined for i n 3 We know that the constraint S (t ) = in Equation (5) implies c = The constraint S n(t n ) in Equation (5) implies n c n nd n = If we define c n =, this constraint fits into the construction here and Equation (16) also holds for i = n 2 Essentially, we are thinking of the spline curve having one additional segment S n (t) = f n + b n (t t n ), a line segment, so that c n = d n = Equation (16) and boundary conditions c = c n = may be solved by setting up a linear system Mc = p (17) the (n 1) (n 1) matrix M is tridiagonal, symmetric, and diagonally dominant: 2( + 1 ) 1 1 2( ) 2 2 2( ) 3 M = n 4 2( n 4 + n 3 ) n 3 n 3 2( n 3 + n 2 ) n 2 n 2 2( n 2 + n ) (18) The (n 1) 1 column vector c has components c 1 through c n and the (n 1) 1 column vector p has components p 1 through p n p i = 3(q i q i ) (19) The matrix M and vector p may be computed and then used as input to a linear system solver that is specialized for tridiagonal systems The solver computes c 1 through c n, and we already know that c = The d-terms are computed from Equation (11), d i = c i+1 c i 3 i (2) and the b-terms are computed from Equation (12), b i = q i i ( ci+1 + 2c i 3 ) (21) The a-terms are already known, a i = f i (22) 4
5 3 Constructing Clamped Splines This section describes how to set up and solve the tridiagonal linear system of equations for the polynomial coefficients when the endpoints are required to have user-specified first-order derivatives Multiply Equation (9) by 3/ i and subtract from it Equation (1) to obtain Increment the index in Equation (23) to obtain 3b i + i c i = 3q i (b i+1 b i ) (23) 3b i+1 + i+1 c i+1 = 3q i+1 (b i+2 b i+1 ) (24) Multiply Equation (23) by i+1, multiply Equation (24) by i and add the two together to obtain 3( i+1 b i + i b i+1 ) + i+1 i (c i+1 + c i ) = 3( i+1 q i + i q i+1 ) i+1 (b i+1 b i ) i (b i+2 b i+1 ) (25) In the left-hand side of Equation (25), substitute Equation (13) to eliminate the c-terms Some algebra simplifies this to i b i+1 + 2( i + i+1 )b i+1 + i+1 b i = 3( i+1 q i + i q i+1 ) (26) which is defined for i n 3 We know that the constraint S (t ) = σ in Equation (6) implies b = σ The constraint S n(t n ) = σ 1 in Equation (6) implies n b n +2 2 nc n +3 3 nd n = σ 1 If we define b n = σ 1, this constraint fits into the construction here and Equation (26) also holds for i = n 2 Essentially, we are thinking of the spline curve having one additional segment S n (t) = f n + b n (t t n ) + c n (t t n ) 2, a parabolic arc, so that 2c n = 2 n c n nd n and d n = Equation (26) and boundary conditions b = σ and b n = σ 1 may be solved by setting up a linear system Nb = r (27) the (n 1) (n 1) matrix N is tridiagonal, symmetric, and diagonally dominant: 2( + 1 ) 2 2( ) 1 3 2( ) 2 N = n 3 2( n 4 + n 3 ) n 4 n 2 2( n 3 + n 2 ) n 3 n 2( n 2 + n ) (28) The (n 1) 1 column vector b has components b 1 through b n and the (n 1) 1 column vector r has components r 1 through r n r 1 = 3( 1 q + q 1 ) 1 σ r i = 3( i q i + i q i ), 2 i n 2 (29) r n = 3( n q n 2 + n 2 q n ) n 2 σ 1 The matrix N and vector r may be computed and then used as input to a linear system solver that is specialized for tridiagonal systems The solver computes b 1 through b n The c-terms are computed from Equation (23), c i = 3q i b i+1 2b i i (3) 5
6 and the d-terms are computed from Equation (1), d i = b i+1 b i 2 i c i 3 2 i (31) The a-terms are already known, a i = f i (32) 4 Constructing Periodic Splines This section describes how to set up and solve the linear system of equations for the polynomial coefficients when the endpoints are required to have the same function value, same first-order derivative, and same second-order derivative 41 Three Samples Consider the simplest case of three samples: (t, f ), (t 1, f 1 ), and (t 2, f ) The two cubic polynomials are S (t) = a + b (t t ) + c (t t ) 2 + d (t t ) 3 and S 1 (t) = a 1 + b 1 (t t 1 ) + c 1 (t t 1 ) 2 + d 1 (t t 1 ) 3 The eight constraints are f = S (t ) = S 1 (t 2 ) S (t ) = S 1(t 2 ) S (t ) = S 1 (t 2 ) f 1 = S (t 1 ) = S 1 (t 1 ) S (t 1 ) = S 1(t 1 ) S (t 1 ) = S 1 (t 1 ) Define = t 1 t and 1 = t 2 t 1 In terms of the polynomial coefficients, the constraints are (33) a + b + c 2 + d 3 = a 1 = f 1 b + 2c + 3d 2 = b 1 2c + 6d = 2c 1 a 1 + b c d = a = f b 1 + 2c d = b (34) 2c 1 + 6d 1 1 = 2c The linear system of equations is a b c d a 1 b 1 c 1 d 1 = f f 1 (35) 6
7 In block matrix form, M L K = f U (36) L M 1 K 1 f 1 U 1 1 M i = i 2 i 3 i i 3 2, L =, K i = i Notice that M i is invertible a i b i c i d i 1, U = The linear system may be solved by row-reducing the augmented block matrix as shown next The matrix I is the 4 4 identity matrix M L f U I M f M L M 1 f 1 U L M 1 f 1 U (37) I I M L f M U M 1 LM L (f 1I + f LM )U M L f M U I (M 1 LM L) (f 1 I + f LM )U I [f M + M 1 LM (f 1 I + f LM I (M 1 LM (f 1 I + f LM (38) = I K I K 1 K 1 = (M 1 LM L) (f 1 I + f LM ), K = f M U + M LK 1 (39) 42 Four Samples Similar to the case of three samples, the constraints may be formulated as a block-matrix system, M L K f U M 1 L K 1 = f 1 U L M 2 K 2 f 2 U (4) 7
8 The linear system may be solved by row-reducing the augmented block matrix, M L f U M 1 L f 1U L M 2 f 2U I M L fm U I M 1 L f1m1 L M 2 f 2U I M L fm U I M 1 L f1m1 U LM L M2 (f2i + flm )U I M L fm U I M 1 L f1m1 M 2 LM LM1 L [f2i + LM (fi + f1lm1 I M L fm U I M 1 L f1m1 I K 2 (41) K 2 = (M 2 LM LM1 L) [f 2 I + LM (f I + f 1 LM1 )]U (42) Back substitution may now be used to compute K and K 1 The first back substitution is I M L f M U I M L f M U I M 1 L f 1M1 U I K 1 I K 2 I K 2 (43) K 1 = f 1 M 1 U + M 1 LK 2 (44) The second back substitution is I M M I K I K 1 I K 1 I K 2 I K 2 (45) K = f M U + M LK 1 (46) Equations (46), (44), and (42) are the polynomial coefficients for the spline 8
9 43 The General Case In general, the constraints may be formulated as a block-matrix system, M L K M 1 L K 1 M n 2 L L M n K n = f U f 1 U f n U (47) whose solution is obtained by row reducing the augmented block matrix M L f U M 1 L f 1 U M n 2 L f n 2 U L M n f n U (48) Inverting the M i matrices in the first n 1 block rows leads to I M L f M U I M1 L f 1M1 U I Mn 2 L f n 2Mn 2 U L M n f n U Forward elimination by the first n 1 block rows leads to I M L f M U I M1 L f 1M1 U I Mn 2 L f n 2Mn 2 U P Q (49) (5) and Q = [ f n I + f LM + f 1 LM = f n U + LM P = M n ( LM LM1 LMn 2) L (51) ( LM1 + + f n 2 LM LM1 LMn 2)] U [ f U + LM1 (f 1 U + ) ] (52) 9
10 The last expression shows that Q may be computed efficiently in a nested manner Inverting P in the last row of the reduced augmented block matrix, I M L f M U I M1 L f 1M1 U I Mn 2 L f n 2Mn 2 U I K n (53) K n = P Q (54) Back substitution may be applied n 1 times to solve for the other coefficients Generally, these are K i = f i M i U + M i LK i+1, n 2 i (55) In the calculations, we need the inverses of the M i matrices These are M i = i i 2 i i 2 i i 2 3 i 1 3 i 2 i 1 i (56) 1
The Laplace Expansion Theorem: Computing the Determinants and Inverses of Matrices
The Laplace Expansion Theorem: Computing the Determinants and Inverses of Matrices David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative
More informationA Relationship Between Minimum Bending Energy and Degree Elevation for Bézier Curves
A Relationship Between Minimum Bending Energy and Degree Elevation for Bézier Curves David Eberly, Geometric Tools, Redmond WA 9852 https://www.geometrictools.com/ This work is licensed under the Creative
More informationDistance Between Ellipses in 2D
Distance Between Ellipses in 2D David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To
More informationLeast-Squares Fitting of Data with Polynomials
Least-Squares Fitting of Data with Polynomials David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International
More informationPerspective Projection of an Ellipse
Perspective Projection of an Ellipse David Eberly, Geometric ools, Redmond WA 98052 https://www.geometrictools.com/ his work is licensed under the Creative Commons Attribution 4.0 International License.
More informationSolving Systems of Polynomial Equations
Solving Systems of Polynomial Equations David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationReconstructing an Ellipsoid from its Perspective Projection onto a Plane
Reconstructing an Ellipsoid from its Perspective Projection onto a Plane David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons
More informationInformation About Ellipses
Information About Ellipses David Eberly, Geometric Tools, Redmond WA 9805 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view
More informationInterpolation of Rigid Motions in 3D
Interpolation of Rigid Motions in 3D David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationIntersection of Ellipsoids
Intersection of Ellipsoids David Eberly, Geometric Tools, Redmond WA 9805 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view
More informationCubic Splines MATH 375. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Cubic Splines
Cubic Splines MATH 375 J. Robert Buchanan Department of Mathematics Fall 2006 Introduction Given data {(x 0, f(x 0 )), (x 1, f(x 1 )),...,(x n, f(x n ))} which we wish to interpolate using a polynomial...
More informationIntersection of Rectangle and Ellipse
Intersection of Rectangle and Ellipse David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationDistance from Line to Rectangle in 3D
Distance from Line to Rectangle in 3D David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationCubic Splines; Bézier Curves
Cubic Splines; Bézier Curves 1 Cubic Splines piecewise approximation with cubic polynomials conditions on the coefficients of the splines 2 Bézier Curves computer-aided design and manufacturing MCS 471
More informationIntersection of Objects with Linear and Angular Velocities using Oriented Bounding Boxes
Intersection of Objects with Linear and Angular Velocities using Oriented Bounding Boxes David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the
More information1 Piecewise Cubic Interpolation
Piecewise Cubic Interpolation Typically the problem with piecewise linear interpolation is the interpolant is not differentiable as the interpolation points (it has a kinks at every interpolation point)
More informationMA 323 Geometric Modelling Course Notes: Day 07 Parabolic Arcs
MA 323 Geometric Modelling Course Notes: Day 07 Parabolic Arcs David L. Finn December 9th, 2004 We now start considering the basic curve elements to be used throughout this course; polynomial curves and
More informationIntersection of Infinite Cylinders
Intersection of Infinite Cylinders David Eberly, Geometric Tools, Redmond WA 980 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationSPLINE INTERPOLATION
Spline Background SPLINE INTERPOLATION Problem: high degree interpolating polynomials often have extra oscillations. Example: Runge function f(x = 1 1+4x 2, x [ 1, 1]. 1 1/(1+4x 2 and P 8 (x and P 16 (x
More informationMatrix Basic Concepts
Matrix Basic Concepts Topics: What is a matrix? Matrix terminology Elements or entries Diagonal entries Address/location of entries Rows and columns Size of a matrix A column matrix; vectors Special types
More informationLearning Module 1 - Basic Algebra Review (Appendix A)
Learning Module 1 - Basic Algebra Review (Appendix A) Element 1 Real Numbers and Operations on Polynomials (A.1, A.2) Use the properties of real numbers and work with subsets of the real numbers Determine
More informationMoving Along a Curve with Specified Speed
Moving Along a Curve with Specified Speed David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International
More informationConverting Between Coordinate Systems
Converting Between Coordinate Systems David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationIntersection of Ellipses
Intersection of Ellipses David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view
More informationLow-Degree Polynomial Roots
Low-Degree Polynomial Roots David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view
More informationCollege Algebra. Course Text. Course Description. Course Objectives. StraighterLine MAT101: College Algebra
College Algebra Course Text Barnett, Raymond A., Michael R. Ziegler, and Karl E. Byleen. College Algebra, 8th edition, McGraw-Hill, 2008, ISBN: 9780072867381 [This text is available as an etextbook at
More informationChapter 7. Tridiagonal linear systems. Solving tridiagonal systems of equations. and subdiagonal. E.g. a 21 a 22 a A =
Chapter 7 Tridiagonal linear systems The solution of linear systems of equations is one of the most important areas of computational mathematics. A complete treatment is impossible here but we will discuss
More informationPiecewise Polynomial Interpolation
Piecewise Polynomial Interpolation 1 Piecewise linear interpolation Suppose we have data point (x k,y k ), k =0, 1,...N. A piecewise linear polynomial that interpolates these points is given by p(x) =p
More informationComputing Orthonormal Sets in 2D, 3D, and 4D
Computing Orthonormal Sets in 2D, 3D, and 4D David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International
More informationGlobal polynomial interpolants suffer from the Runge Phenomenon if the data sites (nodes) are not chosen correctly.
Piecewise polynomial interpolation Global polynomial interpolants suffer from the Runge Phenomenon if the data sites (nodes) are not chosen correctly. In many applications, one does not have the freedom
More informationNumerical Analysis Fall. Gauss Elimination
Numerical Analysis 2015 Fall Gauss Elimination Solving systems m g g m m g x x x k k k k k k k k k 3 2 1 3 2 1 3 3 3 2 3 2 2 2 1 0 0 Graphical Method For small sets of simultaneous equations, graphing
More informationLecture 1 INF-MAT : Chapter 2. Examples of Linear Systems
Lecture 1 INF-MAT 4350 2010: Chapter 2. Examples of Linear Systems Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo August 26, 2010 Notation The set of natural
More informationMTH5112 Linear Algebra I MTH5212 Applied Linear Algebra (2017/2018)
MTH5112 Linear Algebra I MTH5212 Applied Linear Algebra (2017/2018) COURSEWORK 3 SOLUTIONS Exercise ( ) 1. (a) Write A = (a ij ) n n and B = (b ij ) n n. Since A and B are diagonal, we have a ij = 0 and
More informationQuestion: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI?
Section 5. The Characteristic Polynomial Question: Given an n x n matrix A, how do we find its eigenvalues? Idea: Suppose c is an eigenvalue of A, then what is the determinant of A-cI? Property The eigenvalues
More informationLinear Algebra Practice Problems
Math 7, Professor Ramras Linear Algebra Practice Problems () Consider the following system of linear equations in the variables x, y, and z, in which the constants a and b are real numbers. x y + z = a
More informationLinear Algebra I Lecture 8
Linear Algebra I Lecture 8 Xi Chen 1 1 University of Alberta January 25, 2019 Outline 1 2 Gauss-Jordan Elimination Given a system of linear equations f 1 (x 1, x 2,..., x n ) = 0 f 2 (x 1, x 2,..., x n
More informationNotes on Row Reduction
Notes on Row Reduction Francis J. Narcowich Department of Mathematics Texas A&M University September The Row-Reduction Algorithm The row-reduced form of a matrix contains a great deal of information, both
More information1 - Systems of Linear Equations
1 - Systems of Linear Equations 1.1 Introduction to Systems of Linear Equations Almost every problem in linear algebra will involve solving a system of equations. ü LINEAR EQUATIONS IN n VARIABLES We are
More informationand let s calculate the image of some vectors under the transformation T.
Chapter 5 Eigenvalues and Eigenvectors 5. Eigenvalues and Eigenvectors Let T : R n R n be a linear transformation. Then T can be represented by a matrix (the standard matrix), and we can write T ( v) =
More informationEigenvalues and Eigenvectors: An Introduction
Eigenvalues and Eigenvectors: An Introduction The eigenvalue problem is a problem of considerable theoretical interest and wide-ranging application. For example, this problem is crucial in solving systems
More informationEmpirical Models Interpolation Polynomial Models
Mathematical Modeling Lia Vas Empirical Models Interpolation Polynomial Models Lagrange Polynomial. Recall that two points (x 1, y 1 ) and (x 2, y 2 ) determine a unique line y = ax + b passing them (obtained
More informationAn overview of key ideas
An overview of key ideas This is an overview of linear algebra given at the start of a course on the mathematics of engineering. Linear algebra progresses from vectors to matrices to subspaces. Vectors
More informationLECTURES 4/5: SYSTEMS OF LINEAR EQUATIONS
LECTURES 4/5: SYSTEMS OF LINEAR EQUATIONS MA1111: LINEAR ALGEBRA I, MICHAELMAS 2016 1 Linear equations We now switch gears to discuss the topic of solving linear equations, and more interestingly, systems
More informationMath 4377/6308 Advanced Linear Algebra
3.1 Elementary Matrix Math 4377/6308 Advanced Linear Algebra 3.1 Elementary Matrix Operations and Elementary Matrix Jiwen He Department of Mathematics, University of Houston jiwenhe@math.uh.edu math.uh.edu/
More informationDistance to Circles in 3D
Distance to Circles in 3D David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view
More informationPolyhedral Mass Properties (Revisited)
Polyhedral Mass Properties (Revisited) David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationElementary maths for GMT
Elementary maths for GMT Linear Algebra Part 2: Matrices, Elimination and Determinant m n matrices The system of m linear equations in n variables x 1, x 2,, x n a 11 x 1 + a 12 x 2 + + a 1n x n = b 1
More informationAlgebra 2 Matrices. Multiple Choice Identify the choice that best completes the statement or answers the question. 1. Find.
Algebra 2 Matrices Review Multiple Choice Identify the choice that best completes the statement or answers the question. 1. Find. Evaluate the determinant of the matrix. 2. 3. A matrix contains 48 elements.
More informationLecture 12: Solving Systems of Linear Equations by Gaussian Elimination
Lecture 12: Solving Systems of Linear Equations by Gaussian Elimination Winfried Just, Ohio University September 22, 2017 Review: The coefficient matrix Consider a system of m linear equations in n variables.
More informationSomething that can have different values at different times. A variable is usually represented by a letter in algebraic expressions.
Lesson Objectives: Students will be able to define, recognize and use the following terms in the context of polynomials: o Constant o Variable o Monomial o Binomial o Trinomial o Polynomial o Numerical
More informationThe Area of Intersecting Ellipses
The Area of Intersecting Ellipses David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationQuaternion Algebra and Calculus
Quaternion Algebra and Calculus David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To
More informationLeast Squares Fitting of Data by Linear or Quadratic Structures
Least Squares Fitting of Data by Linear or Quadratic Structures David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution
More informationLinear Algebra Section 2.6 : LU Decomposition Section 2.7 : Permutations and transposes Wednesday, February 13th Math 301 Week #4
Linear Algebra Section. : LU Decomposition Section. : Permutations and transposes Wednesday, February 1th Math 01 Week # 1 The LU Decomposition We learned last time that we can factor a invertible matrix
More informationHow to (Almost) Square a Circle
How to (Almost) Square a Circle Moti Ben-Ari Department of Science Teaching Weizmann Institute of Science http://www.weizmann.ac.il/sci-tea/benari/ c 2017 by Moti Ben-Ari. This work is licensed under the
More informationCalculating determinants for larger matrices
Day 26 Calculating determinants for larger matrices We now proceed to define det A for n n matrices A As before, we are looking for a function of A that satisfies the product formula det(ab) = det A det
More informationEigenvalue and Eigenvector Homework
Eigenvalue and Eigenvector Homework Olena Bormashenko November 4, 2 For each of the matrices A below, do the following:. Find the characteristic polynomial of A, and use it to find all the eigenvalues
More informationAn Approximation for the Inverse Square Root Function
An Approximation for the Inverse Square Root Function David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International
More informationPrincipal Curvatures of Surfaces
Principal Curvatures of Surfaces David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationDistance from a Point to an Ellipse, an Ellipsoid, or a Hyperellipsoid
Distance from a Point to an Ellipse, an Ellipsoid, or a Hyperellipsoid David Eberly, Geometric Tools, Redmond WA 9805 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution
More informationMA 138 Calculus 2 with Life Science Applications Matrices (Section 9.2)
MA 38 Calculus 2 with Life Science Applications Matrices (Section 92) Alberto Corso albertocorso@ukyedu Department of Mathematics University of Kentucky Friday, March 3, 207 Identity Matrix and Inverse
More informationSystem of Linear Equations
Math 20F Linear Algebra Lecture 2 1 System of Linear Equations Slide 1 Definition 1 Fix a set of numbers a ij, b i, where i = 1,, m and j = 1,, n A system of m linear equations in n variables x j, is given
More informationDynamic Collision Detection using Oriented Bounding Boxes
Dynamic Collision Detection using Oriented Bounding Boxes David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0
More informationDefinition of Equality of Matrices. Example 1: Equality of Matrices. Consider the four matrices
IT 131: Mathematics for Science Lecture Notes 3 Source: Larson, Edwards, Falvo (2009): Elementary Linear Algebra, Sixth Edition. Matrices 2.1 Operations with Matrices This section and the next introduce
More informationLinear Algebra and Matrix Inversion
Jim Lambers MAT 46/56 Spring Semester 29- Lecture 2 Notes These notes correspond to Section 63 in the text Linear Algebra and Matrix Inversion Vector Spaces and Linear Transformations Matrices are much
More informationChapter 4 Systems of Linear Equations; Matrices
Chapter 4 Systems of Linear Equations; Matrices Section 5 Inverse of a Square Matrix Learning Objectives for Section 4.5 Inverse of a Square Matrix The student will be able to identify identity matrices
More informationGEOMETRY OF MATRICES x 1
GEOMETRY OF MATRICES. SPACES OF VECTORS.. Definition of R n. The space R n consists of all column vectors with n components. The components are real numbers... Representation of Vectors in R n.... R. The
More informationB-Spline Interpolation on Lattices
B-Spline Interpolation on Lattices David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More informationElementary Matrices. MATH 322, Linear Algebra I. J. Robert Buchanan. Spring Department of Mathematics
Elementary Matrices MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Outline Today s discussion will focus on: elementary matrices and their properties, using elementary
More informationLecture 1 Systems of Linear Equations and Matrices
Lecture 1 Systems of Linear Equations and Matrices Math 19620 Outline of Course Linear Equations and Matrices Linear Transformations, Inverses Bases, Linear Independence, Subspaces Abstract Vector Spaces
More informationThree-Dimensional Coordinate Systems. Three-Dimensional Coordinate Systems. Three-Dimensional Coordinate Systems. Three-Dimensional Coordinate Systems
To locate a point in a plane, two numbers are necessary. We know that any point in the plane can be represented as an ordered pair (a, b) of real numbers, where a is the x-coordinate and b is the y-coordinate.
More informationA Review of Matrix Analysis
Matrix Notation Part Matrix Operations Matrices are simply rectangular arrays of quantities Each quantity in the array is called an element of the matrix and an element can be either a numerical value
More informationchapter 12 MORE MATRIX ALGEBRA 12.1 Systems of Linear Equations GOALS
chapter MORE MATRIX ALGEBRA GOALS In Chapter we studied matrix operations and the algebra of sets and logic. We also made note of the strong resemblance of matrix algebra to elementary algebra. The reader
More informationMath 1314 Week #14 Notes
Math 3 Week # Notes Section 5.: A system of equations consists of two or more equations. A solution to a system of equations is a point that satisfies all the equations in the system. In this chapter,
More informationLecture 1 INF-MAT3350/ : Some Tridiagonal Matrix Problems
Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems Tom Lyche University of Oslo Norway Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems p.1/33 Plan for the day 1. Notation
More informationThere is a unique function s(x) that has the required properties. It turns out to also satisfy
Numerical Analysis Grinshpan Natural Cubic Spline Let,, n be given nodes (strictly increasing) and let y,, y n be given values (arbitrary) Our goal is to produce a function s() with the following properties:
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 informationLesson U2.1 Study Guide
Lesson U2.1 Study Guide Sunday, June 3, 2018 2:05 PM Matrix operations, The Inverse of a Matrix and Matrix Factorization Reading: Lay, Sections 2.1, 2.2, 2.3 and 2.5 (about 24 pages). MyMathLab: Lesson
More informationSolving Consistent Linear Systems
Solving Consistent Linear Systems Matrix Notation An augmented matrix of a system consists of the coefficient matrix with an added column containing the constants from the right sides of the equations.
More informationLINEAR SYSTEMS, MATRICES, AND VECTORS
ELEMENTARY LINEAR ALGEBRA WORKBOOK CREATED BY SHANNON MARTIN MYERS LINEAR SYSTEMS, MATRICES, AND VECTORS Now that I ve been teaching Linear Algebra for a few years, I thought it would be great to integrate
More information7.5 Operations with Matrices. Copyright Cengage Learning. All rights reserved.
7.5 Operations with Matrices Copyright Cengage Learning. All rights reserved. What You Should Learn Decide whether two matrices are equal. Add and subtract matrices and multiply matrices by scalars. Multiply
More informationHani Mehrpouyan, California State University, Bakersfield. Signals and Systems
Hani Mehrpouyan, Department of Electrical and Computer Engineering, Lecture 26 (LU Factorization) May 30 th, 2013 The material in these lectures is partly taken from the books: Elementary Numerical Analysis,
More information(Linear equations) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Linear equations) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3
ANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3 ISSUED 24 FEBRUARY 2018 1 Gaussian elimination Let A be an (m n)-matrix Consider the following row operations on A (1) Swap the positions any
More informationMaria Cameron Theoretical foundations. Let. be a partition of the interval [a, b].
Maria Cameron 1 Interpolation by spline functions Spline functions yield smooth interpolation curves that are less likely to exhibit the large oscillations characteristic for high degree polynomials Splines
More informationDerivation of the Kalman Filter
Derivation of the Kalman Filter Kai Borre Danish GPS Center, Denmark Block Matrix Identities The key formulas give the inverse of a 2 by 2 block matrix, assuming T is invertible: T U 1 L M. (1) V W N P
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 informationLinear Algebra. The analysis of many models in the social sciences reduces to the study of systems of equations.
POLI 7 - Mathematical and Statistical Foundations Prof S Saiegh Fall Lecture Notes - Class 4 October 4, Linear Algebra The analysis of many models in the social sciences reduces to the study of systems
More informationCS123 INTRODUCTION TO COMPUTER GRAPHICS. Linear Algebra 1/33
Linear Algebra 1/33 Vectors A vector is a magnitude and a direction Magnitude = v Direction Also known as norm, length Represented by unit vectors (vectors with a length of 1 that point along distinct
More informationMATH10212 Linear Algebra B Homework Week 5
MATH Linear Algebra B Homework Week 5 Students are strongly advised to acquire a copy of the Textbook: D C Lay Linear Algebra its Applications Pearson 6 (or other editions) Normally homework assignments
More informationMath Matrix Theory - Spring 2012
Math 440 - Matrix Theory - Spring 202 HW #2 Solutions Which of the following are true? Why? If not true, give an example to show that If true, give your reasoning (a) Inverse of an elementary matrix is
More informationFinish section 3.6 on Determinants and connections to matrix inverses. Use last week's notes. Then if we have time on Tuesday, begin:
Math 225-4 Week 7 notes Sections 4-43 vector space concepts Tues Feb 2 Finish section 36 on Determinants and connections to matrix inverses Use last week's notes Then if we have time on Tuesday, begin
More informationLecture 1: Systems of linear equations and their solutions
Lecture 1: Systems of linear equations and their solutions Course overview Topics to be covered this semester: Systems of linear equations and Gaussian elimination: Solving linear equations and applications
More informationMatrix notation. A nm : n m : size of the matrix. m : no of columns, n: no of rows. Row matrix n=1 [b 1, b 2, b 3,. b m ] Column matrix m=1
Matrix notation A nm : n m : size of the matrix m : no of columns, n: no of rows Row matrix n=1 [b 1, b 2, b 3,. b m ] Column matrix m=1 n = m square matrix Symmetric matrix Upper triangular matrix: matrix
More informationx n -2.5 Definition A list is a list of objects, where multiplicity is allowed, and order matters. For example, as lists
Vectors, Linear Combinations, and Matrix-Vector Mulitiplication In this section, we introduce vectors, linear combinations, and matrix-vector multiplication The rest of the class will involve vectors,
More informationA FIRST COURSE IN LINEAR ALGEBRA. An Open Text by Ken Kuttler. Matrix Arithmetic
A FIRST COURSE IN LINEAR ALGEBRA An Open Text by Ken Kuttler Matrix Arithmetic Lecture Notes by Karen Seyffarth Adapted by LYRYX SERVICE COURSE SOLUTION Attribution-NonCommercial-ShareAlike (CC BY-NC-SA)
More informationPhysics 6303 Lecture 3 August 27, 2018
Physics 6303 Lecture 3 August 27, 208 LAST TIME: Vector operators, divergence, curl, examples of line integrals and surface integrals, divergence theorem, Stokes theorem, index notation, Kronecker delta,
More informationLS.1 Review of Linear Algebra
LS. LINEAR SYSTEMS LS.1 Review of Linear Algebra In these notes, we will investigate a way of handling a linear system of ODE s directly, instead of using elimination to reduce it to a single higher-order
More informationConceptual Questions for Review
Conceptual Questions for Review Chapter 1 1.1 Which vectors are linear combinations of v = (3, 1) and w = (4, 3)? 1.2 Compare the dot product of v = (3, 1) and w = (4, 3) to the product of their lengths.
More information