arxiv: v1 [math.na] 10 Dec 2017

Transcription

1 On a perturbation theory of Hamitonian systems with periodic coefficients Traoré G. Y. Arouna, Mouhamadou Dosso, Jean-Caude Koua Brou arxiv:7.v [math.na] Dec 7 Abstract A theory of rank k perturbation of sympectic matrices and Hamitonian systems with periodic coefficients ( ) using a base of isotropic subspaces, is presented. After showing that the fundamenta matrix X(t) of the rank k perturbation of Hamitonian system with periodic coefficients and t the rank k perturbation of the fundamenta matrix (X(t)) t of the unperturbed system are the ( ) same, the Jordan canonica form of X(t) is given. Two numerica exampes iustrating this t theory and the consequences of rank k perturbations on the strong stabiity of Hamitonian systems were aso given. Mathematics Subject Cassification : 5E5-5F - 5F5-7H5-9D. Key words : Hamitonian system - sympectic matrix - isotropic subspace - perturbations - strong stabiity Introduction The Hamitonian systems with periodic coefficients are generay derived from physica probems and engineering []. These systems are differentia equations with periodic coefficients that originate from the theory of optima contro [, ] and parametric resonance []. They can be put in the form J dx(t) = H(t)x(t), t R (.) where H(t) R N N is symmetric and P-periodic i.e. H(t + P) = H(t) = (H(t)) T and J is skewsymmetric matrix of R N N. The square matrix X(t) with coumns x (t),x (t),...,x N (t) beonging to fundamenta set of soutions of equation (.), is caed a fundamenta matrix. Considering the foowing matrix system [, Vo., chap. ] J dx(t) X() = I = H(t)X(t), t R + whose matrix soution (X(t)) t R+ satisfies the reationship X(t + np) = X(t)X n (P), t R + and n N. We have the foowing definition (.) Definition. The matrix X(t) satisfying equation (.) is caed the matrizant of equation (.). The vaue at the period P of the matrizant X(t) defined by the initia condition X() = I N, is caed the monodromy matrix and its eigenvaues are the mutipiers of system (.). Université Féix Houphouët-Boigny de Cocody-Abidjan. Laboratoire de Mathématiques Appiquées. UFR Mathématiques et Informatique, BP 58 Abidjan, Côte d Ivoire (traorearounagogbeyves@yahoo.ca) Université Féix Houphouët-Boigny de cocody-abidjan. Laboratoire de Mathématiques Appiquées. UFR de Mathématiques et Informatique, BP 58 Abidjan, Côte d Ivoire (mouhamadou.dosso@univ-fhb.edu.ci or mouhamadoudoss@yahoo.fr) Université Féix Houphouët-Boigny de Cocody-Abidjan. Laboratoire de Mathématiques Appiquées. UFR Mathématiques et Informatique, BP 58 Abidjan, Côte d Ivoire (k brou@hotmai.com)

2 An important property of Hamitonian system with periodic coefficients is that the matrizant X(t) R N N of (.) verifies the identity X(t) T JX(t) = J, (.) i.e. X(t) is J-orthogona or J-sympectic. These matrices were studied in [, 5, 7, 8, 9]. We reca that the spectrums of the sympectic matrices are generay divided into three groups of eigenvaues (see e.g [, 9]) : N eigenvaues outside the unite circe, N = N eigenvaues inside the unite circe and (N N ) eigenvaues on the unite circe. Considering the symmetric matrix [, 7, 8, 9, ] S = (/) ( (JW)+(JW) T), where W is a J-sympectic matrix of R N N and J a skew-symmetric matrix of system (.). S.K. Godunov and Sadkane in [9] gave a cassification of the eigenvaues which ie on the unit circe as foows Definition. An eigenvaue ρ of W on the unit circe is an eigenvaue of red coor or r-eigenvaue (respectivey an eigenvaue of green coor or g-eigenvaue) if (S x, x) > (respectivey (S x, x) < ) for any eigenvector x associated with ρ. However if (S x, x) =, then ρ is of mixed coor. From this definition, we give the foowing theorem [] Theorem. The matrix W is strongy stabe if and ony if, one of the foowing conditions is verified. W has ony r- and/or g- eigenvaues and the quantity shoud not be cose to zero. δ S = min{ e iθ e iθj : e iθ,e iθj are r- and g-eigenvaues of W} (.). P r + P g = I, and P T r S P g =, where P r and P g are the projectors associated respectivey with r-eigenvaues and g-eigenvaues of W, and S r = P T r S P r = S T r and S g = P T gs P g = S T g.. the sequence of averaged matrix ( S (n)) n defined by S(n) = n n k= (W T ) k W k converges to a positive definite symmetric constant matrix S ( ) and the quantity defined in (.) is not cose to zero. Regarding the strong stabiity anaysis of the Hamitonian system with periodic coefficients, we give the foowing theorem (see [,, 5]) Theorem. System (.) is strongy stabe if and ony if, one of the foowing conditions is verified ) If there exists ε > such that any Hamitonian system with P-periodic coefficients of the form and verifying H H P J dx(t) H(t) H(t) < ε is stabe. ) The monodromy matrix X(P) of system (.) is strongy stabe. = H(t)X(t), (.5) Thus the anaysis of the strong stabiity of a Hamitonian system with periodic coefficients is inked to the stabiity of any sma perturbation of the system preserving its structure. Which eads us to study the perturbation of these type of system. In this paper, we are interesting in a type of perturbation caed perturbation of rank k of Hamitonian system with periodic coefficients. A study of rank one perturbations was made in [] from a study of rank one perturbation of sympectic in [, 7]. In our study, we use matrices whose coumns generate Lagrangian invariant subspaces. Thus to understand the

