More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms

Size: px
Start display at page:

Download "More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms"

Transcription

1 More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing, Cina b epartment of Psycology, McGill University, Montréal, Québec, Canada bstract. Necessary and sufficient conditions are derived for a 2-by-2 partitioned matrix to ave {1}-, {1,2}-, {1,3}-, {1,4}-inverses and te Moore-Penrose inverse wit anaciewicz-scur forms. s applications, te anaciewicz-scur forms of {1}-, {1,2}-, {1,3}-, {1,4}-inverses and te Moore-Penrose inverse of a 2-by-2 partitioned Hermitian matrix are also given. Keywords: anaciewicz-scur form; partitioned matrix; generalized inverse; Moore-Penrose inverse; Hermitian matrix; Scur complement; matrix rank metod Matematics Subject Classifications: 1503; Introduction Trougout tis paper, C m n stands for te set of all m n matrices over te field of complex numbers. Te symbols, r() and R() stand for te conjugate transpose, te rank and te range (column space) of a matrix C m n, respectively;, denotes a row block matrix consisting of and. Te Moore-Penrose inverse of C m n, denoted by, is defined to be te unique matrix X C n m satisfying te following four matrix equations (1) X = (2) XX = X, (3) (X) = X, (4) (X) = X. Furter, let E = I m and F = I n stand for te two ortogonal projectors. matrix X is called an {i,..., j}-inverse of, denoted by (i,...,j), if it satisfies te i,..., jt equations. Te collection of all {i,..., j}-inverses of is denoted by { (i,...,j) }. Some frequently used generalized inverses of are (1), (1,2), (1,3) and (1,4). Let M be a 2 2 block matrix M =, (1.1) C were C m n, C m k, C C l n and C l k. If in (1.1) is square and nonsingular, ten M can be decomposed as Im 0 0 Im M = 1 C 1 I l 0 C 1. (1.2) 0 I l Tis decomposition is often called itken block-diagonalization formula in te literature, see Puntanen and Styan 9. Moreover, if bot M and are nonsingular, ten te Scur complement S = C 1 is nonsingular too, and te inverse of M can be written in te following form M 1 Im = I m 0 0 I l 0 S 1 C 1 I l = S 1 C 1 1 S 1 S 1 C 1 S 1. (1.3) Tis well-known formula is called te anaciewicz inversion formula for te inverse of a nonsingular matrix in te literature, see Puntanen and Styan 9, and can be found in most linear algebra books. Te two formulas in (1.2) and (1.3) and teir consequences are widely used in manipulating partitioned matrices and teir operations. Wen bot and M in (1.1) are singular, te s: yongge@mail.sufe.edu.cn, takane@psyc.mcgill.ca 1

