THE Hadamard product of two nonnegative matrices and

Size: px
Start display at page:

Download "THE Hadamard product of two nonnegative matrices and"

Transcription

1 IAENG Internatonal Journal of Appled Mathematcs 46:3 IJAM_46_3_5 Some New Bounds for the Hadamard Product of a Nonsngular M-matrx and Its Inverse Zhengge Huang Lgong Wang and Zhong Xu Abstract Some new convergent sequences of the lower bounds of the mnmum egenvalue for the Hadamard product of a nonsngular M-matrx and an nverse M-matrx are gven Numercal examples show that these bounds could reach the true value of the mnmum egenvalue n some cases These bounds n ths paper mprove some exstng ones Index Terms sequences nonsngular M-matrx Hadamard product mnmum egenvalue lower bounds I INTRODUCTION THE Hadamard product of two nonnegatve matrces and the Fan product of two M-matrces are two specal matrces product In mathematcs many problems can be transformed nto computatonal problems of the Hadamard product and the Fan product such as numercal method for solvng volterra ntegral equatons wth a convoluton ernel ] solvng a system of Wener-Hopf Integral Equatons ] and so on The upper bounds for the spectral radus ρa B) of the Hadamard product of two nonnegatve matrces A and B the lower bounds for the mnmum egenvalue τa B ) of the Hadamard product of a M- matrx A and an nverse of M-matrx B are research hotspots n matrx theory researchng In ths paper by constructng sequence of teratons and combnng wth the slls of nequaltes scalng we wll conduct further research n the lower bounds for the mnmum egenvalue τa B ) of the Hadamard product of a M-matrx A and an nverse of M-matrx B In theory we prove that the sequences obtaned n ths paper are more accurate than the exstng ones Some results of the comparson are also consdered To llustrate our results some examples are gven For a postve nteger n N denotes the set N = n throughout The set of all n n real matrces s denoted by R n n and C n n denotes the set of all n n complex matrces Let Z n denote the set of n n real matrces all of whose off-dagonal entres are nonpostve A matrx A = a ) R n n s called an M-matrx 3] f there exst a nonnegatve matrx B and a nonnegatve real number λ such that A = λi B λ ρb) where I s the dentty matrx ρb) s the spectral radus of the matrx B If λ = ρb) then A s a sngular M-matrx; f λ > ρb) then A s a nonsngular M-matrx If C s Manuscrpt receved August 4 05; revsed October 05 Ths wor was supported by the Natonal Natural Scence Foundatons of Chna No773) Zhengge Huang s wth the Department of Appled Mathematcs Northwestern Polytechncal Unversty X an Shaanx 7007 PR Chna e-mal: ZhenggeHuang@malnwpueducn Lgong Wang e-mal: lgwang@nwpueducn) and Zhong Xu e-mal: zhongxu@nwpueducn) are wth the Department of Appled Mathematcs Northwestern Polytechncal Unversty X an Shaanx 7007 PR Chna a prncpal submatrx of A then C s also a nonsngular M-matrx Denote by M n the set of all n n nonsngular M-matrces Let us denote τa) = mnreλ) : λ σa) and σa) denotes the spectrum of A It nown that 4] τa) = ρa ) s a postve real egenvalue of A M n The Hadamard product of two matrces A = a ) n n and B = b ) n n s defned as the matrx A B = a b ) n n If A and B are nonsngular M-matrces then t s proved n 5] that A B s also a nonsngular M-matrx A matrx A s rreducble 6] 7] 8] f there does not exst a permutaton matrx P such that ) P AP T A A = 0 A where A and A are square matrces For α N denote by α the cardnalty of α and α = N α If α β N then Aα β) s the submatrx of A lyng n the rows ndcated by α and the columns ndcated by β In partcular Aα α) s abbrevated to Aα) Assume that Aα) s nonsngular Then A/α = A/Aα) = Aα ) Aα α)aα)] Aα α ) s called the Schur complement of A respect to Aα) 9] Let A = a ) R n n be a strctly row dagonally domnant matrx for any N t = denote R = a d = R a a + a d s = s = max a s ; a + a m = m = max a m ; a r = a a r = max r ; a + a r t = t = max a t ; a + a t u = u = max a u ; a l = max a t a t a + a t l w = a Advance onlne publcaton: 6 August 06)

2 IAENG Internatonal Journal of Appled Mathematcs 46:3 IJAM_46_3_5 w = max w q = mns t a h = max a q a q a + a q h v 0) = a p t) = a + h t) = max v t) = a + a a p t) a v t ) a a p t) ht) a p t) a p t) = max pt) Let A be a nonsngular M-matrx and A = α ) n n be a doubly stochastc matrx Before presentng our results we revew the exstng results that relate to lower bounds of the mnmum egenvalue for the Hadamard product of a nonsngular M-matrx and an nverse M-matrx as follows It was proved n 0] that 0 < τa A ) Fedler and Marham 5] gave a lower bound on τa A ): τa A ) n and conectured that τa A ) n Chen ] Song ] and Yong 3] have ndependently proved ths conecture L et al n 4] and L et al n 5] obtaned the followng results τa A ) mn N τa A ) mn N a s R + s a t R + t respectvely In 6] Cheng et al mproved the results n 4] and 5] showng that τa A ) mn τa A ) mn N a u R + u Recently Chen n 7] mproved the results n 4] and gave the followng result: a α + a α a α a α ) τa A ) mn ) m a α m a α )] In 8] Cheng et al proposed τa A ) ρ J A ) and n 9] L et al presented the followng result: a α + a α a α a α ) ) w a α w a α )] In the sequel Zhao et al 0] mproved the results n 4] 5] 6] and arrved at τa A a p t) R ) max mn N + a s mn N p t) a p t) s + p t) In ] Zhou et al proved the followng result: If A = a ) B = b ) R n n are nonsngular M-matrces and A = α ) n n then b t τb A ) mn N a b In ths paper we exhbt some new lower bounds for τa A ) and τb A ) These bounds mprove the results n 4] 5] 6] 7] 8] 9] 0] ] The rest of ths paper s organzed as follows In Secton II we propose some notatons and lemmas whch are useful n the followng proofs We propose some new lower bounds for τb A ) and τa A ) n Secton III Secton IV s devoted to some numercal experments to show that the advantages and precse of the new convergent sequences of the lower bounds of the mnmum egenvalue for the Hadamard product of a nonsngular M-matrx and an nverse M-matrx Fnally some concludng remars are gven II PRELIMINARIES In ths secton we start wth some notatons and lemmas that nvolve nequaltes for the entres of A and the strctly dagonally domnant matrx They wll be useful n the proofs Let A = a ) R n n For N t = denote a + a t l g 0) = = w a a + a g t ) f t) = a l t) = max g t) = a + a f t) a a f t) lt) a f t) a f t) = max t) f Lemma 7] Let A = a ) n n be a nonsngular M- matrx B = b ) n n Z n and A B Then B s a nonsngular M-matrx and A B O Lemma 7] Let A = a ) n n α N and assume that Aα) s nonsngular Then det A = det Aα) det A/α Lemma 3 0] If A = a ) R n n s a strctly row dagonally domnant nonsngular M-matrx then for all Advance onlne publcaton: 6 August 06)

