Conformal Mapping among Orthogonal, Symmetric, and Skew-Symmetric Matrices

Size: px
Start display at page:

Download "Conformal Mapping among Orthogonal, Symmetric, and Skew-Symmetric Matrices"

Transcription

1 AAS Conformal Mapping among Orthogonal, Symmetri, and Skew-Symmetri Matries Daniele Mortari Department of Aerospae Engineering, Texas A&M University, College Station, TX Abstrat This paper shows that Cayley Transforms, whih map Orthogonal and Skew- Symmetri matries, may be onsidered the extension to matrix field of the omplex onformal mapping funtion f 1 (z) = 1 z. Then, by using a set 1 + x of real matries whih are, simultaneously, Orthogonal and Symmetri (the Ortho Sym matries), it similarly shows how to extend two omplex onformal mapping funtions (namely, the f (z) = i z i + z, and the f 3(z) = 1 z 1 + z i - here alled the lokwise, and the ounter lokwise funtions), to matrix field. This extension onsists of some new one-to-one mapping relationships between Orthogonal and Symmetri, and between Symmetri and Skew-Symmetri matries. This new relationships omplete the piture of the one-to-one matrix mapping among Orthogonal, Symmetri, and Skew- Symmetri matries. Finally, this paper shows how to map among Orthogonal, Symmetri, and Skew-Symmetri matries, by means of a diret produt. Introdution Cayley Transforms are a beautiful and useful tool to perform one-to-one mapping between Orthogonal C and Skew-Symmetri Q real matries. Cayley Transforms onsists of two formally idential relationships that an be written as { C = (I Q)(I + Q) 1 = (I + Q) 1 (I Q) (Forward) Q = (I C)(I + C) 1 = (I + C) 1 (I C) (Inverse) (1) where I indiates the identity matrix, and C T C = I, and Q = Q T. Assoiate Professor of Aerospae Engineering, Department of Aerospae Engineering, 701 H.R. Bright Building, Room 741A, Texas A&M University, 3141 TAMU, College Station, Texas Tel. (979) , FAX (979) AIAA and AAS Member. 1

2 Several appliations of these equations an be found in attitude dynami and ontrol. In fat, if C represents the diretion osine matrix desribing a given attitude (or rigid rotation), then the mapped Q matrix ontains the three elements of the Gibbs vetor, as the off-diagonal elements, a minimum number parameters to desribe the attitude. Moreover, it is easy to derive the following relationships (See pp of Ref. 1), whih relate the angular veloity ω with the rate of hange of C and Q ω = (I + Q) 1 Q(I Q) 1 = ĊCT = 0 ω 3 ω ω 3 0 ω 1 ω ω 1 0 () Another very important appliation of Cayley Transforms is for real Symmetri S matries. In fat, any Symmetri matrix an be deomposed as S = CΛ S C T = (I Q)(I + Q) 1 Λ S (I Q) 1 (I + Q) (3) that is, a minimal parametrization to desribe Symmetri matries. Equations (1,,3) an be speialized to any n-dimensional Eulidean spae. Finally, Eqs. (1) are extensively applied in all the fields where matrix transformation plays a key role. An important aspet of the Cayley Transforms, easily to demonstrate, is the fat that C and Q share the same omplex eigenvetor matrix W. In fat, let CW = W Λ C be the spetral deomposition of C, then, sine W W = I, we an write Q = (I C)(I + C) 1 = (W IW W Λ C W )(W IW + W Λ C W ) 1 = = W (I Λ C )W W (I + Λ C ) 1 W = W (I Λ C )(I + Λ C ) 1 W = = W Λ Q W where Λ Q = (I Λ C )(I + Λ C ) 1, that in salar form beomes λ Q = 1 λ C 1 + λ C (4) whih shows that the eigenvalues of C and Q are related through a bilinear onformal transformation. Therefore, Cayley Transforms represent the extension to matrix field of the partiular omplex onformal mapping funtion f 1 (z) = 1 z. The mapping 1 + z rule of this funtion, here alled Cayley funtion f 1 (z), is shown in Fig. 1, where the presene of a mirror helps to understand the mapping geometry. The relationship between Cayley mapping and the omplex onformal funtion f 1 (z), arises the question whether it is possible to develop similar matrix mapping with different omplex onformal funtions. In partiular, we will look to develop the Symmetri matries are widely used in mehanis (as representing inertia and stress tensors, or as mass, stiffness, and damping matries), as well as in dynamis and ontrol (as representing ovariane matries in Kalman filter, gain matries in optimal ontrol, and as solution of the matrix Lyapunov and the Riati differential equations).

3 missing similar matrix mapping to omplete the mapping between the three most ommon real n n matries: the Orthogonal C, the Symmetri S, and the Skew- Symmetri Q matries, respetively. To this end, the problem to share the real eigenvetor matrix P of S with the omplex eigenvetor matrix W of C and Q, arises. Spetral deomposition properties of C, Q, and S are summarized in Table 1. Matrix Definition Eigenvalues Eigenvetors Orthogonal C T C = I Unit-Cirle Unitary Skew-Symmetri Q = Q T Imaginary Axis Unitary Symmetri S = S T Real Axis Orthogonal Table 1: Spetral deomposition properties for C, Q, and S matries. The eigenvetor problem will be solved thanks to the Nuleus matrix M, while the eigenvalue problem (Whih are the mapping funtions to map the eigenvalues) will be solved by introduing the Counter Clokwise f (z), and the Clokwise f 3 (z) omplex mapping funtions. However, prior to enter into these arguments, the Ortho-Skew, and the Ortho-Sym matries, will be introdued. The Ortho-Skew I, and the Ortho-Sym R matries The intersetion of the unit irle (eigenvalue lous of the Unitary matries) and the imaginary axis (eigenvalue lous of the Skew-Hermitian matries) yields to matries whih are, simultaneously, Unitary and Skew-Hermitian. These matries an be, therefore, defined by the relationships { I I = I, I = I } = II = I (5) The resulting matrix I is, therefore, Unitary and Skew-Hermitian matrix set, here alled improperly Ortho-Skew, whih has been introdued in Ref. 3, presents very interesting properties. The general expression of these matries, whose eigenvalues are only pure imaginary ±i, is I = i p w k w k i n k=p+1 w k w k (w i w j = δ ij ) (6) where the omplex diretions w k, whih form the unitary eigenvetor matrix W onstitute an Orthogonal set, and where p indiates the number of positive eigenvalues +i (0 p n). Note that, sine the w k are Orthogonal, then if p = n then I = i I, while if p = 0 then I = i I. 3

4 The Ortho-Skew-Hermitian matries are, in general, omplex. However, in even dimensional spaes, it is possible to build real I e matries as follows I e = n/ n/ W k I W k = w k w k that an be simply alled Ortho-Skew. i 0 0 i Based on Eq. (5), subsequent powers of I satisfy the following rule = +I (k = 1, 5, 9, ) I k = I (k =, 6, 10, ) = = I (k = 3, 7, 11, ) = +I (k = 4, 8, 1, ) w k w k (7) Analogously, the intersetion of the unit-radius irle (eigenvalue lous of C) with the real axis (eigenvalue lous of S) yields to the definition of the Ortho-Sym R matries, whih are, simultaneously, Orthogonal and Symmetri. The Ortho-Sym real matries R, whih have been introdued in Ref. 4, are defined by their own definitions (8) R T R = I and R = R T (9) Subsequent powers of these matries obey to the simple rule R k = { = R (k = odd) = I (k = even) (10) The general expression of R, whose eigenvalues are only ±1, is p R = p k p T k n p k p T k (p T i p j = δ ij ) (11) k=p+1 where the real diretions p k form an orthogonal set and where p indiates the number of positive eigenvalues +1 (0 p n). Sine the p k are orthogonal, then if p = n then R = I, while if p = 0 then R = I. These matries, in partiular the Ortho-Sym matries (that an be built assoiated with a given C or S), will allow us to omplete the one-to-one mapping between Orthogonal, Skew-Symmetri, and Symmetri matries. The Nuleus Matrix The nuleus matrix is the omplex matrix 1 i M = 1 i (1) Also the omplex I matries will be here improperly alled Ortho-Skew. 4

5 This unitary matrix (MM = I ) allows to extrat the information of the plane, provided by a omplex-onjugate eigenvetor pair W k = ˆp 1k + iˆp k. ˆp 1k iˆp k = ŵ k. ŵ k, into a real plane defined by a real eigenvetor pair P k = ˆp 1k. ˆp k. In fat 1 i W k M = ˆp 1k + iˆp k. ˆp 1k iˆp k = ˆp 1 i 1k. ˆp k = P k (13) This property allows us to map a omplex unitary eigenvetor matrix W into an assoiated real Orthogonal matrix P. In fat, M M 0 P = P 1. P. p = W 1. W. p.... = W D(M) (14) where D(M) is a diagonal-blok matrix ontaining, as a diagonal elements, M matries assoiated with the n/ omplex eigenvalue pairs, and a one assoiated to the real eigenvalue (if n is odd, only). The nuleus matrix M allows us to build Ortho-Skew and Ortho-Sym matries, I(W ) and R(W ), whih are assoiated with a given unitary eigenvetor matrix W. In fat, for a even dimensional spae, we may write that where and where I(W ) = R(W ) = I = R = n/ n/ +i 0 0 i n/ W k I W k = P k R P T k = and n/ and P k M I M P T k = n/ W k M R M W k = n/ I (R) = M I M = R (I) = M R M = P k I (R) P T k W k R (I) W k The nuleus matrix M has many important properties. First of all, it allows to build all the three Pauli Spin Matries (See pp of Ref. ) σ 1 = = M M T σ 3 = (15) (16) (17) = M T M (18) and σ = M M M T M T M M T = 5 0 i i 0 (19)

6 whih are used in the omplex Cayley-Klein matrix to parameterize rigid rotation, and in quantum mehanis where they represent the angular-momentum matries for spin ( h). The seond property (really unusual), is that M as well as M T and M, all present a power periodiity of 4 whih allows us to write M 4 = (M T ) 4 = (M ) 4 = I (0) M 4l+k = (M T ) 4l+k = (M ) 4l+k = M k = (M T ) k = (M ) k (1) where l and k represent any integer. Conformal Complex Transformations A onformal mapping, also alled a onformal map, onformal transformation, anglepreserving transformation, or biholomorphi map, is a transformation w = f(z) that preserves loal angles. Let us onsider the following interesting omplex onformal mapping funtions (See pp of Ref. ), together with their inverses w = f 1 (z) = 1 z z = f 1 (w) = 1 w 1 + z 1 + w w = f (z) = i z i + z w = f 4 (z) = i z 1 + z z = f 3 (w) = 1 w 1 + w i z = f 4 (w) = i w 1 + w () w = f 5 (z) = 1 z i + z z = f 6 (w) = 1 iw 1 + w Table 1 summarizes the mapping results for these funtions when applied to the real axis (z = ±Φ), the imaginary axis (z = ±iφ), and the unit-irle (z = e ±i Φ ). From Real axis From Imaginary axis From Unit-Cirle f 1 (z) Real Unit-Cirle Imaginary f (z) Unit-Cirle Real Imaginary f 3 (z) Imaginary Unit-Cirle Real f 4 (z) I = R + 1 Cirle at 1, +1/ I = R f 5 (z) Cirle at 1, 1/ I = R 1 I = R f 6 (z) I = R 1 Cirle at +1, 1/ I = R Table 1: Mapping summary of w = R + i I = f(z). 6

7 It is easy to demonstrate that the first four of these omplex mapping funtions obey to the following general rules { { z = f1 (f 1 (z)) z = f (f and (f (z))) (3) z = f 4 (f 4 (z)) z = f 3 (f 3 (f 3 (z))) In the following setions, how to extend the first three omplex onformal funtions, f 1 (z), f (z), and f 3 (z), to matrix mapping, is shown. The analysis for the remaining three funtions, f 4 (z), f 5 (z), and f 6 (z), whih are muh more ompliated, will be the subjet of a future work. The Cayley mapping funtion f 1 (z) = 1 z 1 + z As already stated, Cayley Transforms, given in Eq. (1), an be seen as the extension to matrix field of this mapping funtion. In fat, the eigenvalues of C and Q are related through the omplex trigonometri identity ( ) Φ ( ) Φ i tan = 1 e±iφ e ±iφ = 1 + e ±iφ i tan ( ) (4) Φ i tan The geometry assoiated with this trigonometri identity or with the Cayley mapping funtion f 1 (z), is shown in Fig. 1, where the point 1, 0 works like a mapping refleting point. Φ / Im i tan Φ Φ e iφ Figure 1: f 1 (z) mapping rules Re Q Im U S S 1 Cayley funtion Re Figure : f 1 (z) = 1 z 1 + z funtion Figure summarizes the mapping properties of the f 1 (z) funtion. In partiular, f 1 (z) maps the real axis onto the real axis. This property an also be seen as the 7

8 trigonometri identities os Φ = ( ) Φ 1 tan ( ) Φ ( ) tan Φ 1 + tan = 1 os Φ 1 + os Φ (5) whose equivalent matrix mapping has the form { S1 = (I S )(I + S ) 1 = (I + S ) 1 (I S ) S = (I S 1 )(I + S 1 ) 1 = (I + S 1 ) 1 (I S 1 ) (6) where S 1 and S are two Symmetri matries. These two matries have the same ( ) Φ eigenvetor matrix, and whose eigenvalues are λ k = os Φ k, and tan, respe- tively. Trigonometri/Matrix relationships It is well known that any real square matrix an be deomposed as a sum of a Skew- Symmetri and a Symmetri matries. Speified to an orthogonal matrix C, this deomposition is as follows ( ) ( ) C + C T C C T C = + = S + A (7) where S and A have the eigenvalue pairs os Φ k, and ±i sin Φ k, respetively. ( ) Therefore, it is possible to apply the first of Eq. (6) with S = S =. The C + C T transformed matrix S 1, as it easy to demonstrate, an be set as S = QQ T, whih satisfies the two requested properties of having the same eigenvetor matrix of C, and whose eigenvalue pairs are obtained as Therefore, Eq. (6) beomes ( ) ( ) ( ) Φ Φ Φ i tan ±i tan = tan (8) { C + C T = (I QQ T )(I + QQ T ) 1 = (I + QQ T ) 1 (I QQ T ) QQ T = (I C C T )(I + C + C T ) 1 = (I + C + C T ) 1 (I C C T ) (9) Analogously, the trigonometri property ( ) Φ i tan ±i sin Φ = ( ) (30) Φ 1 + tan 8

9 allows us to write the assoiated matrix identity A = C C T = Q (I + QQ T ) 1 = (I + QQ T ) 1 Q (31) (the sign omes from the fat that, assoiated with Q we have the eigenvalues ( ) Φ i tan, while assoiated with A, we have the eigenvalues ±i sin Φ). Analogously, the definition of the tangent trigonometri funtion, whih an be written as ( ) Φ ±i sin Φ i tan ±i tan Φ = os Φ = ( ) (3) Φ 1 tan allows us to introdue a Skew-Symmetri matrix T that an be seen as a tangent matrix T = (C C T )(C + C T ) 1 = (C + C T ) 1 (C C T ) (33) that an also be expressed or deomposed as T = Q (I QQ T ) 1 = (I QQ T ) 1 Q (34) Analogously, the expression of the tangent half-angle ( ) Φ ±i sin Φ i tan = 1 + os Φ = 1 os Φ ±i sin Φ allows us to introdue the matrix identities { Q = (C C T )(I + C + C T ) 1 = (I + C + C T ) 1 (C C T ) Q = (I C C T )(C C T ) 1 = (C C T ) 1 (I C C T ) (35) (36) Equations (9), (31), (33), (34), and (36) and the assoiated trigonometri relationships given in Eqs. (8), (30), (3), and (35), are, therefore, stritly onneted. Note that, some of the above relationships represent a different way to write the Cayley Transforms. It is lear that matrix properties are substantially more omplex than trigonometri properties. However, the above relationships allows us to see the Orthogonal, Symmetri, and Skew-Symmetri matries as illuminated with an additional dimmer light. This small additional knowledge will ertainly help to understand the nature, and to develop the Matri Trigonometry, of whih the known one (the salar), only onstitutes a small subset. It is easy to demonstrate that the Ortho-Sym matrix R represents the key matrix to map Symmetri onto Skew-Symmetri matries, and vieversa. In fat, the relationship S = R(W ) A = A R(W ) A = R(W ) S = S R(W ) (37) where R(W ) means a Ortho-Sym matrix R whih is built with the eigenvetor matrix W of A, while the eigenvetor P = W D(M) of S is obtained by Eq. (14). 9

10 The Counter Clokwise funtion f (z) = i z i + z = 1 + iz 1 iz To extend the appliation of f (z) to matrix field (see Fig. 3 for its mapping summary), it is possible to demonstrate that, for the Symmetri matries S, the extension beomes { C = (I + RS)(I RS) 1 = (I RS) 1 (I + RS) C T = (I + RS) 1 (I RS) = (I RS)(I + RS) 1 (38) where R = R(W ) is the Ortho-Sym matrix assoiated with the eigenvetor matrix of C. The previous equation an be written also in the form { C = (R S)(R + S) 1 = (R S) 1 (R + S) C T = (R + S)(R S) 1 = (R + S) 1 (39) (R S) these equations implies the trigonometri identities ( ) Φ i tan ( ) = os Φ ± i sin Φ (40) Φ i ± tan Just for example, let us to demonstrate Eq. (39). If, as it happens in Cayley Transforms (that is, the formal struture of the eigenvalue mapping is idential to the matrix mapping), then we need to look for an unknown matrix X suh that C = (X S) (X + S) 1 is Orthogonal. From this equation we an easily derive that the unknown matrix X has the expression X = (I C) 1 (I + C) S. ( ) Φ Setting os Φ, s sin Φ, and t tan, we an write the expression for X in the -D ase 1 1 s 1 + s X s 1 s 1 + t 0 0 t = this matrix, whih is independent from Φ, is Symmetri, and Orthogonal, with eigenvalues ±1. This implies the onlusion that the unknown matrix X is the Ortho-Sym matrix R. Demonstration of Eq. (38) may be done on the same way. The f (z) mapping funtion, here named Counter Clokwise (beause it moves from real axis, to imaginary axis, to unit irle, and bak to real axis), annot be applied to Orthogonal and to Skew-Symmetri matries beause f (z) f (z ). In partiular, this mapping funtion satisfies f (e i Φ )f (e i Φ ) = 1 and f (i Φ)f ( i Φ) = +1 (41) 10

11 Im Counter Clokwise U funtion Im U Clokwise funtion Q Q Q U Q U Re Re S S S S Figure 3: f (z) = i z i + z = 1 + iz 1 iz funtion Figure 4: f 3 (z) = 1 z 1 + z i funtion The Clokwise funtion f 3 (z) = 1 z 1 + z i The mapping funtion f 3 (z) (see Fig. 4 for its mapping summary), an be used to find a way to map Orthogonal to Symmetri. The resulting mapping is S = (I C)(I + C) 1 R = (I + C) 1 (I C)R (4) whose demonstration is easily derived from Eq. (66) whih, in turn, an be seen as the a way to map Skew-Symmetri to Symmetri and vieversa S = QR = RQ and Q = SR = RS (43) The f 3 (z) mapping funtion, here named Clokwise (beause it moves from real axis, to unit-irle, to imaginary axis, and bak to real axis), annot be applied to Symmetri matries beause f 3 (Φ) f 3 ( Φ), and annot be applied to Skew-Symmetri matries beause f 3 (i Φ) f 3 ( i Φ). In partiular, this mapping funtion satisfies f 3 (Φ)f 3 ( Φ) = 1 and f 3 (i Φ)f 3 ( i Φ) = 1 (44) Mapping Summary All the above relationship are summarized below: From Orthogonal to Skew-Symmetri: Q = (I C)(I + C) 1 = (I + C) 1 (I C) (45) 11

12 From Skew-Symmetri to Orthogonal: C = (I Q)(I + Q) 1 = (I + Q) 1 (I Q) (46) From Orthogonal to Symmetri: S = (I C)(I + C) 1 R = (I + C) 1 (I C)R (47) From Symmetri to Orthogonal: C = (R S)(R + S) 1 = (R S) 1 (R + S) (48) From Skew-Symmetri to Symmetri: S = QR = RQR (49) From Symmetri to Skew-Symmetri: Q = SR = RS (50) The above set of equations onstitute the omplete matrix mapping relationships among Orthogonal, Skew-Symmetri, and Symmetri matries. From this set of equations it is possible to demonstrate the following R C R = C T R Q R = Q = Q T R S R = S and I C I = C I Q I = Q = Q T I S I = S (51) whih yield to the relationships R = i W D(I ) W and I = i W D(R ) W (5) where D(I ), and D(R ) are the n n diagonal matries obtained with the blok elements I, and R, respetively. Diret Produt Mapping The purpose of this setion is to find the relationships allowing one-to-one mapping among Orthogonal, Skew-Symmetri, and Symmetri matries, by a diret produt. This implies, for instane, to find out whih is the matrix G suh that S = G C maps a given orthogonal matrix C onto a symmetri matrix S, whih has assigned eigenvalues. Here again, the problem to share the eigenvetors, is solved using the Nuleus matrix. 1

13 In fat, setting k = os Φ k and s k = sin Φ k, the spetral deomposition of a given Orthogonal matrix C an be written as C = whih an also be written as W k Λ k W k + p pt where Λ k = k + is k 0 0 k is k (53) where n p C = W k MM Λ k MM W k + p pt = W k M = P k and M Λ k M = os Φ k sin Φ k P k C k P T k + p p T (54) sin Φ k = C os Φ k (55) k Therefore, the spetral deomposition of real Orthogonal matries an also be written in the full real form C = P k C k P T k + p p T where C k = Similarly, Skew-Symmetri matries an be expanded as A = P k A k P T k where A k = while for the Symmetri matries we have S = P k S k P T k where S k = Thanks to the eigenvetor orthogonality k s k 0 a k a k 0 rk 0 0 r k s k k (56) (57) (58) P T i P j = 0, i j; P T i P I ; P T i p = 0,1 ; p T p = 1. (59) the expressions provided by Eqs. (56-58) allows us to obtain ompat formulations for subsequent produts of Orthogonal, Symmetri, and Skew-Symmetri matries. For instane, using these expansions, the produt between an Orthogonal and a Symmetri matrix beomes C S = np i=1 P i C i P T i n p + p p T P j S j P T j (60) j=1 where all the terms with i j vanish beause of Eq. (59). Thus, C S = P k C k S k P T k = 13 P k F k P T k (61)

14 where F k = C k S k. This result provides us with the tool to develop general matrix mapping by means of a diret produt. In fat, the G matrix that maps the matrix B onto the final matrix F = G B an be built as follows G = F B 1 = P k G k P T k where G k = F k B 1 k (6) This proedure holds provided that the matries B k are not singular. The inverse mapping transformation (from F onto B) is simply G 1 = B F 1 = P k G 1 k P k T where G 1 k = B k Fk 1 (63) that holds if the matries F k are not singular. From the above, the following produt properties A S = S A C S = S C T C A = A C = { C S A = C A S = A C S A S C = S A C = S C A (64) are easily obtained. Speifially to our end, omitting the k index, setting s = sin Φ and = os Φ, and indiating with G BF and G k the matries used to map the B matrix onto the F matrix, we deal with the following nine ases: From Symmetri onto Symmetri G SS = S f S 1 rf 0 0 r f ri 0 0 r i 1 = r f r i I (65) From Skew-Symmetri onto Symmetri G AS = S f A 1 r 0 0 r 0 a a 0 1 = r a R (66) From Orthogonal onto Symmetri G US = S f C 1 r 0 0 r s s 1 = r s s (67) 14

15 From Symmetri onto Skew-Symmetri G SA = A f S 1 0 a a 0 r 0 0 r 1 = a r R (68) From Skew-Symmetri onto Skew-Symmetri G AA = A f A 1 0 ai a i 0 0 a f a f 0 1 = a i a f I (69) From Orthogonal onto Skew-Symmetri G UA = A f C 1 0 a a 0 s s 1 = a s s (70) From Symmetri onto Orthogonal G SU = C f S 1 s s r 0 0 r 1 = 1 r s s (71) From Skew-Symmetri onto Orthogonal G AU = C f A 1 s s 0 a a 0 1 = 1 a s s (7) From Orthogonal onto Orthogonal G UU = C f C 1 1 s 1 s 1 1 where 3 = os(φ 1 Φ ), and s 3 = sin(φ 1 Φ ) s s 1 = 3 s 3 s 3 3 (73) 15

16 Example The above equations, from Eq. (65) through Eq. (73), an be used, for instane, to keep the information of rotation (ontained in C), into Symmetri, Skew-Symmetri forms. An an example, let us to find the matrix mapping an Orthogonal matrix C onto Symmetri and Skew-Symmetri forms, whose eigenvalues are expressed as ( ) ( ) Φk tan Φk. S is obtained using Eq. (67), with r k = tan 4 4 ( ) Φk os Φ G k = G US = tan k sin Φ k (74) 4 sin Φ k os Φ k thus, the matrix m T US = P k G k Pk T (75) is suh that the matrix S = T US C (76) ( ) Φk is Symmetri and has eigenvalues ± tan. 4 Analogously, to obtain the Skew-Symmetri matrix, we should apply Eq. (70), with ( ) Φk a k = tan 4 ( ) Φk sin Φ G k = G UA = tan k os Φ k (77) 4 os Φ k sin Φ k thus, the matrix m T UA = P k G k Pk T (78) is suh that the matrix A = T UA C (79) ( ) Φk is Skew-Symmetri and has eigenvalues ±i tan. 4 Conlusion This paper shows that Cayley Transforms an be onsidered the extension to the matrix field of the omplex onformal mapping funtion f 1 (z). This results has been extended to similar matrix transforms whih map Orthogonal from/onto Symmetri, and Symmetri from/onto Skew-Symmetri matries. This has been ahieved Note that there is an expliit referene to the Modified Rodriguez Parameters (MRP), a minimum non-singular parameter set to desribe Orientation 16

17 using a set of real matries whih are simultaneously Orthogonal, and Symmetri (the Ortho Sym matries), and using two different omplex onformal mapping funtions (namely, the f (z), and the f 3 (z), that are here alled lokwise, and ounter lokwise funtions). Moreover, the lose relationship between some trigonometri identities and some matrix mappings is shown. This paper ends with a desription of a matrix mapping tehnique (always mapping Orthogonal, Symmetri, and Skew-Symmetri matries) that is obtained by diret matrix produt. Aknowledgment I would warmly like to thank Dr John L. Junkins and all the olleagues of the Aerospae Department of Texas A&M University for diretly and indiretly enouraging me to keep going with this researh. Referenes 1 J.L. Junkins and Y. Kim. Introdution to Dynamis and Control of Flexible Strutures. Amerian Institute of Aeronautis and Astronautis, AIAA Eduation Series, Washington, D.C., G.A. Korn and T.M. Korn. Mathematial Handbook for Sientists and Engineers. MGraw-Hill Book Company, 11 Avenue of the Amerias, New York, N.Y. 1000, Seond edition, D. Mortari. On the Rigid Rotation Conept in n-dimensional Spaes. The Journal of the Astronautial Sienes, 49(3):401 40, July-September A.K. Sanyal. Geometrial Transformations in Higher Dimensional Eulidean Spaes. Master s thesis, Department of Aerospae Engineering, Texas A&M University, College Station, TX, May

SURFACE WAVES OF NON-RAYLEIGH TYPE

SURFACE WAVES OF NON-RAYLEIGH TYPE SURFACE WAVES OF NON-RAYLEIGH TYPE by SERGEY V. KUZNETSOV Institute for Problems in Mehanis Prosp. Vernadskogo, 0, Mosow, 75 Russia e-mail: sv@kuznetsov.msk.ru Abstrat. Existene of surfae waves of non-rayleigh

More information

Control Theory association of mathematics and engineering

Control Theory association of mathematics and engineering Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology

More information

Quantum Mechanics: Wheeler: Physics 6210

Quantum Mechanics: Wheeler: Physics 6210 Quantum Mehanis: Wheeler: Physis 60 Problems some modified from Sakurai, hapter. W. S..: The Pauli matries, σ i, are a triple of matries, σ, σ i = σ, σ, σ 3 given by σ = σ = σ 3 = i i Let stand for the

More information

Advanced Computational Fluid Dynamics AA215A Lecture 4

Advanced Computational Fluid Dynamics AA215A Lecture 4 Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas

More information

Hankel Optimal Model Order Reduction 1

Hankel Optimal Model Order Reduction 1 Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both

More information

Aharonov-Bohm effect. Dan Solomon.

Aharonov-Bohm effect. Dan Solomon. Aharonov-Bohm effet. Dan Solomon. In the figure the magneti field is onfined to a solenoid of radius r 0 and is direted in the z- diretion, out of the paper. The solenoid is surrounded by a barrier that

More information

Study of EM waves in Periodic Structures (mathematical details)

Study of EM waves in Periodic Structures (mathematical details) Study of EM waves in Periodi Strutures (mathematial details) Massahusetts Institute of Tehnology 6.635 partial leture notes 1 Introdution: periodi media nomenlature 1. The spae domain is defined by a basis,(a

More information

1 sin 2 r = 1 n 2 sin 2 i

1 sin 2 r = 1 n 2 sin 2 i Physis 505 Fall 005 Homework Assignment #11 Solutions Textbook problems: Ch. 7: 7.3, 7.5, 7.8, 7.16 7.3 Two plane semi-infinite slabs of the same uniform, isotropi, nonpermeable, lossless dieletri with

More information

MODELING MATTER AT NANOSCALES. 4. Introduction to quantum treatments Eigenvectors and eigenvalues of a matrix

MODELING MATTER AT NANOSCALES. 4. Introduction to quantum treatments Eigenvectors and eigenvalues of a matrix MODELING MATTER AT NANOSCALES 4 Introdution to quantum treatments 403 Eigenvetors and eigenvalues of a matrix Simultaneous equations in the variational method The problem of simultaneous equations in the

More information

Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties

Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties Amerian Journal of Applied Sienes 4 (7): 496-501, 007 ISSN 1546-939 007 Siene Publiations Robust Flight ontrol Design for a urn oordination System with Parameter Unertainties 1 Ari Legowo and Hiroshi Okubo

More information

After the completion of this section the student should recall

After the completion of this section the student should recall Chapter I MTH FUNDMENTLS I. Sets, Numbers, Coordinates, Funtions ugust 30, 08 3 I. SETS, NUMERS, COORDINTES, FUNCTIONS Objetives: fter the ompletion of this setion the student should reall - the definition

More information

A 4 4 diagonal matrix Schrödinger equation from relativistic total energy with a 2 2 Lorentz invariant solution.

A 4 4 diagonal matrix Schrödinger equation from relativistic total energy with a 2 2 Lorentz invariant solution. A 4 4 diagonal matrix Shrödinger equation from relativisti total energy with a 2 2 Lorentz invariant solution. Han Geurdes 1 and Koji Nagata 2 1 Geurdes datasiene, 2593 NN, 164, Den Haag, Netherlands E-mail:

More information

COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION

COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION 4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and

More information

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge The Conept of Mass as Interfering Photons, and the Originating Mehanism of Gravitation D.T. Froedge V04 Formerly Auburn University Phys-dtfroedge@glasgow-ky.om Abstrat For most purposes in physis the onept

More information

The gravitational phenomena without the curved spacetime

The gravitational phenomena without the curved spacetime The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,

More information

Nonreversibility of Multiple Unicast Networks

Nonreversibility of Multiple Unicast Networks Nonreversibility of Multiple Uniast Networks Randall Dougherty and Kenneth Zeger September 27, 2005 Abstrat We prove that for any finite direted ayli network, there exists a orresponding multiple uniast

More information

A EUCLIDEAN ALTERNATIVE TO MINKOWSKI SPACETIME DIAGRAM.

A EUCLIDEAN ALTERNATIVE TO MINKOWSKI SPACETIME DIAGRAM. A EUCLIDEAN ALTERNATIVE TO MINKOWSKI SPACETIME DIAGRAM. S. Kanagaraj Eulidean Relativity s.kana.raj@gmail.om (1 August 009) Abstrat By re-interpreting the speial relativity (SR) postulates based on Eulidean

More information

MOLECULAR ORBITAL THEORY- PART I

MOLECULAR ORBITAL THEORY- PART I 5.6 Physial Chemistry Leture #24-25 MOLECULAR ORBITAL THEORY- PART I At this point, we have nearly ompleted our rash-ourse introdution to quantum mehanis and we re finally ready to deal with moleules.

More information

Ortho-Skew and Ortho-Sym Matrix Trigonometry 1

Ortho-Skew and Ortho-Sym Matrix Trigonometry 1 Abstract Ortho-Skew and Ortho-Sym Matrix Trigonometry Daniele Mortari Department of Aerospace Engineering, Texas A&M University, College Station, TX 77843-34 This paper introduces some properties of two

More information

Particle-wave symmetry in Quantum Mechanics And Special Relativity Theory

Particle-wave symmetry in Quantum Mechanics And Special Relativity Theory Partile-wave symmetry in Quantum Mehanis And Speial Relativity Theory Author one: XiaoLin Li,Chongqing,China,hidebrain@hotmail.om Corresponding author: XiaoLin Li, Chongqing,China,hidebrain@hotmail.om

More information

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution.

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution. arxiv:physis/99536v1 [physis.lass-ph] 15 May 1999 Eletromagneti radiation of the travelling spin wave propagating in an antiferromagneti plate. Exat solution. A.A.Zhmudsky November 19, 16 Abstrat The exat

More information

Critical Reflections on the Hafele and Keating Experiment

Critical Reflections on the Hafele and Keating Experiment Critial Refletions on the Hafele and Keating Experiment W.Nawrot In 1971 Hafele and Keating performed their famous experiment whih onfirmed the time dilation predited by SRT by use of marosopi loks. As

More information

). In accordance with the Lorentz transformations for the space-time coordinates of the same event, the space coordinates become

). In accordance with the Lorentz transformations for the space-time coordinates of the same event, the space coordinates become Relativity and quantum mehanis: Jorgensen 1 revisited 1. Introdution Bernhard Rothenstein, Politehnia University of Timisoara, Physis Department, Timisoara, Romania. brothenstein@gmail.om Abstrat. We first

More information

Complexity of Regularization RBF Networks

Complexity of Regularization RBF Networks Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw

More information

A Characterization of Wavelet Convergence in Sobolev Spaces

A Characterization of Wavelet Convergence in Sobolev Spaces A Charaterization of Wavelet Convergene in Sobolev Spaes Mark A. Kon 1 oston University Louise Arakelian Raphael Howard University Dediated to Prof. Robert Carroll on the oasion of his 70th birthday. Abstrat

More information

COMPARISON OF GEOMETRIC FIGURES

COMPARISON OF GEOMETRIC FIGURES COMPARISON OF GEOMETRIC FIGURES Spyros Glenis M.Ed University of Athens, Department of Mathematis, e-mail spyros_glenis@sh.gr Introdution the figures: In Eulid, the geometri equality is based on the apability

More information

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite. Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need

More information

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential

More information

n n=1 (air) n 1 sin 2 r =

n n=1 (air) n 1 sin 2 r = Physis 55 Fall 7 Homework Assignment #11 Solutions Textbook problems: Ch. 7: 7.3, 7.4, 7.6, 7.8 7.3 Two plane semi-infinite slabs of the same uniform, isotropi, nonpermeable, lossless dieletri with index

More information

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six

More information

Sensitivity Analysis in Markov Networks

Sensitivity Analysis in Markov Networks Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores

More information

Lecture 15 (Nov. 1, 2017)

Lecture 15 (Nov. 1, 2017) Leture 5 8.3 Quantum Theor I, Fall 07 74 Leture 5 (Nov., 07 5. Charged Partile in a Uniform Magneti Field Last time, we disussed the quantum mehanis of a harged partile moving in a uniform magneti field

More information

Berry s phase for coherent states of Landau levels

Berry s phase for coherent states of Landau levels Berry s phase for oherent states of Landau levels Wen-Long Yang 1 and Jing-Ling Chen 1, 1 Theoretial Physis Division, Chern Institute of Mathematis, Nankai University, Tianjin 300071, P.R.China Adiabati

More information

Math 220A - Fall 2002 Homework 8 Solutions

Math 220A - Fall 2002 Homework 8 Solutions Math A - Fall Homework 8 Solutions 1. Consider u tt u = x R 3, t > u(x, ) = φ(x) u t (x, ) = ψ(x). Suppose φ, ψ are supported in the annular region a < x < b. (a) Find the time T 1 > suh that u(x, t) is

More information

The Hanging Chain. John McCuan. January 19, 2006

The Hanging Chain. John McCuan. January 19, 2006 The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a

More information

Relative Maxima and Minima sections 4.3

Relative Maxima and Minima sections 4.3 Relative Maxima and Minima setions 4.3 Definition. By a ritial point of a funtion f we mean a point x 0 in the domain at whih either the derivative is zero or it does not exists. So, geometrially, one

More information

Extending LMR for anisotropic unconventional reservoirs

Extending LMR for anisotropic unconventional reservoirs Extending LMR for anisotropi unonventional reservoirs Maro A. Perez Apahe Canada Ltd Summary It has beome inreasingly advantageous to haraterize rok in unonventional reservoirs within an anisotropi framework.

More information

Scalable Positivity Preserving Model Reduction Using Linear Energy Functions

Scalable Positivity Preserving Model Reduction Using Linear Energy Functions Salable Positivity Preserving Model Redution Using Linear Energy Funtions Sootla, Aivar; Rantzer, Anders Published in: IEEE 51st Annual Conferene on Deision and Control (CDC), 2012 DOI: 10.1109/CDC.2012.6427032

More information

arxiv: v2 [cs.dm] 4 May 2018

arxiv: v2 [cs.dm] 4 May 2018 Disrete Morse theory for the ollapsibility of supremum setions Balthazar Bauer INRIA, DIENS, PSL researh, CNRS, Paris, Frane Luas Isenmann LIRMM, Université de Montpellier, CNRS, Montpellier, Frane arxiv:1803.09577v2

More information

Singular Event Detection

Singular Event Detection Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate

More information

On solution of Klein-Gordon equation in scalar and vector potentials

On solution of Klein-Gordon equation in scalar and vector potentials On solution of Klein-Gordon equation in salar and vetor potentials Liu Changshi Physis division, Department of mehanial and eletrial engineering, Jiaxing College, Zhejiang, 314001, P. R. China Lius4976@sohu.om

More information

Integration of the Finite Toda Lattice with Complex-Valued Initial Data

Integration of the Finite Toda Lattice with Complex-Valued Initial Data Integration of the Finite Toda Lattie with Complex-Valued Initial Data Aydin Huseynov* and Gusein Sh Guseinov** *Institute of Mathematis and Mehanis, Azerbaijan National Aademy of Sienes, AZ4 Baku, Azerbaijan

More information

Bäcklund Transformations: Some Old and New Perspectives

Bäcklund Transformations: Some Old and New Perspectives Bäklund Transformations: Some Old and New Perspetives C. J. Papahristou *, A. N. Magoulas ** * Department of Physial Sienes, Helleni Naval Aademy, Piraeus 18539, Greee E-mail: papahristou@snd.edu.gr **

More information

Chapter 9. The excitation process

Chapter 9. The excitation process Chapter 9 The exitation proess qualitative explanation of the formation of negative ion states Ne and He in He-Ne ollisions an be given by using a state orrelation diagram. state orrelation diagram is

More information

Array Design for Superresolution Direction-Finding Algorithms

Array Design for Superresolution Direction-Finding Algorithms Array Design for Superresolution Diretion-Finding Algorithms Naushad Hussein Dowlut BEng, ACGI, AMIEE Athanassios Manikas PhD, DIC, AMIEE, MIEEE Department of Eletrial Eletroni Engineering Imperial College

More information

Directional Coupler. 4-port Network

Directional Coupler. 4-port Network Diretional Coupler 4-port Network 3 4 A diretional oupler is a 4-port network exhibiting: All ports mathed on the referene load (i.e. S =S =S 33 =S 44 =0) Two pair of ports unoupled (i.e. the orresponding

More information

Beams on Elastic Foundation

Beams on Elastic Foundation Professor Terje Haukaas University of British Columbia, Vanouver www.inrisk.ub.a Beams on Elasti Foundation Beams on elasti foundation, suh as that in Figure 1, appear in building foundations, floating

More information

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry 9 Geophysis and Radio-Astronomy: VLBI VeryLongBaseInterferometry VLBI is an interferometry tehnique used in radio astronomy, in whih two or more signals, oming from the same astronomial objet, are reeived

More information

4. (12) Write out an equation for Poynting s theorem in differential form. Explain in words what each term means physically.

4. (12) Write out an equation for Poynting s theorem in differential form. Explain in words what each term means physically. Eletrodynamis I Exam 3 - Part A - Closed Book KSU 205/2/8 Name Eletrodynami Sore = 24 / 24 points Instrutions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try to

More information

Ayan Kumar Bandyopadhyay

Ayan Kumar Bandyopadhyay Charaterization of radiating apertures using Multiple Multipole Method And Modeling and Optimization of a Spiral Antenna for Ground Penetrating Radar Appliations Ayan Kumar Bandyopadhyay FET-IESK, Otto-von-Guerike-University,

More information

Application of the Dyson-type boson mapping for low-lying electron excited states in molecules

Application of the Dyson-type boson mapping for low-lying electron excited states in molecules Prog. Theor. Exp. Phys. 05, 063I0 ( pages DOI: 0.093/ptep/ptv068 Appliation of the Dyson-type boson mapping for low-lying eletron exited states in moleules adao Ohkido, and Makoto Takahashi Teaher-training

More information

Discrete Bessel functions and partial difference equations

Discrete Bessel functions and partial difference equations Disrete Bessel funtions and partial differene equations Antonín Slavík Charles University, Faulty of Mathematis and Physis, Sokolovská 83, 186 75 Praha 8, Czeh Republi E-mail: slavik@karlin.mff.uni.z Abstrat

More information

QUANTUM MECHANICS II PHYS 517. Solutions to Problem Set # 1

QUANTUM MECHANICS II PHYS 517. Solutions to Problem Set # 1 QUANTUM MECHANICS II PHYS 57 Solutions to Problem Set #. The hamiltonian for a lassial harmoni osillator an be written in many different forms, suh as use ω = k/m H = p m + kx H = P + Q hω a. Find a anonial

More information

the following action R of T on T n+1 : for each θ T, R θ : T n+1 T n+1 is defined by stated, we assume that all the curves in this paper are defined

the following action R of T on T n+1 : for each θ T, R θ : T n+1 T n+1 is defined by stated, we assume that all the curves in this paper are defined How should a snake turn on ie: A ase study of the asymptoti isoholonomi problem Jianghai Hu, Slobodan N. Simić, and Shankar Sastry Department of Eletrial Engineering and Computer Sienes University of California

More information

University of Wollongong Department of Economics Working Paper Series 2000

University of Wollongong Department of Economics Working Paper Series 2000 University of Wollongong Department of Eonomis Working Paper Series 000 Rational Non-additive Eating: Cyles, Overweightness, and Underweightness Amnon Levy WP 00-07 RATIONAL NON-ADDICTIVE EATING: CYCLES,

More information

Aircraft CAS Design with Input Saturation Using Dynamic Model Inversion

Aircraft CAS Design with Input Saturation Using Dynamic Model Inversion International Journal of Control, Automation, and Systems Vol., No. 3, September 003 35 Airraft CAS Design with Input Saturation sing Dynami Model Inversion Sangsoo Lim and Byoung Soo Kim Abstrat: This

More information

RESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED

RESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED The Seventh Asia-Paifi Conferene on Wind Engineering, November 8-1, 9, Taipei, Taiwan RESEARCH ON RANDOM FORIER WAVE-NMBER SPECTRM OF FLCTATING WIND SPEED Qi Yan 1, Jie Li 1 Ph D. andidate, Department

More information

MAC Calculus II Summer All you need to know on partial fractions and more

MAC Calculus II Summer All you need to know on partial fractions and more MC -75-Calulus II Summer 00 ll you need to know on partial frations and more What are partial frations? following forms:.... where, α are onstants. Partial frations are frations of one of the + α, ( +

More information

Evaluation of effect of blade internal modes on sensitivity of Advanced LIGO

Evaluation of effect of blade internal modes on sensitivity of Advanced LIGO Evaluation of effet of blade internal modes on sensitivity of Advaned LIGO T0074-00-R Norna A Robertson 5 th Otober 00. Introdution The urrent model used to estimate the isolation ahieved by the quadruple

More information

A NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM

A NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM NETWORK SIMPLEX LGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM Cen Çalışan, Utah Valley University, 800 W. University Parway, Orem, UT 84058, 801-863-6487, en.alisan@uvu.edu BSTRCT The minimum

More information

arxiv:physics/ v4 [physics.gen-ph] 9 Oct 2006

arxiv:physics/ v4 [physics.gen-ph] 9 Oct 2006 The simplest derivation of the Lorentz transformation J.-M. Lévy Laboratoire de Physique Nuléaire et de Hautes Energies, CNRS - IN2P3 - Universités Paris VI et Paris VII, Paris. Email: jmlevy@in2p3.fr

More information

Packing Plane Spanning Trees into a Point Set

Packing Plane Spanning Trees into a Point Set Paking Plane Spanning Trees into a Point Set Ahmad Biniaz Alfredo Garía Abstrat Let P be a set of n points in the plane in general position. We show that at least n/3 plane spanning trees an be paked into

More information

A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE. W. J.

A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE. W. J. A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE W. J. Haering* Senior Projet Engineer General Motors Corporation Warren, Mihigan R.

More information

On Component Order Edge Reliability and the Existence of Uniformly Most Reliable Unicycles

On Component Order Edge Reliability and the Existence of Uniformly Most Reliable Unicycles Daniel Gross, Lakshmi Iswara, L. William Kazmierzak, Kristi Luttrell, John T. Saoman, Charles Suffel On Component Order Edge Reliability and the Existene of Uniformly Most Reliable Uniyles DANIEL GROSS

More information

A Recursive Approach to the Kauffman Bracket

A Recursive Approach to the Kauffman Bracket Applied Mathematis, 204, 5, 2746-2755 Published Online Otober 204 in SiRes http://wwwsirporg/journal/am http://ddoiorg/04236/am20457262 A Reursive Approah to the Kauffman Braet Abdul Rauf Nizami, Mobeen

More information

Sufficient Conditions for a Flexible Manufacturing System to be Deadlocked

Sufficient Conditions for a Flexible Manufacturing System to be Deadlocked Paper 0, INT 0 Suffiient Conditions for a Flexile Manufaturing System to e Deadloked Paul E Deering, PhD Department of Engineering Tehnology and Management Ohio University deering@ohioedu Astrat In reent

More information

CMSC 451: Lecture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017

CMSC 451: Lecture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017 CMSC 451: Leture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017 Reading: Chapt 11 of KT and Set 54 of DPV Set Cover: An important lass of optimization problems involves overing a ertain domain,

More information

LECTURE 2 Geometrical Properties of Rod Cross Sections (Part 2) 1 Moments of Inertia Transformation with Parallel Transfer of Axes.

LECTURE 2 Geometrical Properties of Rod Cross Sections (Part 2) 1 Moments of Inertia Transformation with Parallel Transfer of Axes. V. DEMENKO MECHNCS OF MTERLS 05 LECTURE Geometrial Properties of Rod Cross Setions (Part ) Moments of nertia Transformation with Parallel Transfer of xes. Parallel-xes Theorems S Given: a b = S = 0. z

More information

Examples of Tensors. February 3, 2013

Examples of Tensors. February 3, 2013 Examples of Tensors February 3, 2013 We will develop a number of tensors as we progress, but there are a few that we an desribe immediately. We look at two ases: (1) the spaetime tensor desription of eletromagnetism,

More information

A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS

A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Vietnam Journal of Mehanis, VAST, Vol. 4, No. (), pp. A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Le Thanh Tung Hanoi University of Siene and Tehnology, Vietnam Abstrat. Conventional ship autopilots are

More information

EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION

EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION Journal of Mathematial Sienes: Advanes and Appliations Volume 3, 05, Pages -3 EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION JIAN YANG, XIAOJUAN LU and SHENGQIANG TANG

More information

The Effectiveness of the Linear Hull Effect

The Effectiveness of the Linear Hull Effect The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports

More information

Brazilian Journal of Physics, vol. 29, no. 3, September, Classical and Quantum Mechanics of a Charged Particle

Brazilian Journal of Physics, vol. 29, no. 3, September, Classical and Quantum Mechanics of a Charged Particle Brazilian Journal of Physis, vol. 9, no. 3, September, 1999 51 Classial and Quantum Mehanis of a Charged Partile in Osillating Eletri and Magneti Fields V.L.B. de Jesus, A.P. Guimar~aes, and I.S. Oliveira

More information

Velocity Addition in Space/Time David Barwacz 4/23/

Velocity Addition in Space/Time David Barwacz 4/23/ Veloity Addition in Spae/Time 003 David arwaz 4/3/003 daveb@triton.net http://members.triton.net/daveb Abstrat Using the spae/time geometry developed in the previous paper ( Non-orthogonal Spae- Time geometry,

More information

arxiv:cond-mat/ v1 [cond-mat.str-el] 3 Aug 2006

arxiv:cond-mat/ v1 [cond-mat.str-el] 3 Aug 2006 arxiv:ond-mat/0608083v1 [ond-mat.str-el] 3 Aug 006 Raman sattering for triangular latties spin- 1 Heisenberg antiferromagnets 1. Introdution F. Vernay 1,, T. P. Devereaux 1, and M. J. P. Gingras 1,3 1

More information

Armenian Theory of Special Relativity (Illustrated) Robert Nazaryan 1 and Haik Nazaryan 2

Armenian Theory of Special Relativity (Illustrated) Robert Nazaryan 1 and Haik Nazaryan 2 29606 Robert Nazaryan Haik Nazaryan/ Elixir Nulear & Radiation Phys. 78 (205) 29606-2967 Available online at www.elixirpublishers.om (Elixir International Journal) Nulear Radiation Physis Elixir Nulear

More information

The Unified Geometrical Theory of Fields and Particles

The Unified Geometrical Theory of Fields and Particles Applied Mathematis, 014, 5, 347-351 Published Online February 014 (http://www.sirp.org/journal/am) http://dx.doi.org/10.436/am.014.53036 The Unified Geometrial Theory of Fields and Partiles Amagh Nduka

More information

LECTURE NOTES FOR , FALL 2004

LECTURE NOTES FOR , FALL 2004 LECTURE NOTES FOR 18.155, FALL 2004 83 12. Cone support and wavefront set In disussing the singular support of a tempered distibution above, notie that singsupp(u) = only implies that u C (R n ), not as

More information

University of Groningen

University of Groningen University of Groningen Port Hamiltonian Formulation of Infinite Dimensional Systems II. Boundary Control by Interonnetion Mahelli, Alessandro; van der Shaft, Abraham; Melhiorri, Claudio Published in:

More information

Non-Abelian quantum Hall states and their quasiparticles: From the pattern of zeros to vertex algebra

Non-Abelian quantum Hall states and their quasiparticles: From the pattern of zeros to vertex algebra Non-Abelian quantum Hall states and their quasipartiles: From the pattern of zeros to vertex algebra Yuan-Ming Lu, 1 Xiao-Gang Wen,, Zhenghan Wang, 4 and Ziqiang Wang 1 1 Department of Physis, Boston College,

More information

The Corpuscular Structure of Matter, the Interaction of Material Particles, and Quantum Phenomena as a Consequence of Selfvariations.

The Corpuscular Structure of Matter, the Interaction of Material Particles, and Quantum Phenomena as a Consequence of Selfvariations. The Corpusular Struture of Matter, the Interation of Material Partiles, and Quantum Phenomena as a Consequene of Selfvariations. Emmanuil Manousos APM Institute for the Advanement of Physis and Mathematis,

More information

The Electromagnetic Radiation and Gravity

The Electromagnetic Radiation and Gravity International Journal of Theoretial and Mathematial Physis 016, 6(3): 93-98 DOI: 10.593/j.ijtmp.0160603.01 The Eletromagneti Radiation and Gravity Bratianu Daniel Str. Teiului Nr. 16, Ploiesti, Romania

More information

Bäcklund Transformations and Explicit Solutions of (2+1)-Dimensional Barotropic and Quasi-Geostrophic Potential Vorticity Equation

Bäcklund Transformations and Explicit Solutions of (2+1)-Dimensional Barotropic and Quasi-Geostrophic Potential Vorticity Equation Commun. Theor. Phys. Beijing, China) 53 010) pp. 803 808 Chinese Physial Soiety and IOP Publishing Ltd Vol. 53, No. 5, May 15, 010 Bäklund Transformations and Expliit Solutions of +1)-Dimensional Barotropi

More information

Model-based mixture discriminant analysis an experimental study

Model-based mixture discriminant analysis an experimental study Model-based mixture disriminant analysis an experimental study Zohar Halbe and Mayer Aladjem Department of Eletrial and Computer Engineering, Ben-Gurion University of the Negev P.O.Box 653, Beer-Sheva,

More information

arxiv:cond-mat/ v1 [cond-mat.stat-mech] 16 Aug 2004

arxiv:cond-mat/ v1 [cond-mat.stat-mech] 16 Aug 2004 Computational omplexity and fundamental limitations to fermioni quantum Monte Carlo simulations arxiv:ond-mat/0408370v1 [ond-mat.stat-meh] 16 Aug 2004 Matthias Troyer, 1 Uwe-Jens Wiese 2 1 Theoretishe

More information

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field Four-dimensional equation of motion for visous ompressible substane with regard to the aeleration field, pressure field and dissipation field Sergey G. Fedosin PO box 6488, Sviazeva str. -79, Perm, Russia

More information

arxiv:math/ v1 [math.ca] 27 Nov 2003

arxiv:math/ v1 [math.ca] 27 Nov 2003 arxiv:math/011510v1 [math.ca] 27 Nov 200 Counting Integral Lamé Equations by Means of Dessins d Enfants Sander Dahmen November 27, 200 Abstrat We obtain an expliit formula for the number of Lamé equations

More information

General Equilibrium. What happens to cause a reaction to come to equilibrium?

General Equilibrium. What happens to cause a reaction to come to equilibrium? General Equilibrium Chemial Equilibrium Most hemial reations that are enountered are reversible. In other words, they go fairly easily in either the forward or reverse diretions. The thing to remember

More information

The Laws of Acceleration

The Laws of Acceleration The Laws of Aeleration The Relationships between Time, Veloity, and Rate of Aeleration Copyright 2001 Joseph A. Rybzyk Abstrat Presented is a theory in fundamental theoretial physis that establishes the

More information

A Queueing Model for Call Blending in Call Centers

A Queueing Model for Call Blending in Call Centers A Queueing Model for Call Blending in Call Centers Sandjai Bhulai and Ger Koole Vrije Universiteit Amsterdam Faulty of Sienes De Boelelaan 1081a 1081 HV Amsterdam The Netherlands E-mail: {sbhulai, koole}@s.vu.nl

More information

arxiv: v1 [cs.cc] 13 Feb 2008

arxiv: v1 [cs.cc] 13 Feb 2008 A review of the Statistial Mehanis approah to Random Optimization Problems Fabrizio Altarelli 1,2, Rémi Monasson 2, Guilhem Semerjian 2 and Franeso Zamponi 2 1 Dipartimento di Fisia and CNR, Università

More information

Verka Prolović Chair of Civil Engineering Geotechnics, Faculty of Civil Engineering and Architecture, Niš, R. Serbia

Verka Prolović Chair of Civil Engineering Geotechnics, Faculty of Civil Engineering and Architecture, Niš, R. Serbia 3 r d International Conferene on New Developments in Soil Mehanis and Geotehnial Engineering, 8-30 June 01, Near East University, Niosia, North Cyprus Values of of partial fators for for EC EC 7 7 slope

More information

A FETI Method for a Class of Indefinite or Complex Second- or Fourth-Order Problems

A FETI Method for a Class of Indefinite or Complex Second- or Fourth-Order Problems A FETI Method for a Class of Indefinite or Complex Seond- or Fourth-Order Problems Charbel Farhat 1, Jing Li 2, Mihel Lesoinne 1 and Philippe Avery 1 1 University of Colorado at Boulder, Department of

More information

A Functional Representation of Fuzzy Preferences

A Functional Representation of Fuzzy Preferences Theoretial Eonomis Letters, 017, 7, 13- http://wwwsirporg/journal/tel ISSN Online: 16-086 ISSN Print: 16-078 A Funtional Representation of Fuzzy Preferenes Susheng Wang Department of Eonomis, Hong Kong

More information

The simulation analysis of the bridge rectifier continuous operation in AC circuit

The simulation analysis of the bridge rectifier continuous operation in AC circuit Computer Appliations in Eletrial Engineering Vol. 4 6 DOI 8/j.8-448.6. The simulation analysis of the bridge retifier ontinuous operation in AC iruit Mirosław Wiślik, Paweł Strząbała Kiele University of

More information

Remarks Around Lorentz Transformation

Remarks Around Lorentz Transformation Remarks Around Lorentz Transformation Arm Boris Nima arm.boris@gmail.om Abstrat After diagonalizing the Lorentz Matrix, we find the frame where the Dira equation is one derivation and we alulate the speed

More information

Stability Analysis of Orbital Motions around Uniformly Rotating Irregular Asteroids

Stability Analysis of Orbital Motions around Uniformly Rotating Irregular Asteroids Stability Analysis of Orbital Motions around Uniformly Rotating Irregular Asteroids By Xiyun HoU, ) Daniel J. SCHEERES, ) Xiaosheng XIN, ) Jinglang FENG, ) Jingshi TANG, ) Lin LIU, ) ) Shool of Astronomy

More information

Estimating the probability law of the codelength as a function of the approximation error in image compression

Estimating the probability law of the codelength as a function of the approximation error in image compression Estimating the probability law of the odelength as a funtion of the approximation error in image ompression François Malgouyres Marh 7, 2007 Abstrat After some reolletions on ompression of images using

More information

Wuming Li and Fan Yang

Wuming Li and Fan Yang NEW ZEALAND JOURNAL OF MATHEMATICS Volume 33 (004), 159 164 N DIMENSIONAL MINKOWSKI SPACE AND SPACE TIME ALGEBRA Wuming Li and Fan Yang (Reeived February 003) Abstrat. By using an n Dimensional Minkowski

More information

Nuclear Shell Structure Evolution Theory

Nuclear Shell Structure Evolution Theory Nulear Shell Struture Evolution Theory Zhengda Wang (1) Xiaobin Wang () Xiaodong Zhang () Xiaohun Wang () (1) Institute of Modern physis Chinese Aademy of SienesLan Zhou P. R. China 70000 () Seagate Tehnology

More information