2 two formulas in (1.2) and (1.3) can be extended to generalized inverses of matrices. reasonable extension of (1.3) wit generalized inverses of submatrices in (1.1) is given by N( (i,...,j), S (i,...,j) In (i,...,j) (i,...,j) 0 I m 0 0 I k 0 S (i,...,j) C (i,...,j) I l = (i,...,j) + (i,...,j) S (i,...,j) C (i,...,j) (i,...,j) S (i,...,j) S (i,...,j) C (i,...,j) S (i,...,j), (1.4) were S = C (i,...,j). Eq. (1.4) is called te anaciewicz-scur form induced from M, see aksalary and Styan 1. It can be seen from (1.4) tat te matrix N( (i,...,j), S (i,...,j) ) varies over te coice of (i,...,j) and S (i,...,j). Let {N( (i,...,j), S (i,...,j) )} denote te collection of all N( (i,...,j), S (i,...,j) ). ltoug te rigt-and side of (1.4) is obtained by replacing inverses wit generalized inverses, it is not necessarily an {i,..., j}-inverse of M. In tis case, it is of interest to investigate relations between generalized inverses of M in (1.1) and te matrix N( (i,...,j), S (i,...,j) ) in (1.4), in particular, to derive necessary and sufficient conditions for N( (i,...,j), S (i,...,j) ) to be generalized inverses of M in (1.1). some autors ave investigated te relations between M (i,...,j) and N( (i,...,j), S (i,...,j) ) for some special coices of {i,..., j}. well-known result asserts tat M = N(, S ) R() R(), R(C ) R( ), R(C) R(S) and R( ) R(S ), (1.5) were S = C ; see, e.g., 1. Oter results can be found, e.g., in 2, 4, 5, 8, 15, 16. In a recent paper 16, we considered relations between M (1) and N( (1), S (1) ) troug te matrix rank metod and obtain te following two rank formulas r M (1) N( (1), S (1) ) (1), M (1) { } = max r(m) r() r C,, r(m) r() r, 0, (1.6) max r M (1) N( (1), S (1) ) (1) M (1) = r r r r C, 2r(). (1.7) 0 C C y setting te rigt-and sides of te two rank equalities to zero, we obtain tat (a) Tere exist (1) and S (1) suc tat N( (1), S (1) ) is a {1}-inverse of M if and only if { r(m) r() + r C,, } r() + r (1.8) olds. (b) Te set inclusion {N( (1), S (1) )} {M (1) } olds if and only if 0 r = r() + r C, and r 0 0 = r() +r 0 C C, (1.9) old, or equivalently, 0 C R (E ) R 0 and R R CF old. 2

3 s an extension of te previous investigation, we derive in tis paper some rank formulas for te difference M (i,...,j) N( (i,...,j), S (i,...,j) ) (1.10) for {1,2}-, {1,3}-, {1,4}-inverses and te Moore-Penroses of matrices, and use te rank formulas to caracterize te following relations {N( (i,...,j), S (i,...,j) )} {M (i,...,j) }, (1.11) {N( (i,...,j), S (i,...,j) )} {M (i,...,j) }. (1.12) In order to establis rank equalities associated wit (1.10), we need a variety of rank formulas for partitioned matrices and generalized Scur complements. Lemma 1.1 (7) Let C m n, C m k, C C l n and C l k. Ten r, = r() + r( (1) r() + r( (1) ), (1.13) r = r() + r( C C (1) r(c) + r( C (1) C ), (1.14) C r = r() + r(c) + r (I C 0 m (1) )(I n C (1) C), (1.15) 0 r = r() + r (1) C C C (1) C (1). (1.16) Lemma 1.2 (12, 14) Let M be as given in (1.1). Ten r( C (1) r() + r C, + r + r (1) C max (1) { r( C (1) r C,, r, r C 0 0 r r 0, 0 C C (1.17) } r(), (1.18) r( C (1,2) r + r C, + r() + max{ r (1,2) 1, r 2 }, (1.19) { } max r( C (1,2) r() + r(), r C,, r, r r(), (1.20) (1,2) C r( C (1,3) r + r r 0 0, (1.21) (1,3) C C { } max r( C (1,3) r r(), r, (1.22) (1,3) C r( C (1,4) 0 r C, + r (1,4) C r, (1.23) 0 C { } max r( C (1,4) r C,, r (1,4) C r(), (1.24) r( C r C r(), (1.25) were 0 r 1 = r r r 0 0 C 0 C C 0, r 2 = r() r C r 0. 3

4 Te following lemma is derived from Lemma 1.2 by setting and C to identity matrices in (1.17), (1.19), (1.21) and (1.23). Lemma 1.3 Let C m n and C n m. Ten (1) r( (1) r( ), (1.26) (1,2) r( (1,2) max{ r( ), r() + r() r() r() }, (1.27) (1,3) r( (1,3) r( ), (1.28) (1,4) r( (1,4) r( ). (1.29) 2 Generalized inverses of partitioned matrices wit anaciewicz-scur forms We first give two rank formulas for te difference M (1,2) N( (1,2), S (1,2) ). Teorem 2.1 Let M and N( (1,2), S (1,2) ) be as given in (1.1) and (1.4), respectively, were S = C (1,2). Ten were (1,2), M (1,2) r M (1,2) N( (1,2), S (1,2) ) = max{r 1, r 2, r 3, 0}, (2.1) max r M (1,2) N( (1,2), S (1,2) ) = { s 1, s 2 }, (2.2) (1,2) M (1,2) r 1 = r(m) 2r() r(), r 2 = r(m) r() r C,, r 3 = r(m) r() r, s 1 = r r 0 C C 0 s 2 = r + r C 0 r r C, 2r(), + r(m) 2r() r() r C, r. Hence, (a) 16 Tere exist (1,2) and S (1,2) suc tat N( (1,2), S (1,2) ) is a {1, 2}-inverse of M if and only if { } r(m) 2r() + r(), r() + r, r() + r C,. (2.3) (b) Te set inclusion {N( (1,2), S (1,2) )} {M (1,2) } olds if and only if (1.9) olds, or 0 r(m r + r C, r() and r = r = r() + r() (2.4) C 0 old. Proof It was sown in 16 tat r M (1,2) N( (1,2), S (1,2) ) = r(m) r() r( C (1,2) ), (2.5) M (1,2) 4

5 so tat (1,2),M (1,2) r M (1,2) N( (1,2), S (1,2) ) = r(m) r() max (1,2) r( C (1,2) ), (2.6) max r M (1,2) N( (1,2), S (1,2) ) = r(m) r() r( C (1,2) ). (2.7) (1,2) M (1,2) (1,2) Substituting (1.19) and (1.20) into (2.6) and (2.7) gives (2.1) and (2.2). Setting te rigt-and side of (2.1) to zero leads to r 1 0, r 2 0 and r 3 0, tat is, (2.3) olds. Setting te rigt-and side of (2.2) to zero leads to s 1 0 or s 2 0. It can be seen from (1.7) tat s 1 0 is equivalent to (1.9). Note tat s 2 in (2.2) can be rewritten as a sum of tree parts ( ) ( ) 0 s 2 = + r() r() r C ( r(m) + r() r C, r + r 0 ), r() r() were eac part is nonnegative. In tis case, Setting s 2 = 0 leads to (2.4). Teorem 2.2 Let M and N( (1,3), S (1,3) ) be as given in (1.1) and (1.4). Ten were Hence, (1,3), M (1,3) r M (1,3) N( (1,3), S (1,3) ) = max{ r 1, r 2 }, (2.8) max r M (1,3) N( (1,3), S (1,3) ) = r 1 + r 3, (2.9) (1,3) M (1,3) r 1 = r, + r C, r, C r 2 = r, + r C, r() r, r 3 = r 0 0 r r(). C (a) 16 Tere exist (1,3) and S (1,3) suc tat N( (1,3), S (1,3) ) is a {1, 3}-inverse of M if and only if R() R(), R C R = {0} and r C, r old. (b) Te set inclusion {N( (1,3), S (1,3) )} {M (1,3) } olds if and only if R() R(), R C R = {0} and R(CF ) R(F ) old. Proof Te following formula r M (1,3) N( (1,3), S (1,3) ) = r, + r C, r() r( C (1,3) ) (2.10) M (1,3) was sown in 16. Substituting (1.21) and (1.22) into (2.10) gives (2.8) and (2.9). Setting te rigtand side of (2.8) to zero, we see tat tere exist (1,3) and S (1,3) suc tat N( (1,3), S (1,3) ) {M (1,3) } if and only if r, + r C, =r and r, + r C, r() + r. (2.11) C 5

6 Note tat r, + r C, r() + r C, r C Hence, te first equality in (2.11) is equivalent to R() R() and r = r, + r C,, C and te second inequality in (2.11) is equivalent to r C, r. Setting te rigt-and side of (2.9) to zero and noting tat te two terms on te rigt-and side of (2.9) are nonnegative, we obtain (b). Te following teorem can be sown similarly. Teorem 2.3 Let M and N( (1,4), S (1,4) ) be as given in (1.1) and (1.4). Ten. were (1,4), M (1,4) r M (1,4) N( (1,4), S (1,4) ) = max{ r 1, r 2 }, max r M (1,4) N( (1,4), S (1,4) ) = r 1 + r 3, (1,4) M (1,4) r 1 = r + r r C C, r 2 = r + r r() r C,, C 0 r 3 = r r C, r(). 0 C Hence, (a) 16 Tere exist (1,4) and S (1,4) suc tat N( (1,4), S (1,4) ) is a {1, 4}-inverse of M if and only if R(C ) R( ), R C R = {0} and r r C, old. (b) Te set inclusion {N( (1,4), S (1,4) )} {M (1,4) } olds if and only if R(C ) R( ), R C R = {0} and R( E ) R( E C ) old. special case of (1.4) corresponding to te Moore-Penrose inverse is given by N(, S + S C S S C S, (2.12) were S = C. Te relations between N(, S ) and {i,..., j}-inverse of M are given in te following teorems. Teorem 2.4 Let M and N(, S ) be as given in (1.1) and (2.12), respectively. Ten r M (1) N(, S ) = r M (1,2) N(, S ) = r(m) r M (1) M (1,2) C. (2.13) Hence, te following statements are equivalent: 6

7 (a) N(, S ) is a {1}-inverse of M. (b) N(, S ) is a {1, 2}-inverse of M. (c) r(m r() + r( C ). (d) r(m r C. (e) R C = R(M) and R C = R(M ). Proof It follows from (1.26) and (1.27) tat M (1) r M (1) N(, S ) = r M MN(, S )M, (2.14) M (1,2) r M (1,2) N(, S ) = max{ r M MN(, S )M, It is also easy to verify tat rn(, S ) + r(m) rmn(, S ) rn(, S )M }. (2.15) rmn(, S rn(, S )M = r() + r(s), r M MN(, S )M = r(m) r() r(s), were S = C. Substituting tese two equalities and (1.25) into (2.14) and (2.15) gives (2.13). Te equivalence of (a), (b), (c) and (d) follows from (2.13). Recall a simple fact tat r(p Q r() R( P R( ) and R(Q R(). pplying tis result to (d) gives te equivalence of (d) and (e). Teorem 2.5 Let M and N(, S ) be as given in (1.1) and (2.12), respectively. Ten r M (1,3) N(, S ) = r, + r C, r M (1,3) C. (2.16) Hence, te following statements are equivalent: (a) N(, S ) is a {1, 3}-inverse of M. (b) r C = r, + r C,. (c) r C = r, + r C,, r C, = r C, and r, = r,. (d) R C R = {0}, R() R() and R C, = R C,. Proof It follows from (1.28) tat It is easy to verify tat r M (1,3) N(, S ) = r M MN(, S ) M. (2.17) M (1,3) MN(, S E S C SS Im 0 C, I l M MN(, S ) M = M I m+l MN(, S ) = M E + E S C E S E S C E S. 7

8 Recall tat elementary block matrix operations (EMOs) don t cange te rank of a matrix. Hence we can derive by elementary block matrix operations tat ( ) r M MN(, S ) M = r M E + E S C E S E S C E S ( ) = r M E 0 0 E S E E = r E S C E S 0 = r r() r(s) (by (1.13)) 0 S C 0 0 = r r() r(s) 0 C 0 = r, + r C, r() r( C ). Tus we ave (2.16) by (1.25). lso note tat r C r, + r C, r, + r C,. pplying tis inequality to (b) leads to te equivalence of (b) and (c). Te equivalence of (c) and (d) is obvious. Te following result can be sown similarly. Teorem 2.6 Let M and N(, S ) be as given in (1.1) and (2.12), respectively. Ten r M (1,4) N(, S ) = r + r r M (1,4) C C. Hence, te following statements are equivalent: (a) N(, S ) is a {1, 4}-inverse of M. (b) r C = r + r. C (c) r C = r C + r (d) R C R, r C = r and r C = {0}, R, = R, and R(C ) R( ). = r. 3 Generalized inverses of partitioned Hermitian matrices Let M be an Hermitian matrix, and partition M as M =, (3.1) were = C m m, C m n and = C n n. Te anaciewicz-scur form induced from M is N( (i,...,j), S (i,...,j) (i,...,j) + (i,...,j) S (i,...,j) (i,...,j) (i,...,j) S (i,...,j) S (i,...,j) (i,...,j) S (i,...,j), (3.2) were S = (i,...,j). pplying te results in Section 2 to (3.1) and (3.2) gives te following results. 8

9 Teorem 3.1 Let M and N( (1), S (1) ) be as given in (3.1) and (3.2). Ten Hence, r M (1) N( (1), S (1) ) = max {r(m) r() r,, 0}, (1), M (1) max r M (1) N( (1), S (1) 0 ) = 2r (1) M (1) 0 2r() 2r,. (a) Tere exist (1) and S (1) suc tat N( (1), S (1) ) is a {1}-inverse of M if and only if r(m) r() + r,. (b) Te set inclusion {N( (1), S (1) )} {M (1) } olds if and only if 0 r 0 = r() + r,. Proof It follows from (1.6) and (1.7) by setting C =. Teorem 3.2 Let M and N( (1,2), S (1,2) ) be as given in (3.1) and (3.2), were S = (1,2). Ten were Hence, r M (1,2) N( (1,2), S (1,2) ) = max{r 1, r 2, 0}, (1,2), M (1,2) max r M (1,2) N( (1,2), S (1,2) ) = { s 1, s 2 }, (1,2) M (1,2) r 1 = r(m) 2r() r(), r 2 = r(m) r() r,, 0 s 1 = 2r 0 s 2 = 2r 0 2r, 2r(), + r(m) 2r() r() 2r,. (a) Tere exist (1,2) and S (1,2) suc tat N( (1,2), S (1,2) ) is a {1, 2}-inverse of M if and only if r(m) { 2r() + r(), r() + r, }. (3.3) (b) Te set inclusion {N( (1,2), S (1,2) )} {M (1,2) } olds if and only if (3.3) olds, or r(m 2r r() and r 0 = r() + r(). Proof It follows from Teorem 2.1 by setting C =. Teorem 3.3 Let M and N( (1,3), S (1,3) ) be as given in (3.1) and (3.2). Ten r M (1,3) N( (1,3), S (1,3) ) = max{ r 1, r 2 }, (1,3), M (1,3) max r M (1,3) N( (1,3), S (1,3) ) = r 1 + r 3, (1,3) M (1,3) 9

10 were r 1 = r, + r 2, r r 2 = r, r(), r 3 = r 0 0 r r()., Hence, (a) Tere exist (1,3) and S (1,3) suc tat N( (1,3), S (1,3) ) is {1, 3}-inverse of M if and only if 2 R() R() and R R = {0}. (b) Te set inclusion {N( (1,3), S (1,3) )} {M (1,3) } olds if and only if 2 R() R() and R R = {0}. Proof It follows from Teorem 2.2 by setting C =. special case of (3.2) corresponding to te Moore-Penrose inverse is N(, S + S S S S, (3.4) were S =. Te relations between N(, S ) and {i,..., j}-inverse of M are given in te following teorems. Teorem 3.4 Let M and N(, S ) be as given in (3.1) and (3.4), respectively. Ten r M (1) N(, S ) = r M (1,2) N(, S 3 ) = r(m) r M (1) M (1,2). Hence, te following statements are equivalent: (a) N(, S ) is a {1}-inverse of M. (b) N(, S ) is a {1, 2}-inverse of M. (c) r(m r() + r( ). 3 (d) r(m r. 2 (e) R = R(M). Proof It follows from Teorem 2.4 by setting C =. Teorem 3.5 Let M and N(, S ) be as given in (3.1) and (3.4), respectively. Ten r M (1,3) N(, S ) = r M (1,4) N(, S 3 ) = r M (1,3) M (1,4) + r r Hence, te following statements are equivalent: (a) N(, S ) is a {1, 3}-inverse of M.. 10

11 (b) N(, S ) is a {1, 4}-inverse of M. 3 (c) r = r + r. 3 (d) R R = {0}, R() R() and R, = R,. Proof It follows from Teorems 2.5 and 2.6 by setting C =. ssume te Hermitian matrix in (3.1) is nonnegative definite, tat is, tere exists a matrix U suc tat M = UU. In tis case, R() R() and R( ) R() and S = (1) = for any (1) old. Relations between {1}- and {1, 2}-inverses of M and N( (1), S (1) ) and N( (1,2), S (1,2) ) were considered by Rode 11. pplying Teorems 3.1, 3.2, 3.3 and (1.5) to te nonnegative Hermitian matrix in (3.1) gives te following result. Teorem 3.6 Let M and N( (1), S (1) ) be as given in (3.1) and (3.2), and assume M is nonnegative definite. lso let S =. Ten, (a) 11 Te set inclusion {N( (1), S (1) )} {M (1) } always olds. (b) 11 Te set inclusion {N( (1,2), S (1,2) )} {M (1,2) } always olds. (c) Te following statements are equivalent: (i) Te set inclusion {N( (1,3), S (1,3) )} {M (1,3) } olds. (ii) Te set inclusion {N( (1,4), S (1,4) )} {M (1,4) } olds. (iii) M = N(, S ). (iv) R R = {0}. n {i,..., j}-inverse of a square matrix is said to be an Hermitian {i,..., j}-inverse of and is denoted by (i,...,j), if it Hermitian. It is easy to verify tat any Hermitian matrix always as an Hermitian {i,..., j}-inverse for any give set {i,..., j}. Furter, te Hermitian anaciewicz- Scur form induced from te Hermitian matrix M in (3.1) is defined to be (i,...,j) + (i,...,j) S (i,...,j) (i,...,j) (i,...,j) S (i,...,j) N( (i,...,j), S (i,...,j) S (i,...,j) (i,...,j) S (i,...,j), (3.5) were S = (i,...,j). In order to caracterize relations between M and N( (i,...,j), S (i,...,j) ), we need to know te extremal ranks of (i,...,j) wit respect to (i,...,j), wic now are open problems. 4 Generalized inverses of a bordered Hermitian matrix Setting = 0 in (3.1) gives M =, (4.1) 0 were C m m is nonnegative definite and C m n. Tis matrix occurs widely in various problems in matrix teory, in particular, in regression analysis. Te anaciewicz-scur form induced from M is N( (i,...,j), S (i,...,j) (i,...,j) + (i,...,j) S (i,...,j) (i,...,j) (i,...,j) S (i,...,j) S (i,...,j) (i,...,j) S (i,...,j), (4.2) were S = (i,...,j). pplying te results in Section 2 to (4.1) and (4.2) gives us te following results. 11

12 Teorem 4.1 Let M and N( (1), S (1) ) be as given in (4.1) and (4.2). Ten te following statements are equivalent: (a) Tere exist (1) and S (1) suc tat N( (1), S (1) ) {M (1) }. (b) Te set inclusion {N( (1), S (1) )} {M (1) } olds. (c) Tere exist (1,3) and S (1,3) suc tat N( (1,3), S (1,3) ) {M (1,3) }. (d) Te set inclusion {N( (1,3), S (1,3) )} {M (1,3) } olds. (e) Tere exist (1,4) and S (1,4) suc tat N( (1,4), S (1,4) ) {M (1,4) }. (f) Te set inclusion {N( (1,4), S (1,4) )} {M (1,4) } olds. (g) M = N(, S ). () R() R(). Teorem 4.2 Let M and N( (1), S (1) ) be as given in (4.1) and (4.2). Ten te following statements are equivalent: (a) Tere exist (1,2) and S (1,2) suc tat N( (1,2), S (1,2) ) {M (1,2) }. (b) Te set inclusion {N( (1,2), S (1,2) )} {M (1,2) } olds. (c) r, 2r() r() or R() R(). 5 Concluding remarks If in (1.1) is square and nonsingular, ten te matrix in (1.1) can also be decomposed as Im M = 1 1 C 0 Im 0 0 I l 0 1. (5.1) C I l If bot M and are nonsingular in (1.1), ten te Scur complement T = 1 C is nonsingular, too, and te inverse of M can also be written as M 1 T = 1 T CT CT 1 1. (5.2) y symmetry, anoter type of anaciewicz-scur form induced from M is given by K( (i,...,j), T (i,...,j) T (i,...,j) T (i,...,j) (i,...,j) (i,...,j) CT (i,...,j) (i,...,j) + (i,...,j) CT (i,...,j) (i,...,j), (5.3) were T = (i,...,j) C. ltoug (1.3) and (5.2) are identical, (1.4) and (5.3) are not necessarily te same. pplying te results in te previous sections to (5.3), we can derive various conclusions on relations between M and K( (i,...,j), T (i,...,j) ). Furtermore, it is of interest to give necessary and sufficient conditions for te following equality N( (i,...,j), S (i,...,j) K( (i,...,j), T (i,...,j) ) (5.4) to old, or equalities of submatrices in tem to old. ll te results in tis paper on partitioned matrices and anaciewicz-scur forms induced from te matrices can be used to furter study various problems related to block matrices and teir generalized inverses. cknowledgements. We are grateful to anonymous referees for teir elpful comments and suggestions on an earlier version of tis paper. 12

13 References 1 J.K. aksalary and G.P.H. Styan. Generalized inverses of partitioned matrices in anaciewicz-scur form. Linear lgebra ppl. 354(2002) en-israel. note on partitioned matrices and equations. SIM Rev. 11(1969), en-israel and T.N.E. Greville. Generalized Inverses: Teory and pplications, 2nd Ed., Springer- Verlag, New York, P. imasankaram. On generalized inverses of partitioned matrices. Sankyā Ser. 33(1971), F. urns,. Carlson, E. Haynswort and T. Markam. Generalized inverse formulas using te Scur complement. SIM. J. ppl. Mat. 26(1974), S.L. Campbell and C.. Meyer. Generalized Inverses of Linear Transformations, Corrected reprint of te 1979 original, over Publications, Inc., New York, G. Marsaglia and G.P.H. Styan. Equalities and inequalities for ranks of matrices. Linear and Multilinear lgebra 2(1974), G. Marsaglia and G.P.H. Styan. Rank conditions for generalized inverses of partitioned matrices. Sankyā Ser. 36(1974), S. Puntanen and G.P.H. Styan. Scur complements in statistics and probabiliy. In: Te Scur Complement and Its pplications (F. Zang ed.), Springer, pp , C.R. Rao and S.K. Mitra. Generalized Inverse of Matrices and Its pplications. Wiley, New York, C.. Rode. Generalized inverses of partitioned matrices. SIM J. ppl. Mat. 13(1965), Y. Tian. Upper and lower bounds for ranks of matrix expressions using generalized inverses. Linear lgebra ppl. 355(2002), Y. Tian. Rank equalities for block matrices and teir Moore-Penrose inverses. Houston J. Mat. 30(2004), Y. Tian. More on maximal and imal ranks of Scur complements wit applications. ppl. Mat. Comput. 152(2004), Y. Tian. Eigt expressions for generalized inverses of a bordered matrix. Linear and Multilinear lgebra, in press. 16 Y. Tian and Y. Takane. Scur complements and anaciewicz-scur forms. Electron. J. Linear lgebra 13(2005),

ELA

ELA Electronic Journal of Linear Algebra ISSN 181-81 A publication of te International Linear Algebra Society ttp://mat.tecnion.ac.il/iic/ela RANK AND INERTIA OF SUBMATRICES OF THE MOORE PENROSE INVERSE OF

More information

Rank and inertia of submatrices of the Moore- Penrose inverse of a Hermitian matrix

Rank and inertia of submatrices of the Moore- Penrose inverse of a Hermitian matrix Electronic Journal of Linear Algebra Volume 2 Volume 2 (21) Article 17 21 Rank and inertia of submatrices of te Moore- Penrose inverse of a Hermitian matrix Yongge Tian yongge.tian@gmail.com Follow tis

More information

SPECIAL FORMS OF GENERALIZED INVERSES OF ROW BLOCK MATRICES YONGGE TIAN

SPECIAL FORMS OF GENERALIZED INVERSES OF ROW BLOCK MATRICES YONGGE TIAN Electronic Journal of Linear lgebra ISSN 1081-3810 publication of the International Linear lgebra Society EL SPECIL FORMS OF GENERLIZED INVERSES OF ROW BLOCK MTRICES YONGGE TIN bstract. Given a row block

More information

On V-orthogonal projectors associated with a semi-norm

On V-orthogonal projectors associated with a semi-norm On V-orthogonal projectors associated with a semi-norm Short Title: V-orthogonal projectors Yongge Tian a, Yoshio Takane b a School of Economics, Shanghai University of Finance and Economics, Shanghai

More information

A revisit to a reverse-order law for generalized inverses of a matrix product and its variations

A revisit to a reverse-order law for generalized inverses of a matrix product and its variations A revisit to a reverse-order law for generalized inverses of a matrix product and its variations Yongge Tian CEMA, Central University of Finance and Economics, Beijing 100081, China Abstract. For a pair

More information

Symmetry Labeling of Molecular Energies

Symmetry Labeling of Molecular Energies Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry

More information

On convexity of polynomial paths and generalized majorizations

On convexity of polynomial paths and generalized majorizations On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal

More information

Optimization problems on the rank and inertia of the Hermitian matrix expression A BX (BX) with applications

Optimization problems on the rank and inertia of the Hermitian matrix expression A BX (BX) with applications Optimization problems on the rank and inertia of the Hermitian matrix expression A BX (BX) with applications Yongge Tian China Economics and Management Academy, Central University of Finance and Economics,

More information

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of

More information

Rank and inertia optimizations of two Hermitian quadratic matrix functions subject to restrictions with applications

Rank and inertia optimizations of two Hermitian quadratic matrix functions subject to restrictions with applications Rank and inertia optimizations of two Hermitian quadratic matrix functions subject to restrictions with applications Yongge Tian a, Ying Li b,c a China Economics and Management Academy, Central University

More information

Solutions of a constrained Hermitian matrix-valued function optimization problem with applications

Solutions of a constrained Hermitian matrix-valued function optimization problem with applications Solutions of a constrained Hermitian matrix-valued function optimization problem with applications Yongge Tian CEMA, Central University of Finance and Economics, Beijing 181, China Abstract. Let f(x) =

More information

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these. Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra

More information

Continuity and Differentiability Worksheet

Continuity and Differentiability Worksheet Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;

More information

1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist

1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter

More information

Analytical formulas for calculating the extremal ranks and inertias of A + BXB when X is a fixed-rank Hermitian matrix

Analytical formulas for calculating the extremal ranks and inertias of A + BXB when X is a fixed-rank Hermitian matrix Analytical formulas for calculating the extremal ranks and inertias of A + BXB when X is a fixed-rank Hermitian matrix Yongge Tian CEMA, Central University of Finance and Economics, Beijing 100081, China

More information

MA455 Manifolds Solutions 1 May 2008

MA455 Manifolds Solutions 1 May 2008 MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using

More information

Continuity. Example 1

Continuity. Example 1 Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)

More information

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note

More information

A = h w (1) Error Analysis Physics 141

A = h w (1) Error Analysis Physics 141 Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.

More information

232 Calculus and Structures

232 Calculus and Structures 3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE

More information

ch (for some fixed positive number c) reaching c

ch (for some fixed positive number c) reaching c GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan

More information

Yongge Tian. China Economics and Management Academy, Central University of Finance and Economics, Beijing , China

Yongge Tian. China Economics and Management Academy, Central University of Finance and Economics, Beijing , China On global optimizations of the rank and inertia of the matrix function A 1 B 1 XB 1 subject to a pair of matrix equations B 2 XB 2, B XB = A 2, A Yongge Tian China Economics and Management Academy, Central

More information

Differentiation in higher dimensions

Differentiation in higher dimensions Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends

More information

Function Composition and Chain Rules

Function Composition and Chain Rules Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function

More information

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225 THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:

More information

f a h f a h h lim lim

f a h f a h h lim lim Te Derivative Te derivative of a function f at a (denoted f a) is f a if tis it exists. An alternative way of defining f a is f a x a fa fa fx fa x a Note tat te tangent line to te grap of f at te point

More information

Derivatives of Exponentials

Derivatives of Exponentials mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function

More information

MVT and Rolle s Theorem

MVT and Rolle s Theorem AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state

More information

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example, NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing

More information

INTRODUCTION AND MATHEMATICAL CONCEPTS

INTRODUCTION AND MATHEMATICAL CONCEPTS Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips

More information

Characterization of reducible hexagons and fast decomposition of elementary benzenoid graphs

Characterization of reducible hexagons and fast decomposition of elementary benzenoid graphs Discrete Applied Matematics 156 (2008) 1711 1724 www.elsevier.com/locate/dam Caracterization of reducible exagons and fast decomposition of elementary benzenoid graps Andrej Taranenko, Aleksander Vesel

More information

Math 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0

Math 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0 3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,

More information

Material for Difference Quotient

Material for Difference Quotient Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient

More information

WYSE Academic Challenge 2004 Sectional Mathematics Solution Set

WYSE Academic Challenge 2004 Sectional Mathematics Solution Set WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means

More information

2.11 That s So Derivative

2.11 That s So Derivative 2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point

More information

Generic maximum nullity of a graph

Generic maximum nullity of a graph Generic maximum nullity of a grap Leslie Hogben Bryan Sader Marc 5, 2008 Abstract For a grap G of order n, te maximum nullity of G is defined to be te largest possible nullity over all real symmetric n

More information

A SHORT INTRODUCTION TO BANACH LATTICES AND

A SHORT INTRODUCTION TO BANACH LATTICES AND CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,

More information

Finding and Using Derivative The shortcuts

Finding and Using Derivative The shortcuts Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex

More information

VARIANCE ESTIMATION FOR COMBINED RATIO ESTIMATOR

VARIANCE ESTIMATION FOR COMBINED RATIO ESTIMATOR Sankyā : Te Indian Journal of Statistics 1995, Volume 57, Series B, Pt. 1, pp. 85-92 VARIANCE ESTIMATION FOR COMBINED RATIO ESTIMATOR By SANJAY KUMAR SAXENA Central Soil and Water Conservation Researc

More information

The Column and Row Hilbert Operator Spaces

The Column and Row Hilbert Operator Spaces Te Column and Row Hilbert Operator Spaces Roy M Araiza Department of Matematics Purdue University Abstract Given a Hilbert space H we present te construction and some properties of te column and row Hilbert

More information

Notes on wavefunctions II: momentum wavefunctions

Notes on wavefunctions II: momentum wavefunctions Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles

More information

MTH 119 Pre Calculus I Essex County College Division of Mathematics Sample Review Questions 1 Created April 17, 2007

MTH 119 Pre Calculus I Essex County College Division of Mathematics Sample Review Questions 1 Created April 17, 2007 MTH 9 Pre Calculus I Essex County College Division of Matematics Sample Review Questions Created April 7, 007 At Essex County College you sould be prepared to sow all work clearly and in order, ending

More information

Generalized Principal Pivot Transform

Generalized Principal Pivot Transform Generalized Principal Pivot Transform M. Rajesh Kannan and R. B. Bapat Indian Statistical Institute New Delhi, 110016, India Abstract The generalized principal pivot transform is a generalization of the

More information

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x

More information

Tis article was originally publised in a journal publised by Elsevier, and te attaced copy is provided by Elsevier for te autor s benefit and for te benefit of te autor s institution, for non-commercial

More information

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department

More information

Math 31A Discussion Notes Week 4 October 20 and October 22, 2015

Math 31A Discussion Notes Week 4 October 20 and October 22, 2015 Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes

More information

Tail Conditional Expectations for Extended Exponential Dispersion Models

Tail Conditional Expectations for Extended Exponential Dispersion Models American Researc Journal of Matematics Original Article ISSN 378-704 Volume 1 Issue 4 015 Tail Conditional Expectations for Extended Exponential Dispersion Models Ye (Zoe) Ye Qiang Wu and Don Hong 1 Program

More information

3.4 Worksheet: Proof of the Chain Rule NAME

3.4 Worksheet: Proof of the Chain Rule NAME Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are

More information

64 IX. The Exceptional Lie Algebras

64 IX. The Exceptional Lie Algebras 64 IX. Te Exceptional Lie Algebras IX. Te Exceptional Lie Algebras We ave displayed te four series of classical Lie algebras and teir Dynkin diagrams. How many more simple Lie algebras are tere? Surprisingly,

More information

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra

More information

MAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points

MAT 145. Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points MAT 15 Test #2 Name Solution Guide Type of Calculator Used TI-89 Titanium 100 points Score 100 possible points Use te grap of a function sown ere as you respond to questions 1 to 8. 1. lim f (x) 0 2. lim

More information

Math 312 Lecture Notes Modeling

Math 312 Lecture Notes Modeling Mat 3 Lecture Notes Modeling Warren Weckesser Department of Matematics Colgate University 5 7 January 006 Classifying Matematical Models An Example We consider te following scenario. During a storm, a

More information

A note on the equality of the BLUPs for new observations under two linear models

A note on the equality of the BLUPs for new observations under two linear models ACTA ET COMMENTATIONES UNIVERSITATIS TARTUENSIS DE MATHEMATICA Volume 14, 2010 A note on the equality of the BLUPs for new observations under two linear models Stephen J Haslett and Simo Puntanen Abstract

More information

. If lim. x 2 x 1. f(x+h) f(x)

. If lim. x 2 x 1. f(x+h) f(x) Review of Differential Calculus Wen te value of one variable y is uniquely determined by te value of anoter variable x, ten te relationsip between x and y is described by a function f tat assigns a value

More information

On equality and proportionality of ordinary least squares, weighted least squares and best linear unbiased estimators in the general linear model

On equality and proportionality of ordinary least squares, weighted least squares and best linear unbiased estimators in the general linear model Statistics & Probability Letters 76 (2006) 1265 1272 www.elsevier.com/locate/stapro On equality and proportionality of ordinary least squares, weighted least squares and best linear unbiased estimators

More information

Bob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk

Bob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk Bob Brown Mat 251 Calculus 1 Capter 3, Section 1 Completed 1 Te Tangent Line Problem Te idea of a tangent line first arises in geometry in te context of a circle. But before we jump into a discussion of

More information

Chapter 1 Functions and Graphs. Section 1.5 = = = 4. Check Point Exercises The slope of the line y = 3x+ 1 is 3.

Chapter 1 Functions and Graphs. Section 1.5 = = = 4. Check Point Exercises The slope of the line y = 3x+ 1 is 3. Capter Functions and Graps Section. Ceck Point Exercises. Te slope of te line y x+ is. y y m( x x y ( x ( y ( x+ point-slope y x+ 6 y x+ slope-intercept. a. Write te equation in slope-intercept form: x+

More information

Math Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim

Math Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z

More information

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as

More information

2011 Fermat Contest (Grade 11)

2011 Fermat Contest (Grade 11) Te CENTRE for EDUCATION in MATHEMATICS and COMPUTING 011 Fermat Contest (Grade 11) Tursday, February 4, 011 Solutions 010 Centre for Education in Matematics and Computing 011 Fermat Contest Solutions Page

More information

1. Which one of the following expressions is not equal to all the others? 1 C. 1 D. 25x. 2. Simplify this expression as much as possible.

1. Which one of the following expressions is not equal to all the others? 1 C. 1 D. 25x. 2. Simplify this expression as much as possible. 004 Algebra Pretest answers and scoring Part A. Multiple coice questions. Directions: Circle te letter ( A, B, C, D, or E ) net to te correct answer. points eac, no partial credit. Wic one of te following

More information

Polynomial Functions. Linear Functions. Precalculus: Linear and Quadratic Functions

Polynomial Functions. Linear Functions. Precalculus: Linear and Quadratic Functions Concepts: definition of polynomial functions, linear functions tree representations), transformation of y = x to get y = mx + b, quadratic functions axis of symmetry, vertex, x-intercepts), transformations

More information

Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems

Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for

More information

Introduction to Derivatives

Introduction to Derivatives Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))

More information

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1 Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.

More information

The equalities of ordinary least-squares estimators and best linear unbiased estimators for the restricted linear model

The equalities of ordinary least-squares estimators and best linear unbiased estimators for the restricted linear model The equalities of ordinary least-squares estimators and best linear unbiased estimators for the restricted linear model Yongge Tian a and Douglas P. Wiens b a School of Economics, Shanghai University of

More information

4.2 - Richardson Extrapolation

4.2 - Richardson Extrapolation . - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence

More information

EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS

EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS Statistica Sinica 24 2014, 395-414 doi:ttp://dx.doi.org/10.5705/ss.2012.064 EFFICIENCY OF MODEL-ASSISTED REGRESSION ESTIMATORS IN SAMPLE SURVEYS Jun Sao 1,2 and Seng Wang 3 1 East Cina Normal University,

More information

arxiv: v1 [math.dg] 4 Feb 2015

arxiv: v1 [math.dg] 4 Feb 2015 CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola

More information

Preface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.

Preface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed. Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser

More information

The generalized Schur complement in group inverses and (k + 1)-potent matrices

The generalized Schur complement in group inverses and (k + 1)-potent matrices The generalized Schur complement in group inverses and (k + 1)-potent matrices Julio Benítez Néstor Thome Abstract In this paper, two facts related to the generalized Schur complement are studied. The

More information

Packing polynomials on multidimensional integer sectors

Packing polynomials on multidimensional integer sectors Pacing polynomials on multidimensional integer sectors Luis B Morales IIMAS, Universidad Nacional Autónoma de México, Ciudad de México, 04510, México lbm@unammx Submitted: Jun 3, 015; Accepted: Sep 8,

More information

Exam 1 Review Solutions

Exam 1 Review Solutions Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),

More information

Poisson Equation in Sobolev Spaces

Poisson Equation in Sobolev Spaces Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on

More information

Linear Algebra and its Applications

Linear Algebra and its Applications Linear Algebra and its Applications 433 (2010) 476 482 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa Nonsingularity of the

More information

MATH1151 Calculus Test S1 v2a

MATH1151 Calculus Test S1 v2a MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums

More information

Higher Derivatives. Differentiable Functions

Higher Derivatives. Differentiable Functions Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.

More information

7.1 Using Antiderivatives to find Area

7.1 Using Antiderivatives to find Area 7.1 Using Antiderivatives to find Area Introduction finding te area under te grap of a nonnegative, continuous function f In tis section a formula is obtained for finding te area of te region bounded between

More information

Cointegration in functional autoregressive processes

Cointegration in functional autoregressive processes Dipartimento di Scienze Statistice Sezione di Statistica Economica ed Econometria Massimo Franci Paolo Paruolo Cointegration in functional autoregressive processes DSS Empirical Economics and Econometrics

More information

1 + t5 dt with respect to x. du = 2. dg du = f(u). du dx. dg dx = dg. du du. dg du. dx = 4x3. - page 1 -

1 + t5 dt with respect to x. du = 2. dg du = f(u). du dx. dg dx = dg. du du. dg du. dx = 4x3. - page 1 - Eercise. Find te derivative of g( 3 + t5 dt wit respect to. Solution: Te integrand is f(t + t 5. By FTC, f( + 5. Eercise. Find te derivative of e t2 dt wit respect to. Solution: Te integrand is f(t e t2.

More information

2.8 The Derivative as a Function

2.8 The Derivative as a Function .8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open

More information

3.4 Algebraic Limits. Ex 1) lim. Ex 2)

3.4 Algebraic Limits. Ex 1) lim. Ex 2) Calculus Maimus.4 Algebraic Limits At tis point, you sould be very comfortable finding its bot grapically and numerically wit te elp of your graping calculator. Now it s time to practice finding its witout

More information

1 Lecture 13: The derivative as a function.

1 Lecture 13: The derivative as a function. 1 Lecture 13: Te erivative as a function. 1.1 Outline Definition of te erivative as a function. efinitions of ifferentiability. Power rule, erivative te exponential function Derivative of a sum an a multiple

More information

INTRODUCTION AND MATHEMATICAL CONCEPTS

INTRODUCTION AND MATHEMATICAL CONCEPTS INTODUCTION ND MTHEMTICL CONCEPTS PEVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips of sine,

More information

CHAPTER 2. Unitary similarity and unitary equivalence

CHAPTER 2. Unitary similarity and unitary equivalence CHAPTER 2 Unitary similarity and unitary equivalence 2.0 Introduction yin Capter 1, we made an initial study of similarity of A 2 M n via a general nonsingular matrix S, tat is, te transformation A! S

More information

, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim

, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim Mat 132 Differentiation Formulas Stewart 2.3 So far, we ave seen ow various real-world problems rate of cange and geometric problems tangent lines lead to derivatives. In tis section, we will see ow to

More information

5.1 We will begin this section with the definition of a rational expression. We

5.1 We will begin this section with the definition of a rational expression. We Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go

More information

Pre-Calculus Review Preemptive Strike

Pre-Calculus Review Preemptive Strike Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly

More information

Differential Calculus (The basics) Prepared by Mr. C. Hull

Differential Calculus (The basics) Prepared by Mr. C. Hull Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit

More information

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

Some Results on the Growth Analysis of Entire Functions Using their Maximum Terms and Relative L -orders

Some Results on the Growth Analysis of Entire Functions Using their Maximum Terms and Relative L -orders Journal of Matematical Extension Vol. 10, No. 2, (2016), 59-73 ISSN: 1735-8299 URL: ttp://www.ijmex.com Some Results on te Growt Analysis of Entire Functions Using teir Maximum Terms and Relative L -orders

More information

Rank equalities for idempotent and involutory matrices

Rank equalities for idempotent and involutory matrices Linear Algebra and its Applications 335 (2001) 101 117 www.elsevier.com/locate/laa Rank equalities for idempotent and involutory matrices Yongge Tian a, George P.H. Styan a, b, a Department of Mathematics

More information

Lab 6 Derivatives and Mutant Bacteria

Lab 6 Derivatives and Mutant Bacteria Lab 6 Derivatives and Mutant Bacteria Date: September 27, 20 Assignment Due Date: October 4, 20 Goal: In tis lab you will furter explore te concept of a derivative using R. You will use your knowledge

More information

Section 2.1 The Definition of the Derivative. We are interested in finding the slope of the tangent line at a specific point.

Section 2.1 The Definition of the Derivative. We are interested in finding the slope of the tangent line at a specific point. Popper 6: Review of skills: Find tis difference quotient. f ( x ) f ( x) if f ( x) x Answer coices given in audio on te video. Section.1 Te Definition of te Derivative We are interested in finding te slope

More information

DIGRAPHS FROM POWERS MODULO p

DIGRAPHS FROM POWERS MODULO p DIGRAPHS FROM POWERS MODULO p Caroline Luceta Box 111 GCC, 100 Campus Drive, Grove City PA 1617 USA Eli Miller PO Box 410, Sumneytown, PA 18084 USA Clifford Reiter Department of Matematics, Lafayette College,

More information

LATTICE EXIT MODELS S. GILL WILLIAMSON

LATTICE EXIT MODELS S. GILL WILLIAMSON LATTICE EXIT MODELS S. GILL WILLIAMSON ABSTRACT. We discuss a class of problems wic we call lattice exit models. At one level, tese problems provide undergraduate level exercises in labeling te vertices

More information

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001 Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x

More information

Math 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006

Math 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2

More information

AMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.

AMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error. Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error

More information