3 IAENG Internatonal Journal of Appled Mathematcs 46:3 IJAM_46_3_5 N t = a) > q v 0) p ) v ) p ) v ) p t) v t) 0; b) h 0 h t) 0 Usng the same technque as the proof of Lemma n 0] we can obtan the Lemma 4 Lemma 4 If A = a ) R n n s a strctly row dagonally domnant nonsngular M-matrx then for all N t = a) > t g 0) f ) g ) f ) g ) f t) g t) 0; b) l 0 l t) 0 Lemma 5 0] Let A = a ) R n n be a strctly row dagonally domnant nonsngular M-matrx Then for A = b ) n n N t = 0 we have a + α a a v t) α = p t+) α Lemma 6 Let A = a ) R n n be a strctly row dagonally domnant nonsngular M-matrx Then for A = b ) n n we have α z t+) α z t+) α N t = 0 where = mnp t) f t) zt) = max zt) t = Proof Usng the same technques as the proof of Lemma n ] for N t = 0 we have α a + a a g t) α = f t+) α Then t follows from the above nequalty and Lemma 5 that α mnp t+) f t+) α = z t+) α for N t = 0 Remar By Lemma 3 we fnd that 0 h n vew of q = mns t s and by Lemma 3 for N t = we have a + a q h p t) v 0) = a a + a a = m Furthermore by Lemma 4 0 l then we have a + a t l f t) g 0) = a a + a t = w = u a By Theorem 33 n 4] 5] and the above two nequaltes for N t = we nfer that v 0) m s g 0) = w u t whch results n v 0) m s g 0) = w u t N Moreover nasmuch as = max zt) = maxzt) zt) zt) + zt) n maxp t) pt) pt) + pt) n = max pt) = pt) = max zt) = maxzt) zt) zt) + zt) n maxf t) f t) f t) + f t) n t) = maxf = f t) for N the result of Lemma 6 s sharper than the results of Theorem n 4] Lemma n 5] 6] 7] and Lemma n 9] 0] Lemma 7 Let A = a ) R n n be a strctly row dagonally domnant nonsngular M-matrx Then for A = α ) n n N t = we have a a a a α a a ) where = mnp t) f t) N Proof Let B = A Snce A s a nonsngular M-matrx B 0 Denote A s the submatrx of A obtaned by deletng the -th row and the -th column of A Then A s a nonsngular M-matrx and A daga a a ++ a nn ) = A Thus by Lemma and A beng a nonsngular M-matrx we have A daga a a ++ a nn) ) For N by Lemma 3 we have α = det A det A = det A det A det A/A = By Inequalty ) we deduce that Hence det A/A a a a a det A/A α = det A det A a a a N 3) a Combnng Lemma 7 and AB = I results n n = a α e = a α = a α a α a α = a ) a α α a N 4) a Advance onlne publcaton: 6 August 06)

4 IAENG Internatonal Journal of Appled Mathematcs 46:3 IJAM_46_3_5 By Inequaltes 3) and 4) the result s obtaned Remar Accordng to Remar for N t = we have s m t u w = g 0) pt) whch mples that a a a a a a a a m a a s a Moreover t s easy to get a a p t) a a g 0) = a a w a a u a a t a a a a N So the bounds of Lemma 7 are sharper than the ones of Theorem 5 n 4] Lemma 3 n 5] 6] 7] and Lemma 3 n 0] Lemma 8 ] If A s a doubly stochastc matrx then Ae = e A T e = e where e = ) T Lemma 9 7] Let A = a ) R n n be a nonsngular M-matrx Then 0 τa) a for all N Lemma 0 3] Let A = a ) C n n and x x x n be postve real numbers Then all the egenvalues of A le n the regon n z C : z a z a = x ) a x x ) a x Lemma 7] Let A Z n A s a nonsngular M-matrx f and only f all ts leadng prncpal mnors are postve III MAIN RESULTS In ths secton we exhbt some new lower bounds for τb A ) and τa A ) whch mprove the ones n 4] 5] 6] 7] 8] 9] 0] ] Theorem 3 Let A = a ) B = b ) R n n be nonsngular M-matrces and A = α ) n n Then for t = τb A ) mn α b + α b α b α b ) ) b α b α )] = t 5) Proof It s evdent that the result 5) holds wth equalty for n = Below we assume that n Snce A s a nonsngular M-matrx there exsts a postve dagonal matrx D such that D AD s a strctly row dagonally domnant nonsngular M-matrx and τb A ) = τd B A )D) = τb D AD) ) Therefore for convenence and wthout loss of generalty we assume that A s a strctly row dagonally domnant matrx Now let us dstngush two cases: ) Frst we assume that A and B are rreducble matrces Then for N t = wthout loss of generalty we assume that = max zt) zt) + zt) n = maxp t) pt) pt) + pt) n By the defnton of p t) obvously 0 < zt) < for N t = Let τb A ) = λ By Lemma 0 and Lemma 6 there exsts a par ) of postve ntegers wth such that λ α b λ α b ) α b = ) b α ) b α ) α b ) b α b α ) 6) Snce A and B are nonsngular M-matrces B A s a nonsngular M-matrx as well By Lemma 9 we have 0 λ α b for all N It follows from Inequalty 6) that λ α b )λ α b ) ) b α b α ) 7) Thus Inequalty 7) s equvalent to τb A ) α b + α b α b α b ) that s τb A ) mn ) b α b α )] α b + α b α b α b ) ) b α b α )] ) Now assume that one of A and B s reducble By Lemma we now that all the leadng prncpal mnors of A and B are postve If we denote by P = p ) the n n permutaton matrx wth p = p 3 = = p n n = p n = the remanng p beng zero Then for any chosen postve real number ε suffcently small such that all the leadng prncpal mnors of A εp and B εp are postve t follows that A εp and B εp are rreducble nonsngular M-matrx Now we substtute A εp for A B εp for B respectvely n the prevous case and then lettng ε 0 8) Advance onlne publcaton: 6 August 06)

