A FORMULA FOR COMPUTING INTEGER POWERS FOR ONE TYPE OF TRIDIAGONAL MATRIX
|
|
- Marsha Fox
- 5 years ago
- Views:
Transcription
1 Hacettepe Journal of Mathematcs and Statstcs Volume 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 :0 :0 bstract In ths paper we derve the general expresson of the r th power r N for one type of trdagonal matrx Keywords: Trdagonal matrx Matrx power Egenvalues Chebyshev polynomal 000 MS Classfcaton: C C E Introducton In the present paper we derve a general expresson for the rth power for one type of trdagonal matrx r N and N denotes the set of natural numbers In [4] Rmas derved a general expresson for the lth power for one type of symmetrc trdagonal matrces of even order General expressons for the rth power of a matrx are obtaned by usng the equalty r = PJ r P [] J s the Jordan form of and P the transformng matrx We need the egenvalues and egenvectors of the matrx to compute the matrces J and P The egenvalues of are the roots x nk = cos kπ k = n n + of the nth degree Chebyshev polynomal U nx = of the second knd snn + arccosx x snarccosx Selcuk Unversty Scence Faculty Department of Mathematcs 450 Campus Konya Turkey E-mal: H Kıyak kyak8@hotmalcom I G Gürses rematca@hotmalcom F Yılmaz fylmaz@selcukedutr D Bozkurt dbozkurt@selcukedutr Correspondng uthor
2 35 K Kıyak I G Gürses F Yılmaz D Bozkurt Man results Let be the followng n-square trdagonal matrx = 0 0 n n 0 The egenvalues of are the roots of the characterstc equaton λe = 0 and also the roots of U nx the nth degree Chebyshev polynomal of the second knd [3] In the present paper we nvestgate nteger powers of ths matrx Intally we wll scrutnze the nteger powers for even orders Let us wrte α D nα = λe = α α and α 3 nα = α α α = λ R By and 3 we obtan the followng two results 4 D nα = nα 5 nα = α n α n α α = α α = α and 0α = We obtan 4 by the nducton prncple For n = D α = α = α
3 Computng Integer Powers 353 Let D k α = k α for n = k We have to show that D k+ α = k+ α for n = k + Let α D k+ α = α k+ k+ If we expand the determnant accordng to frst row then we have α 0 α 7 D k+ α = α + α α and f we expand the frst determnant n 7 accordng to last row and the second determnant accordng to frst column then we obtan α D k+ α = α α α + α 0 α α α }} D k α= k α } } k α
4 354 K Kıyak I G Gürses F Yılmaz D Bozkurt Then α D k+ α = α α k α k α α } } k α = α[α k α k α] k α }} k+ α = α k+ α k α = k+ α By solvng the dfference equaton 5 we obtan [] α 8 nα = U n It s known that the roots of the polynomal U nx are [3] 9 x nk = cos kπ k = n n + whch are ncluded n the nterval [ ] Takng and 7 nto account we fnd the egenvalues of the matrx to be λ k = cos kπ k = n n + Snce the multplcty of each of the egenvalues λ k s and applyng the relaton λ k = λ n k+ k = n we wrte down the Jordan form of the matrx as J = dag λ n λ n λ n λ n + λ n + λ n λ n λ n Let us solve the homogeneous lnear equatons λ E x = 0 λ s the th egenvalue of the matrx n By elementary row operatons the coeffcent matrx of the system s λ λ λ λ λ D n λ
5 Computng Integer Powers 355 Snce D n λ = 0 rankλ E = n Then x + λ x x 3 = 0 x λ x 3 x 4 = 0 x n λ x n x n = 0 x n + λ x n = 0 Now we wll scrutnze the solutons of the above system of lnear equatons accordng to n and k λ Let x n = = U 0 and for any n and k defne n k or mod 4 a = n k 0 or 3 mod 4 n or mod 4 e = n 0 or 3 mod 4 We consder the followng two cases: Case Let n be an even number If k s an odd number then λ x n = λ = U x n = λ λ = U x k+ = aλ n k + cλ n k n k x k+ = aλ n k + bλ n k 3 + dλ n k 5 x k = aλ n k + a bλ n k λ λ = au n k λ + = au n k + c dλ n k n k + λ λ = au n k x = eλ n en λ n n λ = eun λ b = an k c = an k 3 and d s the sum of the coeffcents of the terms λ n k 5 n x k 3 and λ n k n x k
6 35 K Kıyak I G Gürses F Yılmaz D Bozkurt If k s an even number then λ x n = λ = U x n = λ λ = U x k+3 = aλ n k n k x k+ = aλ n k + cλ n k 4 x k+ = aλ n k + bλ n k 3 λ = au n k 3 λ + = au n k λ + dλ n k n k + + n λ = eun λ λ λ = au n k x k = aλ n k + b aλ n k + c + dλ n k 4 λ + = au n k x = eλ n en λ n 3 b = an k c = an k 3 and d s the dfference of the coeffcents of the terms λ n k 5 x k+ Case Let n be an odd number If k s an odd number x n = λ = U λ x n = λ λ = U x k+3 = aλ n k n k x k+ = aλ n k + cλ n k 4 x k+ = aλ n k + bλ n k 3 λ = au n k 3 λ + = au n k λ + dλ n k n k + = eu n λ n x k+3 and λ n k λ λ = au n k x k = aλ n k + b aλ n k + c + dλ n k 4 λ + = au n k x = eλ n en λ n 3 b = an k c = an k 3 n and d s the dfference of the coeffcents of the terms λ n k 5 x k+ n x k+3 and λ n k n
7 Computng Integer Powers 357 If k s an even number then λ x n = λ = U x n = λ λ = U x k+ = aλ n k + cλ n k n k x k+ = aλ n k + bλ n k 3 + dλ n k 5 x k = aλ n k + a bλ n k λ λ = au n k λ + = au n k + c dλ n k n k + λ λ = au n k x = eλ n en λ n 3 b = an k c = an k 3 + = eu n λ and d s the sum of the coeffcents of the terms λ n k 5 Takng P = x λ nto account P = [ P P P n ] n x k 3 and λ n k n x k Usng ths expresson we can wrte down the transformng matrx P as eu n λ eu n λ eu n λn au n k λ au n k λ au n k λn P = U 3 λ U 3 λ U 3 λn U λ U λ U λn U λ U λ U λn U 0 λ U 0 λ U 0 λn Denotng the th row of P by R = n from we have 3 R = [ t U n λ t U n λ t U n λn ] t = n or mod 4 n 0 or 3 mod 4
8 358 K Kıyak I G Gürses F Yılmaz D Bozkurt Denotng jth column of the nverse matrx P by ρ j we obtan m j n+ n j m j U λ n+ n j 4 ρ j = m 3 j U λ3 n+ n j m n j U λn n+ n j m j = n j or mod 4 n j 0 or 3 mod 4 Takng expresson 4 nto account we obtan the matrx P as follows: e n+ n n+ 3 n+ P e U λ n+ n U λ n+ 3 U λ n+ = e n U λn n+ n n U λn n n+ 3 U λn n+ n+ n+ n n+ n+ 0 U λ U λ n+ 0 U λn n n+ U λn 0 Usng the equalty r = PJ r P we derve r } j = PJ r P } j = RJr ρ j = tmj n + n k= λ r k4 λ λk ku n λ r m j n+ n j λ r m j U λ n+ n j = R λ r 3m 3 j U λ3 n+ n j λ r nm n j U λn n+ n j U n j λk λ are the egenvalues of the matrx U k x s the kth degree Chebyshev polynomal of the second knd j = n and we equate ths matrx wth ts unque element as usual
9 Computng Integer Powers 359 Now we consder the case when the order n s odd Workng as before we obtan α 5 D nα = λe = α α and α nα = α α for odd orders α = λ R lso the equaltes and 3 are vald for odd orders Snce the multplcty of all the egenvalues λ k s and applyng the relaton λ k = λ n k+ k = n and λ = 0 we can wrte down the Jordan form of the matrx as: 7 J = dag λ n λ n λ n λ n+3 0 λ n+3 λ n λ n λ n Denotng jth column of the nverse matrx P by ρ j we obtan 8 ρ j = λ m + j n+ λ m j λ m + j n+ λ m +3 j n+ +n n+ n j U λ n j U λ3 n j U λn n j m j = n j or mod 4 n j 0 or 3 mod 4
10 30 K Kıyak I G Gürses F Yılmaz D Bozkurt Takng the expresson 8 nto account we obtan the matrx P as follows: P e = λ e + n+ λ + n+ λ e +n n+ n U λ n λ + n+ λ + n+ U λn λ n +n n+ 3 U λ 3 λ + n+ λ + n+ U λ U λn λ +n 3 n+ U λn Usng the equalty r = PJ r P we derve r } j = PJ r P } j = RJr ρ j = r t m j n + n k= λ + n+ λ + n+ U λ λ +n n+ = R λ λ r nm j λ r kλ λk +kun λ + n+ λ + n+ U λn λ +n n+ λ λ r m + j n+ λ λ r m + j n+ λ λ r 3m +3 j n+ +n n+ U n j λk n j U λ n j U λ3 n j U λn n j 0 U λ 0 U λn 0 λ are the egenvalues of the matrx U k x s the kth degree Chebyshev polynomal of the second knd j = n and we equate the matrx wth ts unque element 3 Numercal examples Takng nto account the derved expressons for n even one can see the effectveness of the formula For n = 4 we obtan; J = dag λ 4 λ 3 λ 3 λ 4 = dag a b b a a = cos π 5 b = cos π 5
11 Computng Integer Powers 3 If r s even r} = a r 4 a a 3 a + b r 4 b b 3 b r} = r} = 0 r} = r } 3 3 = a r+ 4 a a 3 a + b r+ 4 b b 3 b r} = r} = r} = a r 4 a a + b r 4 b b r} = r} = r} = r} = 4 4 a r 4 a a + b r 4 b b r} = 33 a r+ 4 a + b r+ 4 b r} = r} = r} = 44 a r 4 a + b r 4 b For example for r = 4; = If r s odd r} = 0 r} = r} = [ a r 4 a a 3 aa + b r 4 b b 3 bb ] r} = r} = r} = r} = 4 4 a r 4 a a 3 a + b r 4 b b 3 b r} = 0 r} = r} = 3 3 a r+ 4 a a + b r+ 4 b b r} = r} = r} = 0 33 r} = r} = a r+ 4 a + b r+ 4 b r} = 0 44 For example for r = =
12 3 K Kıyak I G Gürses F Yılmaz D Bozkurt We wll now see the effectveness of the formula for odd orders For example f n = 5 we obtan If r s even J = dag λ 5 λ 4 0 λ 4 λ 5 = dag c d 0 d c c = 3 d = r } = c r d + c d r = 3 r + 3 r } = r} = 0 r } 3 = r} 3 = 3 cr d = 3 r r } 4 = r} 4 = 0 r } = r} = 5 5 c r d c d r = 3 r 3 r } = cr+ d + c d r = 3 r+ + 3 r } 3 = r} 3 = 0 r } = r} = 4 4 c d r c r+ d = 3 3 r+ r } = r} = r } 33 = 3 ar b = 3 cr = c r = 3 r r } 34 = r} 43 = 0 r } 35 = r} 53 = 3 cr d = c r = 3 r r } = 44 c r+ d + c d r = 3 r+ + 3 } r = r} = } r = 55 c r d + c d r = 3 r + 3 For example for r = 4 we have = If r s odd r } = 0 r } = r} = c r+ d + c d r = 3 r+ + 3 r } = r} = r } = r} = 4 4 c r+ d c d r = 3 r+ 3 } r = r} = r } = 0
13 Computng Integer Powers 33 r } 3 = r} 3 = 3 cr+ d = 3 r r } 4 = r} 4 = 0 r } = r} = 5 5 c d r c r+ d = 3 3 r+ } r = 0 33 r } 34 = r} 43 = 3 cr+ d = c r = 3 r r } 35 = r} 53 = 0 r } 44 = 0 } r = r} = c r+ d + c d r = 3 r+ + 3 r } 55 = 0 For example for r = 3 we have: = cknowledgement We are grateful to Jonas Rmas for comments and contrbutons References [] garwal R P Dfference Equatons and Inequltes Marcel Dekker New York 99 [] Horn R and Johnson Ch Matrx nalyss Cambrdge Unversty Press Cambrdge 985 [3] Mason J C and Handscomb D C Chebyshev Polynomals CRC Press Washngton 003 [4] Rmas J On computng of arbtrary postve nteger powers for one type of symmetrc trdagonal matrces of even order-i ppled Mathematcs and Computaton [5] Rmas J On computng of arbtrary postve nteger powers for one type of symmetrc trdagonal matrces of odd order-i ppled Mathematcs and Computaton
arxiv: v1 [math.co] 12 Sep 2014
arxv:1409.3707v1 [math.co] 12 Sep 2014 On the bnomal sums of Horadam sequence Nazmye Ylmaz and Necat Taskara Department of Mathematcs, Scence Faculty, Selcuk Unversty, 42075, Campus, Konya, Turkey March
More informationA property of the elementary symmetric functions
Calcolo manuscrpt No. (wll be nserted by the edtor) A property of the elementary symmetrc functons A. Esnberg, G. Fedele Dp. Elettronca Informatca e Sstemstca, Unverstà degl Stud della Calabra, 87036,
More informationRandić 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 informationThe Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices
Internatonal Mathematcal Forum, Vol 11, 2016, no 11, 513-520 HIKARI Ltd, wwwm-hkarcom http://dxdoorg/1012988/mf20166442 The Jacobsthal and Jacobsthal-Lucas Numbers va Square Roots of Matrces Saadet Arslan
More informationConvexity 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 informationThe 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 informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationFormulas 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 informationPerron 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 informationOn 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 informatione - 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 informationOn 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= = = (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 informationDifference 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 informationMath 217 Fall 2013 Homework 2 Solutions
Math 17 Fall 013 Homework Solutons Due Thursday Sept. 6, 013 5pm Ths homework conssts of 6 problems of 5 ponts each. The total s 30. You need to fully justfy your answer prove that your functon ndeed has
More information1 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 information2.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 informationSL 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 informationCOMPARISON 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 informationLectures - 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 informationHyper-Sums of Powers of Integers and the Akiyama-Tanigawa Matrix
6 Journal of Integer Sequences, Vol 8 (00), Artcle 0 Hyper-Sums of Powers of Integers and the Ayama-Tangawa Matrx Yoshnar Inaba Toba Senor Hgh School Nshujo, Mnam-u Kyoto 60-89 Japan nava@yoto-benejp Abstract
More informationA 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 informationSalmon: 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 informationYong 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 informationThe Degrees of Nilpotency of Nilpotent Derivations on the Ring of Matrices
Internatonal Mathematcal Forum, Vol. 6, 2011, no. 15, 713-721 The Degrees of Nlpotency of Nlpotent Dervatons on the Rng of Matrces Homera Pajoohesh Department of of Mathematcs Medgar Evers College of CUNY
More informationThe 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 informationA 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 informationBy Samuel Schechter. (1) H = i-î -1 with 1 ^ i,j g n. [a, b,j
On tne Inverson of Certan Matrces By Samuel Schechter 1. Introducton. Let a, a2,, an ; 61, 62,, 6 be 2n dstnct, but otherwse arbtrary, complex numbers. For the matrx H, of order n, (1) H = -î -1 wth 1
More informationThe 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 informationChapter 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 informationSome 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 informationA 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 informationHowever, since P is a symmetric idempotent matrix, of P are either 0 or 1 [Eigen-values
Fall 007 Soluton to Mdterm Examnaton STAT 7 Dr. Goel. [0 ponts] For the general lnear model = X + ε, wth uncorrelated errors havng mean zero and varance σ, suppose that the desgn matrx X s not necessarly
More informationNOTES 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 informationRemarks 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 informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased
More informationDeterminants 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 informationLECTURE 9 CANONICAL CORRELATION ANALYSIS
LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of
More informationNorms, 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 informationTHE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens
THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of
More informationLinear 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 informationApplied 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 informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationBinomial transforms of the modified k-fibonacci-like sequence
Internatonal Journal of Mathematcs and Computer Scence, 14(2019, no. 1, 47 59 M CS Bnomal transforms of the modfed k-fbonacc-lke sequence Youngwoo Kwon Department of mathematcs Korea Unversty Seoul, Republc
More informationTHE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD
Journal of Appled Mathematcs and Computatonal Mechancs 06, 5(), 7-36 www.amcm.pcz.pl p-iss 99-9965 DOI: 0.75/jamcm.06..4 e-iss 353-0588 THE STURM-LIOUVILLE EIGEVALUE PROBLEM - A UMERICAL SOLUTIO USIG THE
More informationCME 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 informationInexact 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 informationNorm 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 informationAPPENDIX 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 informationBallot Paths Avoiding Depth Zero Patterns
Ballot Paths Avodng Depth Zero Patterns Henrch Nederhausen and Shaun Sullvan Florda Atlantc Unversty, Boca Raton, Florda nederha@fauedu, ssull21@fauedu 1 Introducton In a paper by Sapounaks, Tasoulas,
More informationThe binomial transforms of the generalized (s, t )-Jacobsthal matrix sequence
Int. J. Adv. Appl. Math. and Mech. 6(3 (2019 14 20 (ISSN: 2347-2529 Journal homepage: www.jaamm.com IJAAMM Internatonal Journal of Advances n Appled Mathematcs and Mechancs The bnomal transforms of the
More informationAGC Introduction
. Introducton AGC 3 The prmary controller response to a load/generaton mbalance results n generaton adjustment so as to mantan load/generaton balance. However, due to droop, t also results n a non-zero
More informationEvaluation of a family of binomial determinants
Electronc Journal of Lnear Algebra Volume 30 Volume 30 2015 Artcle 22 2015 Evaluaton of a famly of bnomal determnants Charles Helou Pennsylvana State Unversty, cxh22@psuedu James A Sellers Pennsylvana
More informationSolving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions
ISSN 746-7659 England UK Journal of Informaton and Computng Scence Vol. 9 No. 3 4 pp. 69-8 Solvng Fractonal Nonlnear Fredholm Integro-dfferental Equatons va Hybrd of Ratonalzed Haar Functons Yadollah Ordokhan
More informationTHE RING AND ALGEBRA OF INTUITIONISTIC SETS
Hacettepe Journal of Mathematcs and Statstcs Volume 401 2011, 21 26 THE RING AND ALGEBRA OF INTUITIONISTIC SETS Alattn Ural Receved 01:08 :2009 : Accepted 19 :03 :2010 Abstract The am of ths study s to
More informationCOMPLEX NUMBERS AND QUADRATIC EQUATIONS
COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not
More informationThe 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 information763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.
7636S ADVANCED QUANTUM MECHANICS Soluton Set 1 Sprng 013 1 Warm-up Show that the egenvalues of a Hermtan operator  are real and that the egenkets correspondng to dfferent egenvalues are orthogonal (b)
More informationGeneral 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 informationMATH Sensitivity of Eigenvalue Problems
MATH 537- Senstvty of Egenvalue Problems Prelmnares Let A be an n n matrx, and let λ be an egenvalue of A, correspondngly there are vectors x and y such that Ax = λx and y H A = λy H Then x s called A
More informationCausal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim
Causal Damonds M. Aghl, L. Bombell, B. Plgrm Introducton The correcton to volume of a causal nterval due to curvature of spacetme has been done by Myrhem [] and recently by Gbbons & Solodukhn [] and later
More informationHongyi 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 informationQuantum Mechanics I - Session 4
Quantum Mechancs I - Sesson 4 Aprl 3, 05 Contents Operators Change of Bass 4 3 Egenvectors and Egenvalues 5 3. Denton....................................... 5 3. Rotaton n D....................................
More informationMaejo International Journal of Science and Technology
Maejo Int. J. Sc. Technol. () - Full Paper Maejo Internatonal Journal of Scence and Technology ISSN - Avalable onlne at www.mjst.mju.ac.th Fourth-order method for sngularly perturbed sngular boundary value
More informationSome 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 informationSingular 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 informationComputing π with Bouncing Balls
by Let the mass of two balls be M and m, where M = (6 n )m for n N. The larger ball s rolled towards the lghter ball, whch ear a wall and ntally at rest. Fnd the number of collsons between the two balls
More informationU.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 information5 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 informationTHE GENERALIZED (s, t)-fibonacci AND FIBONACCI MATRIX SEQUENCES
TJMM 7 205), No 2, 37-48 THE GENERALIZED s, t)-fibonacci AND FIBONACCI MATRIX SEQUENCES AHMET İPEK, KAMIL ARI, AND RAMAZAN TÜRKMEN Abstract In ths paper, we study the generalzatons of the s, t)-fbonacc
More informationMATH 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 informationAn efficient algorithm for multivariate Maclaurin Newton transformation
Annales UMCS Informatca AI VIII, 2 2008) 5 14 DOI: 10.2478/v10065-008-0020-6 An effcent algorthm for multvarate Maclaurn Newton transformaton Joanna Kapusta Insttute of Mathematcs and Computer Scence,
More informationarxiv: v6 [math.nt] 23 Aug 2016
A NOTE ON ODD PERFECT NUMBERS JOSE ARNALDO B. DRIS AND FLORIAN LUCA arxv:03.437v6 [math.nt] 23 Aug 206 Abstract. In ths note, we show that f N s an odd perfect number and q α s some prme power exactly
More informationAppendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis
A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems
More informationSemi-Analytical Solution of the 1D Helmholtz Equation, Obtained from Inversion of Symmetric Tridiagonal Matrix
Journal of Electromagnetc nalyss and pplcatons, 04, 6, 45-438 Publshed Onlne December 04 n ScRes. http://www.scrp.org/ournal/emaa http://d.do.org/0.436/emaa.04.64044 Sem-nalytcal Soluton of the D Helmholtz
More informationarxiv: v2 [quant-ph] 29 Jun 2018
Herarchy of Spn Operators, Quantum Gates, Entanglement, Tensor Product and Egenvalues Wll-Hans Steeb and Yorck Hardy arxv:59.7955v [quant-ph] 9 Jun 8 Internatonal School for Scentfc Computng, Unversty
More informationρ 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 informationFirst 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 informationFoundations 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 informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton
More informationAsymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation
Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton
More informationSolution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method
Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222
More informationA summation on Bernoulli numbers
Journal of Number Theory 111 (005 37 391 www.elsever.com/locate/jnt A summaton on Bernoull numbers Kwang-Wu Chen Department of Mathematcs and Computer Scence Educaton, Tape Muncpal Teachers College, No.
More informationA generalization of a trace inequality for positive definite matrices
A generalzaton of a trace nequalty for postve defnte matrces Elena Veronca Belmega, Marc Jungers, Samson Lasaulce To cte ths verson: Elena Veronca Belmega, Marc Jungers, Samson Lasaulce. A generalzaton
More informationChat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980
MT07: Multvarate Statstcal Methods Mke Tso: emal mke.tso@manchester.ac.uk Webpage for notes: http://www.maths.manchester.ac.uk/~mkt/new_teachng.htm. Introducton to multvarate data. Books Chat eld, C. and
More informationLecture 5 Decoding Binary BCH Codes
Lecture 5 Decodng Bnary BCH Codes In ths class, we wll ntroduce dfferent methods for decodng BCH codes 51 Decodng the [15, 7, 5] 2 -BCH Code Consder the [15, 7, 5] 2 -code C we ntroduced n the last lecture
More informationExample: (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 informationNumber of cases Number of factors Number of covariates Number of levels of factor i. Value of the dependent variable for case k
ANOVA Model and Matrx Computatons Notaton The followng notaton s used throughout ths chapter unless otherwse stated: N F CN Y Z j w W Number of cases Number of factors Number of covarates Number of levels
More informationThe Synchronous 8th-Order Differential Attack on 12 Rounds of the Block Cipher HyRAL
The Synchronous 8th-Order Dfferental Attack on 12 Rounds of the Block Cpher HyRAL Yasutaka Igarash, Sej Fukushma, and Tomohro Hachno Kagoshma Unversty, Kagoshma, Japan Emal: {garash, fukushma, hachno}@eee.kagoshma-u.ac.jp
More informationLecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2
P470 Lecture 6/7 (February 10/1, 014) DIRAC EQUATION The non-relatvstc Schrödnger equaton was obtaned by notng that the Hamltonan H = P (1) m can be transformed nto an operator form wth the substtutons
More informationVector 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 informationHomework 1 Lie Algebras
Homework 1 Le Algebras Joshua Ruter February 9, 018 Proposton 0.1 Problem 1.7a). Let A be a K-algebra, wth char K. Then A s alternatve f and only f the folowng two lnear) denttes hold for all a, b, y A.
More informationExercises. 18 Algorithms
18 Algorthms Exercses 0.1. In each of the followng stuatons, ndcate whether f = O(g), or f = Ω(g), or both (n whch case f = Θ(g)). f(n) g(n) (a) n 100 n 200 (b) n 1/2 n 2/3 (c) 100n + log n n + (log n)
More informationProblem Solving in Math (Math 43900) Fall 2013
Problem Solvng n Math (Math 43900) Fall 2013 Week four (September 17) solutons Instructor: Davd Galvn 1. Let a and b be two nteger for whch a b s dvsble by 3. Prove that a 3 b 3 s dvsble by 9. Soluton:
More information18.781: Solution to Practice Questions for Final Exam
18.781: Soluton to Practce Questons for Fnal Exam 1. Fnd three solutons n postve ntegers of x 6y = 1 by frst calculatng the contnued fracton expanson of 6. Soluton: We have 1 6=[, ] 6 6+ =[, ] 1 =[,, ]=[,,
More informationComplex Numbers. x = B B 2 4AC 2A. or x = x = 2 ± 4 4 (1) (5) 2 (1)
Complex Numbers If you have not yet encountered complex numbers, you wll soon do so n the process of solvng quadratc equatons. The general quadratc equaton Ax + Bx + C 0 has solutons x B + B 4AC A For
More informationA 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 informationSTAT 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 informationMEM 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