Algebraic Attacks on Summation Generators

Size: px
Start display at page:

Download "Algebraic Attacks on Summation Generators"

Transcription

1 Algebraic Aacks on Summaion Generaors Dong Hoon Lee, Jaeheon Kim, Jin Hong, Jae Woo Han, and Dukjae Moon Naional Securiy Research Insiue 161 Gajeong-dong, Yuseong-gu, Daejeon, , Korea Absrac. We apply he algebraic aacks on sream ciphers wih memories o he summaion generaor. For a summaion generaor ha uses n LFSRs, an algebraic equaion relaing he key sream bis and LFSR oupu bis can be made o be of degree less han or equal o 2 log 2 n, using log 2 n + 1 consecuive key sream bis. This is much lower han he upper bound given by previous general resuls. We also show ha he echniques of [5] can be applied o summaion generaors using 2 k LFSRs o reduce he effecive degree of he algebraic equaion. Keywords: sream ciphers, algebraic aacks, summaion generaors 1 Inroducion Among recen developmens on sream ciphers, he algebraic aack has gahered much aenion. In his aack, an algebraic equaion relaing he iniial key bis and he oupu key sream bis is se up and solved hrough linearizaion echniques. Algebraic aack was firs applied o block ciphers and public key cryposysems [7, 2]. And is firs successful applicaion o sream cipher was done on Toyocryp [3]. As he mehod was soon exended o LILI-128 [6], i gahered much aenion. Sream ciphers ha uilize memory were firs hough o be much more resisan o hese aacks, bu soon i was shown ha even hese cases were subjec o algebraic aacks [1, 4]. The summaion generaor proposed by Ruepel [10] is a nonlinear combiner wih memory. I is known ha he generaor produces sequences whose period and correlaion immuniy are maximum, and whose linear complexiy is conjecured o be close o he period. Hence i serves as a good building block for sream ciphers. However, a correlaion aack on he summaion generaor ha uses wo LFSRs was presened in [9], even hough i is also saed in [9] ha his aack is no plausible if here are more han wo LFSRs in use. Anoher well known aack on summaion generaor is given by [8] and poins in he opposie direcion. I uses FCSRs o simulae he summaion generaor and indicaes ha for a fixed iniial key size, breaking hem ino oo many LFSRs will add o is weakness. In his work, we sudy he summaion generaor wih he algebraic aack in mind. We show ha for a summaion generaor ha uses n LFSRs, an algebraic equaion relaing he key sream bis and LFSR oupu bis can be made o be of degree less han or equal o 2 log 2 n, using log 2 n + 1 consecuive key

2 sream bis. This is much lower han he upper bound on he degree of algebraic equaions ha is guaraneed by he general works [1, 4]. We also show ha he echniques of [5] can be applied o summaion generaors using 2 k LFSRs o reduce he effecive degree of he equaion furher. The impaien reader may confer o Tables 3 and 4, appearing in laer secions, for a quick look on he improvemens we have given o undersanding he acual srengh of algebraic aacks on summaion generaors. Acually, Table 3 should be a good reference in view of algebraic aacks for anyone considering he use of a summaion generaor. 2 Summaion generaors and symmeric polynomials We consider a summaion generaor ha uses n binary LFSRs. The oupu of he j-h LFSR a ime will be denoed by x j. Since we are dealing wih binary values, we will no be using powers of hese erms, hence he superscrip should no cause any confusion. The curren carry value from he previous sage is denoed by c. Noe ha he carry can be expressed in k = log 2 n bis. We shall le c = (c k 1,..., c 1, c 0) be he binary expression of he carry value. The binary oupu z and he carry value for he nex sage from he summaion generaor is given by z = x 1 x 2 x n c 0, (1) = (x 1 + x x n + c )/2. (2) Le us denoe by Si, he i-h elemenary symmeric polynomial in he variables {x 1,..., x n}. We view hem as boolean funcions raher han as polynomials. Explicily, hey are S 0 = 1, S 1 = n j=1x j, S 2 = 1 j1 <j 2 nx j 1 x j 2,. S n = x 1x 2 x n. We shall also call hese by he name elemenary symmeric boolean funcions. For each 0 a, b n, le m a,b denoe he value of Sa under he condiion j x j = b. This condiion saes ha b of he n variable x j are equal o 1 and ha he ohers are equal o 0. Since Sa is symmeric, i makes sense o evaluae Sa under his condiion. I is also clear ha he values (m a,b ) n b=0 compleely deermines he value of Sa a an arbirary inpu. To calculae m a,b, one has only o coun how many of he monomials conained in Sa is nonzero. Hence we have ( ) b m a,b (mod 2). (3) a