5 IAENG Internatonal Journal of Appled Mathematcs 46:3 IJAM_46_3_5 the result follows by contnuty Ths completes our proof of Theorem 3 Remar 3 Wthout loss of generalty for t = assume that α b b α α b Thus Inequalty 9) can be wrtten as b α α b α b + b α 9) b α 0) From Inequaltes 5) 0) and by Lemma 0 and Lemma 3 n 0] t follows that α b + α b α b α b ) ) b α )] b α α b + α b ) b α = α b + α b α b α b ) + 4 b α b α ) + 4 α b α b + ) b α ) )] b α α b α b = α b + α b α b α b + = α b b b α = b a Ths means that mn α b + α b ) b α b mn N a b t mn N a b b b α ) ] b ) b α b α b ) mn N b α )] α b α )] b p t) a b So the lower bound n Theorem 3 s better than the lower bounds n Theorem n 0] and Theorem 48 n ] Moreover snce mn mn w for N t = we have α b + α b ) b α α b α b ) b α )] α b + α b α b α b ) ) w b α w b α )] Thus the lower bound n Theorem 3 s better than that n Theorem n 9] Theorem 3 The sequence t t = obtaned from Theorem 3 s monotone ncreasng wth an upper bound τb A ) and consequently s convergent Proof By Lemma 3 and Lemma 4 we have p t) p t+) 0 and f t) f t+) 0 t = So by the defntons of p t) sequence p t) snce and f t) = mnf t) = mnf t) pt) pt) so by the defnton of t = we have and f t) t s easy to see that are monotone decreasng and for t = we have mnf t+) = max p t+) = z t+) zt) = max N zt) zt) zt) + zt) n max N zt+) z t+) zt+) for N for + zt+) n = z t+) whch mples the sequence s monotone decreasng sequence Then t s a monotoncally ncreasng sequence Hence the sequence s convergent Let A = a ) R n n be a nonsngular M-matrx By Lemma 9 we now that f A s a doubly stochastc matrx then A T e = e Ae = e that s a = + a = + a So A s a strctly dagonally domnant matrx by row and by column If B = A accordng to Theorem 3 the followng corollary s establshed Corollary 3 Let A = a ) n n be a nonsngular M- matrx and A = α ) n n be doubly stochastc Then for t = τa A ) mn α a + α a α a α a ) ) a α a α )] = Γ t ) Remar 3 Accordng to Remar for N N and t = we have s m u w = g 0) pt) s m u w = g 0) p t) Advance onlne publcaton: 6 August 06)

6 IAENG Internatonal Journal of Appled Mathematcs 46:3 IJAM_46_3_5 whch results n τa A ) mn mn a α + a α a α a α ) τa A ) mn ) a α a α + a α ) m a α m a α )] a α a α ) a α )] and smlar to the proof of Remar 3 and by Lemma 3 n 6] we have a α + a α a α a α ) ) a α mn a α + a α ) w a α w mn N a u R + u a α )] a α a α ) Furthermore by Theorem 33 n 6] we have a u R mn N + a t R mn u N + t a α )] So the lower bound n Corollary 3 s an mprovement on the lower bounds n Theorem 3 n 7] Corollary 3 n 9] and Theorem 3 n 5] 6] In addton by Theorem 33 n 7] we have mn a b + a b ) m a b m mn N a s R + s a b a b ) a b )] Hence the result of Corollary 3 s sharper than that of Theorem 3 n 4] Theorem 33 Let A = a ) B = b ) R n n be nonsngular M-matrces and A = α ) n n Then for t = τb A ) mn α b + α b s α α b α b ) b ) s α b )] = Ω t ) Proof It s not dffcult to verfy that the result holds wth equalty for n = We next assume that n Snce A s a nonsngular M-matrx there exsts a postve dagonal matrx D such that D AD s a strctly row dagonally domnant nonsngular M-matrx and τb A ) = τd B A )D) = τb D AD) ) So for convenence and wthout loss of generalty we assume that A s a strctly row dagonally domnant matrx ) If A and B are rreducble matrces Then for any N we derve a + a d 0 < s = max < Let τb A ) = λ By Lemma 0 and Lemma 6 there exsts a par ) of postve ntegers wth for t = such that λ α b λ α b s s = s α b α b α a ) s ) s ) b α b ) s α b α ) b ) 3) Snce A and B are nonsngular M-matrces B A s a nonsngular M-matrx Havng n mnd that 0 λ α b for all N from Inequalty 3) we see that whch yelds that λ α b )λ α b ) s α b ) s α τb A ) α b + α b α b α b ) s α and therefore τb A ) b ) s α mnα b + α b b s α b ) 4) b )] α b α b ) ) s α b )] 5) If one of A and B s reducble smlar to the proof of Theorem 3 ) the result s obtaned Remar 33 Smlar to the proof of Remar 3 and by Advance onlne publcaton: 6 August 06)

7 IAENG Internatonal Journal of Appled Mathematcs 46:3 IJAM_46_3_5 Lemma 3 n 5] for t = s gven by τb A ) mn α b + α b s α b s mn N α b α b ) b ) s α a b p t) s b )] Therefore the bound of Theorem 33 s sharper than that n 0] Usng the same method as the proof of Theorem 3 the followng theorem can be deduced Theorem 34 The sequence Ω t t = obtaned from Theorem 33 s monotone ncreasng wth an upper bound τb A ) and consequently s convergent Tang B = A and usng Theorem 33 we can get the followng corollary Corollary 3 Let A = a ) n n be a nonsngular M- matrx and let A = α ) n n be doubly stochastc Then for t = τa A ) mn α a + α a s α a ) s α α a α a ) a )] = Υ t Remar 34 Accordng to Remar 3 and by Lemma 4 n 0] for t = we can conclude that the bound of Corollary 3 s sharper than that n Corollary 4 n 0] Remar 35 The sequence Υ t t = obtaned from Corollary 3 s monotone ncreasng wth an upper bound τa A ) and consequently s convergent Let Φ t = max t Ω t Combnng Theorems 3 and 33 Theorem 35 s easly obtaned Theorem 35 Let A = a ) B = b ) R n n be nonsngular M-matrces and A = α ) n n Then for t = 6) τb A ) Φ t 7) Let Ψ t = maxγ t Υ t By Corollares 3 and 3 Theorem 36 s easly derved Theorem 36 Let A = a ) n n be a nonsngular M- matrx and let A = α ) n n be doubly stochastc Then for t = τa A ) Ψ t 8) IV NUMERICAL EXAMPLES Example 4 Consder the followng two nonsngular M- matrces 9]: A = B = Numercal results are gven n Table I and Table II for the total number of teratons T = 5 In fact τb A ) = We defne RES = T a) τb A ) where T a) stands for the result derved by usng sequence a when the number of teraton s T In Table I for gven matrces A and B we lst the lower bounds of τb A ) calculated by Theorem 48 n ] Theorem n 9] Theorem n 9] Theorem n 0] and Theorem 3 n ths paper For the sequences obtaned by Theorem n 0] and Theorem 3 n ths paper the results are dsplayed for every step of teraton In Table II we present the lower bounds of τb A ) calculated by Theorem 3 n 0] and Theorem 33 n ths paper respectvely For the sequences obtaned by Theorem 3 n 0] and Theorem 33 n ths paper the results are dsplayed for every step of teraton Numercal results n Table I and Table II show that: a) Lower bounds obtaned from Theorem 3 and Theorem 33 are greater than the ones n Theorem n 9] Theorem n 0] Theorem 48 n ] and Theorem 3 n 0] b) Sequences obtaned from Theorems 3 and 33 both are monotone ncreasng c) The sequences obtaned from Theorems 3 and 33 are closer to the true value of τb A ) generally when the number of teratons T s ncreasng d) The sequences obtaned from Theorems 3 and 33 are much closer to the true value of τb A ) than those obtaned from Theorems and 3 n 0] respectvely In Fgure and Fgure we dsplay the RES generated by Theorem n 0] Theorem 3 and Theorem 3 n 0] Theorem 33 respectvely aganst number of teratons for T = 00 From these two fgures we note that the four Advance onlne publcaton: 6 August 06)

