Computing The Determinants By Reducing The Orders By Four

Similar documents
Chapter 2. Determinants

a a a a a a a a a a a a a a a a a a a a a a a a In this section, we introduce a general formula for computing determinants.

Chapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY

MATRICES AND VECTORS SPACE

CHAPTER 2d. MATRICES

Chapter 5 Determinants

On a Method to Compute the Determinant of a 4 4 Matrix

INTRODUCTION TO LINEAR ALGEBRA

Matrices 13: determinant properties and rules continued

Introduction To Matrices MCV 4UI Assignment #1

Determinants Chapter 3

Numerical Linear Algebra Assignment 008

Introduction to Determinants. Remarks. Remarks. The determinant applies in the case of square matrices

HW3, Math 307. CSUF. Spring 2007.

ECON 331 Lecture Notes: Ch 4 and Ch 5

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Engineering Analysis ENG 3420 Fall Dan C. Marinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00

Algebra Of Matrices & Determinants

DETERMINANTS. All Mathematical truths are relative and conditional. C.P. STEINMETZ

Matrices and Determinants

Matrices. Elementary Matrix Theory. Definition of a Matrix. Matrix Elements:

Here we study square linear systems and properties of their coefficient matrices as they relate to the solution set of the linear system.

SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics

Geometric Sequences. Geometric Sequence a sequence whose consecutive terms have a common ratio.

The Algebra (al-jabr) of Matrices

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

Contents. Outline. Structured Rank Matrices Lecture 2: The theorem Proofs Examples related to structured ranks References. Structure Transport

308K. 1 Section 3.2. Zelaya Eufemia. 1. Example 1: Multiplication of Matrices: X Y Z R S R S X Y Z. By associativity we have to choices:

DETERMINANTS. All Mathematical truths are relative and conditional. C.P. STEINMETZ

New Expansion and Infinite Series

Matrix Algebra. Matrix Addition, Scalar Multiplication and Transposition. Linear Algebra I 24

Results on Planar Near Rings

Multivariate problems and matrix algebra

Matrix Solution to Linear Equations and Markov Chains

A modern modification of Gjonbalaj-Salihu cornice determinant, transformation to semi-diagonal determinant

Matrix Eigenvalues and Eigenvectors September 13, 2017

Elements of Matrix Algebra

N 0 completions on partial matrices

Matrix & Vector Basic Linear Algebra & Calculus

Lecture Note 9: Orthogonal Reduction

Things to Memorize: A Partial List. January 27, 2017

A basic logarithmic inequality, and the logarithmic mean

MathCity.org Merging man and maths

Theoretical foundations of Gaussian quadrature

A Matrix Algebra Primer

Elementary Linear Algebra

MATHEMATICS FOR MANAGEMENT BBMP1103

Quadratic Forms. Quadratic Forms

DonnishJournals

Arithmetic Mean Derivative Based Midpoint Rule

Definite Integrals. The area under a curve can be approximated by adding up the areas of rectangles = 1 1 +

AN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION

Chapter 14. Matrix Representations of Linear Transformations

Solution to Fredholm Fuzzy Integral Equations with Degenerate Kernel

Operations with Polynomials

Some estimates on the Hermite-Hadamard inequality through quasi-convex functions

September 13 Homework Solutions

MATRIX DEFINITION A matrix is any doubly subscripted array of elements arranged in rows and columns.

set is not closed under matrix [ multiplication, ] and does not form a group.

Path product and inverse M-matrices