3 Now, consider he (n + 1) (n + 1) marix M = (m a,b ). The n = 4 case is given in Table 1, as an example. Denoe by M, he n n marix obained by x j S S S S S Table 1. The marix M for n = 4 removing he las row and column from M. When we wan o make explici he number of variable used in defining M and M, we shall wrie M(n) and M (n). Lemma 1. For powers of 2, he marix M saisfies ( M (2 k+1 M ) = (2 k ) M (2 k ) ) 0 M (2 k. ) Proof. Le us firs divide M ino four pars and wrie ( ) M II I =. III IV I is clear from (3) ha he second quadran, par II, is M (2 k ), as claimed. The equaion also shows ha he lower riangular par of M is zero. Hence he hird quadran is filled wih zero, as claimed. We now show ha I, II, and IV are idenical. Once more, referring o (3), i suffices o show ( ) b a ( 2 k ) + b a ( 2 k ) + b 2 k + a for all 0 a, b < 2 k. To his end, we may easily check ha (1 + x) 2k +b = (1 + x) 2k (1 + x) b (1 + x 2k )(1 + x) b (mod 2) = (1 + x) b + x 2k (1 + x) b. (mod 2) (4) The coefficien of x a ha would appear in he expansion of he lef hand side (1 + x) 2k +b is he middle erm of (4). Since a < 2 k, he x a erm in he righ hand side may appear only in (1 + x) b and is equal o he firs erm of (4). This shows he firs equaliy. Similarly, comparison of he coefficiens of x 2k +a shows he equaliy beween he firs and he las erms of (4). This complees he proof.

4 This lemma allows one o wrie he marix M (2 k ) explicily for any given k. Then, owing o (3), he marix M(n) for any n may be obained as a submarix of a big enough M (2 k ). For example, Table 1 is a submarix of M (8). The basic heory on symmeric polynomials ells us ha he se consising of producs of elemenary symmeric polynomials forms a basis for he space of symmeric polynomials. In case of symmeric boolean funcions, he following lemma can easily be seen o be rue. Lemma 2. Any boolean funcion ha is symmeric in is variables, may be wrien as a linear combinaion of he elemenary symmeric boolean funcions. 3 Summaion generaor on 4 LFSRs We fix n = 4 hroughou his secion and derive an algebraic equaion of degree 4 ha is saisfied by he summaion generaor. The variables of he final equaion will consis of he LFSR oupu bis and he key sream bis bu will no conain any carry bis. Since we are dealing wih 4 = 2 2 LFSRs, we have k = 2, and i suffices o use 2 bis in expressing he carry value. I is clear ha he carry value should be symmeric wih respec o he order of he LFSRs in use. Recalling Lemma 2, his implies ha, when he curren carry value c is fixed, he carry bis 0 and 1 may be expressed as linear combinaions of he elemenary symmeric boolean funcions Si. In he general case when he curren carry bis are no fixed, hey will be linear combinaions of he Si wih boolean funcions of he carry bis c 0 and c 1 used as coefficiens. For each value of c and LFSR inpus possible, we explicily calculaed 0 and 1. We hen used Table 1 o expressed hem using he elemenary symmeric boolean funcions. The resul is given in Table 2. c S2 S1 S2 S0 S2 S0 S1 S2 1 S4 S3 S4 S2 S4 S1 S2 S3 S4 Table 2. The nex carry bis in relaion o he curren carry value Lemma 3. For a summaion generaor on 4 LFSRs, he following expresses he nex sage carry bis as funcions of he curren LFSR oupu bis and curren carry bis. Proof. Using Table 2, we may wrie 0 = S2 c 0S1 c 1 (5) 1 = S4 c 0S3 c 1S2 c 0c 1S1. (6) 0 = {(1 c 0)(1 c 1)S 2} {c 0(1 c 1)(S 1 S 2)} {(1 c 0)c 1(S 0 S 2)} {c 0c 1(S 0 S 1 S 2)}

5 and 1 = {(1 c 0)(1 c 1)S 4} {c 0(1 c 1)(S 3 S 4)} {(1 c 0)c 1(S 2 S 4)} {c 0c 1(S 1 S 2 S 3 S 4)}. Simplificaion of hese equaions gives he claimed saemens. Lemma 4. For n = 4, he following expresses he carry bis of he summaion generaor as polynomials in he LFSR oupu bis and he key sream bis. c 0 = S 1 z. (7) c 1 = S 2 (1 z )S 1 S +1 1 z +1. (8) The firs carry bi c 0 is linear and he second carry bi c 1 LFSR oupu bis. is of degree 2 in he Proof. The firs equaion follows immediaely from (1). I is a linear funcion on he variables x j. Subsiuing his and is shif c+1 0 ino (5) gives c 1 = S 2 (S 1 z )S 1 S +1 1 z +1. Now, since we are dealing wih boolean funcions, we have (S1 )2 = S1 second equaion follows. and he Finally, wih his lemma, we may remove all occurrence of he carry bis from (6) o obain he following proposiion. Proposiion 1. The following algebraic equaion holds rue for a summaion generaor on 4 LFSRs. 0 = S4 (1 z )S3 S2S 1 +1 (1 z +1 )S2 S2 +1 (1 z )S1S 1 +1 (1 z )(1 z +1 )S1 (1 z +1 )S1 +1 S1 +2 z +2. I is of degree 4 in he oupu bis of he LFSRs and uses 3 consecuive key sream bis. While simplifying he subsiuion of (7), (8), and shif of (8) ino (6), we have used he equaliy S 1S 3 = S 3 of boolean funcions. 4 Summaion generaor on n = 2 k LFSRs Le us denoe by Fi +1 (n), he (symmeric) polynomial ha expresses he nex sage carry bi i in erms of LFSR oupus x j and curren carry bis c j. Since we shall be dealing wih polynomials on differen number of variables, we