8 IAENG Internatonal Journal of Appled Mathematcs 46:3 IJAM_46_3_5 TABLE I THE LOWER BOUNDS OF τb A ) Method t τb A ) Method t τb A ) Theorem 48 n ] 0007 Theorem n 9] Theorem n 9] 07 Theorem n 0] 03 Theorem sequences of τb A ) are convergent whle the RES of the sequence obtaned from Theorem 3 are much less than that obtaned from Theorem n 0] For the sequences obtaned from Theorem 33 and Theorem 3 n 0] we can get the same concluson TABLE II THE LOWER BOUNDS OF τb A ) Method t τb A ) Method t τb A ) Theorem 3 n 0] Theorem Example 4 0] Let A = By Ae = e A T e = e we now that A s a strctly dagonally domnant by row and column Base on A Z n t s easy to see that A s a nonsngular M-matrx and A s doubly stochastc Numercal results are gven n Table III and Table IV for the total number of teratons T = 0 In fact τa A ) = In Table III for gven matrx A we lst the lower bounds of τa A ) calculated by Theorem 3 n 4] 5] 6] Corollary 5 n 8] Theorem 3 n 7] and Corollary 3 n 9] In Table IV we show the lower bounds of τa A ) calculated by Theorem 5 n 0] and Theorem 36 n ths paper respectvely For the sequences obtaned by Theorem 5 n 0] and Theorem 36 n ths paper the results are dsplayed for every step of teraton Numercal results n Table III and Table IV show that: a) Lower bounds obtaned from Theorem 36 are greater than those n Theorem 3 n 4] 5] 6] Corollary 5 n 8] Theorem 3 n 7] Corollary 3 n 9] and Theorem 5 n 0] b) Sequence obtaned from Theorem 36 s monotone ncreasng c) Sequence obtaned from Theorem 36 s convergent to the value whch s close to the true value of τa A ) d) The sequence obtaned from Theorem 36 s much closer to the true value of τa A ) than that obtaned from Theorem 5 n 0] and the sequence obtaned from Theorem 36 approxmates effectvely to the true value of τa A ) so we can estmate τa A ) by Theorem 36 In Fgure 3 we present the RES generated by Theorem 5 n 0] and Theorem 36 respectvely aganst number of teratons for T = 0 where RES s defned as Example 4 From ths fgure we fnd that the two sequences of τa A ) are convergent whle the RES of the sequence obtaned from Theorem 36 are much less than that obtaned from Theorem 5 n 0] TABLE III THE LOWER BOUNDS OF τa A ) Method τa A ) Theorem 3 n 4] 059 Theorem 3 n 5] 045 Theorem 3 n 6] 0447 Corollary 5 n 8] 040 Theorem 3 n 7] 0473 Corollary 3 n 9] Example 43 0] Let A = a ) R n n where a = a = = a nn = a = a 3 = = a n n = a n = and a = 0 elsewhere It s easy to see that A s a nonsngular M-matrx and A s doubly stochastc The results obtaned from Theorem 5 n 0] and Theorem 36 for n = 5 and T = 0 are lsted n Table 5 where T s defned n Example 4 In fact τa A ) = In Table V we show the lower bounds of τa A ) calculated by Theorem 5 n 0] and Theorem 36 n ths paper respectvely For the sequences obtaned by Theorem 5 n 0] and Theorem 36 n ths paper the results are dsplayed for every step of teraton Numercal results n Table V show that the lower bound obtaned from Theorem 36 could reach the true value of τa A ) n some cases Moreover applyng Theorem 36 s faster to reach the convergence value than Theorem 5 n 0] From Table 5 t can be seen that applyng Theorem 36 reaches the true value of τa A ) after only one teraton but applyng Theorem 5 n 0] needs 9 teratons for n = 5 when they reach the true value of τa A ) Advance onlne publcaton: 6 August 06)