3 rank k perturbation theory of Hamitonian systems with periodic coefficients, we give some basic properties of the isotropic subspaces in section. in section the theory of rank k perturbations of sympectic matrices is proposed. Section expains the concept of rank k perturbation of Hamitonian systems with periodic coefficients. In section 5, we anayse the Jordan canonica form of matrizant of rank k perturbation of.. In section, We give some numerica exampes which iustrate our theoretica resuts. Finay, we make some concuding remarks in section 7. Throughout the paper, we use the foowing notation: The identity and zero matrices of order k are respectivey denoted by I k, and k or just I and when the order is cear from the context. And by the symbos A and U T we denote the -norm of the matrix A and the transposed matrix (or vector) U respectivey. Some basic notions on some types of subspaces Start by basic notions on the Lagragian and isotropic subspaces.. Lagrangian subspaces These subspaces are defined as foow [7] Definition. Let J R N N be either skew-symmetric and invertibe (or in the compex case ony, Hermitian and invertibe, respectivey). A subspace L of C N is caed J-Lagrangian if it has the dimension N and < Jx,y >=, x,y L. or in the case J Hermitian if < x,y > =, x,y L where the standard biinear and sesquiinear forms are defined as foow < x,y >= < x,y > = N j= N j= x j y j, x = [x,...,x N ] T, y = [y,...,y N ] T R N, x j y j, x = [x,...,x N ] T, y = [y,...,y N ] T C N, Speciay, a subspace L is caed Lagrangian subspace if and ony if there exists a matrix L whose coumns generating L satisfies rank(l) = N and L J N L =. Consider the foowing definition Definition [. A] matrix H C N N is caed Hamitonian if JH = (JH) is Hermitian, where I J = N and the superscript denotes the conjugate transpose. I N The foowing emma gives the ink between the J-sympectic and the J-Hamitonian matrices via the Caey transform, see e.g., [, 5, 7] C (M) = (I M) (I +M), C (N) = (I +N) (I N) for M,N R N N with and not beonging to the spectrum of M and N respectivey. Lemma. (Caey Transform) Let W R N N be J-sympectic. (i) If W has not of eigenvaues, then the matrix A = C (W) = (I W) (I +W) is J-Hamitonian and ± are not eigenvaues of A. Moreover, we have W = C (A) = (A I)(A+I). (ii) If W has not of eigenvaues, then the matrix B = C (W) = (I+W) (I W) is J-Hamitonian and ± are not eigenvaues of B. Moreover, we have W = C (B) = (I B)(B +I).

4 The foowing proposition gives us a reations between the Lagrangian subspaces and the sympectic matrices (see in [, 7]) [ ] Proposition.. Let W C N N IN be a sympectic matrix. Then the coumns of W span N a Lagrangian subspace. Moreover, if the coumns of a matrix L( C[ N N ]) span a Lagrangian subspace, then there exists a sympectic matrix W IN such that range W = range(l). N. Let H C N N be a Hamitonian matrix. There exists a Lagragian ( [ invariant ]) subspace L of H if IN and ony if there exists a sympectic matrix W such that range W = L and we have the Hamitonian bock trianguar form. Isotropic subspaces W HW = [ ] R D R. The isotropic subspaces of certain types of matrices are usuay of interest in appications [, 9]. Definition. A subspace X R N is caed isotropic if X J N X. A maxima isotropic subspace is caed Lagrangian. We coect some properties on the isotropic subspaces in the theorem beow Proposition. to N.. Let X be an isotropic subspace. Then the dimension of X is ess than or equa. A isotropic subspace is contained in a Lagrangian subspaces.. Let S = [S S ] R N N be a sympectic matrix with S i R N N, i =, ; then the coumns of S and S span isotropic subspaces. Reca us two usefu emmas on the isotropic subspace [] Lemma. Let X S R N be a subspace that is invariant under a Hamitonian matrix S which has a its eigenvaues associated with X S satisfying R(λ) <. Then X S is isotropic. The beow emma gives a ink between the invariant isotropic subspaces and the ( existence of ) the orthogona sympectic matrices i.e. the matrix U which has the representation U = U U, U U U,U R n,n []. Lemma. Let S R n n be a skew-hamitonian matrix and X R n k (k n) with orthogona coumns. Then the coumns of X span an isotropic invariant subspace of S if and ony if there exists an orthogona sympectic matrix U = [X,Z,J T X,J T Z] with some Z R n (n k) so that U T SU = k n k k n k k n k k n k A A G G A G T G A T H A T A T We can buid isotropic subspaces from the methods of Kryov subspace. Reca that the Kryov subspaces are of the form K m K(A,v) = span { v,a,av,a v,...,a m v } where A R n,m and v R m. The Kryov subspace methods are : the Hermitian or skew-hermitian Lanczos agorithm and Arnodi s method and its variations. We give the foowing proposition which contains some properties of these subspaces (see [8, p. ])