6 shall wrie Sj (n) o denoe he elemenary symmeric boolean funcion on n variables. As an example, we saw in he previous secion ha by F0 +1 (2 2 ) = S2(2 2 ) c 0S1(2 2 ) c 1, (9) F1 +1 (2 2 ) = S4(2 2 ) c 0S3(2 2 ) c 1S2(2 2 ) c 0c 1S1(2 2 ). (10) Le us suppose ha for some n = 2 k, he polynomials Fi +1 (2 k ) are given F +1 i (2 k ) = j f i,j S j(2 k ). (11) Here, each coefficien f i,j is a boolean funcion defined on he curren carry bis c 0, c 1,..., c k 1. For example, we see from (10) ha f 1,4 = 1, f 1,3 = c 0, f 1,2 = c 1, f 1,1 = c 0c 1, f 1,0 = 0, for k = 2. The following proposiion will allow us o inducively calculae all F +1 i (2 k ) for any i and k. Proposiion 2. Suppose equaion (11) holds for some k. Then we have F +1 i (2 k+1 ) = j f i,j S j(2 k+1 ), for i < k 1, (12) F +1 k 1 (2k+1 ) = ( j f k 1,j S j(2 k+1 ) ) c k, (13) F +1 k (2 k+1 ) = ( j f k 1,j S j+2 k (2 k+1 ) ) c k( j f k 1,j S j(2 k+1 ) ). (14) The proof of his proposiion may be found in he appendices. As an immediae applicaion of his proposiion o (9) and (10), we may wrie F0 +1 (2 3 ) = S2(2 3 ) c 0S1(2 3 ) c 1, (15) F1 +1 (2 3 ) = S4(2 3 ) c 0S3(2 3 ) c 1S2(2 3 ) c 0c 1S1(2 3 ) c 2, (16) F2 +1 (2 3 ) = S8(2 3 ) c 0S7(2 3 ) c 1S6(2 3 ) c 0c 1S5(2 3 ) c 2S4(2 3 ) c 0c 2S3(2 3 ) c 1c 2S2(2 3 ) c 0c 1c 2S1(2 3 ). (17) Now, le us briefly recall he process we wen hrough in Secion 3 in obaining he degree 4 equaion of Proposiion 1. We sared ou wih hree equaions. z = S 1 c 0. (18) 0 = S2 c 0S1 c 1. (19) 1 = S4 c 0S3 c 1(S2 c 0S1). (20) We used (18) o wrie c 0 as a degree 1 equaion ha involves jus he key sream bi and he LFSR oupus. This was subsiued in he nex equaion (19) o wrie c 1 as a degree 2 equaion of he same kind. Finally, hese expressions for c 0 and c 1 were subsiued in he las equaion o obain he degree 4 equaion ha involves only key sream bis and LFSR oupu bis.

7 Wha would happen if we waned o do he same for n = 2 3. We would sar wih he following se of equaions. z = S 1 c 0. (21) 0 = S2 c 0S1 c 1. (22) 1 = S4 c 0S3 c 1S2 c 0c 1S1 c 2. (23) 2 = S8 c 0S7 c 1S6 c 0c 1S5 c 2(S4 c 0S3 c 1S2 c 0c 1S1). (24) Noice ha he firs wo equaions here are idenical o (18) and (19), as saed by (12) of Proposiion 2. Hence, as before, c 0 and c 1 will be wrien as degree 1 and 2 polynomials. The main par of (23) is idenical o (20), as saed by (13). They only differ in ha c 2 appears a he end of (23). Since (20) is of degree 4, equaion (23) gives a degree 4 expression for c 2. Now, as given by (14), he righ hand side of (24) may be broken ino wo big erms of degree (less han or equal o) 8. The firs erm is a degree 4 equaion shifed by degree 4 and he second erm is a produc of wo degree 4 equaions. Finally, he lef hand side of (24) is of degree 4. Hence, subsiuion of c 0, c 1 and c 2 ino (24) gives a degree 8 polynomial connecing various LFSR oupu bis and key sream bis. One can easily see ha he above argumen is general enough o be seen as he inducion sep needed in proving he following heorem. Theorem 1. Consider a summaion generaor on n = 2 k LFSRs. There exiss an algebraic equaion connecing LFSR oupu bis and k + 1 consecuive key sream bis in such a way ha i is of degree 2 k in he LFSR oupu bis. 5 The general case The following is an easy corollary o Theorem 1. Theorem 2. Consider a summaion generaor of n LFSRs. We shall le k = log 2 n. There exiss an algebraic equaion connecing LFSR oupu bis and k + 1 consecuive key sream bis in such a way ha i is of degree less han or equal o 2 k in he LFSR oupu bis. Proof. We may model a summaion generaor on n LFSRs as a summaion generaor on 2 k LFSRs wih (2 k n) of he LFSRs se o zero. Hence, our claim follows from Theorem 1. For small n s, we explicily calculaed he algebraic equaions. Table 3 compares he upper bounds on he degree of he algebraic equaion claimed by various mehods. We believe his able is big enough o cover any pracically usable summaion generaor and should serve as a good reference for anyone implemening a summaion generaor and considering is immuniy o algebraic aacks.