9 IAENG Internatonal Journal of Appled Mathematcs 46:3 IJAM_46_3_5 TABLE IV THE LOWER BOUNDS OF τa A ) for τb A ) and τa A ) s also a further wor Method t τa A ) Method t τa A ) Theorem 5 n 0] Theorem TABLE V THE LOWER UPPER OF τa A ) N=5 Method t τa A ) Method t τa A ) Theorem 5 n 0] 0905 Theorem V CONCLUSIONS In ths paper base on the constructng sequence of teratons we have establshed new convergent sequences Φ t and Ψ t t = n Theorem 35 and Theorem 36 whch are more accurate than the exstng ones n 4] 5] 6] 7] 8] 9] 0] ] to approxmate τb A ) and τa A ) respectvely Numercal results gven n Secton IV Tables I-IV) show that the results obtaned by the new convergent sequences Φ t and Ψ t t = are more sharper than the those obtaned by the lower bounds n 4] 5] 6] 7] 8] 9] 0] ] Numercal results also present the feasblty and effectveness of the new convergent sequences when they are used to estmate τb A ) and τa A ) Inasmuch as the calculaton n experments are nvolve n controllng accuracy so there an nterestng problem s how to accurately these bounds can be computed whch may mprove the accuracy of the of these sequences At present t s very dffcult for us to gve the error analyss We wll contnue to study ths problem n the future In addton although the effcency convergent sequences for Φ t and Ψ t are gven they are almost useless for the mplementaton of these new sequences for Φ t and Ψ t when the sze n of matrx s large snce the calculatons of the nverse of the matrces s very expensve How to fnd easer calculated convergent sequences for Φ t and Ψ t s stll a tough tas whch should be further studed n the further Fnally we can fnd that the sequences obtaned by Theorem 3 Theorem 33 Theorem 35 and Theorem 36 are sharper than the exstng ones by theory analyzng and numercal examples also verfy that whereas the sequences are less precse sometmes we can see that the result 006 n Table I obtaned by Theorem 3 and the result 305 n Table II obtaned by Theorem 33 they are not close to the true value Hence fndng more accurate sequences ACKNOWLEDGMENTS Ths wor s supported by the Natonal Natural Scence Foundatons of Chna No773) REFERENCES ] C Q Yang and J H Hou Numercal method for solvng volterra ntegral equatons wth a convoluton ernel IAENG Internatonal Journal of Appled Mathematcs vol 4 no 4 pp ] A T P Naafabad and D Z Kucerovsy On solutons of a system of wener-hopf ntegral equatons IAENG Internatonal Journal of Appled Mathematcs vol 44 no pp ] A Berman and R J Plemmons Nonnegatve Matrces n the Mathematcal Scences SIAM Phladelpha 979 4] R A Horn and C R Johnson Topcs n Matrx Analyss Cambrdge Unversty Press 99 5] M Fedler and T L Marham An nequalty for the Hadamard product of an M-matrx and ts nverse Lnear Algebra Appl vol 0 pp ] A J L Improvng AOR teratve methods for rreducble L- matrces Engneerng Letters vol 9 no pp ] Z Xu Q Lu K Y Zhang and X H An Theory and Applcatons of H-matrces Scence Press 03 8] A J L A New Precondtoned AOR Iteratve Method and Comparson Theorems for Lnear Systems IAENG Internatonal Journal of Appled Mathematcs vol 4 no 3 pp ] J Zhang J Z Lu and G Tu The mproved dsc theorems for the Schur complements of dagonally domnant matrces J Inequal Appl vol 03 no 03 0] M Fedler T L Marham and M Neumann A trace nequalty for M-matrx and the symmetrzablty of a real matrx by a postve dagonal matrx Lnear Algebra Appl vol 7 pp ] S C Chen A lower bound for the mnmum egenvalue of the Hadamard product of matrces Lnear Algebra Appl vol 378 pp ] Y Z Song On an nequalty for the Hadamard product of an M- matrx and ts nverse Lnear Algebra Appl vol 305 pp ] X R Yong Proof of a conecture of Fedler and Marham Lnear Algebra Appl vol 30 pp ] H B L T Z Huang S Q Shen and H L Lower bounds for the mnmum egenvalue of Hadamard product of an M-matrx and ts nverse Lnear Algebra Appl vol 40 pp ] Y T L F B Chen and D F Wang New lower bounds on egenvalue of the Hadamard product of an M-matrx and ts nverse Lnear Algebra Appl vol 430 pp ] G H Cheng Q Tan and Z D Wang Some nequaltes for the mnmum egenvalue of the Hadamard product of an M-matrx and ts nverse J Inequal Appl vol 03 no ] F B Chen New nequaltes for the Hadamard product of an M- matrx and ts nverse J Inequal Appl vol 05 no ] D M Zhou G L Chen G X Wu and X Y Zhang Some nequaltes for the Hadamard product of an M-matrx and ts nverse J Inequal Appl vol 03 no ] Y T L F Wang C Q L and J X Zhao Some new bounds for the mnmum egenvalue of the Hadamard product of an M-matrx and an nverse M-matrx J Inequal Appl vol 03 no ] J X Zhao F Wang and C L Sang Some nequaltes for the mnmum egenvalues of the Hadamard product of an M-matrx and an nverse M-matrx J Inequal Appl vol 05 no 9 05 ] D M Zhou G L Chen G X Wu and X Y Zhang On some new bounds for egenvalues of the Hadamard product and the Fan product of matrces Lnear Algebra Appl vol 438 pp ] X Y Yong On a conecture of Fedler and Maham Lnear Algebra Appl vol 88 pp ] R A Horn and C R Johnson Matrx Analyss Cambrdge Unversty Press 985 Advance onlne publcaton: 6 August 06)

10 RES RES RES IAENG Internatonal Journal of Appled Mathematcs 46:3 IJAM_46_3_ Theorem of 9] * Theorem Number of teratonst) Fg The teratons curves of two sequences obtaned from Theorem n 0] and Theorem 3 respectvely wth T = Theorem 3 n 9] 3 * Theorem Number of teratons) Fg The teratons curves of two sequences obtaned from Theorem 3 n 0] and Theorem 33 respectvely wth T = Theorem 5 n 9] * Theorem Number of teratons Fg 3 The teratons curves of two sequences obtaned from Theorem 5 n 0] and Theorem 36 respectvely wth T = 0 Advance onlne publcaton: 6 August 06)

The lower and upper bounds on Perron root of nonnegative irreducible matrices

The lower and upper bounds on Perron root of nonnegative irreducible matrices Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College

More information

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Perron Vectors of an Irreducible Nonnegative Interval Matrix Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Linear Algebra and its Applications

Linear Algebra and its Applications Lnear Algebra and ts Applcatons 4 (00) 5 56 Contents lsts avalable at ScenceDrect Lnear Algebra and ts Applcatons journal homepage: wwwelsevercom/locate/laa Notes on Hlbert and Cauchy matrces Mroslav Fedler

More information

Discrete Mathematics. Laplacian spectral characterization of some graphs obtained by product operation

Discrete Mathematics. Laplacian spectral characterization of some graphs obtained by product operation Dscrete Mathematcs 31 (01) 1591 1595 Contents lsts avalable at ScVerse ScenceDrect Dscrete Mathematcs journal homepage: www.elsever.com/locate/dsc Laplacan spectral characterzaton of some graphs obtaned

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Bounds for eigenvalues of nonsingular H-tensor

Bounds for eigenvalues of nonsingular H-tensor Electronc Journal of Lnear Algebra Volume 9 Specal volume for Proceedngs of the Internatonal Conference on Lnear Algebra and ts Applcatons dedcated to Professor Ravndra B. Bapat Artcle 05 Bounds for egenvalues

More information

Bounds for Spectral Radius of Various Matrices Associated With Graphs