4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX. be a real symmetric matrix. ; (where we choose θ π for.

New data structures to reduce data size and search time

Lecture 14: Quadrature

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Section 6.1 INTRO to LAPLACE TRANSFORMS

MTH 5102 Linear Algebra Practice Exam 1 - Solutions Feb. 9, 2016

Integral points on the rational curve

Pavel Rytí. November 22, 2011 Discrete Math Seminar - Simon Fraser University

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

1 Linear Least Squares

Lecture Solution of a System of Linear Equation

arxiv: v2 [math.nt] 2 Feb 2015

Acta Universitatis Carolinae. Mathematica et Physica

The Regulated and Riemann Integrals

MTH3101 Spring 2017 HW Assignment 6: Chap. 5: Sec. 65, #6-8; Sec. 68, #5, 7; Sec. 72, #8; Sec. 73, #5, 6. The due date for this assignment is 4/06/17.

Summary: Method of Separation of Variables

Tutorial Worksheet. 1. Find all solutions to the linear system by following the given steps. x + 2y + 3z = 2 2x + 3y + z = 4.

Rudimentary Matrix Algebra

Modified Crout s method for an LU decomposition of an interval matrix

Each term is formed by adding a constant to the previous term. Geometric progression

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3

DEFINITION OF ASSOCIATIVE OR DIRECT PRODUCT AND ROTATION OF VECTORS

Journal of Inequalities in Pure and Applied Mathematics

Bernoulli Numbers Jeff Morton

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

Pre-Session Review. Part 1: Basic Algebra; Linear Functions and Graphs

WHEN IS A FUNCTION NOT FLAT? 1. Introduction. {e 1 0, x = 0. f(x) =

Self-similarity and symmetries of Pascal s triangles and simplices mod p

Module 6: LINEAR TRANSFORMATIONS

Higher Checklist (Unit 3) Higher Checklist (Unit 3) Vectors

ON A CONVEXITY PROPERTY. 1. Introduction Most general class of convex functions is defined by the inequality

1 Techniques of Integration

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

SCHOOL OF ENGINEERING & BUILT ENVIRONMENT

Lecture 19: Continuous Least Squares Approximation

Operations with Matrices

and that at t = 0 the object is at position 5. Find the position of the object at t = 2.

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

A Bernstein polynomial approach for solution of nonlinear integral equations

Journal of Inequalities in Pure and Applied Mathematics

Transcription:

Applied Mthemtics E-Notes, 10(2010), 151-158 c ISSN 1607-2510 Avilble free t mirror sites of http://wwwmthnthuedutw/ men/ Computing The Determinnts By Reducing The Orders By Four Qefsere Gjonblj, Armend Slihu Received 2 August 2009 Abstrct We present new method to compute the determinnts of n n(n 5) mtrices by reducing their sizes by four To prove our results we use the so-clled cornice determinnts, ie, squre determinnts of order n (n 5) where, with the exception of the first nd lst entries, the entries of the 2nd row nd (n 1)- th row, s well the 2nd column nd (n 1)-th column re ll zero The method introduced here hs the dvntge of reducing the size of determinnt by four, nd thus enbling their quicker nd esier computtion 1 Introduction Let A be n n n mtrix A = 11 12 1n 21 22 2n n1 n2 nn DEFINITION 1 A determinnt of order n, or size n n, (see [5], [9], [10], [11]) is the sum 11 12 1n 21 22 2n D = det(a) = A = = ε j1,j 2,,j n j1 j2 jn, S n n1 n2 nn rnging over the symmetric permuttion group S n, where { +1, ifj 1, j 2,, j n is n even permuttion ε j1,j 2,,j n = 1, ifj 1, j 2,, j n is n odd permution Mthemtics Subject Clssifictions: 15A15, 11C20, 65F40 Deprtment of Mthemtics, Fculty of Electricl nd Computer Engineering, University of Prishtin, Pristin 10000, Kosovo Deprtment of Telecommuniction, Fculty of Electricl nd Computer Engineering, University of Prishtin, Pristin 10000, Kosovo 151

152 A New Method for Computing the Determinnts 11 Properties Chrcterizing Determinnts Let A nd B be ny n n mtrices 1 If A is tringulr mtrix, ie ij = 0 whenever i > j or, lterntively, whenever i < j, then det(a) = 11 22 nn 2 If B results from A by interchnging two rows or columns, then det(b) = det(a) 3 If B results from A by multiplying one row or column with number c, then det(b) = c det(a) 4 If B results from A by dding multiple of one row to nother row, or multiple of one column to nother column, then det(b) = det(a) These four properties (see eg [3], [7], [8]) cn be used to compute the determinnt of ny mtrix, using Gussin elimintion This is n lgorithm tht trnsforms ny given mtrix to tringulr mtrix, only by using the opertions from the lst three items bove Since the effect of these opertions on the determinnt cn be trced, the determinnt of the originl mtrix is known, once Gussin elimintion is performed It is lso possible to expnd determinnt long row or column using Lplce s formul, which is efficient for reltively smll mtrices To do this long the row i, sy, we write n n det(a) = ij C ij = ij ( 1) i+j M ij, j=1 where the C ij represent the mtrix cofctors, ie, C ij is ( 1) i+j times the minor M ij, which is the determinnt of the mtrix tht results from A by removing the i-th row nd the j-th column, nd n is the size of the mtrix 5 If I is identity mtrix, ie ij = j=1 { 1, ifi = j 0, ifi j then the determinnt of the identity mtrix is one, det(i) = 1 6 Multiplying mtrix by number r ffects the determinnt s follows: det(ra) = r n det(a) 7 If A contins zero row (or column), then det(a) = 0 8 If A contins two identicl or proportionl rows (or columns), then det(a) = 0 9 A mtrix nd its trnspose hve the sme determinnt: det(a T ) = det(a), PROPOSITION 1 Let B be obtined from A by one of the following elementry row (column) opertions:

Q Gjonblj nd A Slihu 153 1) Two rows (or columns) of A re switched, or 2) A row (or column) of A is multiplied by number α, or 3) A multiplier of row (or column) of A is dded to nother row (or column) Then we hve, respectively, det(b) = det(a), or det(b) = α det(a), or det(b) = det(a) 12 Chio s Condenstion Method Chio s condenstion is method for evluting n n n determinnt in terms of (n 1) (n 1) determinnts; see [1], [4]: A = 11 12 1n 21 22 2n n1 n2 nn = 1 n 2 11 11 12 21 22 11 12 31 32 11 12 n1 n2 11 13 21 23 11 13 31 33 11 13 n1 n3 11 21 11 31 11 n1 1n 2n 1n 3n 1n nn 13 Dodgson Condenstion Method Dodgson s condenstion method computes determinnts of size n n by expressing them in terms of those of size (n 1) (n 1), nd then expresses the ltter in terms of determinnts of size (n 2) (n 2), nd so on (see [2]) 2 A New Method Let A n n be n n n determinnt: 11 12 13 1n 21 22 23 2n A n n = n1 n2 n3 nn (1) Let us now nlyze determinnts of size n n (n 5), where, with the exception of the first nd lst entries, the entries of the 2nd row nd (n 1)-th row, s well the 2nd column nd (n 1)-th column re ll zero We cll such determinnts cornice