8 n [1, 4] Thm explici calc Table 3. Degree bounds on algebraic equaions for summaion generaors 6 Reducing he degree furher for he n = 2 k case For he case when n = 2 k, we may reduce he degree of he algebraic equaion a lile bi furher, assuming ha we have access o consecuive key sream bis. We will apply he aack described in [5] due o Courois. Following he noaion of [2, 5], we classify mulivariae equaions ha relae key bis k i and oupu bis z j ino ypes given by heir degrees of k i and z j. We say ha a polynomial is of ype k d z f if all of is monomial erms are of he form k i1 k id z j1 z jf. Capial leers will be used in a similar manner o denoe ypes of equaions ha may conain lower degree monomials also. For example, K 2 = k 2 k 1 and KZ = kz k z 1. In [5], double-decker equaions (DDE) of degree (d, e, f) are defined o be any mulivariae equaion of ype K d K e Z f. Cases where d > e are of ineres from he aacker s poin of view. The following proposiion saes ha he equaions describing he summaion generaor on 2 k inpu LFSRs, given by Theorem 1, are DDEs. Theorem 3. A double-decker equaion of degree (2 k, 2 k 1, 2 k 1 ) ha relaes iniial key bis and oupu sream bis exiss for he summaion generaor on 2 k LFSRs. Proof. Le us, once more, recall he n = 2 2 case. The following is a very simple illusraion of he process we wen hrough in obaining he degree 4 equaion. (18) Z 1 = K 1 c 0 c 0 is of ype K 1 Z 1 (25) (19), (25) K 1 Z 1 = K 2 (K 1 Z 1 )K 1 c 1 c 1 is of ype K 2 K 1 Z 1 (26) (20), (26) K 2 K 1 Z 1 = (K 2+2 K 1+2 Z 1 ) (K 2 K 1 Z 1 )(K 2 K 1 Z 1 ) relaion of ype K 4 K 3 Z 2 (27) Le us nex consider he n = 2 3 case also. Since changing he number of variables, i.e., replacing S j (22 ) wih S j (23 ) does no change he degree of hese

9 equaions, he degree 2 3 equaion is obained as follows. (18) = (21), (25) c 0 is of ype K 1 Z 1 (28) (19) = (22), (26) c 1 is of ype K 2 K 1 Z 1 (29) (20) (23), (27) c 2 is of ype K 4 K 3 Z 2 (30) (24), (30) K 4 K 3 Z 2 = (K 4+4 K 3+4 Z 2 ) (K 4 K 3 Z 2 )(K 4 K 3 Z 2 ) relaion of ype K 8 K 7 Z 4 (31) This shows ha he general case may be proved by inducion. To prove he inducion sep, i suffices o show ha K 2k K 2k 1 Z 2k 1 = (K 2k +2 k K 2k 1+2 k Z 2k 1 ) (K 2k K 2k 1 Z 2k 1 )(K 2k K 2k 1 Z 2k 1 ) gives a DDE of ype K 2k+1 K 2k+1 1 Z 2k. Checking he validiy of his saemen is rivial. And his complees he proof. We may assume ha all periods of he LFSRs used in he summaion generaor are relaively prime. Then he summaion generaor saisfies he requiremen of he aack described in [5], if we have access o consecuive key sream bis. The following seps may be aken o reduce he complexiy of he algebraic aack. 1. Compue an algebraic equaion explicily using he formulae given in Proposiion Wrie he resuling DDE in he following form. L (k) = R (k, z). Here, L (k) is he sum of all monomials of ype K d appearing in he equaion and R (k, z) is he sum of all oher monomials of ype K e Z f. 3. Fix an arbirary nonrivial iniial key k and compue he value L (k ) for a sequence of lengh 2 ( ) m d. 4. Using he Berlekamp-Massey algorihm, find a linear relaion α = (α ) such ha α L (k ) = 0. We noe ha seps 1,2,3,4 are independen of he iniial key k, hence we can pre-compue he relaion α. 5. We have obained an algebraic equaion of degree e. I is given by α R (k, z) = Apply he general algebraic aack given in [6] o he above equaion.

