Linear convergence of an algorithm for computing the largest eigenvalue of a nonnegative tensor

Size: px
Start display at page:

Download "Linear convergence of an algorithm for computing the largest eigenvalue of a nonnegative tensor"

Transcription

1 NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer Lnear Algebra Al 22; 9:83 84 Publshed onlne 8 October 2 n Wley Onlne Lbrary (wleyonlnelbrarycom) OI: 2/nla822 Lnear convergence of an algorthm for comutng the largest egenvalue of a nonnegatve tensor Lng Zhang and Lqun Q 2, *, eartment of Mathematcal Scences, Tsnghua Unversty, Bejng 84, Chna 2 eartment of Aled Mathematcs, The Hong Kong Polytechnc Unversty, Hong Kong SUMMARY An teratve method for fndng the largest egenvalue of a nonnegatve tensor was roosed by Ng, Q, and Zhou n 29 In ths aer, we establsh an exlct lnear convergence rate of the Ng Q Zhou method for essentally ostve tensors Numercal results are gven to demonstrate lnear convergence of the Ng Q Zhou algorthm for essentally ostve tensors Coyrght 2 John Wley & Sons, Ltd Receved 25 Setember 2; Revsed 7 August 2; Acceted 28 August 2 KEY WORS: nonnegatve tensor; essentally ostve tensor; egenvalue; teratve method; convergence INTROUCTION Consder an m-order n-dmensonal tensoraconsstng of n m entres n the real feld <: A a m, a m 2 <, 6, : : :, m 6 n An m-order n-dmensonal tensor A s called nonnegatve (or, resectvely, ostve) f a m > (or, resectvely, a m > ) A tensoras called reducble, f there exsts a nonemty roer ndex subset J f, 2, : : :, ng such that a m, 8 2 J, 8 2, : : :, m 62 J If A s not reducble, then we say that A s rreducble Ths defnton was used n [ 6] In [7, 8], ths roerty s called ndecomosable To an n-dmensonal column vector x x I x 2 I : : : I x n /, real or comlex, and any comlex number, we defne n-dmensonal column vectors Ax and x Œ : Ax a 2 ::: m x 2 : : : x m A, x Œ W xį 2,:::, m 66n 66n Let C be the comlex feld A ar, x/ 2 C C n nfg/ s called an egenvalue egenvector ar of A, f they satsfy: Ax x Œ () *Corresondence to: Lqun Q, eartment of Aled Mathematcs, The Hong Kong Polytechnc Unversty, Hong Kong E-mal: maqlq@olyueduhk Coyrght 2 John Wley & Sons, Ltd

2 LINEAR CONVERGENCE FOR FINING THE LARGEST EIGENVALUE 83 Ths defnton was ntroduced by Q [9] when m s even and A s symmetrc and extended to the general case n [] Indeendently, Lm [] gave such a defnton but restrcted x and to be real Unlke matrces, the egenvalue roblem for tensors s nonlnear Recently, the largest egenvalue roblem for a nonnegatve tensor has attracted much attenton because t has many mortant alcatons such as multlnear agerank [], hyergrahs [2], hgher-order Markov chans [4, 3], and ostve defnteness of a multvarate form [3] Chang, Pearson, and Zhang [] generalzed the Perron Frobenus theorem from nonnegatve matrces to rreducble nonnegatve tensors Yang and Yang [6] generalzed the weak Perron Frobenus theorem to general nonnegatve tensors Bulò and Pelllo [2] gave new bounds on the clque number of grahs on the bass of sectral hyergrah theory The calculaton of these new bounds reles on fndng the largest egenvalue of a f, g nonnegatve tensor Ng, Q, and Zhou [4] roosed an teratve method for fndng the largest egenvalue of an rreducble nonnegatve tensor, whch s an extenson of the Collatz method [4] for calculatng the sectral radus of an rreducble nonnegatve matrx Pearson [5] ntroduced the noton of essentally ostve tensors and roved that the unque ostve egenvalue s real geometrcally smle when the tensor s essentally ostve wth even order Here, real geometrcally smle means that the corresondng real egenvector s unque u to a scalng constant Lu, Zhou, and Ibrahm [3] modfed the Ng Q Zhou method such that the modfed algorthm s always convergent for fndng the largest egenvalue of an rreducble nonnegatve tensor Chang, Pearson, and Zhang [2] ntroduced rmtve tensors An essentally ostve tensor s a rmtve tensor, and a rmtve tensor s an rreducble nonnegatve tensor but not vce versa Chang, Pearson, and Zhang [2] establshed convergence of the Ng Q Zhou method for rmtve tensors Fredland, Gaubert, and Han [7] onted out that the Perron Frobenus theorem for nonnegatve tensors has a very close lnk wth the Perron Frobenus theorem for homogeneous monotone mas, ntated by Nussbaum [5] and further studed by Gaubert and Gunawardena [8] Frendland, Gaubert, and Han [7] ntroduced weakly rreducble nonnegatve tensors and establshed the Perron Frobenus theorem for them enote < n C fx 2 <n W x > g and < n CC fx 2 <n W x > g A ma F W < n C! <n C s called a homogeneous monotone ma f Ftx/ tfx/ for any x 2 < n C and any ostve number t and f Fx/ 6 Fy/ for any x, y 2 < n C wth x 6 y For an m-order n-dmensonal tensor A, defne F W < n C! <n C by Fx/ Ax Œ for any x 2 < n C Then F s a homogeneous monotone ma Hence, n < C < n C, egenvalues and egenvectors of nonnegatve tensors, dscussed here, fall n the class of egenvalues and egenvectors of homogeneous monotone mas Egenvalues and egenvectors of nonnegatve tensors are defned ncc n nfg/ In Corollary 43 of [7], Frendland, Gaubert, and Han showed how to extend the Perron Frobenus result for a nonnegatve tensor from < C < n C toc Cn nfg/ In ths aer, we establsh an exlct lnear convergence rate of the Ng Q Zhou method for essentally ostve tensors Such a result has not aeared n the lterature of egenvalues of homogeneous monotone mas and nonnegatve tensors Ths aer s organzed as follows In Secton 2, we recall some relmnary results In Secton 3, we rewrte the Ng Q Zhou method n the srt of the earler work of Hall and Porschng [6, 7] An exlct lnear convergence rate of the Ng Q Zhou method s establshed for essentally ostve tensors, wth a secfed startng ont, n Secton 4 Fnally, n Secton 5, we reort some numercal results Note that Chang, Pearson, and Zhang [2] showed that lnear convergence rate dd not hold for some rmtve but not essentally ostve tensors for the Ng Q Zhou method Ths shows that our result s shar 2 PRELIMINARIES Frst, we state the Perron Frobenus theorem for nonnegatve tensors gven n [, Theore4] and the mnmax theorem for rreducble nonnegatve tensors gven n [, Theorem 42] These two results, as we stated n the ntroducton, may be derved from the earler results on homogeneous monotone mas Coyrght 2 John Wley & Sons, Ltd Numer Lnear Algebra Al 22; 9:83 84 OI: 2/nla