Bounds for Spectral Radius of Various Matrices Associated With Graphs 45 5 Vol.45, No.5 016 9 AVANCES IN MATHEMATICS (CHINA) Sep., 016 o: 10.11845/sxjz.015015b Bouns for Spectral Raus of Varous Matrces Assocate Wth Graphs CUI Shuyu 1, TIAN Guxan, (1. Xngzh College, Zhejang

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information

Deriving the X-Z Identity from Auxiliary Space Method

Deriving the X-Z Identity from Auxiliary Space Method Dervng the X-Z Identty from Auxlary Space Method Long Chen Department of Mathematcs, Unversty of Calforna at Irvne, Irvne, CA 92697 chenlong@math.uc.edu 1 Iteratve Methods In ths paper we dscuss teratve

More information

Vector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence.

Vector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence. Vector Norms Chapter 7 Iteratve Technques n Matrx Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematcs Unversty of Calforna, Berkeley Math 128B Numercal Analyss Defnton A vector norm

More information

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

General viscosity iterative method for a sequence of quasi-nonexpansive mappings Avalable onlne at www.tjnsa.com J. Nonlnear Sc. Appl. 9 (2016), 5672 5682 Research Artcle General vscosty teratve method for a sequence of quas-nonexpansve mappngs Cuje Zhang, Ynan Wang College of Scence,

More information

e - c o m p a n i o n

e - c o m p a n i o n OPERATIONS RESEARCH http://dxdoorg/0287/opre007ec e - c o m p a n o n ONLY AVAILABLE IN ELECTRONIC FORM 202 INFORMS Electronc Companon Generalzed Quantty Competton for Multple Products and Loss of Effcency

More information

On the spectral norm of r-circulant matrices with the Pell and Pell-Lucas numbers

On the spectral norm of r-circulant matrices with the Pell and Pell-Lucas numbers Türkmen and Gökbaş Journal of Inequaltes and Applcatons (06) 06:65 DOI 086/s3660-06-0997-0 R E S E A R C H Open Access On the spectral norm of r-crculant matrces wth the Pell and Pell-Lucas numbers Ramazan

More information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b Int J Contemp Math Scences, Vol 3, 28, no 17, 819-827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,

More information

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16 STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus

More information

for Linear Systems With Strictly Diagonally Dominant Matrix

for Linear Systems With Strictly Diagonally Dominant Matrix MATHEMATICS OF COMPUTATION, VOLUME 35, NUMBER 152 OCTOBER 1980, PAGES 1269-1273 On an Accelerated Overrelaxaton Iteratve Method for Lnear Systems Wth Strctly Dagonally Domnant Matrx By M. Madalena Martns*

More information

MATH Homework #2

MATH Homework #2 MATH609-601 Homework #2 September 27, 2012 1. Problems Ths contans a set of possble solutons to all problems of HW-2. Be vglant snce typos are possble (and nevtable). (1) Problem 1 (20 pts) For a matrx

More information

CSCE 790S Background Results

CSCE 790S Background Results CSCE 790S Background Results Stephen A. Fenner September 8, 011 Abstract These results are background to the course CSCE 790S/CSCE 790B, Quantum Computaton and Informaton (Sprng 007 and Fall 011). Each

More information

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k) ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

5 The Rational Canonical Form

5 The Rational Canonical Form 5 The Ratonal Canoncal Form Here p s a monc rreducble factor of the mnmum polynomal m T and s not necessarly of degree one Let F p denote the feld constructed earler n the course, consstng of all matrces

More information

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13 CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,

More information

MEM 255 Introduction to Control Systems Review: Basics of Linear Algebra

MEM 255 Introduction to Control Systems Review: Basics of Linear Algebra MEM 255 Introducton to Control Systems Revew: Bascs of Lnear Algebra Harry G. Kwatny Department of Mechancal Engneerng & Mechancs Drexel Unversty Outlne Vectors Matrces MATLAB Advanced Topcs Vectors A

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

Maximizing the number of nonnegative subsets

Maximizing the number of nonnegative subsets Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum

More information

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

Graph Reconstruction by Permutations

Graph Reconstruction by Permutations Graph Reconstructon by Permutatons Perre Ille and Wllam Kocay* Insttut de Mathémathques de Lumny CNRS UMR 6206 163 avenue de Lumny, Case 907 13288 Marselle Cedex 9, France e-mal: lle@ml.unv-mrs.fr Computer

More information

Inexact Newton Methods for Inverse Eigenvalue Problems

Inexact Newton Methods for Inverse Eigenvalue Problems Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.

More information

Cauchy-Schwarz Inequalities Associated with Positive Semidenite. Matrices. Roger A. Horn and Roy Mathias. October 9, 2000

Cauchy-Schwarz Inequalities Associated with Positive Semidenite. Matrices. Roger A. Horn and Roy Mathias. October 9, 2000 Cauchy-Schwarz Inequaltes ssocated wth Postve Semdente Matrces Roger. Horn and Roy Mathas Department of Mathematcal Scences The Johns Hopkns Unversty, altmore, Maryland 228 October 9, 2 Lnear lgebra and

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

More information

HADAMARD PRODUCT VERSIONS OF THE CHEBYSHEV AND KANTOROVICH INEQUALITIES

HADAMARD PRODUCT VERSIONS OF THE CHEBYSHEV AND KANTOROVICH INEQUALITIES HADAMARD PRODUCT VERSIONS OF THE CHEBYSHEV AND KANTOROVICH INEQUALITIES JAGJIT SINGH MATHARU AND JASPAL SINGH AUJLA Department of Mathematcs Natonal Insttute of Technology Jalandhar 144011, Punab, INDIA

More information

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS Journal of Mathematcal Scences: Advances and Applcatons Volume 25, 2014, Pages 1-12 A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS JIA JI, WEN ZHANG and XIAOFEI QI Department of Mathematcs

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

ON THE GENERAL ALGEBRAIC INVERSE EIGENVALUE PROBLEMS. 1. Introduction