10 Example 1. The DDE for n = 2 2 is given as follows: S4 S3 S1 +1 S 2 S2 +1 S2 S1S 1 +1 = S1 +2 S1 +1 S1 z S3 z +1 S2 z S1S +1 1 z +1 S1 +1 z z +1 S1. z +1 S1 z S1 z +2 We ake 4 LFSRs defined by he following characerisic polynomials. L 1 : x 3 + x + 1 L 2 : x 5 + x L 3 : x 7 + x + 1 L 4 : x 11 + x Since ( ) , 000, we compue 30,000 bis from he lef hand side of he above equaion for some arbirary nonrivial key bis k. Afer applying he Berlekamp-Massey algorihm, we found ha he sequence has a linear relaion of lengh 3892 < ( ) Wih a consecuive key sream of lengh of he order 6520 = (3 1) + ( ( ) ), one will be able o find he 26 bi iniial key k. Table 4 gives a simple comparison of he daa and compuaional complexiies needed for aacks on summaion generaors. Le w be he Gaussian eliminaion exponen. We shall use w = log 2 7, as given by he Srassen algorihm. Le he summaion generaor wih 2 k inpu LFSRs use an m-bi iniial key. generaor size (m, 2 k ) (128, 2 2 ) (256, 2 2 ) (256, 2 3 ) ( ) m daa k 1 (k+1) [1, 4] ) w compuaion ( m 2 k 1 (k+1) daa T = 2 m 2 k +k [8] compuaion T 2 log 2 T log 2 log 2 T ( daa m ) Thm 1 ( k compuaion m ) w ( k daa m ) k 1 Thm 3 ( compuaion m ) w k 1 Table 4. Complexiy comparison of aacks on summaion generaors We remark ha for [8] and Theorem 3, he key sream needs o be consecuive. For [1, 4] and Theorem 1, he key sream need only be parially consecuive, i.e., we need groups of k + 1 consecuive bis, bu hese groups may be far apar from each oher. Hence, a sraighforward comparison of daa complexiy migh no be fair. Also, he values for [8] have been calculaed assuming ha he LFSRs in use have been chosen well, so ha heir 2-adic span is maximal.

11 7 Conclusion We have applied he general resuls of [1, 4] and [5] on sream ciphers wih memories o he summaion generaor. Our resuls show ha he degree of algebraic equaion obainable and he complexiy of he aack applicable are much lower han given by he general resuls. For a summaion generaor ha uses n LFSRs, he algebraic equaion relaing he key sream bis and LFSR oupu bis can be made o be of degree less han or equal o 2 log 2 n, using log 2 n + 1 consecuive key sream bis. Under cerain condiions, for he n = 2 k case, he effecive degree may furher be reduced by 1. And for small n s we have summarized he degrees of he explici equaions in Table 3. The able should be aken ino accoun by anyone using a summaion generaor. References 1. F. Armknech and M. Krause, Algebraic aacks on combiners wih memory, Advances in Crypology - Crypo 2003, LNCS 2729, Springer-Verlag, pp , N. Courois, The securiy of Hidden Field Equaions (HFE), CT-RSA 2001, LNCS 2020, Springer-Verlag, pp , N. Courois, Higher order correlaion aacks, XL algorihm and Crypanalysis of Toyocryp, ICISC 2002, LNCS 2587, Springer-Verlag, pp , N. Courois, Algebraic aacks on combiners wih memory and several oupus, E-prin archive, 2003/ N. Courois, Fas algebraic aack on sream ciphers wih linear feedback, Advances in Crypology - Crypo 2003, LNCS 2729, Springer-Verlag, pp , N. Courois and W. Meier, Algebraic aacks on sream ciphers wih linear feedback, Advances in Crypology - Eurocryp 2003, LNCS 2656, Springer-Verlag, pp , N. Courois and J. Pieprzyk, Crypanalysis of block ciphers wih overdefined sysems of equaions, Asiacryp 2002, LNCS 2501, Springer-Verlag, pp , A. Klapper and M. Goresky, Crypanalysis based on 2-adic raional approximaion, Advances in Crypology - Crypo 95, LNCS 963, Springer-Verlag, pp , W. Meier and O. Saffelbach, Correlaion Properies of Combiners wih Memory in Sream Cipher, Journal of Crypology, vol.5, pp.67 86, R.A. Rueppel, Correlaion immuniy and he summaion generaor, Advances in Crypology - Crypo 85, LNCS 219, Springer-Verlag, pp , A Proof of Proposiion 2, equaion (12) To prove (12), we shall evaluae is righ hand side a some arbirary inpu value and show ha i equals he nex carry bi i. Bu before we do his, le us make some observaions. Noe ha we may ake i ( j x j + c )/2 i+1 (mod 2) (32) as he definiion of he carry bis. We may evaluae he righ hand side of (11) a j x j = 0 and equae i wih he evaluaion of (32) a he same poin o shows c i+1 = f i,0. Now, from he discussions in Secion 2 on he marix M, we know ha a j x j = 2k, we have S0 (2k ) = S (2 k ) = 1 and all oher S 2 k j (2k ) = 0. Once more, evaluaing (11) and (32) a j x j = 2k wih his in mind shows

12 c i+1 = f i,0 f i,2 k. Since we already know c i+1 = f i,0, his implies f i,2 k = 0, or equivalenly, ha he erm S does no appear in he linear sum (11). 2 k We shall now evaluae he righ hand side of (12) a some fixed LFSR oupu values (x 1,..., x ) and carry value c. Se 2 k+1 x = j x j, x = remainder of x divided by 2 k, c = remainder of c divided by 2 k. From he above observaion, we know ha he righ hand side of (12) conains some of he erms S0 (2k+1 ),..., S 2 k 1 (2k+1 ), bu does no conain any of he erms S j (2 k+1 ) for j 2 k. We also know from discussions of Secion 2 ha he evaluaion of each Sj (2k+1 ) a x for j < 2 k is equal o is evaluaion a x. Noe also ha he erm c does no appear as inpu o any of he coefficiens f 2 k i,j in he righ hand side of (12). Hence, RHS of (12) a x and c = RHS of (12) a x and c = RHS of (11) a x and c = (x + c )/2 i+1 (mod 2) = (x + c )/2 i+1 (mod 2) = i a x and c. The condiion i k 2 has been used in he fourh equaliy. And he las equaliy is jus (32). We have compleed he proof ha (12) is a valid expression for he nex carry bis. B Proof of Proposiion 2, equaion (13) The proof for (13) is very similar o ha of (12) and hence we shall be very brief. We ask he readers o read Appendix A before reading his secion. The carry bi is given by k 1 ( j x j + c )/2 k (mod 2). Evaluaion of (11) a j x j = 0 and 2k shows ha f k 1,0 = 0 and f k 1,2 k = 1, hence we now always have S as a linear erm, wih coefficien equal o 1, in 2 k he sums (11) and (13). The emporary values x, x, and c may be defined as before. From he discussions of Secion 2, one may wrie ( S 2 k (2 k+1 ) a x ) = x/2 k = ( S 2 k (2 k+1 ) a x ) x/2 k (mod 2) and ( S j (2 k+1 ) a x ) = ( S j(2 k+1 ) a x )

13 for j < 2 k. Also noe ha none of he erms S j (2 k+1 ), for j > 2 k, appears in he sum (13) and ha he erm c k visible in (13) is is only use in (13). Hence, RHS of (13) a x and c = ( RHS of (13) a x and c ) c k = ( RHS of (13) a x and c ) c k x/2k (mod 2) = ( RHS of (11) a x and c ) c k x/2k (mod 2) = (x + c )/2 k c k x/2k (mod 2) = (x + c )/2 k (mod 2) = k 1 a x and c. This complees he proof. C Proof of Proposiion 2, equaion (14) Define x = j x j. Careful reading of Appendices A and B shows ha he firs erm of (14) simplifies o { j f k 1,j S j+2 (2 k+1 0 for 0 x 2 k, ) = k (x 2 k + c )/2 k c k for 2 k < x 2 k+1 (33), and ha he second erm saisfies c k( j f k 1,j S j(2 k+1 ) ) = c k( (x + c )/2 k c k). (34) I now suffices o compare he sum of hese wo values wih k (x + c )/2 k+1 (mod 2). (35) Define x = x 2 k when x > 2 k and define c = c 2 k when c k = 1. Case 1) x 2 k, c k = 0. This is he mos easy case. Case 2) x 2 k, c k = 1. Case 3) x > 2 k, c k = 0. (33) (34) = 0 = (x + c )/2 k+1. (33) (34) = (x + c )/2 k 1 = (x + c )/2 k = {(x + c ) + 2 k }/2 k+1 = (x + c)/2 k+1. (33) (34) = (x + c )/2 k = {(x + c ) + 2 k }/2 k+1 = (x + c )/2 k+1.