5 Proposition.. The Kryov subspace K m is the subspace of a vectors in C n which can be written as x = p(a)v, where p is a poynomia of degree ess than or equa to m.. Let m be the degree of the minima poynomia of v. Then K m is invariant under A and K m = K m for a m m.. The Kryov subspace K m is of dimension m if and ony if the grade of v with respect to A is arger than m. Thus any Kryov process constructed from a skew-hamitonian matrix automaticay produces an isotropic subspace. Hence the foowing proposition (see [9, p. 99]) Proposition. Let S R N N be a skew-hamitonian matrix and u R N be an arbitrary nonzero vector. Then the Kryov subspace K j (S,u) is isotropic for a j. Rank k perturbation of sympectic matrices Consider a sympectic matrix W and a J-Lagrangian subspace L of dimension N. Let u,,u k be k vectors of L, where k N. Setting U = [u ;...;u k ], and considering the matrix we have the foowing proposition Proposition. the matrix W is J-sympectic. Proof We have the foowing inequaities W = ( I +UU T J ) W, W T J W = W T (I JUU T )(J +JUU T J)W = W T J +JUU } T J {{ JUU T J } JU (U T JU) U T J W }{{} = = = W T JW = J. The foowing proposition is a set of resuts deduced from []. Proposition. Consider the matrix Ĩ = (I +UUT J). Then ) Ĩ is J-sympectic. ) Ĩ = I UU T J. ) ) dim (ker(ĩ I) = N k, where k is the rank of U. ) σ(ĩ), where σ(ĩ) is the spectrum of Ĩ. Proof The proof is easiy deduced from those of []. From the foregoing, we give the foowing definition Definition. Let W be a sympectic matrix. We ca rank k perturbation of W, any matrix of the form W = (I +UU T J)W, (.) where U is a matrix of rank k whose coumns beong in a J-Lagrangian subspace. 5