3 832 L ZHANG AN L QI Theorem 2 If A s an rreducble nonnegatve tensor of order m and dmenson n, then there exst > and x 2 < n CC such that Ax x Œ Moreover, f s an egenvalue wth a nonnegatve egenvector, then If s an egenvalue of A, then jj 6 Theorem 22 Assume that A s an rreducble nonnegatve tensor of order m and dmenson n Then mn max x2< n 66n CC Ax x max mn x2< n 66n CC Ax / x, where s the unque ostve egenvalue corresondng to a ostve egenvector On the bass of Theorem 22, the Ng Q Zhou method resented n [4] works as follows Choose x / 2 < n CC and let y/ A x / For k,, 2, : : :, comute x kc/ yk/ Œ, yk/ Œ y kc/ A x kc/, kc mn x kc/ > y kc/ x kc/, kc max x kc/ > y kc/ x kc/ (2) It s shown n [4] that the obtaned sequences f k g and f k g converge to some numbers and, resectvely, and we have 6 6, where s the largest egenvalue of A, defned n Theorem 2 If, that s, the ga s zero, then both the sequences f k g and f k g converge to However, a ostve ga may haen, whch can be seen n an examle gven n [4] That examle s rreducble but not rmtve Chang, Pearson, and Zhang [2] establshed convergence of the Ng Q Zhou method for rmtve tensors and gave an examle, whch s rmtve but not essentally ostve, such that lnear convergence fals for ths examle wth the Ng Q Zhou method We thus ntend to establsh an exlct lnear convergence rate of the Ng Q Zhou method for essentally ostve tensors For ths urose, we need the followng defnton gven n [5] efnton 2 A nonnegatve m-order n-dmensonal tensor A s essentally ostve f Ax 2 < n CC for any nonzero x 2 < n C By Theorem 32 and efnton 2 n [5], t s easy to obtan the followng result Prooston 2 A nonnegatve m-order n-dmensonal tensor A s essentally ostve f and only f a j :::j > for, j 2 f, 2, : : :, ng For the remander of ths aer, denote I f 2, : : :, m / j 2, : : :, m 2 f, : : :, ngg Coyrght 2 John Wley & Sons, Ltd Numer Lnear Algebra Al 22; 9:83 84 OI: 2/nla

4 LINEAR CONVERGENCE FOR FINING THE LARGEST EIGENVALUE ALGORITHM We frst gve ncluson bounds for the largest egenvalue of an rreducble nonnegatve tensor The followng lemma was gven n [, Lemma 22] Lemma 3 If a nonnegatve tensor A of order m and dmenson n s rreducble, then R W 2,:::, m a 2 ::: m >,, 2, : : :, n On the bass of Theorem 22 and Lemma 3, we obtan the followng ncluson bounds Prooston 3 Let A be an rreducble nonnegatve tensor of order m and dmenson n, and let be the largest egenvalue of A Then mn M 6 6 max M, 66n 66n where M R 2,:::, m m a 2 ::: m R m 2 R m,, 2, : : :, n efne an n-dmensonal column vector R W R / 66n, where R s defned n Lemma 3 By Lemma 3, the vector R Π2 < n CC Takng x RΠnto the equaltes n Theorem 22, we mmedately get the lower and uer bounds Let For, 2, : : :, n, we have R R R 2,:::, m R max 66n R, a 2 ::: m 6 M 6 R R R mn 66n R (3) 2,:::, m a 2 ::: m R Ths shows that the lower and uer bounds gven n Prooston 3 are better than the bounds n [6, Lemma 56] On the bass of Prooston 3 and n the srt of the earler work of Hall and Porschng [6,7], we rewrte the Ng Q Zhou method as follows Algorthm 3 Ste Let A / A and S / np 2,:::, m " > be a suffcently small number and set k W Ste Set A kc/ a kc/ ::: m, S kc/ a 2 ::: m for, : : :, n Let the accurate tolerance 2,:::, m a kc/ 2 ::: m Coyrght 2 John Wley & Sons, Ltd Numer Lnear Algebra Al 22; 9:83 84 OI: 2/nla