154 A New Method for Computing the Determinnts determinnts nd write them down s C n n (n 5): C n n = 11 12 13 1,n 2 1,n 1 1n 21 0 0 0 0 2n 31 0 33 3,n 2 0 3n n 2,1 0 n 2,3 n 2,n 2 0 n 2,n n 1,1 0 0 0 0 n 1,n n1 n2 n3 n,n 2 n,n 1 nn (2) By using elementry row nd column opertions nd bsed on determinnts properties, considerble number of determinnts of size n n cn be trnsformed into cornice determinnts We will prove tht ny cornice determinnt of size n (n 5) cn be computed by reducing its size by four, thus trnsforming it to n (n 4) (n 4) determinnt For exmple, determinnt of size 5 5 is trnsformed to number, one of size 6 6 to determinnt of size 2 2, etc THEOREM 1 Every cornice determinnt C n n of size n n (n 5) cn be computed by reducing the order of the determinnt by four: C n n = ( 12 21 n,n 1 n 1,n 12 2n n,n 1 n 1,1 21 n2 1,n 1 n 1,n + 1,n 1 2n n2 n 1,1 ) C (n 4) (n 4), (3) where C (n 4) (n 4) = 33 3,n 2 n 2,3 n 2,n 2 We illustrte this using the following picture: n n 1,1 11 21 31 n1 n2 n n3 n 2, n 2 1, n 1 Cn! n "!!! "!!! 0 0 2,1 12 0 0 13 33 2,3 1, n 2 0 0 3, n 2 0 0 0 0 n, n 2 0 0 n, n 1 1n 2n 3n n 2, n n 1, n n, n PROOF Let n = 5 We will prove tht Theorem 1 holds for C 5 5 cornice determinnts Bsed on Lplce s formul, the determinnt C 5 5 cn be expnded long

