State Estimation with ARMarkov Models
|
|
- Kelley Hutchinson
- 6 years ago
- Views:
Transcription
1 Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University, New York, NY Minh Q. Phan 2 Princeton University, Princeton, NJ Richard W. Longman 3 Columbia University, New York, NY Abstract The ARMarkov models were originally develoed for adative neural control, and later for redictive control, and state-sace identification. Recently, an interaction matrix formulation has been develoed that exlains the internal structure of the ARMarkov models and their connection to the state-sace reresentation. Using the interaction matrix formulation, we show in this aer how a state estimator can be identified directly from inut-outut data. The conventional aroach is to design such a state estimator from knowledge of the lant, and the difficult-to-obtain rocess and measurement noise statistics. A numerical examle comares the identified state estimator with an otimal Kalman filter derived with erfect knowledge of the lant and noise statistics. 1 Graduate Student, Deartment of Mechanical Engineering. 2 Assistant Professor, Deartment of Mechanical and Aerosace Engineering, Dissertation Advisor. 3 Professor, Deartment of Mechanical Engineering, Dissertation Co-Advisor.
2 1. Introduction State estimation is an imortant element of modern control theory. Given a known model of the system under the influence of rocess and measurement noise with known statistics secified in terms of their covariances, it is well known that the Kalman filter is an otimal state estimator in the sense that its state estimation error is minimized. In ractice, it is difficult to design such an otimal estimator because neither the system nor the noise statistics can be known exactly. From the oint of view of system identification, information about the system and the noise statistics are embedded in a sufficiently long set of inut-outut data. Thus it would be advantageous to be able to obtain such an estimator directly from inut-outut data without having to identify the system and the noise statistics searately. This is the roblem of observer identification. Recently, a class of models known as ARMarkov models has been develoed in the context of adative neural control, Ref. 1. The term ARMarkov refers to Auto- Regressive model with exlicit Markov arameter coefficients. ARMarkov models form a bridge between the common ARX model (Auto-Regressive model with exogenous inuts) where the Markov arameters are imlicit, and the non-auto-regressive ulse resonse model where every coefficient is a Markov arameter. Later the ARMarkov models are used for state-sace system identification, Refs In articular, it was found that when the ARMarkov models are used to identify the system Hankel matrix, the true or effective order of the system can be detected more effectively than with an ARX model. This issue has been investigated extensively in Ref. 4. In fact, with an ARX 2
3 model it is also ossible to identify a state-sace model together with an observer gain as shown in Refs. 5 and 6. But with this technique it is not always ossible to determine the order of the state sace realization by Hankel singular value truncation alone and a searate ost-identification model reduction rocedure must be used. With ARMarkov models, we have the oortunity to identify state estimators with true or effective orders without having to invoke a searate model reduction ste as normally required. This is one motivation for the resent aer. It is clear that for state estimation, efficient detection of the dimension of the effective state sace model is imortant because it is comutationally a burden to have a state estimator with unnecessarily large dimensions. The ARMarkov model used in Ref. 4 is based on the develoment of the interaction matrix in Ref. 7 that can be exlained in terms of a generalization of the wellknown Cayley-Hamilton theorem. In all of these develoments the role of the interaction matrix has been to justify the structures of various inut-outut models and the relationshi among the coefficients of these models, but there has been no need to identify the interaction matrix itself from inut-outut data. The second motivation of this aer is an investigation of the role of this interaction matrix in the context of state estimation as oosed to treating it as a convenient mathematical construct. In this work, we establish that fact that identifying a state estimator amounts to identifying this interaction matrix. Another interesting asect of this formulation is that this interaction matrix based state estimator has a non-standard form, different from the usual form of a Luenberger observer or a Kalman filter. However, it will be shown that this new form is very convenient both from the oint of view of state estimator design as well as from the oint of view of its 3
4 identification from inut-outut data. This oint makes u the third motivation for this work. In this aer we will quickly derive such a state estimator using the interaction matrix. We then show how this state estimator can be identified from inut-outut data. This identification will first be derived in the deterministic (noise-free) setting. Then a stochastic analysis will be carried out that shows why the calculations involved in the deterministic case are indeed justified in the stochastic case. Following the theoretical justification, a numerical examle illustrates how an identified state estimator using the technique develoed here comares to that of an otimal Kalman filter designed with erfect knowledge of the system and erfect knowledge of the noise statistics. In articular, we show that the outut residuals obtained with this identified state estimator indeed match the residuals of the otimal Kalman filter. 2. State Estimation by Interaction Matrix In the following we briefly derive a state estimator via an interaction matrix. Consider an n-th order, r-inut, m-outut discrete-time model of a system in state-sace format xk ( + 1) = Axk ( ) + Buk ( ) y( k) = Cx( k) + Du( k) (1) 4
5 By reeated substitution, we have for some 0, x( k + ) = A x( k) + Cu ( k) y ( k) = Ox( k) + Tu ( k) (2) where u ( k ) and y ( k ) are defined as column vectors of inut and outut data going stes into the future starting with u( k) and yk ( ), resectively, u uk ( ) uk ( + 1) ( k) = M uk ( + 1) yk ( ) yk ( + 1), y( k) = M yk ( + 1) (3) For a sufficiently large, C in Eq. (2) is an n r controllability matrix, O is a m n observability matrix, T is a m r Toelitz matrix of the system Markov arameters, C = [ A 1 B, K, AB, B], O = C CA M CA 1 D 0 0 L 0 CB D O O M, T = CAB CB D O 0 M O O O 0 CA 2 B K CAB CB D (4) As long as m n, it is guaranteed for an observable system that an interaction matrix M exists such that A + MO = 0. The existence of M ensures that for k 0 an exression 5
6 for xk ( + ) exists where the state variable is comletely eliminated from the right hand side of Eq. (2), x(k + ) = A x(k) + Cu (k) = A x(k) + Cu (k) + M[ Ox(k) + T u (k)] My (k) = ( A + MO)x(k) + ( C + MT )u (k) My (k) = ( C + MT )u (k) My (k) (5) Shifting the time indices back by time stes, we have for k an exression that relates the current state of the system in terms of ast inut and ast outut measurements, x(k) = α i u(k i) + β i y(k i) (6) i=1 i=1 where [ α, K, α 2, α 1 ]= C+ MT, and [ β, K, β 2, β 1 ]= M. Note that these formulas are alicable to both the single-inut single-outut and multile-inut multile-outut cases. Equation (6) is a state estimator, but not in the standard form of a Luenberger observer or a Kalman filter. In fact, this non-standard form is quite convenient from both ersectives of design and identification. If an observable state-sace model of the system is known, to design this state estimator, one simly forms the controllability matrix C, the Toelitz matrix T, comute an interaction matrix M from A + MO = 0, 6
7 m n. When is chosen such that m > n then M is no longer unique, in which case the solution M = A + O where O + denotes the seudo-inverse of O, roduces a minimum-norm solution for M, whose elements are recisely the outut gains of this state estimator. We have focused on state estimation, but outut estimation can be similarly derived. It is simly, ( ) + y( k) = C C+ MT u ( k ) CMy ( k ) Du( k) (7) + and the minimum norm solution for CM is CA O, whose elements are the outut gains for the outut estimator is CM = CA + O. On the other hand, if the system is unknown but inut-outut data is available, then the arameter combinations C + MT and M can be comuted for the state estimator. This is shown in the next section. 3. Identification of State Estimator with ARMarkov Models We must first derive an inut-outut exression that involves C + MT and M by combining Eq. (5) with the outut exression in Eq. (2), y ( k) = Ox( k) + Tu ( k) ( ) + = OC+ MTu ( k ) OMy ( k ) Tu ( k) (8) 7
8 Equation (8) is a family of ARMarkov models. The first member of the family is an ARX model for yk ( ), which is a secial case of an ARMarkov model. All remaining members are exressions for yk ( + 1 ),, yk ( + 1 ) which are true ARMarkov models. Each ARMarkov model is different from another in that it has an increasing number of Markov arameters aearing exlicitly as coefficients. We need not one but such a family of ARMarkov models to solve for C + MT and M. Furthermore, we assume to have only outut measurements and not the full state, therefore the above inut-outut exression does not have C + MT and M aear exlicitly as coefficients, but the combinations O( C+ MT ) and OM (and T ). Since O is not known, the identification of C + MT and M from inut-outut data is in fact a non-linear roblem. Fortunately, an exact solution can be found without any kind of iterations as shown below. For simlicity define A = OC+ ( MT ), B = OM, then A, B, and T can be identified from inut-outut data as T T [ A T B ]= YV ( VV ) + (9) where u ( 2 k ) combines u ( k) with u ( k) + for convenience, [ ] Y = y ( ) y ( + 1) L y ( + l ) (10) u u u V = 2 ( 0) 2 ( 1) L 2 ( l) y( 0) y( 1) L y ( l) (11) 8
9 To find ( C + MT ) and M from A and B we must first find O. Taking advantage of the internal structure of the coefficients A, B as revealed by the interaction matrix, the combination H = OC can be comuted from H = A + BT (12) The observability matrix O is obtained by a singular value decomosition of H = U Σ V T n n n, where n is the order of the system, O = U n Σ 1 2 n / (13) Since the state-sace reresentation is uncertain u to a similarity transformation of the state vector, we generally have O in a different set of coordinates from O but they are related to each other by a similarity transformation, O = OT. The identified arameter combinations, however, are invariant with resect to such a transformation, ( ) 1 1 A = OC ( + MT)= OT( T C+ T MT)= OC+ MT 1 B = OM = OTT M = OM (14) 9
10 Hence, as long as is chosen to be sufficiently large such that the observability matrix O or O has rank n, the needed arameters ( C + M T) and M to construct a state estimator can be found from C + T T [( MT) M]= ( O O) 1 O [ A B ] (15) We thus see that it is not necessary to roduce a realization of A, B, C in the above stes of extracting a state estimator from inut-outut data. However, this ste can be easily done as well, and would be necessary if we desire to ut the state estimator in a different set of coordinates than that chosen by the above realization, such as the modal coordinates. Any realization algorithm can be used for that urose. Here we review the realization rovided by the Eigensystem Realization Algorithm (ERA), Ref. 8. The rocedure calls for the extraction of two Hankel matrices H( 0 ) and H( 1 ) from H, H( 0) = CB L M M M n1 CA B n2 n1+ n2 CA B L CA B, H( 1) = n1 + 1 CAB L CA B M M M n2 + 1 n1+ n2 + 1 CA B L CA B (16) The matrix H has the Markov arameters of increasing order going from left to right, whereas H( 0 ) and H( 1 ) are tyically defined with the higher order Markov arameters going from right to left. To maintain the standard notation, a trivial rearrangement of the Markov arameters is needed in forming H( 0 ) and H( 1 ). A s-th order state-sace model T is A = 1/ 2 Σ U VΣ 1/ 2 H () 1, B is the first r columns of Σ 1 / 2 T s V s, C is the first m rows of s s s s 10
11 U s Σ 1 2 s ( ) /, and s min rn ( + 1), mn ( + 1 ). The matrix U s and V s are made u of s left and 1 2 right singular vectors of H( 0 ), resectively. The diagonal matrix Σ s is made u of s corresonding singular values of H( 0 ). With erfect data H( 0 ) has exactly n ositive singular values (all remaining singular values are identically zero), where n is the true minimum order of the system, s = n. Otherwise, the user can secify the order of the state-sace model by s, the number of Hankel singular values to retain. The articular set of coordinates of the realization in O has a secial roerty that it becomes internally balanced when is large and the system is stable, Ref. 9. A realization is said to be internally balanced if and only if the controllability and observability grammians are equal to each other, and both equal to Σ n. As mentioned, the realization can also be ut in another user-secified set of coordinates. We simly need to comute the corresonding observability matrix O in that set of coordinates and use Eq. (15) to obtain the corresonding state estimator. Stochastic Analysis In the revious section the identification of the state estimator is derived in the deterministic setting. Now we will consider the situation where data is corruted by rocess and measurement noise, and show how the same deterministic calculations are justified in the stochastic case. Consider the case where rocess and measurement noise are resent in Eq. (1). The corresonding version of Eq. (2) is 11
12 x( k + ) = A x( k) + Cu ( k) + v( k) y ( k) = Ox( k) + Tu ( k) + w ( k) (17) Although Eq. (17) is -ste ahead, i.e., it relates xk ( + ) to xk ( ) with samling interval t, inut uk ( ), and outut yk ( ), it can be thought of as one-ste ahead with samling interval t, inut u with some gain K, ( k), and outut y ( k ). Hence it admits an estimator of the form + x ( k) = x ( k) + K[ y ( k) y ( k) x ]= ( k) + Ke( k) + x ( k) = A x ( k ) + C u ( k ) y ( k) = Ox ( k) + Tu ( k) (18) The quality of the estimation certainly deends on the gain K. For a choice of K, ek ( ) is the corresonding outut residual defined to be the difference between estimated outut yk ( ) and measured outut yk ( ), ek ( ) = yk ( ) yk ( ). The exression for the state in Eq. (18) can be written as x ( k) = A ( x ( k ) + K( y ( k ) y( k ) ))+ Cu( k ) ( ) + = A x ( k ) + A K Ox ( k ) + Tu ( k ) y ( k ) Cu ( k ) ( ) ( ) = A + A KO x ( k ) + C+ AKT u( k ) AKy( k ) (19) The above equation is somewhat subtle for the following reason. Due to the resence of the term ( A + A KO ), it is a -ste ahead state estimator, i.e., it estimates x ( ) from 12
13 ( ) x ( 0 ), then x ( 2 ) from x ( ), and so on. However, if A + A KO is zero (which will be justified later), then this equation can be used to rovide state estimation every single ste beginning with x ( ), then x ( + 1 ), x ( + 2 ), and so on. When combined with the outut exression we have y ( k) = OAx ƒ ( k ) + O C+ A KT u ( k ) OA Ky ( k ) T u ( k) e( k) (20) ( ) + + where A ƒ = A + A KO. From a given set of inut-outut data of sufficient length we can form the data matrices Y and V as defined in Eqs. (10) and (11), then Y = [ A T B ] V + O AX ƒ + E (21) ( ) [ ] where A = OC+ AK T, B = OA K, X = x ( 0) x () 1 K x () l, and e() 1 ( ) e() 1 ( + 1) K e() 1 ( + l) e e + e + ( 2) ( ) ( 2) ( 1) K ( 2) ( l) E = [ e ( ) e ( + 1) K e ( + l) ]= M M K M e ( ) e ( + ) e ( + ) ( ) ( ) 1 K ( ) l (22) The residual matrix E has the following interretation when A ƒ = 0. The first residual sequence e ( () k ) 1, k =, + 1,.., + l is associated with the ARX model derived from the first m rows of Eq. (20). The second residual sequence e ( ( ) k ) 2 is associated with the first ARMarkov model derived from the second m rows of Eq. (20). Similarly, the 13
14 remaining residual sequences e ( j) ( k), j = 3 through, are associated with the remaining ARMarkov models. Because not a single but a family of models are being used in the inut-outut ma of Eq. (20), we do not have a single but rather a family of residuals. Pre-multilying Eq. (21) by V T and re-arranging it yields, T T T T YV [ A T B ] VV = OAX ƒ V + EV (23) Now let us imose conditions on K so that the observer in Eq. (19) ossesses desirable roerties. If K is chosen such that ƒ A = 0 then by choosing T T [ A T B ]= YV ( VV ) + (24) the left hand side of Eq. (23) vanishes. Referring back to Eq. (21), with A ƒ = 0, the solution given in Eq. (24) is exactly the one that minimizes the Euclidean norm of E which is the sum of the squares of the residuals for the entire data record. Furthermore, from Eq. (23), this solution also results in EV T = 0, which can be written exlicitly as e e e ut ut ut ut () 1 ( ) () 1 ( + 1) K () 1 ( + l) ( 0) K ( 1) ( ) K ( 2 1) e( 2) ( ) e( 2) ( + 1) K e( 2) ( + l) ut() 1 K ut( ) ut( + 1) K ut( 2) M M K M = 0 M K M M K M e( ) ( ) e( ) ( + 1) e( ) ( + ) K l ut ut + ut + ut ( l) K ( 1 l) ( l) K ( l) (25) 14
15 e e e T T () 1 ( ) () 1 ( + 1) K () 1 ( + l) y ( 0) K y ( 1) e e e T T ( 2) ( ) ( 2) ( + 1) K ( 2) ( + l) y () 1 y ( K 1) = 0 M M K M M K M e( ) ( ) e( ) ( + 1) e( ) ( + ) K l T T y ( l) K y ( 1 + l) (26) Let us now examine the imlications of Eqs. (25) and (26) for the j-th member of the residual sequences e ( j) ( k). When the length of the data record tends to infinity, the above equations imly for a stationary random rocess, + l T 1 T E{ e( j) ( k) u ( k i) }= lim e( j) ( k) u ( k i) = 0, i = + 1,..., 0, 1, 2,.., (27) l l k= + l T 1 T E{ e( j) ( k) y ( k i) }= lim e( j) ( k) y ( k i) = 0, i = 1, 2,.., (28) l l k= where E {}. denotes the exectation oerator. Thus, if the data record is sufficiently long, the identified family of ARMarkov models has its combined residuals minimized, and in articular, Eqs. (27) and (28) state that the residual for each member of the family becomes uncorrelated with inut and outut data. Recall these results are obtained while imosing the condition on K such that A ƒ = A + A KO = 0. We note here that this condition can also be satisfied if is large and the system is stable so that A 0. By referring back to the deterministic formulation, it is clear that the interaction matrix M lays the same role as AK at every ste of the derivation, including the condition A + MO = 0. Hence, the deterministic calculations are indeed justified in the stochastic case. 15
16 Illustration Consider a chain of three masses connected by srings and damers with force inut to the first mass and osition measurements of the last two masses. The statesace matrices for this dynamical system are 03 1 I A = M K 1 M C, B = , C = [ ], D = 0 The state vector is made u of ositions and velocities of the three masses in the following order, x x, x, x, x«, x«, x«t, and the mass, stiffness, and daming matrices are = [ ] m1 0 0 M = 0 m m 3 c1 + c2 c2 0 k1 + k2 k2 0, C = c c + c c , K = k k + k k c c k k 3 3 where m1 = 05. Kg, m2 = m3 = 1Kg, k1 = k2 = k3 = 10 N m, c1 = c2 = c3 = 035. Nsec m. The samling interval is 0.1 sec. The system is excited by random inut, and the outut data is corruted by significant rocess and measurement noise. In modal coordinates, each of the modal states is corruted by about 3-5% rocess noise (measured in terms of standard deviation ratios of noise to signal), and the oututs are corruted by about 15% measurement noise. Because this is a simulation, we can actually comute the noise 16
17 statistics given in terms of their covariances. From exact knowledge of the system and the comuted rocess and measurement noise covariances, a Kalman filter is designed. The Kalman filter reresents the best or otimal estimation that can be achieved for the given system with the known noise statistics. Next we use the above set of noise corruted inut-outut data to identify a state estimator with the rocedure described in this aer, and this is done without knowledge of the system and without knowledge of the embedded rocess and measurement noise statistics. As mentioned in the introduction, ARMarkov models are effective in caturing the true or effective system order. Order determination is achieved by examining the quality of the identified state sace model in reroducing the identification data for various model order selection by Hankel singular value truncation. This is shown in Fig. 1 which indicates that the system order is six, which is indeed the case. Figure 2 shows the actual measured (noise-corruted) oututs together with an overlay of the results of the otimal Kalman filter estimation and the identified 6-th order state estimator. Let us examine the first outut. The jagged curve is the measured noisecorruted outut. The smooth curve reresents the otimal filtering by the Kalman filter. Note that result obtained with the identified state estimator closely follows the Kalman filter result. The same attern is observed for the second outut. Figure 3 shows a comarison of the residuals itself, for each of the two oututs. Recall that the Kalman filter results are derived with exact knowledge of the system and noise statistics, whereas 17
18 the identified state estimator is derived from inut-outut data alone. Figure 4 shows the auto-correlation of the Kalman filter residual and our identified state estimator residual. In addition to comaring filtered oututs, we comare the filtered modal states and this is shown in Fig. 5. Since the Kalman filter minimizes the state estimation error, it is interesting to see how the identified state estimator comares to this otimal result. Keeing in mind that in the resence of noise, identification can never erfectly extract the system model and noise statistics with finite data records but with increasing data length and, imrovement in the identification should be exected. Recall that increasing is beneficial because it hels making the residuals more and more uncorrelated with identification data as shown in the theoretical section. Indeed, it is shown in Table 1 that the norm of the state estimation error (and of outut rediction error) of the identified state estimator aroaches that of the Kalman filter with increasing data length and. These illustrations indicate that the roosed state estimator identified from inut-outut data does indeed aroach the otimal Kalman filter designed with erfect knowledge of the system and erfect knowledge of rocess and measurement noise statistics. Conclusions In this aer we have shown how ARMarkov models can be used to identify a state estimator from inut-outut data. This work extends revious develoment of ARMarkov models for system identification and control alications. Being able to identify a state estimator from data is significant in view of the standard aroach of 18
19 designing such a state estimator from an assumed model of the system and the noise statistics. A key ingredient that makes this ossible is the recent derivation of ARMarkov models through the use of an interaction matrix. In our revious work, the role of the interaction matrix is to justify the existence of various inut-outut models and to establish various relationshi among the identified coefficients. For control and system identification roblems, there is no need to recover the interaction matrix itself. In this work, we roceed one ste further by actually recovering the interaction matrix, and exlaining it in the context of a state estimator. This state estimator is not in the standard form of a Luenberger observer or a Kalman filter but in this new form, this state estimator is easy to design from a known model, and if the model is not known, it can be identified from inut-outut data. References 1. Hyland, D.C., Adative Neural Control (ANC) Architecture - a Tutorial, Proceedings of the Industry, Government, and University Forum on Adative Neural Control for Aerosace Structural Systems, Harris Cor., Melbourne, FL, Akers, J.C., and Bernstein, D.S., ARMARKOV Least-Squares Identification, Proceedings of the American Control Conference, Albuquerque, NM, 1997,
20 3. Akers, J.C., and Bernstein, D.S., Time-Domain Identification Using ARMARKOV/Toelitz Models, Proceedings of the American Control Conference, Albuquerque, NM, 1997, Lim, R.K., Phan, M.Q., and Longman, R.W., State-Sace System Identification with Identified Hankel Matrix, Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3045, Princeton University, Set Juang, J.-N., Phan, M., Horta, L.G., and Longman, R.W., Identification of Observer/Kalman Filter Markov Parameters: Theory and Exeriments, Journal of Guidance, Control, and Dynamics, Vol. 16, No. 2, 1993, Juang, J.-N., Alied System Identification, Prentice-Hall, Englewood Cliffs, NJ, 1994, Phan, M.Q., Lim, R.K., and Longman, R.W., Unifying Inut-Outut and State- Sace Persectives of Predictive Control, Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3044, Princeton University, Set Juang, J.-N., and Paa, R.S., An Eigensystem Realization Algorithm for Model Parameter Identification and Model Reduction, Journal of Guidance, Control, and Dynamics, Vol. 8, No.5, 1985, Juang, J.-N., and Lew, J.-S., Integration of System Identification and Robust Controller Designs for Flexible Structures in Sace, Proceedings of the AIAA Guidance, Control, and Navigation Conference, Portland, OR, 1990,
21 1.8 Prediction error Selected state sace model order Fig. 1. Order determination 21
22 0.05 ARMarkov Kalman Measured Outut ARMarkov Kalman Measured Outut Time (sec.) Fig 2. Measured and estimated oututs by ARMarkov and Kalman filter. 22
23 Residual outut ARMarkov Kalman Residual outut Time (sec.) Fig. 3. ARMarkov and Kalman filter outut residuals. 23
24 Auto-correlation (residual 1) Auto-correlation (residual 2) ARMarkov Kalman Number of time shifts Fig. 4. Auto-correlation of outut residuals. 24
25 4 ARMarkov Kalman Actual State State 4 0 State Time (sec.) Fig. 5. True and estimated modal states by ARMarkov and Kalman filter. 25
26 Data Length State Residual ARMarkov Kalman Outut Residual ARMarkov Kalman Table 1. Comarison of state and outut residuals by ARMarkov and Kalman filter. 26
Observer/Kalman Filter Time Varying System Identification
Observer/Kalman Filter Time Varying System Identification Manoranjan Majji Texas A&M University, College Station, Texas, USA Jer-Nan Juang 2 National Cheng Kung University, Tainan, Taiwan and John L. Junins
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationA Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression
Journal of Modern Alied Statistical Methods Volume Issue Article 7 --03 A Comarison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Ghadban Khalaf King Khalid University, Saudi
More informationDistributed Rule-Based Inference in the Presence of Redundant Information
istribution Statement : roved for ublic release; distribution is unlimited. istributed Rule-ased Inference in the Presence of Redundant Information June 8, 004 William J. Farrell III Lockheed Martin dvanced
More informationGeneration of Linear Models using Simulation Results
4. IMACS-Symosium MATHMOD, Wien, 5..003,. 436-443 Generation of Linear Models using Simulation Results Georg Otte, Sven Reitz, Joachim Haase Fraunhofer Institute for Integrated Circuits, Branch Lab Design
More informationAn All-Interaction Matrix Approach to Linear and Bilinear System Identification
An All-Interaction Matrix Approach to Linear and Bilinear System Identification Minh Q. Phan, Francesco Vicario, Richard W. Longman, and Raimondo Betti Abstract This paper is a brief introduction to the
More informationUncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning
TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More informationA Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition
A Qualitative Event-based Aroach to Multile Fault Diagnosis in Continuous Systems using Structural Model Decomosition Matthew J. Daigle a,,, Anibal Bregon b,, Xenofon Koutsoukos c, Gautam Biswas c, Belarmino
More informationOn split sample and randomized confidence intervals for binomial proportions
On slit samle and randomized confidence intervals for binomial roortions Måns Thulin Deartment of Mathematics, Usala University arxiv:1402.6536v1 [stat.me] 26 Feb 2014 Abstract Slit samle methods have
More informationNumerical Linear Algebra
Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and
More informationarxiv: v1 [physics.data-an] 26 Oct 2012
Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch
More informationRadial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationScaling Multiple Point Statistics for Non-Stationary Geostatistical Modeling
Scaling Multile Point Statistics or Non-Stationary Geostatistical Modeling Julián M. Ortiz, Steven Lyster and Clayton V. Deutsch Centre or Comutational Geostatistics Deartment o Civil & Environmental Engineering
More informationProbability Estimates for Multi-class Classification by Pairwise Coupling
Probability Estimates for Multi-class Classification by Pairwise Couling Ting-Fan Wu Chih-Jen Lin Deartment of Comuter Science National Taiwan University Taiei 06, Taiwan Ruby C. Weng Deartment of Statistics
More informationGeneralized Coiflets: A New Family of Orthonormal Wavelets
Generalized Coiflets A New Family of Orthonormal Wavelets Dong Wei, Alan C Bovik, and Brian L Evans Laboratory for Image and Video Engineering Deartment of Electrical and Comuter Engineering The University
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationResearch Article Controllability of Linear Discrete-Time Systems with Both Delayed States and Delayed Inputs
Abstract and Alied Analysis Volume 203 Article ID 97546 5 ages htt://dxdoiorg/055/203/97546 Research Article Controllability of Linear Discrete-Time Systems with Both Delayed States and Delayed Inuts Hong
More informationAbstract. 1 Introduction
Relationship between state-spaceand input-output models via observer Markov parameters M.Q. Phan,' R. W. Longman* ^Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, NJ
More informationPaper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation
Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional
More informationUse of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek
Use of Transformations and the Reeated Statement in PROC GLM in SAS Ed Stanek Introduction We describe how the Reeated Statement in PROC GLM in SAS transforms the data to rovide tests of hyotheses of interest.
More informationKeywords: pile, liquefaction, lateral spreading, analysis ABSTRACT
Key arameters in seudo-static analysis of iles in liquefying sand Misko Cubrinovski Deartment of Civil Engineering, University of Canterbury, Christchurch 814, New Zealand Keywords: ile, liquefaction,
More informationMultiplicative group law on the folium of Descartes
Multilicative grou law on the folium of Descartes Steluţa Pricoie and Constantin Udrişte Abstract. The folium of Descartes is still studied and understood today. Not only did it rovide for the roof of
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationDETERMINISTIC STOCHASTIC SUBSPACE IDENTIFICATION FOR BRIDGES
DEERMINISIC SOCHASIC SBSPACE IDENIFICAION FOR BRIDGES H. hai, V. DeBrunner, L. S. DeBrunner, J. P. Havlicek 2, K. Mish 3, K. Ford 2, A. Medda Deartment o Electrical and Comuter Engineering Florida State
More informationDETC2003/DAC AN EFFICIENT ALGORITHM FOR CONSTRUCTING OPTIMAL DESIGN OF COMPUTER EXPERIMENTS
Proceedings of DETC 03 ASME 003 Design Engineering Technical Conferences and Comuters and Information in Engineering Conference Chicago, Illinois USA, Setember -6, 003 DETC003/DAC-48760 AN EFFICIENT ALGORITHM
More informationCHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules
CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is
More informationThe non-stochastic multi-armed bandit problem
Submitted for journal ublication. The non-stochastic multi-armed bandit roblem Peter Auer Institute for Theoretical Comuter Science Graz University of Technology A-8010 Graz (Austria) auer@igi.tu-graz.ac.at
More informationGeneralized analysis method of engine suspensions based on bond graph modeling and feedback control theory
Generalized analysis method of ine susensions based on bond grah modeling and feedback ntrol theory P.Y. RICHARD *, C. RAMU-ERMENT **, J. BUION *, X. MOREAU **, M. LE FOL *** *Equie Automatique des ystèmes
More informationGeneral Linear Model Introduction, Classes of Linear models and Estimation
Stat 740 General Linear Model Introduction, Classes of Linear models and Estimation An aim of scientific enquiry: To describe or to discover relationshis among events (variables) in the controlled (laboratory)
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationNotes on Instrumental Variables Methods
Notes on Instrumental Variables Methods Michele Pellizzari IGIER-Bocconi, IZA and frdb 1 The Instrumental Variable Estimator Instrumental variable estimation is the classical solution to the roblem of
More informationA Recursive Block Incomplete Factorization. Preconditioner for Adaptive Filtering Problem
Alied Mathematical Sciences, Vol. 7, 03, no. 63, 3-3 HIKARI Ltd, www.m-hiari.com A Recursive Bloc Incomlete Factorization Preconditioner for Adative Filtering Problem Shazia Javed School of Mathematical
More informationFrequency-Weighted Robust Fault Reconstruction Using a Sliding Mode Observer
Frequency-Weighted Robust Fault Reconstruction Using a Sliding Mode Observer C.P. an + F. Crusca # M. Aldeen * + School of Engineering, Monash University Malaysia, 2 Jalan Kolej, Bandar Sunway, 4650 Petaling,
More informationA Bound on the Error of Cross Validation Using the Approximation and Estimation Rates, with Consequences for the Training-Test Split
A Bound on the Error of Cross Validation Using the Aroximation and Estimation Rates, with Consequences for the Training-Test Slit Michael Kearns AT&T Bell Laboratories Murray Hill, NJ 7974 mkearns@research.att.com
More informationUsing the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process
Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process P. Mantalos a1, K. Mattheou b, A. Karagrigoriou b a.deartment of Statistics University of Lund
More information2-D Analysis for Iterative Learning Controller for Discrete-Time Systems With Variable Initial Conditions Yong FANG 1, and Tommy W. S.
-D Analysis for Iterative Learning Controller for Discrete-ime Systems With Variable Initial Conditions Yong FANG, and ommy W. S. Chow Abstract In this aer, an iterative learning controller alying to linear
More informationDeveloping A Deterioration Probabilistic Model for Rail Wear
International Journal of Traffic and Transortation Engineering 2012, 1(2): 13-18 DOI: 10.5923/j.ijtte.20120102.02 Develoing A Deterioration Probabilistic Model for Rail Wear Jabbar-Ali Zakeri *, Shahrbanoo
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationDEPARTMENT OF ECONOMICS ISSN DISCUSSION PAPER 20/07 TWO NEW EXPONENTIAL FAMILIES OF LORENZ CURVES
DEPARTMENT OF ECONOMICS ISSN 1441-549 DISCUSSION PAPER /7 TWO NEW EXPONENTIAL FAMILIES OF LORENZ CURVES ZuXiang Wang * & Russell Smyth ABSTRACT We resent two new Lorenz curve families by using the basic
More informationSTABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS
STABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS Massimo Vaccarini Sauro Longhi M. Reza Katebi D.I.I.G.A., Università Politecnica delle Marche, Ancona, Italy
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationA SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites
More informationDesign of NARMA L-2 Control of Nonlinear Inverted Pendulum
International Research Journal of Alied and Basic Sciences 016 Available online at www.irjabs.com ISSN 51-838X / Vol, 10 (6): 679-684 Science Exlorer Publications Design of NARMA L- Control of Nonlinear
More informationA generalization of Amdahl's law and relative conditions of parallelism
A generalization of Amdahl's law and relative conditions of arallelism Author: Gianluca Argentini, New Technologies and Models, Riello Grou, Legnago (VR), Italy. E-mail: gianluca.argentini@riellogrou.com
More informationCSE 599d - Quantum Computing When Quantum Computers Fall Apart
CSE 599d - Quantum Comuting When Quantum Comuters Fall Aart Dave Bacon Deartment of Comuter Science & Engineering, University of Washington In this lecture we are going to begin discussing what haens to
More informationON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES CONTROL PROBLEMS
Electronic Transactions on Numerical Analysis. Volume 44,. 53 72, 25. Coyright c 25,. ISSN 68 963. ETNA ON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES
More informationTowards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK
Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)
More informationMULTI-CHANNEL PARAMETRIC ESTIMATOR FAST BLOCK MATRIX INVERSES
MULTI-CANNEL ARAMETRIC ESTIMATOR FAST BLOCK MATRIX INVERSES S Lawrence Marle Jr School of Electrical Engineering and Comuter Science Oregon State University Corvallis, OR 97331 Marle@eecsoregonstateedu
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationAutomatic Generation and Integration of Equations of Motion for Linked Mechanical Systems
Automatic Generation and Integration of Equations of Motion for Linked Mechanical Systems D. Todd Griffith a, John L. Junkins a, and James D. Turner b a Deartment of Aerosace Engineering, Texas A&M University,
More informationElliptic Curves and Cryptography
Ellitic Curves and Crytograhy Background in Ellitic Curves We'll now turn to the fascinating theory of ellitic curves. For simlicity, we'll restrict our discussion to ellitic curves over Z, where is a
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationResearch of power plant parameter based on the Principal Component Analysis method
Research of ower lant arameter based on the Princial Comonent Analysis method Yang Yang *a, Di Zhang b a b School of Engineering, Bohai University, Liaoning Jinzhou, 3; Liaoning Datang international Jinzhou
More informationSolved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.
Solved Problems Solved Problems P Solve the three simle classification roblems shown in Figure P by drawing a decision boundary Find weight and bias values that result in single-neuron ercetrons with the
More informationOn Doob s Maximal Inequality for Brownian Motion
Stochastic Process. Al. Vol. 69, No., 997, (-5) Research Reort No. 337, 995, Det. Theoret. Statist. Aarhus On Doob s Maximal Inequality for Brownian Motion S. E. GRAVERSEN and G. PESKIR If B = (B t ) t
More informationVIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES
Journal of Sound and Vibration (998) 22(5), 78 85 VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES Acoustics and Dynamics Laboratory, Deartment of Mechanical Engineering, The
More informationObservers for Bilinear State-Space Models by Interaction Matrices
Observers for Bilinear State-Space Models by Interaction Matrices Minh Q. Phan, Francesco Vicario, Richard W. Longman, and Raimondo Betti Abstract This paper formulates a bilinear observer for a bilinear
More informationMetrics Performance Evaluation: Application to Face Recognition
Metrics Performance Evaluation: Alication to Face Recognition Naser Zaeri, Abeer AlSadeq, and Abdallah Cherri Electrical Engineering Det., Kuwait University, P.O. Box 5969, Safat 6, Kuwait {zaery, abeer,
More informationSystem Reliability Estimation and Confidence Regions from Subsystem and Full System Tests
009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract
More informationy p 2 p 1 y p Flexible System p 2 y p2 p1 y p u, y 1 -u, y 2 Component Breakdown z 1 w 1 1 y p 1 P 1 P 2 w 2 z 2
ROBUSTNESS OF FLEXIBLE SYSTEMS WITH COMPONENT-LEVEL UNCERTAINTIES Peiman G. Maghami Λ NASA Langley Research Center, Hamton, VA 368 Robustness of flexible systems in the resence of model uncertainties at
More informationPeriod-two cycles in a feedforward layered neural network model with symmetric sequence processing
PHYSICAL REVIEW E 75, 4197 27 Period-two cycles in a feedforward layered neural network model with symmetric sequence rocessing F. L. Metz and W. K. Theumann Instituto de Física, Universidade Federal do
More informationSession 5: Review of Classical Astrodynamics
Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve
More informationOn a Markov Game with Incomplete Information
On a Markov Game with Incomlete Information Johannes Hörner, Dinah Rosenberg y, Eilon Solan z and Nicolas Vieille x{ January 24, 26 Abstract We consider an examle of a Markov game with lack of information
More informationThe Recursive Fitting of Multivariate. Complex Subset ARX Models
lied Mathematical Sciences, Vol. 1, 2007, no. 23, 1129-1143 The Recursive Fitting of Multivariate Comlex Subset RX Models Jack Penm School of Finance and lied Statistics NU College of Business & conomics
More information216 S. Chandrasearan and I.C.F. Isen Our results dier from those of Sun [14] in two asects: we assume that comuted eigenvalues or singular values are
Numer. Math. 68: 215{223 (1994) Numerische Mathemati c Sringer-Verlag 1994 Electronic Edition Bacward errors for eigenvalue and singular value decomositions? S. Chandrasearan??, I.C.F. Isen??? Deartment
More informationIntroduction to MVC. least common denominator of all non-identical-zero minors of all order of G(s). Example: The minor of order 2: 1 2 ( s 1)
Introduction to MVC Definition---Proerness and strictly roerness A system G(s) is roer if all its elements { gij ( s)} are roer, and strictly roer if all its elements are strictly roer. Definition---Causal
More informationGenetic Algorithms, Selection Schemes, and the Varying Eects of Noise. IlliGAL Report No November Department of General Engineering
Genetic Algorithms, Selection Schemes, and the Varying Eects of Noise Brad L. Miller Det. of Comuter Science University of Illinois at Urbana-Chamaign David E. Goldberg Det. of General Engineering University
More informationDr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur
Analysis of Variance and Design of Exeriment-I MODULE II LECTURE -4 GENERAL LINEAR HPOTHESIS AND ANALSIS OF VARIANCE Dr. Shalabh Deartment of Mathematics and Statistics Indian Institute of Technology Kanur
More informationAn Improved Calibration Method for a Chopped Pyrgeometer
96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum
More informationRUN-TO-RUN CONTROL AND PERFORMANCE MONITORING OF OVERLAY IN SEMICONDUCTOR MANUFACTURING. 3 Department of Chemical Engineering
Coyright 2002 IFAC 15th Triennial World Congress, Barcelona, Sain RUN-TO-RUN CONTROL AND PERFORMANCE MONITORING OF OVERLAY IN SEMICONDUCTOR MANUFACTURING C.A. Bode 1, B.S. Ko 2, and T.F. Edgar 3 1 Advanced
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationAlgorithms for Air Traffic Flow Management under Stochastic Environments
Algorithms for Air Traffic Flow Management under Stochastic Environments Arnab Nilim and Laurent El Ghaoui Abstract A major ortion of the delay in the Air Traffic Management Systems (ATMS) in US arises
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More informationYixi Shi. Jose Blanchet. IEOR Department Columbia University New York, NY 10027, USA. IEOR Department Columbia University New York, NY 10027, USA
Proceedings of the 2011 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelsach, K. P. White, and M. Fu, eds. EFFICIENT RARE EVENT SIMULATION FOR HEAVY-TAILED SYSTEMS VIA CROSS ENTROPY Jose
More informationPROFIT MAXIMIZATION. π = p y Σ n i=1 w i x i (2)
PROFIT MAXIMIZATION DEFINITION OF A NEOCLASSICAL FIRM A neoclassical firm is an organization that controls the transformation of inuts (resources it owns or urchases into oututs or roducts (valued roducts
More informationElliptic Curves Spring 2015 Problem Set #1 Due: 02/13/2015
18.783 Ellitic Curves Sring 2015 Problem Set #1 Due: 02/13/2015 Descrition These roblems are related to the material covered in Lectures 1-2. Some of them require the use of Sage, and you will need to
More informationChater Matrix Norms and Singular Value Decomosition Introduction In this lecture, we introduce the notion of a norm for matrices The singular value de
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A Dahleh George Verghese Deartment of Electrical Engineering and Comuter Science Massachuasetts Institute of Technology c Chater Matrix Norms
More informationCharacterizing the Behavior of a Probabilistic CMOS Switch Through Analytical Models and Its Verification Through Simulations
Characterizing the Behavior of a Probabilistic CMOS Switch Through Analytical Models and Its Verification Through Simulations PINAR KORKMAZ, BILGE E. S. AKGUL and KRISHNA V. PALEM Georgia Institute of
More informationA New Asymmetric Interaction Ridge (AIR) Regression Method
A New Asymmetric Interaction Ridge (AIR) Regression Method by Kristofer Månsson, Ghazi Shukur, and Pär Sölander The Swedish Retail Institute, HUI Research, Stockholm, Sweden. Deartment of Economics and
More informationOn Fractional Predictive PID Controller Design Method Emmanuel Edet*. Reza Katebi.**
On Fractional Predictive PID Controller Design Method Emmanuel Edet*. Reza Katebi.** * echnology and Innovation Centre, Level 4, Deartment of Electronic and Electrical Engineering, University of Strathclyde,
More informationElementary Analysis in Q p
Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some
More informationBrownian Motion and Random Prime Factorization
Brownian Motion and Random Prime Factorization Kendrick Tang June 4, 202 Contents Introduction 2 2 Brownian Motion 2 2. Develoing Brownian Motion.................... 2 2.. Measure Saces and Borel Sigma-Algebras.........
More informationChapter 3. GMM: Selected Topics
Chater 3. GMM: Selected oics Contents Otimal Instruments. he issue of interest..............................2 Otimal Instruments under the i:i:d: assumtion..............2. he basic result............................2.2
More informationUncertainty Modeling with Interval Type-2 Fuzzy Logic Systems in Mobile Robotics
Uncertainty Modeling with Interval Tye-2 Fuzzy Logic Systems in Mobile Robotics Ondrej Linda, Student Member, IEEE, Milos Manic, Senior Member, IEEE bstract Interval Tye-2 Fuzzy Logic Systems (IT2 FLSs)
More informationShadow Computing: An Energy-Aware Fault Tolerant Computing Model
Shadow Comuting: An Energy-Aware Fault Tolerant Comuting Model Bryan Mills, Taieb Znati, Rami Melhem Deartment of Comuter Science University of Pittsburgh (bmills, znati, melhem)@cs.itt.edu Index Terms
More informationLPC methods are the most widely used in. recognition, speaker recognition and verification
Digital Seech Processing Lecture 3 Linear Predictive Coding (LPC)- Introduction LPC Methods LPC methods are the most widely used in seech coding, seech synthesis, seech recognition, seaker recognition
More informationarxiv:cond-mat/ v2 25 Sep 2002
Energy fluctuations at the multicritical oint in two-dimensional sin glasses arxiv:cond-mat/0207694 v2 25 Se 2002 1. Introduction Hidetoshi Nishimori, Cyril Falvo and Yukiyasu Ozeki Deartment of Physics,
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More informationImproved Capacity Bounds for the Binary Energy Harvesting Channel
Imroved Caacity Bounds for the Binary Energy Harvesting Channel Kaya Tutuncuoglu 1, Omur Ozel 2, Aylin Yener 1, and Sennur Ulukus 2 1 Deartment of Electrical Engineering, The Pennsylvania State University,
More informationEnsemble Forecasting the Number of New Car Registrations
Ensemble Forecasting the Number of New Car Registrations SJOERT FLEURKE Radiocommunications Agency Netherlands Emmasingel 1, 9700 AL, Groningen THE NETHERLANDS sjoert.fleurke@agentschatelecom.nl htt://www.agentschatelecom.nl
More informationNeural network models for river flow forecasting
Neural network models for river flow forecasting Nguyen Tan Danh Huynh Ngoc Phien* and Ashim Das Guta Asian Institute of Technology PO Box 4 KlongLuang 220 Pathumthani Thailand Abstract In this study back
More informationCSC165H, Mathematical expression and reasoning for computer science week 12
CSC165H, Mathematical exression and reasoning for comuter science week 1 nd December 005 Gary Baumgartner and Danny Hea hea@cs.toronto.edu SF4306A 416-978-5899 htt//www.cs.toronto.edu/~hea/165/s005/index.shtml
More informationUncorrelated Multilinear Discriminant Analysis with Regularization and Aggregation for Tensor Object Recognition
Uncorrelated Multilinear Discriminant Analysis with Regularization and Aggregation for Tensor Object Recognition Haiing Lu, K.N. Plataniotis and A.N. Venetsanooulos The Edward S. Rogers Sr. Deartment of
More informationPositivity, local smoothing and Harnack inequalities for very fast diffusion equations
Positivity, local smoothing and Harnack inequalities for very fast diffusion equations Dedicated to Luis Caffarelli for his ucoming 60 th birthday Matteo Bonforte a, b and Juan Luis Vázquez a, c Abstract
More informationIterative Methods for Designing Orthogonal and Biorthogonal Two-channel FIR Filter Banks with Regularities
R. Bregović and T. Saramäi, Iterative methods for designing orthogonal and biorthogonal two-channel FIR filter bans with regularities, Proc. Of Int. Worsho on Sectral Transforms and Logic Design for Future
More informationAggregate Prediction With. the Aggregation Bias
100 Aggregate Prediction With Disaggregate Models: Behavior of the Aggregation Bias Uzi Landau, Transortation Research nstitute, Technion-srael nstitute of Technology, Haifa Disaggregate travel demand
More informationx and y suer from two tyes of additive noise [], [3] Uncertainties e x, e y, where the only rior knowledge is their boundedness and zero mean Gaussian
A New Estimator for Mixed Stochastic and Set Theoretic Uncertainty Models Alied to Mobile Robot Localization Uwe D. Hanebeck Joachim Horn Institute of Automatic Control Engineering Siemens AG, Cororate
More information