5 834 L ZHANG AN L QI Ste 2 Let for, : : :, n, where a kc/ 2 ::: m 2 ::: m 2 m k max S kc/, k mn S kc/ 66n 66n (4) If k k < ", sto Outut the maxmal egenvalue 2 k C k / Otherwse, set k W k C and go to Ste Algorthm 3 s well defned By Lemma 3, for, : : :, n, S / R > and there exsts at least one subndex array 2, : : :, m / 2 I such that a 2 ::: m > Hence, a / 2 ::: m > and S / > By nducton, we have 2 ::: m > and > for k 2, 3, : : : We now rove that Algorthm 3 s actually the Ng Q Zhou method wth a secfed startng ont Prooston 32 Let A be an rreducble nonnegatve tensor of order m and dmenson n Then Algorthm 3 s just the Ng Q Zhou method wth the startng ont x / I I : : : I / 2 < n For k,, : : :, defne an n-dmensonal column vector ky x k/ A S j / j 66n BecauseAs rreducble, by Lemma 3, x k/ 2 < n CC for k,, : : : Set y k/ W xk/ ΠBy the defnton ofa k/ n Ste, we obtan for 6 6 n, A y k/ y k/ x k/ 2,:::, m, k,, : : : (5) xk/ Π2,:::, m a / 2 ::: m 2,:::, m a 2 ::: m x k/ 2 x k/ m Qk j S j / 2 2 ::: m Qk j S j / m Q k j S j / 2 m S kc/, (6) whch concludes that k s just kc n (2) for k,, 2, : : : Ths also holds for the sequence f k g Hence, we conclude the statement of ths rooston Coyrght 2 John Wley & Sons, Ltd Numer Lnear Algebra Al 22; 9:83 84 OI: 2/nla

6 LINEAR CONVERGENCE FOR FINING THE LARGEST EIGENVALUE 835 Numercal results reorted n [3,4] mly that Algorthm 3 s effcent In artcular, Lu, Zhou, and Ibrahm [3] aled a modfcaton of Algorthm 3 to study the ostve defnteness of a multvarate form Testng ostve defnteness of a multvarate form s an mortant roblem n the stablty study of nonlnear autonomous systems va Lyaunov s drect method n automatc control Researchers n automatc control studed the condtons of such ostve defnteness ntensvely [9] For n > 3 and m > 4, ths roblem s a hard roblem n mathematcs There are only a few methods to answer the queston, and these methods are comutatonally exensve when n > 3 [3] Numercal results n [3] show the method of usng the largest egenvalue of nonnegatve tensors s effectve n some cases 4 LINEAR CONVERGENCE By Prooston 32, t suffces to establsh an exlct lnear convergence rate of Algorthm 3 n the case of essentally ostve tensors By some straghtforward comutatons, we mmedately obtan the followng two roostons Prooston 4 For k,, 2, : : :, S kc/ 2,:::, m For, 2, : : :, n, by (4), we have S kc/ Prooston 42 For k,, 2, : : :, 2,:::, m 2 ::: m a kc/ 2 ::: m 2 m 2,:::, m 2 ::: m ::: a :::,, 2, : : :, n j :::j ak/ j ::: a j :::j a j :::, For, 2, : : :, n, t follows from (4) that S / ::: ak ::: k / k / S S k /,, 2, : : :, n 2 m, j 2 f, 2, : : :, ng For, j 2 f, 2, : : :, ng, accordng to (4), by a drect comutaton, k / a ::: a ::: j :::j ak/ j ::: ak / j :::j a k / j :::j k / S j k / S k / k / S j ::: k / S j a ak / j ::: a j :::j a j ::: Coyrght 2 John Wley & Sons, Ltd Numer Lnear Algebra Al 22; 9:83 84 OI: 2/nla

7 836 L ZHANG AN L QI LetAbe an rreducble nonnegatve tensor of order m dmenson n By Theorem 2, there exsts a ostve egenvalue of A, whch s the largest n modulus among all the egenvalues of A Algorthm 3 yelds two sequences of lower and uer bounds for ths largest egenvalue Theorem 4 Let A be an rreducble nonnegatve tensor of order m and dmenson n and be ts largest egenvalue Assume that f k g and f k g are two sequences generated by Algorthm 3 Then R k k R For k,, : : :, let y k/ be defned by (5) Takng the vector y k/ nto the two equaltes n Theorem 22, we obtan from (6), We now rove for any k >, k mn S kc/ 6 6 max S kc/ k, k,, : : : 66n 66n k 6 kc and kc 6 k We assume, wthout loss of generalty, that kc S kc2/ f, 2, : : :, ng We have by Prooston 4, and kc S kc2/ q where, q 2 kc S kc2/ S kc/ 2,:::, m > mn S kc/ 66n k a kc/ 2 ::: m S kc/ 2 S kc/ m Smlarly, we can rove that kc 6 k Ths, together wth (3), comletes our roof Theorem 4 ndcates that f k g and f k g converge, and the sequences f k k g s nonnegatve and monotoncally decreasng Hence, f k k g has a lmt We now show that n the case of essentally ostve tensors, f k k g lnearly converges to zero wth an exlct convergence rate Ths establshes an exlct lnear converge rate of the Ng Q Zhou method for essentally ostve tensors Theorem 42 Let A be a nonnegatve tensor of order m and dmenson n If A s essentally ostve, then k k 6 k k, k, 2, : : :, (7) where W ˇ 2, /, (8) R ˇ W mn,j 2f,2,:::,ng a j :::j, R max 66n R Coyrght 2 John Wley & Sons, Ltd Numer Lnear Algebra Al 22; 9:83 84 OI: 2/nla

8 LINEAR CONVERGENCE FOR FINING THE LARGEST EIGENVALUE 837 and R a 2 ::: m 2,:::, m By Prooston 2, the nonnegatve tensor A s rreducble Wthout loss of generalty, assume that k S kc/ and k kc/ Then by Prooston 4, we have k k S kc/ kc/ 2 ::: m S k/ 2,:::, m! q 2 ::: m k/ 2 m (9) efne Ik/ ( 2, : : :, m / 2 I ˇ 2 ::: m ) > ak/ q 2 ::: m k/ By the defnton of, for,, n, we have 2,:::, m /2Ik/ 2 ::: m C 2 ::: m 2,:::, m /2InIk/ Lettng and q and combnng these two equaltes, we have 2 ::: m! q 2 ::: m S k/ k/ 2,:::, m /2Ik/ 2,:::, m /2InIk/ 2 ::: m! q 2 ::: m () k/ Combnng (9), (), and the defntons of k and k, by Theorem 4, we obtan k k 2,:::, m /2Ik/ C 2,:::, m /2InIk/ 6 k 2,:::, m /2Ik/ 2 ::: m 2 ::: m 2 ::: m k k k 2,:::, m q 2 ::: m k/ 2 m! q 2 ::: m k/ 2 m! q 2 ::: m q 2 ::: m 2,:::, m /2InIk/ 6 k k C 2 6 k k k C 2 R C k q 2 ::: m q 2 ::: m C 2,:::, m /2InIk/! 2 ::: m q 2 ::: m 2,:::, m /2Ik/ k/ AA! q 2 ::: m k/, () Coyrght 2 John Wley & Sons, Ltd Numer Lnear Algebra Al 22; 9:83 84 OI: 2/nla