Q Gjonblj nd A Slihu 155 the first row: C 5 5 is equl to = 11 + 13 + 15 11 12 13 14 15 21 0 0 0 25 31 0 33 0 35 41 0 0 0 45 51 52 53 54 55 0 0 0 25 0 33 0 35 0 0 0 45 52 53 54 55 21 0 0 25 31 0 0 35 41 0 0 45 51 52 54 55 21 0 0 0 31 0 33 0 41 0 0 0 51 52 53 54 21 0 0 25 + ( 1) 1+2 12 31 33 0 35 41 0 0 45 51 53 54 55 21 0 0 25 + +( 1) 1+4 14 31 0 33 35 41 0 0 45 51 52 53 55 Bsed on property 8 from the first section, the first, third nd fifth determinnts re zero, nd hence 33 0 35 C 5 5 = 12 ( 1) 1+1 21 0 0 45 53 54 55 0 33 35 14 ( 1) 1+1 21 0 0 45 52 53 55 + ( 1)1+4 25 + ( 1)1+4 25 31 33 0 41 0 0 51 52 54 31 0 33 41 0 0 51 52 53 = 12 [ 21 ( 33 45 54 ) 25 ( 33 41 54 )] 14 [ 21 ( 33 45 52 ) 25 ( 33 41 52 )] = 12 21 33 45 54 12 25 33 41 54 14 21 33 45 52 + 14 25 33 41 52 = ( 12 21 54 45 12 25 54 41 21 52 14 45 + 41 52 25 14 ) 33, which hs the desired form Next, we prove tht the theorem holds for n 6 Bsed on Lplce s formul, the determinnt C n n cn be expnded long the second column to yield:

156 A New Method for Computing the Determinnts = ( 1) 1+2 12 +( 1) n+2 n2 11 12 13 1,n 2 1,n 1 1n 21 0 0 0 0 2n 31 0 33 3,n 2 0 3n n 2,1 0 n 2,3 n 2,n 2 0 n 2,n n 1,1 0 0 0 0 n 1,n n1 n2 n3 n,n 2 n,n 1 nn 21 0 0 0 2n 31 33 3,n 2 0 3n n 2,1 n 2,3 n 2,n 2 0 n 2,n n 1,1 0 0 0 n 1,n n1 n3 n,n 2 n,n 1 nn 11 13 1,n 2 1,n 1 1n 21 0 0 0 2n 31 33 3,n 2 0 3n n 2,1 n 2,3 n 2,n 2 0 n 2,n n 1,1 0 0 0 n 1,n (4) The determinnts in (4) re expnded long the (n 1)-th column so tht C n n is equl to 21 0 0 2n 31 33 3,n 2 3n ( 1) 1+2 12 ( 1) n 1+n 2 n,n 1 n 2,1 n 2,3 n 2,n 2 n 2,n n 1,1 0 0 n 1,n 21 0 0 2n 31 33 3,n 2 3n +( 1) n+2 n2 ( 1) 1+n 2 1,n 1 n 2,1 n 2,3 n 2,n 2 n 2,n n 1,1 0 0 n 1,n 21 0 0 2n = ( ( 1) 2n 12 n,n 1 + ( 1) 2n+1 ) 31 33 3,n 2 3n n2 1,n 1 n 2,1 n 2,3 n 2,n 2 n 2,n n 1,1 0 0 n 1,n