14 Case 4) x > 2 k, c k = 1. (33) (34) = (x + c )/2 k (x + c )/2 k = (x + c )/2 k ( (x + c )/2 k 1) = 1 = (x + c )/2 k+1. This complees he proof.

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

Bernoulli numbers. Francesco Chiatti, Matteo Pintonello. December 5, 2016

Bernoulli numbers. Francesco Chiatti, Matteo Pintonello. December 5, 2016 UNIVERSITÁ DEGLI STUDI DI PADOVA, DIPARTIMENTO DI MATEMATICA TULLIO LEVI-CIVITA Bernoulli numbers Francesco Chiai, Maeo Pinonello December 5, 206 During las lessons we have proved he Las Ferma Theorem

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

Echocardiography Project and Finite Fourier Series

Echocardiography Project and Finite Fourier Series Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every

More information

Matlab and Python programming: how to get started

Matlab and Python programming: how to get started Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,

More information

Linear Cryptanalysis

Linear Cryptanalysis Linear Crypanalysis T-79.550 Crypology Lecure 5 February 6, 008 Kaisa Nyberg Linear Crypanalysis /36 SPN A Small Example Linear Crypanalysis /36 Linear Approximaion of S-boxes Linear Crypanalysis 3/36

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

FITTING EQUATIONS TO DATA

FITTING EQUATIONS TO DATA TANTON S TAKE ON FITTING EQUATIONS TO DATA CURRICULUM TIDBITS FOR THE MATHEMATICS CLASSROOM MAY 013 Sandard algebra courses have sudens fi linear and eponenial funcions o wo daa poins, and quadraic funcions

More information

Chapter 7: Solving Trig Equations

Chapter 7: Solving Trig Equations Haberman MTH Secion I: The Trigonomeric Funcions Chaper 7: Solving Trig Equaions Le s sar by solving a couple of equaions ha involve he sine funcion EXAMPLE a: Solve he equaion sin( ) The inverse funcions

More information

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n Lecure 3 - Kövari-Sós-Turán Theorem Jacques Versraëe jacques@ucsd.edu We jus finished he Erdős-Sone Theorem, and ex(n, F ) ( /(χ(f ) )) ( n 2). So we have asympoics when χ(f ) 3 bu no when χ(f ) = 2 i.e.

More information

Math 315: Linear Algebra Solutions to Assignment 6

Math 315: Linear Algebra Solutions to Assignment 6 Mah 35: Linear Algebra s o Assignmen 6 # Which of he following ses of vecors are bases for R 2? {2,, 3, }, {4,, 7, 8}, {,,, 3}, {3, 9, 4, 2}. Explain your answer. To generae he whole R 2, wo linearly independen

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

More information

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4. PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence

More information

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

THE MATRIX-TREE THEOREM

THE MATRIX-TREE THEOREM THE MATRIX-TREE THEOREM 1 The Marix-Tree Theorem. The Marix-Tree Theorem is a formula for he number of spanning rees of a graph in erms of he deerminan of a cerain marix. We begin wih he necessary graph-heoreical

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

CHAPTER 12 DIRECT CURRENT CIRCUITS

CHAPTER 12 DIRECT CURRENT CIRCUITS CHAPTER 12 DIRECT CURRENT CIUITS DIRECT CURRENT CIUITS 257 12.1 RESISTORS IN SERIES AND IN PARALLEL When wo resisors are conneced ogeher as shown in Figure 12.1 we said ha hey are conneced in series. As

More information

LAB # 2 - Equilibrium (static)