9 838 L ZHANG AN L QI where 2,:::, m /2InIk/ 2 ::: m, 2 2,:::, m /2Ik/ Because A s essentally ostve, we have from Prooston 2, Because A s rreducble, from (3) and (2), we have q 2 ::: m a j :::j >,, j 2 f, 2, : : :, ng (2) < ˇ < R, (3) where ˇ s defned n the statement of the theorem We now consder the subndex array, : : :, / 2 I If, : : :, / 2 InIk/, then the summaton must nclude a k/ ::: If, : : :, / 2 Ik/, there are two ossbltes: () q, : : :, q/ 2 Ik/ In ths case, the summaton 2 must nclude q q:::q (2) q, : : :, q/ 2 InIk/ In ths case, the summaton must nclude q:::q, and the summaton 2 must nclude q ::: From the dscusson, by Prooston 42, we can obtan n o C 2 > mn a k/ :::, ak/ q q:::q, ak/ q:::q C ak/ q ::: > mn a :::, a q q:::q, 2 a q:::q a q ::: Take > ˇ (4) W ˇ R Then < < follows from (3) Combnng () and (4), we obtan (7) for k, 2, : : : Theorem 42 acheved the target of ths aer, namely, to establsh an exlct lnear convergence rate of the Ng Q Zhou method for essentally ostve tensors In the followng corollary, we gve an exlct asymtotc estmate on the number of teratons Corollary 4 Let and R, R be defned by (8) and (3), resectvely Let " > be the suffcently small number gven n Algorthm 3 IfAs essentally ostve, then Algorthm 3 termnates n at most teratons wth 2 3 " log K R R 6 log / K K < " By (7) n Theorem 42, we have for k, 2, : : :, It follows from (5) and 2, / that log log K/ K log / < log / 7 C (5) k k 6 kr R/ (6) " R R log / " log, R R Coyrght 2 John Wley & Sons, Ltd Numer Lnear Algebra Al 22; 9:83 84 OI: 2/nla

10 LINEAR CONVERGENCE FOR FINING THE LARGEST EIGENVALUE 839 whch yelds K < " R R Ths, together wth (6), mles K K 6 KR R/ < " Ths comletes the roof 5 NUMERICAL EPERIMENTS To demonstrate the lnear convergence of the Ng Q Zhou algorthm (Algorthm 3 ) for essentally ostve tensors, we made numercal exerments on some numercal examles wth the sto rule " 8 We consder the followng fve classes of nonnegatve tensors of order m 3 and dmenson n A W A jj j for, j, 2, : : :, n, and zero elsewhere A2 W A3 W A2 jj A3 jj j for, j, 2, : : :, n, and zero elsewhere for, j, 2, : : :, n, and zero elsewhere Table I Numercal results for Algorthm 3 Tensor menson Eg No Iter Rato (Max, Mn) Rato A , 79, 79 (46, 79) ,, (5, ) e 4, 36e 4, 36e 4 (367e 4, 36e 4) A , 346, 346 (3327, 2736) , 22, 22 (2448, 86) , 76, 76 (2246, 466) A , 667, 667 (78, 472) , 323, 323 (54, 323) , 64, 64 (38, 64) A , 3796, 38 (5635, 3333) , 2588, 2858 (4823, 3) , 36, 33 (447, 877) 5 x n=5 2 5 n=3 5 n= Fgure The curves of rato fora Coyrght 2 John Wley & Sons, Ltd Numer Lnear Algebra Al 22; 9:83 84 OI: 2/nla

11 84 L ZHANG AN L QI A4 W A4 jj for j 2, : : :, n, A4 jj for 2, : : :, n, A5 W j, : : :, n C, and zero elsewhere A5 22 4, A5 2, and zero elsewhere Note that A, A2, and A3 are all essentally ostve A4 s rmtve but not essentally ostve A5 s rreducble but not rmtve and not essentally ostve We aly Algorthm 3 to fnd the 35 3 n= n=3 n= Fgure 2 The curves of rato fora2 8 6 n= n=3 n= Fgure 3 The curves of rato fora n=5 n= n= Fgure 4 The curves of rato fora4 Coyrght 2 John Wley & Sons, Ltd Numer Lnear Algebra Al 22; 9:83 84 OI: 2/nla