Q Gjonblj nd A Slihu 157 The determinnt in the right hnd side is expnded long the first row to yield 33 3,n 2 3n ( 12 n,n 1 n2 1,n 1 ) ( 1)1+1 21 + n 2,3 n 2,n 2 n 2,n 0 0 n 1,n 31 33 3,n 2 +( 1) 1+n 2 2n n 2,1 n 2,3 n 2,n 2 (5) n 1,1 0 0 The determinnts in (5) re expnded long the lst row to yield ( 12 n,n 1 n2 1,n 1 ) ( ) 21 ( 1) n 3+n 3 n 1,n + ( 1) n 1 2n ( 1) n 2 n 1,1 33 34 3,n 2 43 44 4,n 2 n 2,3 n 2,4 n 2,n 2 = ( 12 n,n 1 n2 1,n 1 )( 21 n 1,n 2n n 1,1 ) 33 34 3,n 2 43 44 4,n 2 n 2,3 n 2,4 n 2,n 2 = ( 12 21 n,n 1 n 1,n 12 2n n,n 1 n 1,1 21 n2 1,n 1 n 1,n + 1,n 1 2n n2 n 1,1 ) C (n 4) (n 4), which proves tht Theorem 1 holds We remrk tht if mtrix A of size n n (n 5) cn be trnsformed into cornice determinnt C n n, then it cn be computed using the formul (3) To illustrte the usefulness of our method in computing cornice determinnts, we give n exmple: 7 1 3 11 9 6 8 0 0 0 0 2 C 6 6 = 5 0 5 5 0 1 12 0 3 7 0 2 3 0 0 0 0 11 10 9 2 1 5 7 = [8 1 5 11 1 2 3 5 8 9 9 11 + 9 2 3 9] 5 5 3 7 = (440 30 7128 + 486) 50 = 311600

158 A New Method for Computing the Determinnts References [1] F Chió, Mémoire sur les fonctions connues sous le nom de résultntes ou de déterminnts, Turin: E Pons, 1853 [2] C L Dodgson, Condenstion of Determinnts, Being New nd Brief Method for Computing their Arithmetic Vlues, Proc Roy Soc Ser A, 15(1866), 150 155 [3] H Eves, An Introduction to the History of Mthemtics, pges 405 nd 493, Sunders College Publishing, 1990 [4] H Eves, Chio s Expnsion, 36 in Elementry Mtrix Theory, New York: Dover, (1996), 129 136 [5] E Hmiti, Mtemtik 1, Universiteti i Prishtinës: Fkulteti Elektroteknik, Prishtinë, (2000), 163 164 [6] P H Hnus, An Elementry Tretise on the Theory of Determinnts, A textbook for colleges, Ithc, New York: Cornell University Librry, Boston, Ginn nd Compny (1886), 13, 14, 18 [7] J R Bunch nd J E Hopcroft, Tringulr fctoriztion nd inversion by fst mtrix multipliction, Mthemtics of Computtion, 28 (1974), 231 236 [8] Pierre-Simon (de) Lplce, Expnsion of determinnts in terms of minors, Reserches sur le clcul intégrl et sur le systéme du monde, Histoire de l Acdémie Royle des Sciences (Pris), seconde prtie, pges 267 376 (1772) [9] S Brnrd nd J M Child, Higher Algebr, London Mcmilln LTD New York, ST Mrtin s Press (1959), 131 [10] R F Scott, The theory of determinnts nd their pplictions, Ithc, New York: Cornell University Librry, Cmbridge: University Press, (1904), 3 5 [11] W L Ferrr, Algebr, A Text-Book of Determinnts, Mtrices, nd Algebric Forms, Second edition, Fellow nd tutor of Hertford college Oxford, (1957), 7