6 Remark. The matrix W can be put in the form W = (I + k u j u T j J)W. j= More speciay, this remark shows that any rank k perturbation of W is k rank one perturbations of the sympectic matrix W. We have k ( I +uj u T j J ) k W = I + u j u T j J W. j= Consider a sympectic matrix of function (X(t)) t R ; we can consider for exampe the soution of Hamitonian system (.) which are J-sympectic. We have the foowing definition Definition. We ca rank k perturbation of X(t) any function matrix of the form j= X(t) = (I +UU T J)X(t), (.) where rank(u) = k and the coumns of U beong in a J-Lagrangian subspace. Remark. Since the function matrix (X(t)) t R is J-sympectic, its rank k perturbation wi be J sympectic. From definition. and remark., we can introduce the theory of rank k perturbation of Hamitonian system with periodic coefficients. Rank k perturbation of Hamitonian system with periodic coefficients Let U be a constant matrix of rank k such that its coumns beong in a J-Lagrangian subspace and (X(t)) t be the fundamenta soution of (.). We have the foowing proposition Proposition. Consider the foowing perturbed Hamitonian system where J d X(t) = [H(t)+E(t)] X(t) (.) E(t) = (JUU T H(t)) T +JUU T H(t)+(UU T J) T H(t)(UU T J). Then X(t) = (I +UU T J)X(t) is a soution of system (.). Proof By derivation of X(t), we obtain: J d X(t) Hence system (.) where =J(I +UU T J)J J dx(t) =J(I +UU T J)J H(t)X(t), according form system(.) =[H(t)+JUU T H(t)](I +UU T J) X(t) = H(t)+(JUU T H(t)) T +JUU T H(t)+(UU T J) T H(t)(UU T J) }{{} X(t) E(t) E(t) = (JUU T H(t)) T +JUU T H(t)+(UU T J) T H(t)(UU T J). (.) We can easiy check that E(t) is symmetric and periodic i.e. E(t) T = E(t) and E(t+P) = E(t) for a t R +. The foowing coroary gives us a simpified form of system (.), with X() = I.

7 Coroary. Equation (.) can be put in the form J d X(t) = ( I UU T J ) T ( H(t) I UU T J ) X(t), t R+, X() = I +UU T J (.) Proof Deveoping (I UU T J) T H(t)(I UU T J), we see that (I UU T J) T H(t)(I UU T J) = H(t)+ and X() = (I +UU T J)X() = I +UU T J. We give the foowing coroary (J T UU T H(t)) T +J T UU T H(t)+(UU T J) T H(t)(UU T J) }{{} E(t) Coroary. Any soution ( X(t)) t of the perturbed system (.) of system (.), is of the form X(t) = (I +UU T J)X(t), where (X(t)) t is the fundamenta soution of system (.). Proof From proposition., if X(t) is the soution of (.), then the perturbed matrix X(t) = (I+UU T J)X(t) is the soution of (.). Reciprocay, for any soution X(t) of (.), et X(t) = (I UU T J) X(t) where U is the matrix defined in system (.). Then X(t) = (I +UU T J)X(t). Repacing X(t) in (.), we get J(I +UU T J) dx(t) = (I UU T J) T H(t)X(t) (I +UU T J) T J(I +UU T dx(t) J) = H(t)X(t) }{{} =J J dx(t) = H(t)X(t) and X() = (I UU T J) X() = (I UU T J)(I + UU T J) = I. Consequenty, X(t) is the soution of (.). Remark. Basing on remark., system (.) can be written as beow J d X(t) ( = I ) T ( k j= u ju T j J H(t) I ) k j= u ju T j J X(t) X() = (I + k j= u ju T j J) (.) where each vector (u j ) j k R N beongs in a same J-Lagrangian subspace. We can immediatey see that the rank k perturbation of (.) can be interpreted as k rank one perturbations of (.). In fact, since I UU T J = I k u j u T j J = j= k ( I uj u T j J ), j= 7

8 we easiy see that system (.) can be put in the foowing form J d X(t) = X() = k j= ( k j=( I uj u T j J)) T ( k H(t) j=( I uj u T j J)) X(t) (.5) ( I +uj u T j J) which is the same as the beow system, for a p {,,...,k } : J d X(t) ( k ) T ( = j=p+ (I u ju T j J) k ) H (p) (t) j=p+ (I u ju T j J) X(t) ( k X() = j=p+ (I +u (k+p j+)u T (k+p j+) )X J) (p) () where T p p H (p) (t) = (I u j u T j J) H(t) (I u j u T j J) and X (p) () = j= j= p j= (.) (I +u (p j+) u T (p j+) J ). Now, et us interest to the Jordan canonica form of the soution ( X(t)) t of the perturbed system (.) of (.) in foowing section. 5 Jordan canonica form of ( X(t)) t> Theorem 5. Let J C N N be skew-symmetric and nonsinguar matrix, (X(t)) t> fondamenta soution of system (.) and λ(t) C an eigenvaue of X(t) for a t >. Assume that X(t) has the Jordan canonica form m(t) J n (λ(t)) J n (λ(t)) J nm(t) (λ(t)) J(t), j= j= where n > > n m(t) with m : R N a function of index such that the agebraic mutipicities is a(t) = n + + m(t) n m(t) and J(t) with σ(j(t)) C\{λ(t)} contains a Jordan bocks associated with eigenvaues different from λ(t). Furthermore, et B(t) = UU T JX(t) where U C N k is such that its coumns generate a Lagrangian subspace. () If t >, λ(t) {,}, then genericay with respect to the components of U, the matrix X(t)+ B(t) has the Jordan canonica form j= k J n (λ(t)) j= j= m(t) J n (λ(t)) J nm(t) (λ(t)) J(t), if k i k i i+ m(t) J ni (λ(t)) J ni+ (λ(t)) J nm(t) (λ(t)) J(t), if j= j= j= j= i k = s +k i s= with k i i and i >. where J(t) contains a the Jordan bocks of X(t)+B(t) associated with eigenvaues different from λ(t). () If t >, verifying λ(t ) {+,}, we have 8

9 i (a) if k = s +k i with n,n,...,n i are even and k i i, then genericay with respect to the s= components of U, the matrix X(t )+B(t ) has the Jordan canonica form i k i i+ m(t) J ni (λ(t )) J ni+ (λ(t )) J nm(t) (λ(t )) J(t ), j= j= where J(t ) contains a the Jordan bocks of X(t )+B(t ) associated with eigenvaues different from λ(t ). i (b) if k = s +k i with k i i and n i is odd, then i is even and genericay with respect s= to the components of U, the matrix X(t )+B(t ) has the Jordan canonica form i i m(t) J ni+(λ(t )) J ni (λ(t )) J nm(t) (λ(t )) J(t ), j= where J(t ) contains a the Jordan bocks of X(t )+B(t ) associated with eigenvaues different from λ(t ). Proof we reca that the rank k perturbation X(t) + B(t) of X(t) can be put on the form of k rank one perturbation (X(t)) t> by k ( X(t) = I +uk j+ u T k j+ J) X(t) j= where each vector u j are the coumns of the matrix U.. If λ(t) {,}, t, For k, we have (see [, Theorem ] ) : X = ( I +u u) T X(t) has the foowing Jordan canonica form m(t) J n (λ(t)) J n (λ(t)) J nm(t) (λ(t)) J (t), j= j= where J (t) contains a the Jordan bocks of X (t) associated with eigenvaues different from λ(t). X = ( ( j=( I +u j+ u j+) T = I +u u) T ( I +u u) T X(t) has the foowing Jordan canonica form m(t) J n (λ(t)) J n (λ(t)) J nm(t) (λ(t)) J (t), j= j= where J (t) contains a the Jordan bocks of X (t) associated with eigenvaues different from λ(t). X(t) = X k = ( ) k j= I +u k j+ u T k j+ has the foowing Jordan canonica form k m(t) J n (λ(t)) J n (λ(t)) J nm(t) (λ(t)) J k (t), j= j= where J(t) = J k (t) contains a the Jordan bocks of X(t) associated with eigenvaues different from λ(t). j= j= j= j= j= 9

10 For k = i s= s +k i with k i i ; if i =, then k = +k i. We have k ( X(t) = I +uk j+ u T k j+). j= k j=k + ( I +uk j+ u T k j+) X(t) } {{ } X where X (t) is rank one perturbations of X(t) ; then the sympectic matrix X (t) therefore has the foowing Jordan canonica form m(t) J n (λ(t)) J n (λ(t)) J nm(t) (λ(t)) J (t) j= j= using [, Theorem ]. On the other hand X(t) is k i rank one perturbations of X (t) with k < ; it therefore has the foowing Jordan form k i m(t) J n (λ(t)) J n (λ(t)) J nm(t) (λ(t)) J k (t); j= j= if i >. Putting α(i) = i s= s, we have k α(i) ( X(t) = I +uk j+ u T k j+). j= k j=k α(i)+ j= j= ( I +uk j+ u T k j+) X(t) } {{ } X α(i) where X α(i) is α(i) rank one perturbations of X(t). Using [, Theorem ], the sympectic matrix X α(i) has the foowing Jordan form i i+ m(t) J ni (λ(t)) J ni+ (λ(t)) J nm(t) (λ(t)) J α(i) (t) j= j= where J α(i) (t) contains a the Jordanbocks of X α(i) associatedwith eigenvaues different from λ(t). On the other hand X(t) is k i rank one perturbations of X (t) with k < ; it therefore has the foowing Jordan form i k i i+ m(t) J ni (λ(t)) J ni+ (λ(t)) J nm(t) (λ(t)) J k (t) j= j= where J(t) = J k (t) contains a the Jordan bocks of X(t) associated with eigenvaues different from λ(t).. Consider that there exists t > verifying λ(t ) {, }. i if k = s +k i with n,n,...,n i are even and k i i, then using [, Theorem, (a)], s= we have : the sympectic matrix X(t), k rank one perturbations of X(t), has the foowing canonica Jordan form i k i i+ m(t) J ni (λ(t)) J ni+ (λ(t)) J nm(t) (λ(t)) J(t) j= j= where J(t) contains a the Jordan bocks of X(t) associated with eigenvaues different from λ(t). j= j= j=

11 i if k = s +k i with k i i and n i is odd, then we have s= for i =, k = k and n is odd. According to (b) of [, Theorem ], is even and we have k ( X(t) = I +uk ju T ) k jj X(t) j= and step by step we have X (t) = ( I +u u T J) X(t) has the foowing canonica Jordan form n (λ(t m(t) + )) J ni (λ(t )) J nm(t) (λ(t )) J (t ), J () j= where J (t ) contains a the Jordan bocks of X (t ) associated with eigenvaues different from λ(t ). X (t) = ( I +u u T J)( I +u u T J) X(t) has the foowing canonica Jordan form j= m(t) J ni (λ(t )) J nm(t) (λ(t )) J (t ), j= using (a) of [, Theorem ] because n + is even. X [ (t) = j=( I +u j u T j J)] X(t) has the foowing canonica Jordan form n (λ(t + )) J () j= j= m(t) J n (λ(t )) J nm(t) (λ(t )) J (t ), where J(t ) contains a the Jordan bocks of X (t ) associated with eigenvaues different from λ(t ) using (b) of [, Theorem ]. X(t) [ k )] = (I +u k ju T k j J X(t) has the foowing canonica Jordan form J (k) j= n (λ(t + )) k j= j= m(t) J n (λ(t )) J nm(t) (λ(t )) J k (t ), where J(t ) = J k (t ) contains a the Jordan bocks of X(t ) associated with eigenvaues different from λ(t ). for i =, k = +k and n odd and we have k ( k X(t) = I +uk j+ u T ( k j+). I +uk j+ u T k j+) X(t) j= j=k + j= } {{ } X if n is even, then using (a) [, Theorem 9], X has the foowing Jordan canonica form m(t) J n (λ(t )) J nm(t) (λ(t )) J (t ), j= j=

12 and using the preceding point X(t) has the foowing Jordan canonica form k J (k) n (λ(t m(t) + )) J n (λ(t )) J nm(t) (λ(t )) J k (t ), j= where J(t ) = J k (t ) contains a the Jordan bocks of X(t ) associated with eigenvaues different from λ(t ). if n is odd then according (b) of [, Theorem ], is even and we deduct that X aso has the foowing form m(t) J n (λ(t )) J nm(t) (λ(t )) J (t ), j= by successivey appying rank one perturbations. Using again the previous point, we deduct that X(t) has the Jordan canonica form k J (k) n (λ(t m(t) + )) J n (λ(t )) J nm(t) (λ(t )) J k (t ), j= where J(t ) = J k (t ) contains a the Jordan bocks of X(t ) associated with eigenvaues different from λ(t ). for i >, n i is odd. Whether n, n... n i are even or odd, using successivey (a) and (b) of [, Theorem ], we deduct that X α(i) aso has the foowing Jordan canonica form i m(t) J ni (λ(t )) J nm(t) (λ(t )) J α(i ) (t ), (5.) j= j= i by successivey appying α(i) = s rank one perturbations. Since n i is odd, we affirm, s= using (b) of[, Theorem ], that i is even. To end, using the precedingpoint, we deduct that X(t) has the canonica Jordan form i k i J (ki) n (λ(t m(t) i+ )) J ni (λ(t )) J nm(t) (λ(t )) J k (t ), j= where J(t ) = J k (t ) contains a the Jordan bocks of X(t ) associated with eigenvaues different from λ(t ). j= j= j= j= i Remark 5. In point () of [, Theorem ], if k = s +k i with k i i and n i is odd, then i is even and generay with respect to the components of U, the rank k perturbation X(t ) = X(t )+B(t ) of X(t ), has the canonica Jordan form i k i m(t ) J ni (λ(t )) J (λ(t nm(t ) )) J k (t ), j= where J(t ) = J k (t ) contains a the Jordan bocks of X(t ) associated with eigenvaues different from λ(t ). j= s=

13 From points (a) and (b) of [, Theorem ], we deduce the foowing coroary i Coroary 5. Suppose there exists t > such that λ(t ) {, }. If k = s +k i and n i is even with k i i, then genericay with respect to the components of U, the matrix X(t ) + B(t ) has the Jordan canonica form i k i i+ m(t ) J ni (λ(t )) J ni+ (λ(t )) J (λ(t nm(t ) )) J(t ), (5.) j= j= where J(t ) contains a the Jordan bocks of X(t )+B(t ) associated with eigenvaues different from λ(t ). Proof j= s= If i =, we k = k and n is even. Thus, according to (a) of Theorem 5., Xk (t ) has the Jordan canonica form k J n (λ(t )) j= j= m(t ) J n (λ(t )) J (λ(t nm(t ) )) J(t ), where J(t ) contains a the Jordan bocks of X(t )+B(t ) associated with eigenvaues different from λ(t ). If i =, we have k = +k (with k ) and n is even. Thus X k (t ) = (k ) j= (I +u k j+ u T k j+j) j= j= (I +u j+u T j+j)x(t ) } {{ } X (t ) if n is even then according to (a) of Theorem 5., X (t ) has the Jordan canonica form m(t ) J n (λ(t )) J (λ(t nm(t ) )) J (t ), (5.) j= j= where J (t ) contains a the Jordan bocks of X (t ) associated with eigenvaues different from λ(t ). However X k (t ) is k rank one perturbations of J (t ) ; thus according to (a) of Theorem 5., the Jordan canonica form of X k (t ) is given by (5.). if n is odd then according to (b) of Theorem 5., the Jordan canonica form of X (t ) is given by 5.. Moreover n being even and X k being k rank one perturbations of X (t ), we obtain that the Jordan canonica form of X k (t ) is given by (5.) using (a) of Theorem 5.. i If i > then we have k = s +k i with k i and n i even. Thus X k (t ) = (k α(i )) j= s= α(i ) (I +u k j+ u T k j+ J) (I +u α(i ) j+ u T α(i ) j+ J)X(t ) j= } {{ } X α(i ) (t )

14 where α(i) = i s. j= From (a) and (b) of Theorem 5. and the above, the Jordan canonica form of Xα(i ) (t ) is given by (5.). Thus appying (a) of Theorem 5. to sympectic matrix X α(i ) (t ), we obtain the Jordan canonica form (5.) of X k (t ). Agorithm and numerica exampes We start to reca the foowing two rotation matrices [, ] I j cos(θ) sin(θ) G j,j+n = I N sin(θ) cos(θ) I N j, j N, for some θ [ π, π [ and the direct sum of two identica N N Househoder matrices (H j H j )(v, β) = I N βvv T, I N βvv T where v is a vector of ength N with its first j eements equa to zero and β a scaar satisfying β(βv T v ) =. The symbo denotes the direct sum of matrices. From these matrices, we propose Agorithm. which is the synthesis of Agorithms, and of []. This Agorithm determines a basis of an isotropic subspace from a random matrix. Agorithm. (Computation of isotropic subspace) Imput : A R N k, with N k. Output : U R isotropic subspace. (a) Q = I N (b) for j =,...,k Let x = Ae j Determine v R N and β R such that the ast N j eements of x = (H j H j )(v,β)x are zero Determine θ [ π, π [ such that (N +j)th eement of x = G j,j+k(θ) is zero Dertemine ω R N and γ R such that the (j +)th to the kth eements of x = (H j H j )(ω,γ)x are zero compute E j (x) = (H j H j )(v,β)g j,j+k (θ)(h j H j )(w,γ). Put A = E T j (x)a and Q = QE j(x) (c) U = span(q(:, : k))

15 In the foowing exampes we show that any rank k perturbation of the soution of (.) is the soution of (.). The software used for cacuating and potting the curves of the exampes beow is MATLAB 7.9.(R9b). Exampe. Consider the system of differentia equations (see [, Vo.,P.]) d η q +p η +[ǫη cos(γt)+(δcos(γt)+csin(γt))η ] = q d η +p η +[gη sin(5γt)] = d η q +p η +[(δcos(γt)+csin(γt))η +gη sin(5γt)] = which can be written down as d η +P(t)η = (.) with η p+ǫ cos(γt) δ cos(γt)+c sin(γt) q q qq η η = q and P(t) = p gsin(5γt) q qq ) η δ cos(γt)+c sin(γt) gsin(5γt) p qq qq q q Putting ( ) η(t) X(t) = dη(t), J = We get a canonica Hamitonian system J dx(t) ( ) I, and H(t) = P(t) I I. (.) = H(t)X(t), X() = I (.) where H(t) = H(t+ π 7 ) = H T (t). In this exampe, we take γ = 7, q = q = q =, p =, p =, p =, a = g = ǫ, b = δ and c =. From a random matrix A R, we deduce a matrix U R k of rank k whose coumns generate an isotropic subspace using agorithm. Consider the perturbed system (.) of (.). We show that the rank k =, perturbation of the fundamenta soution of (.) is the soution of perturbed system (.). For that, consider = X(t) X (t), t where X(t) is the soution of (.), and X (t) = (I+UU T J)X(t). We show by numerica exampes that is very cose to zero, t [, π 7 ] for ǫ = and δ =, consider the random matrix A =. Using Agorithm to matrix A, we obtain the matrix V = whose coumns span an isotropic subspace. 5

16 Let s take U = V (:, : ). In Figure, we consider the matrix U of rank which permits to perturb system (.) by the matrices U, U, U, and U. We remark that a the π e ]. figures are so that.5. This proves that X(t) X (t), t [, 7 ε= and δ= with U x ε= and δ= with U x ε= and δ= with U.5 x ε= and δ= with U x Figure : Comparison of two soutions (Exampe ) However, unperturbed system (.) is strongy stabe. We remark that the rank perturbed system (.) of (.) is unstabe for the matrix U of rank and remains strongy stabe for a matrix taken in { U, U, U }. Tabe gives the different norms of projectors, the quantity δs and a convergence iustration of S (n). Tabe : Checking of the (strong) stabiity of (.) by the approachs defined in [, 5] (Exampe ) (n) ks k δs tr(p ) kp P k tr(p ) kp P k tr(pr ) kpr Pr k tr(pg ) kpg Pg k kpr + Pg I k U U U U Tabe justifies the existence of a neighborhood in which any rank perturbation of the system remains strongy stabe. In Figure, we consider U = V to perturb system (.). We can see that π e ]. for a the figures. This shows that X(t) X (t), t [, 7 In this exampe, the unperturbed system is strongy stabe for a U taken in { V, V, V } and not stabe when U = V. This is iustrated in Tabe which gives the norms of different projectors, the quantity δ and a convergence iustration of S (n). The second Tabe justifies the existence of a neighborhood in which any rank perturbation of the system remains strongy stabe.

17 ε= and δ= with V x ε= and δ= with V x ε= and δ= with V x ε= and δ= with V x Figure : Comparison of two soutions (Exampe ) Tabe : Checking of the (strong) stabiity of (.) by the approachs defined in [, 5] (Exampe ) (n) ks k δs tr(p ) kp P k tr(p ) kp P k tr(pr ) kpr Pr k tr(pg ) kpg Pg k kpr + Pg I k U U U U for ǫ = 5 and δ =, we consider the random matrix A =. Using Agorithm. to matrix A, we obtain the matrix V = of rank whose coumns span an isotropic subspace. Let s take U = V (:, : ). Figure is obtained for vaues of any matrix of rank taken in {U, U, U, U }. We remark that a the figures of Figure, verify π e.. This shows that X(t) X (t), t [, ]. 7 In this exampe, the unperturbed system is unstabe, and the rank perturbation systems remain unstabe for any matrix of rank taken in {U, U, U, U }. This is 7

18 . x. ε=5 and δ= with U. x. ε=5 and δ= with U x 7 ε=5 and δ= with U x ε=5 and δ= with U Figure : Comparison of two soutions. iustrated in foowing Tabe Tabe : Checking of the (strong) stabiity of (.) by the dichotomy approach (Exampe ) U U U U S (n) δ S tr(p ) P P tr(p ) P P Thus there doesn t exist of a neighborhood of the unperturbed system in which any rank perturbation of the system is stabe. Taking U = V, Figure shows that X(t) X (t), t [, π 7 ] for any matrix U of rank taken in {U, U, U, U}. In the first two subfigures of Figure, we see that., whie in the other subfigures, we note that 8. However Tabe shows that the perturbedsystem is not stabe for anymatrix taken in {U, U, U, U,O Tabe : Checking of the (strong) stabiity of (.) by the dichotomy approach (Exampe ) U U U U S (n) δ S tr(p ) P P tr(p ) P P

19 . x. ε=5 and δ= with V. x. ε=5 and δ= with V x 7 ε=5 and δ= with V x 5 ε=5 and δ= with V Figure : Comparison of two soutions. Exampe. Consider the foowing differentia system: d β +(+acos(7t))β +bβ cos(t) = d β +(a+bsin(t))β +aβ sin(5t) =, (.) d β +β +bβ cos(t)+aβ sin(5t) = where a R and b R are rea parameters. Let β = β β x = β ( β dβ System (.) can be written as a Hamitonian of the form (.) with T = π 7 and +acos(7t) bcos(t) ( ) P(t) = a+bsin(t) asin(5t) P(t) and H(t) =. b cos(t) a sin(5t) I We show that the rank k =, perturbation of the fundamenta soution of (.) is the soution of its rank k =, perturbation system. Consider ). = X(t) X (t), t [, π 7 ] where X (t) = (I + UU T J)X(t) and X(t)) t [, π 7 ] is the soution of the rank k =, perturbation Hamitonian system (.) of (.). the foowing figures represent the norm of the difference between X (t) and X(t). 9

20 for a = and b =, consider the random matrix A = Appying Agorithm. to matrix A, we get the foowing matrix V = of rank whose coumns generate an isotropic subspace. Considering the matrix U = V (:, : ), we get Figure 5 perturbing system (.) by matrices taken in {U, U, U, U }. e In Figure 5, we note that a the figures verify.5. This shows that X(t) ]. X (t), t [, π 7 a= and b= with U x a= and b= with x U a= and b= with U x a= and b= with U x Figure 5: Comparison of two soutions (Exampe ) In this first exampe, the unperturbed system is strongy stabe and the rank perturbation of the system is aso strongy stabe for any matrix of rank beonging to { U, U, U } and is unstabe for any matrix of rank with U. This discussion is summaries in Tabe 5

21 Tabe 5: Checking of the (strong) stabiity of (.) by the approachs defined in [, 5] (Exampe ) U (n) ks k, (n = ) δs tr(p ) kp P k tr(p ) kp P k tr(pr ) kpr Pr k tr(pg ) kpg Pg k kpr + Pg I k U U U U This justifies the existence of a neighborhood of the unperturbed system in which any random rank perturbation of the system remains strongy stabe. Let s take U = V ; Figure shows that <.5, t [, e This shows that X(t) X (t). for a the figures. a= and b= with V a= and b= with V x π 7 ] x a= and b= with V x a= and b= with V x Figure : Comparison of two soutions (Exampe ) In this case, the unperturbed system (.) is strongy stabe for any random matrix U of rank beonging to {V, V, V, V }. This is iustrated in Tabe

22 Tabe : Checking of the (strong) stabiity of (.) by the approachs defined in [, 5] (Exampe ) V V V V V S (n), (n = ) δ S tr(p ) P P tr(p ) P P tr(p r) P r P r tr(p g) P g P g P r +P g I This justifies the existence of a neighborhood of the unperturbed system in which any random rank perturbation of the system remains strongy stabe. For a = 8.95 and b =, consider the random matrix A = Appying the agorithm. to the matrix A, we have the foowing random matrix V = of rank whose coumns generate an isotropic subspace. Let s take U = V(:, : ). The foowing Figure shows that X(t) π X (t), t [, 7 ], for any matrix of rank beonging to {U, U, U, U}. Thus in Figure 7, we can observe that π, t [, 7 ] for a figures.

23 a=8.95 and b= with U a=8.95 and b= with U x x a=8.95 and b= with x U.5 a=8.95 and b= with x U Figure 7: Comparison of two soutions (Exampe ) In this case, the unperturbed system is unstabe and its rank perturbation systems remain unstabe for any matrix of rank taken in {U, U, U, U }. This is iustrated in Tabe 7 Tabe 7: Checking of the (strong) stabiity of (.) by the approachs defined in [, 5] (Exampe ) (n) ks k, (n = ) tr(p ) kp P k tr(p ) kp P k U U U U U This justifies the existence of a neighborhood of the unperturbed system in which any rank perturbation of the system remains unstabe. In this atter exampe, we consider U = V to perturb system (.). Figure 8 is obtained for vaue of any random matrix U of rank taken in {U, U, U, U }. We can see that π e 5.5, t [, π 7 ] for a figures. Hence, we have X(t) X (t), t [, 7 ].

24 a=8.95 and b= with V a=8.95 and b= with V x 5 x a=8.95 and b= with x V.5 a=8.95 and b= with x V Figure 8: Comparison of two soutions (Exampe ) However the foowing Tabe 8 shows that the perturbed system is not stabe for any random matrix U of rank taken in {U, U, U, U, O, }. Tabe 8: Checking of the (strong) stabiity of (.) by the approachs defined in [, 5] (Exampe ) (n) ks k, (n = ) tr(p ) kp P k tr(p ) kp P k 7 U U U U U Concuding remarks In this research work, after defining a rank k perturbation theory of a Hamitonian system with periodic coefficients with k, we showed that the soution of its rank k perturbation is the same as the rank k perturbation of the soution of unperturbed system. Then we anayzed Jordan canonica form of the soution of the unperturbed system when it is subjected to a rank k perturbation. This anaysis is a generaization of that made by M. Dosso, et a. in [] in the case of a rank one pertubation of Hamitonian system with periodic coefficients. Finay we proposed numerica exampes which confirm this theory. However, these exampes use an agorithm that randomy constructs an isotropic subspace basis. From these numerica exampes we notice that when a system is strongy stabe (respectivey unstabe), there exists a neighborhood in which any rank k perturbation of the system in this neighborhood remains strongy stabe (respectivey unstabe) In future work, we wi compare the zone of stabiity (strong) of the Hamitonian systems with periodic coefficients and their rank k perturbations. Then it woud be boring to find a ink between any random perturbation and rank k perturbation of Hamitonian system with periodic coefficients. References [] C. Brezinski, Computationa Aspects of Linear Contro, Kuwer Academic Pubishers,.

25 [] M.Dosso, T.G.Y. Arouna., and J.C.Koua. Brou, On rank one perturbations of Hamitonian system with periodic coefficients. Wseas Transations on Mathematics. Voume 5,, Pages 5-5. [] M. Dosso and M. Sadkane. On the strong stabiity of sympectic matrices. Numerica Linear Agebra with Appications () (), -9. [] M.Dosso and N. Couibay, Sympectic matrices and strong stabiity of Hamitonian systems with periodic coefficients. Journa of Mathematica Sciences : Advances and Appications. Vo. 8(), pp 5-8. [5] M. Dosso, N. Couibay and L. Samassi, Strong stabiity of sympectic matrices using a spectra dichotomy method, Far East Journa Appied Mathematics. Vo. 79, No,, pp 7-. [] G. Freiing, V. Mehrmann, and H. Xu. Existence, uniqueness and parametrization of Lagrangian invariant subspaces. SIAM J. Matrix Ana. App., :59,. [7] Godunov SK, Verification of boundedness for the powers of sympectic matrices with the hep of averaging. Siber. Math. J. 99 ; : [8] S.K. Godunov, Stabiity of iterations of sympectic transformations, Siberian Math. J., 5- (989). [9] S.K. Godunov, M. Sadkane, Numerica determination of a canonica form of a sympectic matrix, Siberian Math. J., 9-7 () [] S.K. Godunov, M. Sadkane, Spectra anaysis of sympectic matrices with appication to the theory of parametric resonance, SIAM J. Matrix Ana. App. 8, 8-9 () [] I. Gohberg, P. Lancaster, and L. Rodman. Indefinite Linear Agebra and Appications. Birkhuser, Base, 5. [] B. Hassibi, A. H. Sayed, T. Kaiath, Indefinite-Quadratic Estimation and Contro, SIAM, Phiadephia, PA, 999. [] D. Kressner, Perturbation bounds for isotropic invariant subspaces of skew-hamitonian matrices. SIAM J. Matrix Anaysis Appications (): 97-9, 5. [] D. Kressner, Numerica Methods for Genera and Structured Eigenvaue Probems. Lecture Notes in Computationa Science and Engineering, Springer 5. ISBN , pp. I-XIV, -. [5] P. Lancaster and L. Rodman. The Agebraic Riccati Equation. Oxford University Press, Oxford, 995. [] C. Meh, V. Mehrmann, A. C. M. Ran and L. Rodman. Eigenvaues Perturbation theory of structured matrices under generic structured rank one perturbations ; Sympectic, othogona and unitary matrices. BIT, 5(), [7] C. Meh, V. Mehrmann, A. C. M. Ran and L. Rodman. Perturbation anaysis of Lagrangian invariant subspaces of sympectic matrices. Linear and Mutiinear Agebra, 57:-8, 9. [8] Y. Saad, Numerica methods for arge eugenvaue probems. SIAM, nd ed.,. [9] D. S. Watkins, The matrix eigenvaue probem. GR and Kryov Subspace Methods,SIAM, Phiadephia, 7. [] V.A. Yakubovich, V.M. Starzhinskii, Linear differentia equations with periodic coefficients, Vo. &., Wiey, New York (975) [] YAN Qing-you. The properties of a kind og random sympectic matrices. Appied mathematics and Mechanics. Vo, No 5, May. 5