12 LINEAR CONVERGENCE FOR FINING THE LARGEST EIGENVALUE 84 largest egenvalues of all the fve tensors wth dfferent dmensons n 5, 3, 6 We summarze the numercal results n the followng table, where NoIter denotes the number of teratons, Eg denotes the largest egenvalue, and Rato denotes the rato of kc kc to k k at the last three teratons, (Max, Mn) Rato denotes the maxmum rato and the mnmum rato We also draw the curve of the convergence rato to descrbe the lnear convergence n the followng fgures For tensor A5, the Ng Q Zhou algorthm s dvergent From Table I and Fgures 3, we see that the Ng Q Zhou algorthm s lnearly convergent for tensors A, A2, and A3 From Table I and Fgure 4, we see that the Ng Q Zhou algorthm s also convergent fora4, but t s not lnearly convergent Ths echoes the results n [2] Hence, numercal examles show that the Ng Q Zhou algorthm s lnearly convergent for essentally ostve tensors, and ths result s shar ACKNOWLEGEMENTS We would lke to thank edtor and referees for ther comments We also thank r Y u and Prof Zhengha Huang for ther hel n erformng the numercal exerments The work of Lng Zhang was suorted by the Natonal Natural Scence Foundaton of Chna (Grant No 873), whereas that of Lqun Q was suorted by the Hong Kong Research Grant Councl REFERENCES Chang KC, Pearson K, Zhang T Perron Frobenus theorem for nonnegatve tensors Communcatons n Mathematcal Scences 28; 6: Chang KC, Pearson K, Zhang T Prmtvty the convergence of the NZQ method, and the largest egenvalue for nonnegatve tensors SIAM Journal on Matrx Analyss and Alcatons 2; 32: Lu Y, Zhou G, Ibrahm NF An always convergent algorthm for the largest egenvalue of an rreducble nonnegatve tensor Journal of Comutatonal and Aled Mathematcs 2; 235: Ng M, Q L, Zhou G Fndng the largest egenvalue of a nonnegatve tensor SIAM Journal on Matrx Analyss and Alcatons 29; 3: Pearson K Essentally ostve tensors Internatonal Journal of Algebra 2; 4: Yang YN, Yang QZ Further results for Perron Frobenus theorem for nonnegatve tensors SIAM Journal on Matrx Analyss and Alcatons 2; 3: Fredland S, Gaubert S, Han L Perron Frobenus theorem for nonnegatve multlnear forms and extensons To aear n: Lnear Algebra and Its Alcatons 8 Gaubert S, Gunawardena J The Perron Frobenus theorem for homogeneous, monotone functons Transactons of the Amercan Mathematcal Socety 24; 356: Q L Egenvalues of a real suersymmetrc tensor Journal of Symbolc Comutaton 25; 4: Lm LH Sngular values and egenvalues of tensors: a varatonal aroach In Proceedngs of the IEEE Internatonal Worksho on Comutatonal Advances n Mult-Sensor Addatve Processng (CAMSAP 5), Vol IEEE Comuter Socety Press: Pscataway, NJ, 25; Lm LH Multlnear agerank: measurng hgher order connectvty n lnked objects, July 25 The Internet: Today and Tomorrow 2 Bulò SR, Pelllo M New bounds on the clque number of grahs based on sectral hyergrah theory In Learnng and Intellgent Otmzaton, Stützle T (ed) Srnger Verlag: Berln, 29; L W, Ng M Exstence and unqueness of statonary robablty vector of a transton robablty tensor March 2 eartment of Mathematcs, The Hong Kong Batst Unversty 4 Collatz L Enschlessungssatz für de charakterstschen Zahlen von Matrzen Mathematk Zetschrft 942; 48: Nussbaum R Convexty and log convexty for the sectral radus Lnear Algebra and Its Alcatons 986; 73: Hall C, Porschng T Comutng the maxmal egenvalue and egenvector of a ostve matrx SIAM Journal on Numercal Analyss 968; 5: Hall C, Porschng T Comutng the maxmal egenvalue and egenvector of a nonnegatve rreducble matrx SIAM Journal on Numercal Analyss 968; 5: Coyrght 2 John Wley & Sons, Ltd Numer Lnear Algebra Al 22; 9:83 84 OI: 2/nla

On the Connectedness of the Solution Set for the Weak Vector Variational Inequality 1

On the Connectedness of the Solution Set for the Weak Vector Variational Inequality 1 Journal of Mathematcal Analyss and Alcatons 260, 15 2001 do:10.1006jmaa.2000.7389, avalable onlne at htt:.dealbrary.com on On the Connectedness of the Soluton Set for the Weak Vector Varatonal Inequalty

More information

The dominant eigenvalue of an essentially nonnegative tensor

The dominant eigenvalue of an essentially nonnegative tensor NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Lnear Algebra Appl. 2013; 00:1 13 Publshed onlne n Wley InterScence (www.nterscence.wley.com). The domnant egenvalue of an essentally nonnegatve tensor

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

Supplementary Material for Spectral Clustering based on the graph p-laplacian

Supplementary Material for Spectral Clustering based on the graph p-laplacian Sulementary Materal for Sectral Clusterng based on the grah -Lalacan Thomas Bühler and Matthas Hen Saarland Unversty, Saarbrücken, Germany {tb,hen}@csun-sbde May 009 Corrected verson, June 00 Abstract

More information

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

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

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

THE Hadamard product of two nonnegative matrices and

THE Hadamard product of two nonnegative matrices and 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

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

Advanced Topics in Optimization. Piecewise Linear Approximation of a Nonlinear Function

Advanced Topics in Optimization. Piecewise Linear Approximation of a Nonlinear Function Advanced Tocs n Otmzaton Pecewse Lnear Aroxmaton of a Nonlnear Functon Otmzaton Methods: M8L Introducton and Objectves Introducton There exsts no general algorthm for nonlnear rogrammng due to ts rregular

More information

PARTIAL QUOTIENTS AND DISTRIBUTION OF SEQUENCES. Department of Mathematics University of California Riverside, CA

PARTIAL QUOTIENTS AND DISTRIBUTION OF SEQUENCES. Department of Mathematics University of California Riverside, CA PARTIAL QUOTIETS AD DISTRIBUTIO OF SEQUECES 1 Me-Chu Chang Deartment of Mathematcs Unversty of Calforna Rversde, CA 92521 mcc@math.ucr.edu Abstract. In ths aer we establsh average bounds on the artal quotents

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

NECESSARY AND SUFFICIENT CONDITIONS FOR ALMOST REGULARITY OF UNIFORM BIRKHOFF INTERPOLATION SCHEMES. by Nicolae Crainic

NECESSARY AND SUFFICIENT CONDITIONS FOR ALMOST REGULARITY OF UNIFORM BIRKHOFF INTERPOLATION SCHEMES. by Nicolae Crainic NECESSARY AND SUFFICIENT CONDITIONS FOR ALMOST REGULARITY OF UNIFORM BIRKHOFF INTERPOLATION SCHEMES by Ncolae Cranc Abstract: In ths artcle usng a combnaton of the necessary and suffcent condtons for the

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

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

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