ON THE GENERAL ALGEBRAIC INVERSE EIGENVALUE PROBLEMS. 1. Introduction Journal of Computatonal Mathematcs, Vol.22, No.4, 2004, 567 580. ON THE GENERAL ALGEBRAIC INVERSE EIGENVALUE PROBLEMS Yu-ha Zhang (Department of Mathematcs, Shong Unversty, Jnan 250100, Chna) (ICMSEC,

More information

1 GSW Iterative Techniques for y = Ax

1 GSW Iterative Techniques for y = Ax 1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn

More information

Singular Value Decomposition: Theory and Applications

Singular Value Decomposition: Theory and Applications Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real

More information

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity Int. Journal of Math. Analyss, Vol. 6, 212, no. 22, 195-114 Unqueness of Weak Solutons to the 3D Gnzburg- Landau Model for Superconductvty Jshan Fan Department of Appled Mathematcs Nanjng Forestry Unversty

More information

Matrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD

Matrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD Matrx Approxmaton va Samplng, Subspace Embeddng Lecturer: Anup Rao Scrbe: Rashth Sharma, Peng Zhang 0/01/016 1 Solvng Lnear Systems Usng SVD Two applcatons of SVD have been covered so far. Today we loo

More information

arxiv: v1 [math.sp] 3 Nov 2011

arxiv: v1 [math.sp] 3 Nov 2011 ON SOME PROPERTIES OF NONNEGATIVE WEAKLY IRREDUCIBLE TENSORS YUNING YANG AND QINGZHI YANG arxv:1111.0713v1 [math.sp] 3 Nov 2011 Abstract. In ths paper, we manly focus on how to generalze some conclusons

More information

High resolution entropy stable scheme for shallow water equations

High resolution entropy stable scheme for shallow water equations Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal

More information

On a direct solver for linear least squares problems

On a direct solver for linear least squares problems ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear

More information

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars

More information

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system. Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

The L(2, 1)-Labeling on -Product of Graphs

The L(2, 1)-Labeling on -Product of Graphs Annals of Pure and Appled Mathematcs Vol 0, No, 05, 9-39 ISSN: 79-087X (P, 79-0888(onlne Publshed on 7 Aprl 05 wwwresearchmathscorg Annals of The L(, -Labelng on -Product of Graphs P Pradhan and Kamesh

More information

A new construction of 3-separable matrices via an improved decoding of Macula s construction

A new construction of 3-separable matrices via an improved decoding of Macula s construction Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula

More information

Finding Dense Subgraphs in G(n, 1/2)

Finding Dense Subgraphs in G(n, 1/2) Fndng Dense Subgraphs n Gn, 1/ Atsh Das Sarma 1, Amt Deshpande, and Rav Kannan 1 Georga Insttute of Technology,atsh@cc.gatech.edu Mcrosoft Research-Bangalore,amtdesh,annan@mcrosoft.com Abstract. Fndng

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Valuated Binary Tree: A New Approach in Study of Integers

Valuated Binary Tree: A New Approach in Study of Integers Internatonal Journal of Scentfc Innovatve Mathematcal Research (IJSIMR) Volume 4, Issue 3, March 6, PP 63-67 ISS 347-37X (Prnt) & ISS 347-34 (Onlne) wwwarcournalsorg Valuated Bnary Tree: A ew Approach

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

A FORMULA FOR COMPUTING INTEGER POWERS FOR ONE TYPE OF TRIDIAGONAL MATRIX

A FORMULA FOR COMPUTING INTEGER POWERS FOR ONE TYPE OF TRIDIAGONAL MATRIX Hacettepe Journal of Mathematcs and Statstcs Volume 393 0 35 33 FORMUL FOR COMPUTING INTEGER POWERS FOR ONE TYPE OF TRIDIGONL MTRIX H Kıyak I Gürses F Yılmaz and D Bozkurt Receved :08 :009 : ccepted 5

More information

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information

Applied Mathematics Letters

Applied Mathematics Letters Appled Matheatcs Letters 2 (2) 46 5 Contents lsts avalable at ScenceDrect Appled Matheatcs Letters journal hoepage: wwwelseverco/locate/al Calculaton of coeffcents of a cardnal B-splne Gradr V Mlovanovć

More information

Dynamic Systems on Graphs

Dynamic Systems on Graphs Prepared by F.L. Lews Updated: Saturday, February 06, 200 Dynamc Systems on Graphs Control Graphs and Consensus A network s a set of nodes that collaborates to acheve what each cannot acheve alone. A network,

More information

THIS paper considers the following generalized linear

THIS paper considers the following generalized linear Engneerng Letters, 25:3, EL_25_3_06 A New Branch and Reduce Approach for Solvng Generalzed Lnear Fractonal Programmng Yong-Hong Zhang, and Chun-Feng Wang Abstract In ths paper, for solvng generalzed lnear

More information

Algebraic Solutions to Multidimensional Minimax Location Problems with Chebyshev Distance

Algebraic Solutions to Multidimensional Minimax Location Problems with Chebyshev Distance Recent Researches n Appled and Computatonal Mathematcs Algebrac Solutons to Multdmensonal Mnmax Locaton Problems wth Chebyshev Dstance NIKOLAI KRIVULIN Faculty of Mathematcs and Mechancs St Petersburg

More information

Google PageRank with Stochastic Matrix

Google PageRank with Stochastic Matrix Google PageRank wth Stochastc Matrx Md. Sharq, Puranjt Sanyal, Samk Mtra (M.Sc. Applcatons of Mathematcs) Dscrete Tme Markov Chan Let S be a countable set (usually S s a subset of Z or Z d or R or R d

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Overlapping additive and multiplicative Schwarz iterations for H -matrices

Overlapping additive and multiplicative Schwarz iterations for H -matrices Lnear Algebra and ts Applcatons 393 (2004) 91 105 www.elsever.com/locate/laa Overlappng addtve and multplcatve Schwarz teratons for H -matrces Rafael Bru a,1, Francsco Pedroche a, Danel B. Szyld b,,2 a

More information

Norms, Condition Numbers, Eigenvalues and Eigenvectors

Norms, Condition Numbers, Eigenvalues and Eigenvectors Norms, Condton Numbers, Egenvalues and Egenvectors 1 Norms A norm s a measure of the sze of a matrx or a vector For vectors the common norms are: N a 2 = ( x 2 1/2 the Eucldean Norm (1a b 1 = =1 N x (1b

More information

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014)

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014) 0-80: Advanced Optmzaton and Randomzed Methods Lecture : Convex functons (Jan 5, 04) Lecturer: Suvrt Sra Addr: Carnege Mellon Unversty, Sprng 04 Scrbes: Avnava Dubey, Ahmed Hefny Dsclamer: These notes

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

A Simple Research of Divisor Graphs

A Simple Research of Divisor Graphs The 29th Workshop on Combnatoral Mathematcs and Computaton Theory A Smple Research o Dvsor Graphs Yu-png Tsao General Educaton Center Chna Unversty o Technology Tape Tawan yp-tsao@cuteedutw Tape Tawan

More information

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS These are nformal notes whch cover some of the materal whch s not n the course book. The man purpose s to gve a number of nontrval examples

More information

Eigenvalues of Random Graphs

Eigenvalues of Random Graphs Spectral Graph Theory Lecture 2 Egenvalues of Random Graphs Danel A. Spelman November 4, 202 2. Introducton In ths lecture, we consder a random graph on n vertces n whch each edge s chosen to be n the

More information

Remarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence

Remarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence Remarks on the Propertes of a Quas-Fbonacc-lke Polynomal Sequence Brce Merwne LIU Brooklyn Ilan Wenschelbaum Wesleyan Unversty Abstract Consder the Quas-Fbonacc-lke Polynomal Sequence gven by F 0 = 1,

More information

NOTES ON SIMPLIFICATION OF MATRICES

NOTES ON SIMPLIFICATION OF MATRICES NOTES ON SIMPLIFICATION OF MATRICES JONATHAN LUK These notes dscuss how to smplfy an (n n) matrx In partcular, we expand on some of the materal from the textbook (wth some repetton) Part of the exposton

More information

On the correction of the h-index for career length

On the correction of the h-index for career length 1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

arxiv: v2 [math.ca] 24 Sep 2010

arxiv: v2 [math.ca] 24 Sep 2010 A Note on the Weghted Harmonc-Geometrc-Arthmetc Means Inequaltes arxv:0900948v2 [mathca] 24 Sep 200 Gérard Maze, Urs Wagner e-mal: {gmaze,uwagner}@mathuzhch Mathematcs Insttute Unversty of Zürch Wnterthurerstr

More information

SUCCESSIVE MINIMA AND LATTICE POINTS (AFTER HENK, GILLET AND SOULÉ) M(B) := # ( B Z N)

SUCCESSIVE MINIMA AND LATTICE POINTS (AFTER HENK, GILLET AND SOULÉ) M(B) := # ( B Z N) SUCCESSIVE MINIMA AND LATTICE POINTS (AFTER HENK, GILLET AND SOULÉ) S.BOUCKSOM Abstract. The goal of ths note s to present a remarably smple proof, due to Hen, of a result prevously obtaned by Gllet-Soulé,

More information

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41, The greatest common dvsor of two ntegers a and b (not both zero) s the largest nteger whch s a common factor of both a and b. We denote ths number by gcd(a, b), or smply (a, b) when there s no confuson

More information

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS BOUNDEDNESS OF THE IESZ TANSFOM WITH MATIX A WEIGHTS Introducton Let L = L ( n, be the functon space wth norm (ˆ f L = f(x C dx d < For a d d matrx valued functon W : wth W (x postve sem-defnte for all

More information

Min Cut, Fast Cut, Polynomial Identities

Min Cut, Fast Cut, Polynomial Identities Randomzed Algorthms, Summer 016 Mn Cut, Fast Cut, Polynomal Identtes Instructor: Thomas Kesselhem and Kurt Mehlhorn 1 Mn Cuts n Graphs Lecture (5 pages) Throughout ths secton, G = (V, E) s a mult-graph.

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

The Expectation-Maximization Algorithm

The Expectation-Maximization Algorithm The Expectaton-Maxmaton Algorthm Charles Elan elan@cs.ucsd.edu November 16, 2007 Ths chapter explans the EM algorthm at multple levels of generalty. Secton 1 gves the standard hgh-level verson of the algorthm.

More information

On C 0 multi-contractions having a regular dilation

On C 0 multi-contractions having a regular dilation SUDIA MAHEMAICA 170 (3) (2005) On C 0 mult-contractons havng a regular dlaton by Dan Popovc (mşoara) Abstract. Commutng mult-contractons of class C 0 and havng a regular sometrc dlaton are studed. We prove

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

System in Weibull Distribution

System in Weibull Distribution Internatonal Matheatcal Foru 4 9 no. 9 94-95 Relablty Equvalence Factors of a Seres-Parallel Syste n Webull Dstrbuton M. A. El-Dacese Matheatcs Departent Faculty of Scence Tanta Unversty Tanta Egypt eldacese@yahoo.co

More information

Research Article A Generalized Sum-Difference Inequality and Applications to Partial Difference Equations

Research Article A Generalized Sum-Difference Inequality and Applications to Partial Difference Equations Hndaw Publshng Corporaton Advances n Dfference Equatons Volume 008, Artcle ID 695495, pages do:0.55/008/695495 Research Artcle A Generalzed Sum-Dfference Inequalty and Applcatons to Partal Dfference Equatons

More information

form, and they present results of tests comparng the new algorthms wth other methods. Recently, Olschowka & Neumaer [7] ntroduced another dea for choo

form, and they present results of tests comparng the new algorthms wth other methods. Recently, Olschowka & Neumaer [7] ntroduced another dea for choo Scalng and structural condton numbers Arnold Neumaer Insttut fur Mathematk, Unverstat Wen Strudlhofgasse 4, A-1090 Wen, Austra emal: neum@cma.unve.ac.at revsed, August 1996 Abstract. We ntroduce structural

More information

Joint Statistical Meetings - Biopharmaceutical Section

Joint Statistical Meetings - Biopharmaceutical Section Iteratve Ch-Square Test for Equvalence of Multple Treatment Groups Te-Hua Ng*, U.S. Food and Drug Admnstraton 1401 Rockvlle Pke, #200S, HFM-217, Rockvlle, MD 20852-1448 Key Words: Equvalence Testng; Actve

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

More information

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning Journal of Machne Learnng Research 00-9 Submtted /0; Publshed 7/ Erratum: A Generalzed Path Integral Control Approach to Renforcement Learnng Evangelos ATheodorou Jonas Buchl Stefan Schaal Department of

More information

First day August 1, Problems and Solutions

First day August 1, Problems and Solutions FOURTH INTERNATIONAL COMPETITION FOR UNIVERSITY STUDENTS IN MATHEMATICS July 30 August 4, 997, Plovdv, BULGARIA Frst day August, 997 Problems and Solutons Problem. Let {ε n } n= be a sequence of postve

More information

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso Supplement: Proofs and Techncal Detals for The Soluton Path of the Generalzed Lasso Ryan J. Tbshran Jonathan Taylor In ths document we gve supplementary detals to the paper The Soluton Path of the Generalzed

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

Research Article Global Sufficient Optimality Conditions for a Special Cubic Minimization Problem

Research Article Global Sufficient Optimality Conditions for a Special Cubic Minimization Problem Mathematcal Problems n Engneerng Volume 2012, Artcle ID 871741, 16 pages do:10.1155/2012/871741 Research Artcle Global Suffcent Optmalty Condtons for a Specal Cubc Mnmzaton Problem Xaome Zhang, 1 Yanjun

More information

Determinants Containing Powers of Generalized Fibonacci Numbers

Determinants Containing Powers of Generalized Fibonacci Numbers 1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 19 (2016), Artcle 1671 Determnants Contanng Powers of Generalzed Fbonacc Numbers Aram Tangboonduangjt and Thotsaporn Thanatpanonda Mahdol Unversty Internatonal

More information

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS) Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998

More information

TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES

TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES SVANTE JANSON Abstract. We gve explct bounds for the tal probabltes for sums of ndependent geometrc or exponental varables, possbly wth dfferent

More information