LAB # 2 - Equilibrium (static) AB # - Equilibrium (saic) Inroducion Isaac Newon's conribuion o physics was o recognize ha despie he seeming compleiy of he Unierse, he moion of is pars is guided by surprisingly simple aws. Newon's inspiraion

More information

Testing for a Single Factor Model in the Multivariate State Space Framework

Testing for a Single Factor Model in the Multivariate State Space Framework esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

More information

14 Autoregressive Moving Average Models

14 Autoregressive Moving Average Models 14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class

More information

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1 SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision

More information

Notes for Lecture 17-18

Notes for Lecture 17-18 U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

More Digital Logic. t p output. Low-to-high and high-to-low transitions could have different t p. V in (t)

More Digital Logic. t p output. Low-to-high and high-to-low transitions could have different t p. V in (t) EECS 4 Spring 23 Lecure 2 EECS 4 Spring 23 Lecure 2 More igial Logic Gae delay and signal propagaion Clocked circui elemens (flip-flop) Wriing a word o memory Simplifying digial circuis: Karnaugh maps

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

Cryptanalysis of Reduced NORX initially discussed at ASK 2016

Cryptanalysis of Reduced NORX initially discussed at ASK 2016 Crypanalysis of Reduced NORX iniiy discussed a ASK 2016 Nasour Bagheri (SRTTU, Iran), Tao Huang (NTU, Singaopre), Keing Jia (Tsinghua Univ., Sae Key Lab, China), Florian Mendel (TUGraz, Ausria), Yu Sasaki

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

THE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2.