A NOTE ON THE DISCRETE FOURIER RESTRICTION PROBLEM

A NOTE ON THE DISCRETE FOURIER RESTRICTION PROBLEM A NOTE ON THE DISCRETE FOURIER RESTRICTION PROBLEM XUDONG LAI AND YONG DING arxv:171001481v1 [mathap] 4 Oct 017 Abstract In ths aer we establsh a general dscrete Fourer restrcton theorem As an alcaton

More information

2-Adic Complexity of a Sequence Obtained from a Periodic Binary Sequence by Either Inserting or Deleting k Symbols within One Period

2-Adic Complexity of a Sequence Obtained from a Periodic Binary Sequence by Either Inserting or Deleting k Symbols within One Period -Adc Comlexty of a Seuence Obtaned from a Perodc Bnary Seuence by Ether Insertng or Deletng Symbols wthn One Perod ZHAO Lu, WEN Qao-yan (State Key Laboratory of Networng and Swtchng echnology, Bejng Unversty

More information

Research Article The Point Zoro Symmetric Single-Step Procedure for Simultaneous Estimation of Polynomial Zeros

Research Article The Point Zoro Symmetric Single-Step Procedure for Simultaneous Estimation of Polynomial Zeros Aled Mathematcs Volume 2012, Artcle ID 709832, 11 ages do:10.1155/2012/709832 Research Artcle The Pont Zoro Symmetrc Sngle-Ste Procedure for Smultaneous Estmaton of Polynomal Zeros Mansor Mons, 1 Nasruddn

More information

Fuzzy approach to solve multi-objective capacitated transportation problem

Fuzzy approach to solve multi-objective capacitated transportation problem Internatonal Journal of Bonformatcs Research, ISSN: 0975 087, Volume, Issue, 00, -0-4 Fuzzy aroach to solve mult-objectve caactated transortaton roblem Lohgaonkar M. H. and Bajaj V. H.* * Deartment of

More information

Managing Capacity Through Reward Programs. on-line companion page. Byung-Do Kim Seoul National University College of Business Administration

Managing Capacity Through Reward Programs. on-line companion page. Byung-Do Kim Seoul National University College of Business Administration Managng Caacty Through eward Programs on-lne comanon age Byung-Do Km Seoul Natonal Unversty College of Busness Admnstraton Mengze Sh Unversty of Toronto otman School of Management Toronto ON M5S E6 Canada

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

Algorithms for factoring

Algorithms for factoring CSA E0 235: Crytograhy Arl 9,2015 Instructor: Arta Patra Algorthms for factorng Submtted by: Jay Oza, Nranjan Sngh Introducton Factorsaton of large ntegers has been a wdely studed toc manly because of

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

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

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper Games of Threats Elon Kohlberg Abraham Neyman Workng Paper 18-023 Games of Threats Elon Kohlberg Harvard Busness School Abraham Neyman The Hebrew Unversty of Jerusalem Workng Paper 18-023 Copyrght 2017

More information

SMARANDACHE-GALOIS FIELDS

SMARANDACHE-GALOIS FIELDS SMARANDACHE-GALOIS FIELDS W. B. Vasantha Kandasamy Deartment of Mathematcs Indan Insttute of Technology, Madras Chenna - 600 036, Inda. E-mal: vasantak@md3.vsnl.net.n Abstract: In ths aer we study the

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

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

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

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

An application of generalized Tsalli s-havrda-charvat entropy in coding theory through a generalization of Kraft inequality

An application of generalized Tsalli s-havrda-charvat entropy in coding theory through a generalization of Kraft inequality Internatonal Journal of Statstcs and Aled Mathematcs 206; (4): 0-05 ISS: 2456-452 Maths 206; (4): 0-05 206 Stats & Maths wwwmathsjournalcom Receved: 0-09-206 Acceted: 02-0-206 Maharsh Markendeshwar Unversty,

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Exerments-I MODULE III LECTURE - 2 EXPERIMENTAL DESIGN MODELS Dr. Shalabh Deartment of Mathematcs and Statstcs Indan Insttute of Technology Kanur 2 We consder the models

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

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

ρ some λ THE INVERSE POWER METHOD (or INVERSE ITERATION) , for , or (more usually) to

ρ some λ THE INVERSE POWER METHOD (or INVERSE ITERATION) , for , or (more usually) to THE INVERSE POWER METHOD (or INVERSE ITERATION) -- applcaton of the Power method to A some fxed constant ρ (whch s called a shft), x λ ρ If the egenpars of A are { ( λ, x ) } ( ), or (more usually) to,

More information

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1 Problem Set 4 Suggested Solutons Problem (A) The market demand functon s the soluton to the followng utlty-maxmzaton roblem (UMP): The Lagrangean: ( x, x, x ) = + max U x, x, x x x x st.. x + x + x y x,

More information

On Finite Rank Perturbation of Diagonalizable Operators

On Finite Rank Perturbation of Diagonalizable Operators Functonal Analyss, Approxmaton and Computaton 6 (1) (2014), 49 53 Publshed by Faculty of Scences and Mathematcs, Unversty of Nš, Serba Avalable at: http://wwwpmfnacrs/faac On Fnte Rank Perturbaton of Dagonalzable

More information

Power law and dimension of the maximum value for belief distribution with the max Deng entropy

Power law and dimension of the maximum value for belief distribution with the max Deng entropy Power law and dmenson of the maxmum value for belef dstrbuton wth the max Deng entropy Bngy Kang a, a College of Informaton Engneerng, Northwest A&F Unversty, Yanglng, Shaanx, 712100, Chna. Abstract Deng

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

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

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

The control of nonlinear systems is difficult because

The control of nonlinear systems is difficult because LECURE NOES A Lnear atrx Inequalty Aroach for the Control of Uncertan Fuzzy Systems By H.K. Lam, Frank H.F. Leung, Peter K.S. am he control of nonlnear systems s dffcult because no systematc mathematcal

More information

A new eigenvalue inclusion set for tensors with its applications

A new eigenvalue inclusion set for tensors with its applications COMPUTATIONAL SCIENCE RESEARCH ARTICLE A new egenvalue ncluson set for tensors wth ts applcatons Cal Sang 1 and Janxng Zhao 1 * Receved: 30 Deceber 2016 Accepted: 12 Aprl 2017 Frst Publshed: 20 Aprl 2017

More information

An Improved multiple fractal algorithm

An Improved multiple fractal algorithm Advanced Scence and Technology Letters Vol.31 (MulGraB 213), pp.184-188 http://dx.do.org/1.1427/astl.213.31.41 An Improved multple fractal algorthm Yun Ln, Xaochu Xu, Jnfeng Pang College of Informaton

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

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

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

Some congruences related to harmonic numbers and the terms of the second order sequences

Some congruences related to harmonic numbers and the terms of the second order sequences Mathematca Moravca Vol. 0: 06, 3 37 Some congruences related to harmonc numbers the terms of the second order sequences Neşe Ömür Sbel Koaral Abstract. In ths aer, wth hels of some combnatoral denttes,

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

TOWARDS THE GLOBAL SOLUTION OF THE MAXIMAL CORRELATION PROBLEM

TOWARDS THE GLOBAL SOLUTION OF THE MAXIMAL CORRELATION PROBLEM TOWARDS THE GLOBAL SOLUTION OF THE MAXIMAL CORRELATION PROBLEM LEI-HONG ZHANG, LI-ZHI LIAO, AND LI-MING SUN Abstract. The maxmal correlaton problem (MCP) amng at optmzng correlaton between sets of varables

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

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

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

On the Global Linear Convergence of the ADMM with Multi-Block Variables

On the Global Linear Convergence of the ADMM with Multi-Block Variables On the Global Lnear Convergence of the ADMM wth Mult-Block Varables Tany Ln Shqan Ma Shuzhong Zhang May 31, 01 Abstract The alternatng drecton method of multplers ADMM has been wdely used for solvng structured

More information

Report on Image warping

Report on Image warping Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.

More information

MOLECULAR TOPOLOGICAL INDEX OF C 4 C 8 (R) AND C 4 C 8 (S) NANOTORUS

MOLECULAR TOPOLOGICAL INDEX OF C 4 C 8 (R) AND C 4 C 8 (S) NANOTORUS Dgest Journal of Nanomaterals an Bostructures ol No December 9 69-698 MOLECULAR TOPOLOICAL INDEX OF C C 8 R AND C C 8 S NANOTORUS ABBAS HEYDARI Deartment of Mathematcs Islamc Aza Unverst of Ara Ara 85-567

More information

Norm Bounds for a Transformed Activity Level. Vector in Sraffian Systems: A Dual Exercise

Norm Bounds for a Transformed Activity Level. Vector in Sraffian Systems: A Dual Exercise ppled Mathematcal Scences, Vol. 4, 200, no. 60, 2955-296 Norm Bounds for a ransformed ctvty Level Vector n Sraffan Systems: Dual Exercse Nkolaos Rodousaks Department of Publc dmnstraton, Panteon Unversty

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

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

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

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

6. Hamilton s Equations

6. Hamilton s Equations 6. Hamlton s Equatons Mchael Fowler A Dynamcal System s Path n Confguraton Sace and n State Sace The story so far: For a mechancal system wth n degrees of freedom, the satal confguraton at some nstant

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Malaya J. Mat. 2(1)(2014) 49 60

Malaya J. Mat. 2(1)(2014) 49 60 Malaya J. Mat. 2(1)(2014) 49 60 Functonal equaton orgnatng from sum of hgher owers of arthmetc rogresson usng dfference oerator s stable n Banach sace: drect and fxed ont methods M. Arunumar a, and G.

More information

Received on October 24, 2011 / Revised on November 17, 2011

Received on October 24, 2011 / Revised on November 17, 2011 Journal of ath-for-industry Vol 4 202A-2 5 5 On some roertes of a dscrete hungry Lota-Volterra system of multlcatve tye Yosue Hama Ao Fuuda Yusau Yamamoto asash Iwasa Emo Ishwata and Yoshmasa Naamura Receved

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

CHAPTER 4 MAX-MIN AVERAGE COMPOSITION METHOD FOR DECISION MAKING USING INTUITIONISTIC FUZZY SETS

CHAPTER 4 MAX-MIN AVERAGE COMPOSITION METHOD FOR DECISION MAKING USING INTUITIONISTIC FUZZY SETS 56 CHAPER 4 MAX-MIN AVERAGE COMPOSIION MEHOD FOR DECISION MAKING USING INUIIONISIC FUZZY SES 4.1 INRODUCION Intutonstc fuzz max-mn average composton method s proposed to construct the decson makng for

More information

On Graphs with Same Distance Distribution

On Graphs with Same Distance Distribution Appled Mathematcs, 07, 8, 799-807 http://wwwscrporg/journal/am ISSN Onlne: 5-7393 ISSN Prnt: 5-7385 On Graphs wth Same Dstance Dstrbuton Xulang Qu, Xaofeng Guo,3 Chengy Unversty College, Jme Unversty,

More information

The Matrix Analogs of Firey s Extension of Minkowski Inequality and of Firey s Extension of Brunn-Minkowski Inequality

The Matrix Analogs of Firey s Extension of Minkowski Inequality and of Firey s Extension of Brunn-Minkowski Inequality Proceedngs of the 9th WSEAS Internatonal Conference on Aled Mathematcs, Istanbul, Turkey, May 27-29, 2006 (489-494) The Matrx Analogs of Frey s Extenson of Mnkowsk Inequalty and of Frey s Extenson of Brunn-Mnkowsk

More information

Complete weight enumerators of two classes of linear codes

Complete weight enumerators of two classes of linear codes Comlete weght enumerators of two classes of lnear codes Quyan Wang, Fe L, Kelan Dng and Dongda Ln 1 Abstract arxv:1512.7341v1 [cs.it] 23 Dec 215 Recently, lnear codes wth few weghts have been constructed

More information

The internal structure of natural numbers and one method for the definition of large prime numbers

The internal structure of natural numbers and one method for the definition of large prime numbers The nternal structure of natural numbers and one method for the defnton of large prme numbers Emmanul Manousos APM Insttute for the Advancement of Physcs and Mathematcs 3 Poulou str. 53 Athens Greece Abstract

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

Fuzzy Set Approach to Solve Multi-objective Linear plus Fractional Programming Problem

Fuzzy Set Approach to Solve Multi-objective Linear plus Fractional Programming Problem Internatonal Journal of Oeratons Research Vol.8, o. 3, 5-3 () Internatonal Journal of Oeratons Research Fuzzy Set Aroach to Solve Mult-objectve Lnear lus Fractonal Programmng Problem Sanjay Jan Kalash

More information

Lecture 9: Converse of Shannon s Capacity Theorem

Lecture 9: Converse of Shannon s Capacity Theorem Error Correctng Codes: Combnatorcs, Algorthms and Alcatons (Fall 2007) Lecture 9: Converse of Shannon s Caacty Theorem Setember 17, 2007 Lecturer: Atr Rudra Scrbe: Thanh-Nhan Nguyen & Atr Rudra In the

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

New Kronecker product decompositions and its applications

New Kronecker product decompositions and its applications RESEARCH INVENTY: Internatonal Journal of Engneerng and Scence ISBN: 319-6483, ISSN: 78-471, Vol. 1, Issue 11 (December 01), PP 5-30 www.researchnventy.com New Kronecer roduct decomostons and ts alcatons

More information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

SL n (F ) Equals its Own Derived Group

SL n (F ) Equals its Own Derived Group Internatonal Journal of Algebra, Vol. 2, 2008, no. 12, 585-594 SL n (F ) Equals ts Own Derved Group Jorge Macel BMCC-The Cty Unversty of New York, CUNY 199 Chambers street, New York, NY 10007, USA macel@cms.nyu.edu

More information

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several

More information

Lecture 10: May 6, 2013

Lecture 10: May 6, 2013 TTIC/CMSC 31150 Mathematcal Toolkt Sprng 013 Madhur Tulsan Lecture 10: May 6, 013 Scrbe: Wenje Luo In today s lecture, we manly talked about random walk on graphs and ntroduce the concept of graph expander,

More information

The Degree Distribution of Random Birth-and-Death Network with Network Size Decline

The Degree Distribution of Random Birth-and-Death Network with Network Size Decline The Degree Dstrbuton of Random Brth-and-Death etwork wth etwork Sze Declne Xaojun Zhang *, Hulan Yang School of Mathematcal Scences, Unversty of Electronc Scence and Technology of Chna, Chengdu 673, P.R.

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

Non-Ideality Through Fugacity and Activity

Non-Ideality Through Fugacity and Activity Non-Idealty Through Fugacty and Actvty S. Patel Deartment of Chemstry and Bochemstry, Unversty of Delaware, Newark, Delaware 19716, USA Corresondng author. E-mal: saatel@udel.edu 1 I. FUGACITY In ths dscusson,

More information

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function A Local Varatonal Problem of Second Order for a Class of Optmal Control Problems wth Nonsmooth Objectve Functon Alexander P. Afanasev Insttute for Informaton Transmsson Problems, Russan Academy of Scences,

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

A Note on \Modules, Comodules, and Cotensor Products over Frobenius Algebras"

A Note on \Modules, Comodules, and Cotensor Products over Frobenius Algebras Chn. Ann. Math. 27B(4), 2006, 419{424 DOI: 10.1007/s11401-005-0025-z Chnese Annals of Mathematcs, Seres B c The Edtoral Oce of CAM and Sprnger-Verlag Berln Hedelberg 2006 A Note on \Modules, Comodules,

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

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

The Second Anti-Mathima on Game Theory

The Second Anti-Mathima on Game Theory The Second Ant-Mathma on Game Theory Ath. Kehagas December 1 2006 1 Introducton In ths note we wll examne the noton of game equlbrum for three types of games 1. 2-player 2-acton zero-sum games 2. 2-player

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results. Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson

More information

A p-adic PERRON-FROBENIUS THEOREM

A p-adic PERRON-FROBENIUS THEOREM A p-adic PERRON-FROBENIUS THEOREM ROBERT COSTA AND PATRICK DYNES Advsor: Clayton Petsche Oregon State Unversty Abstract We prove a result for square matrces over the p-adc numbers akn to the Perron-Frobenus

More information

MAT 578 Functional Analysis

MAT 578 Functional Analysis MAT 578 Functonal Analyss John Qugg Fall 2008 Locally convex spaces revsed September 6, 2008 Ths secton establshes the fundamental propertes of locally convex spaces. Acknowledgment: although I wrote these

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

Short running title: A generating function approach A GENERATING FUNCTION APPROACH TO COUNTING THEOREMS FOR SQUARE-FREE POLYNOMIALS AND MAXIMAL TORI

Short running title: A generating function approach A GENERATING FUNCTION APPROACH TO COUNTING THEOREMS FOR SQUARE-FREE POLYNOMIALS AND MAXIMAL TORI Short runnng ttle: A generatng functon approach A GENERATING FUNCTION APPROACH TO COUNTING THEOREMS FOR SQUARE-FREE POLYNOMIALS AND MAXIMAL TORI JASON FULMAN Abstract. A recent paper of Church, Ellenberg,

More information

Convergence of random processes

Convergence of random processes DS-GA 12 Lecture notes 6 Fall 216 Convergence of random processes 1 Introducton In these notes we study convergence of dscrete random processes. Ths allows to characterze phenomena such as the law of large

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

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information