THE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2. THE BERNOULLI NUMBERS The Bernoulli numbers are defined here by he exponenial generaing funcion ( e The firs one is easy o compue: (2 and (3 B 0 lim 0 e lim, 0 e ( d B lim 0 d e +e e lim 0 (e 2 lim 0 2(e

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff Laplace ransfom: -ranslaion rule 8.03, Haynes Miller and Jeremy Orloff Inroducory example Consider he sysem ẋ + 3x = f(, where f is he inpu and x he response. We know is uni impulse response is 0 for

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

Spring Ammar Abu-Hudrouss Islamic University Gaza

Spring Ammar Abu-Hudrouss Islamic University Gaza Chaper 7 Reed-Solomon Code Spring 9 Ammar Abu-Hudrouss Islamic Universiy Gaza ١ Inroducion A Reed Solomon code is a special case of a BCH code in which he lengh of he code is one less han he size of he

More information

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion

More information

Some Ramsey results for the n-cube

Some Ramsey results for the n-cube Some Ramsey resuls for he n-cube Ron Graham Universiy of California, San Diego Jozsef Solymosi Universiy of Briish Columbia, Vancouver, Canada Absrac In his noe we esablish a Ramsey-ype resul for cerain

More information

Second Order Linear Differential Equations

Second Order Linear Differential Equations Second Order Linear Differenial Equaions Second order linear equaions wih consan coefficiens; Fundamenal soluions; Wronskian; Exisence and Uniqueness of soluions; he characerisic equaion; soluions of homogeneous

More information

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX J Korean Mah Soc 45 008, No, pp 479 49 THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX Gwang-yeon Lee and Seong-Hoon Cho Reprined from he Journal of he

More information

GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT

GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT Inerna J Mah & Mah Sci Vol 4, No 7 000) 48 49 S0670000970 Hindawi Publishing Corp GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT RUMEN L MISHKOV Received

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

Georey E. Hinton. University oftoronto.   Technical Report CRG-TR February 22, Abstract Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical

More information

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively: XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

BU Macro BU Macro Fall 2008, Lecture 4

BU Macro BU Macro Fall 2008, Lecture 4 Dynamic Programming BU Macro 2008 Lecure 4 1 Ouline 1. Cerainy opimizaion problem used o illusrae: a. Resricions on exogenous variables b. Value funcion c. Policy funcion d. The Bellman equaion and an

More information

BBP-type formulas, in general bases, for arctangents of real numbers

BBP-type formulas, in general bases, for arctangents of real numbers Noes on Number Theory and Discree Mahemaics Vol. 19, 13, No. 3, 33 54 BBP-ype formulas, in general bases, for arcangens of real numbers Kunle Adegoke 1 and Olawanle Layeni 2 1 Deparmen of Physics, Obafemi

More information

Cryptanalysis of RAKAPOSHI Stream Cipher

Cryptanalysis of RAKAPOSHI Stream Cipher Crypanalysis of RAKAPOSHI Sream Cipher Lin Ding, Jie Guan Zhengzhou Informaion Science and Technology Insiue, China E-mail: dinglin_cipher@63.com; guanie7@63.com Absrac. RAKAPOSHI is a hardware oriened

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu ON EQUATIONS WITH SETS AS UNKNOWNS BY PAUL ERDŐS AND S. ULAM DEPARTMENT OF MATHEMATICS, UNIVERSITY OF COLORADO, BOULDER Communicaed May 27, 1968 We shall presen here a number of resuls in se heory concerning

More information

Orthogonal Rational Functions, Associated Rational Functions And Functions Of The Second Kind

Orthogonal Rational Functions, Associated Rational Functions And Functions Of The Second Kind Proceedings of he World Congress on Engineering 2008 Vol II Orhogonal Raional Funcions, Associaed Raional Funcions And Funcions Of The Second Kind Karl Deckers and Adhemar Bulheel Absrac Consider he sequence

More information

Analyze patterns and relationships. 3. Generate two numerical patterns using AC

Analyze patterns and relationships. 3. Generate two numerical patterns using AC envision ah 2.0 5h Grade ah Curriculum Quarer 1 Quarer 2 Quarer 3 Quarer 4 andards: =ajor =upporing =Addiional Firs 30 Day 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 andards: Operaions and Algebraic Thinking

More information

3.1 More on model selection

3.1 More on model selection 3. More on Model selecion 3. Comparing models AIC, BIC, Adjused R squared. 3. Over Fiing problem. 3.3 Sample spliing. 3. More on model selecion crieria Ofen afer model fiing you are lef wih a handful of

More information

= ( ) ) or a system of differential equations with continuous parametrization (T = R

= ( ) ) or a system of differential equations with continuous parametrization (T = R XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17 EES 16A Designing Informaion Devices and Sysems I Spring 019 Lecure Noes Noe 17 17.1 apaciive ouchscreen In he las noe, we saw ha a capacior consiss of wo pieces on conducive maerial separaed by a nonconducive

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

Section 4.4 Logarithmic Properties

Section 4.4 Logarithmic Properties Secion. Logarihmic Properies 59 Secion. Logarihmic Properies In he previous secion, we derived wo imporan properies of arihms, which allowed us o solve some asic eponenial and arihmic equaions. Properies

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

KEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow KEY Mah 334 Miderm III Fall 28 secions and 3 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering

More information

Math Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later.

Math Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later. Mah 50-004 Week 5: Secion 7.4, mass-spring sysems. These are noes for Monday. There will also be course review noes for Tuesday, posed laer. Mon Apr 3 7.4 mass-spring sysems. Announcemens: Warm up exercise:

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a

More information

Logic in computer science

Logic in computer science Logic in compuer science Logic plays an imporan role in compuer science Logic is ofen called he calculus of compuer science Logic plays a similar role in compuer science o ha played by calculus in he physical

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

Computer-Aided Analysis of Electronic Circuits Course Notes 3

Computer-Aided Analysis of Electronic Circuits Course Notes 3 Gheorghe Asachi Technical Universiy of Iasi Faculy of Elecronics, Telecommunicaions and Informaion Technologies Compuer-Aided Analysis of Elecronic Circuis Course Noes 3 Bachelor: Telecommunicaion Technologies

More information

Generalized Chebyshev polynomials

Generalized Chebyshev polynomials Generalized Chebyshev polynomials Clemene Cesarano Faculy of Engineering, Inernaional Telemaic Universiy UNINETTUNO Corso Viorio Emanuele II, 39 86 Roma, Ialy email: c.cesarano@unineunouniversiy.ne ABSTRACT

More information

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method Journal of Applied Mahemaics & Bioinformaics, vol., no., 01, 1-14 ISSN: 179-660 (prin), 179-699 (online) Scienpress Ld, 01 Improved Approimae Soluions for Nonlinear Evoluions Equaions in Mahemaical Physics

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

ON THE DEGREES OF RATIONAL KNOTS

ON THE DEGREES OF RATIONAL KNOTS ON THE DEGREES OF RATIONAL KNOTS DONOVAN MCFERON, ALEXANDRA ZUSER Absrac. In his paper, we explore he issue of minimizing he degrees on raional knos. We se a bound on hese degrees using Bézou s heorem,

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED

0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED 0.1 MAXIMUM LIKELIHOOD ESTIMATIO EXPLAIED Maximum likelihood esimaion is a bes-fi saisical mehod for he esimaion of he values of he parameers of a sysem, based on a se of observaions of a random variable

More information

THE NEW HEURISTIC GUESS AND DETERMINE ATTACK ON SNOW 2.0 STREAM CIPHER

THE NEW HEURISTIC GUESS AND DETERMINE ATTACK ON SNOW 2.0 STREAM CIPHER THE NEW HEURISTIC GUESS AND DETERMINE ATTACK ON SNOW 0 STREAM CIPHER Mohammad Sadegh Nemai Nia 1, Ali Payandeh 1, Faculy of Informaion, Telecommunicaion and Securiy Technologies, Malek-e- Ashar Universiy

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3. Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies

More information

Homework sheet Exercises done during the lecture of March 12, 2014

Homework sheet Exercises done during the lecture of March 12, 2014 EXERCISE SESSION 2A FOR THE COURSE GÉOMÉTRIE EUCLIDIENNE, NON EUCLIDIENNE ET PROJECTIVE MATTEO TOMMASINI Homework shee 3-4 - Exercises done during he lecure of March 2, 204 Exercise 2 Is i rue ha he parameerized

More information

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013 IMPLICI AND INVERSE FUNCION HEOREMS PAUL SCHRIMPF 1 OCOBER 25, 213 UNIVERSIY OF BRIISH COLUMBIA ECONOMICS 526 We have exensively sudied how o solve sysems of linear equaions. We know how o check wheher

More information

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin ACE 56 Fall 005 Lecure 4: Simple Linear Regression Model: Specificaion and Esimaion by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Simple Regression: Economic and Saisical Model

More information

Explore 2 Proving the Vertical Angles Theorem

Explore 2 Proving the Vertical Angles Theorem Explore 2 Proving he Verical Angles Theorem The conjecure from he Explore abou verical angles can be proven so i can be saed as a heorem. The Verical Angles Theorem If wo angles are verical angles, hen

More information