Trajectory Estimation for Tactical Ballistic Missiles in Terminal Phase Using On-line Input Estimator
|
|
- Angel Palmer
- 5 years ago
- Views:
Transcription
1 Proc. Natl. Sci. Counc. ROC(A) Vol. 23, No. 5, pp Trajectory Estimation for Tactical Ballistic Missiles in Terminal Phase Using On-line Input Estimator SOU-CHEN LEE, YU-CHAO HUANG, AND CHENG-YU LIU Department of System Engineering Chung Cheng Institute of Technology Taoyuan, Taiwan, R.O.C. (Receied Noember 20, 1998; Accepted March 11, 1999) ABSTRACT This work presents a noel on-line estimation approach to estimation of the trajectory of a tactical ballistic missile with simple maneuering in the terminal phase. This proposed approach incorporates the extended Kalman filter and an innoatie recursie input estimator into the hypothetical testing scheme. Relatiely simple steering is accomplished by the ballistic missile and is considered to be an abrupt lateral acceleration in the dynamic equation. Also formulated herein is a regression equation between the obsered alue of the residual sequence and the theoretical residual sequence of the extended Kalman filter with no input. In addition, a recursie least-squares estimator is proided not only to extract the magnitude of the input, but also to proide a testing criterion for detecting the onset and presence of the input. Numerical simulation has demonstrated the superior capabilities of this input approach. Moreoer, trajectory estimates of the two ariations of the ballistic coefficient can erify the robustness and computational stability within ery short time. Key Words: anti tactical ballistic missile, trajectory estimation, extended Kalman filter, input estimation I. Introduction The effectieness of the Tactical Ballistic Missile (TBM) during the Persian Gulf War had a startling effect globally. Expertise deeloped during the War led to the aailability of the Anti Tactical Ballistic Missile (ATBM) for use in conentional air defense missiles and, ultimately, to upgraded performance. The modern TBM was designed to maneuer in the terminal phase to escape tracking and interception. Therefore, on-line trajectory estimation of the TBM in the regular fire control sequence for the ATBM is highly desired in radar tracking, direct hit guidance, and early warning systems. Normally, trajectory estimation of a flying ehicle is inestigated primarily in post analysis to identify states and key parameters in aailable flight data measured using radar, satellites, and on-board sensors. For a conentional ATBM air defense missile, engagement is limited in the re-entry phase since radar is the only measurement instrument used to sense the TBM. Thus, reconstructing the states and parameters of a reentry ehicle is relatiely difficult, particularly during maneuers. Trajectory estimation of a maneuering TBM has receied little attention. Chang et al. (1977) introduced a maneuering re-entry ehicle filter with a newly defined augmented state ector. In the filter, the augmented state ector contains position, elocity, and corresponding parameters for drag and maneuer forces; the extended Kalman filter (EKF) has been proposed the task of estimation. The re-entry ehicle filter performance, howeer, is degraded if maneuering forces are absent. The related parameter estimation also relies heaily on the model and is assumed to be a constant Gaussian Marko process in the state equation. Abutaleb (1985) presented a Pontryagin filter based on the Pontryagin minimum principle, in which the unknown maneuering forces are treated as a control used to drie a ehicle s dynamic so that it will follow the noisy measured trajectory. Neertheless, the Pontryagin filter is sensitie to parameter ariation and initial conditions. The input estimation technique has been successfully applied in many engineering areas, including the tracking problem, initial leeling of strapdown inertial naigation systems, and inerse heat conduction problems. In TBM, maneuering is considered to be an extra lateral acceleration input which rapidly changes a missile s position, elocity, and heading. Therefore, deeloping an input estimation approach to identify this lateral acceleration is highly desired. In this study, we present a noel adaptie filter to proide accurate and fast trajectory estimation with low sensitiity to 644
2 Trajectory Estimation for TB Missiles an unknown key parameter, i.e., the ballistic coefficient. The proposed adaptie filter consists of the EKF and an innoatie input estimator with the input detection criterion. The input estimator proposed herein can continuously estimate all possible lateral accelerations and substitute the estimated results into the EKF in order to identify the TBM s position and elocity if the criterion is satisfied. The impact of the ballistic coefficient uncertainty can be treated as an another output drien by extra accelerations, which are accompanied by lateral accelerations that are to be estimated. Simulation results indicate that the proposed filter leads to much improement compared with the EKF without an input estimator and the augmented EKF, whose state ector contains the position, elocity and lateral acceleration. Simulation results further demonstrate that the proposed algorithm proides robustness to the ballistic coefficient. Due to the nature of data reconstruction, recommended trajectory estimation can proide further insight into TBM system behaior and possibly facilitate a more general model design and analysis inoling the ATBM strategy. The rest of this paper is organized as follows. The next section describes a maneuering TBM model in a 2 dimensional plane, where an equation of motion is introduced and the extended Kalman filter is applied to the estimation problem. In Section III, we present adaptie filtering with a noel input estimation scheme for lateral acceleration of TBM. This section also discusses the related test criterion for input detection. The extended Kalman filter is also combined with the input estimator. Next, Section IV presents two simulation examples which demonstrate the effectieness of the proposed filter and its robustness to the ballistic coefficient. Concluding remarks are finally made in the last section. II. Dynamic Model Consider a 2 dimensional TBM in the re-entry phase and in a flat, nonrotating earth with constant graity model, as illustrated in Fig. 1. The distance that the TBM traels in this phase is shorter than the distance traeled before re-entry. The setting for the TBM is a point mass and constant weight along a ballistic trajectory, where three types of significant forces act on the TBM, among which graity and drag are two types. The adanced tactical ballistic missile usually steers at the right time in the terminal plase so as to escape radar tracking and missile interception. This steering may rapidly change the heading, elocity and acceleration of the TBM. The third force is introduced in a plane which is perpendicular to the elocity ector if the TBM performs a maneuer. The force is induced by a maneuer and may be represented by components along X R and Y R, which are relatie to the acceleration terms in the equations of motion. The TBM equations with a maneuer in radar coordinates (O R,X R,Y R ) centered at a radar site can be represented as x = Dg W cosγ + a x = 2 2β gcosγ + a x, (1) y = Dg W sinγ g + a y = 2 2β gsinγ g + a y, (2) where D and W denote the drag and weight of the missile, respectiely, represents air density and, a x and a y are lateral accelerations generated by a maneuer along X R and Y R, respectiely. The flight path angle γ and ballistic coefficient β are defined as γ = tan 1 ( y x ), β = Fig. 1. Tactical ballistic missile geometry. W SC D0, where S and C D0 represent the reference area and zerolift drag coefficient, respectiely. The air density is originally a function of altitude and should be considered in Eqs. (1) and (2) because the altitude dramatically changes in flight at speeds oer 100 km. The conentionally used approximation model for air density related to altitude is = e y/30000 =0.0034e y/22000 Let the states ector be for y<30000 ft, for y ft. X=[x 1 x 2 x 3 x 4 ] T =[x y x y ] T. (3) 645
3 S.C. Lee et al. The nonlinear state equation can be deried as X =F(X,t)+ϕ u+iζ. (4) In Eq. (4), ζ denotes the process noise a ector with a ariance of Q, I stands for the identity matrix, u is the lateral acceleration term of the TBM, F(X, t)= x 3 x 4, ϕ = 2β (x x 2 4 )gcosγ 2β (x x 2 4 )gsinγ g u=[0 0 a x a y ] T =[0 0 u 3 u 4 ] T For a conentional air defense missile system, ground radar is the major instrument used to detect a TBM; it proides the position, elocity, and een the acceleration of tracked targets. The measurement equation is Z=HX+ε, (5) where ε represents the measurement noise ector and H denotes the identity matrix. Equations (4) and (5) form the dynamic equations for the TBM with a maneuer after re-entry. Once all the states at a specific time instance are precisely known, the trajectory is reconstructed using a fixed-point smoother and n-step ahead predictor. The EKF is a conentionally used state estimator for nonlinear dynamic equations and is applied straightforwardedly herein. Let ε denote white and be normally distributed with zero mean and a ariance of R. The predicted and updated state ectors from t=n t to t= (n+1) t under input u n at t=n t are gien by X n+1/n = φ n X n/n + ϕu n, (6) X n+1/n+1 = X n+1/n + K n+1 (Z n+1 HX n+1/n ), (7) where t is the sampling period, Z n+1 denotes measurements at t=(n+1) t, and the transition matrix as φ n = I + F(X, t) X X = X n/n t. A linearization of F(X,t) can be easily obtained, F(X, t) X X = X n/n = f 32 = = 60000β g(x x 4 2 ) cosγ 44000β g(x x 4 2 ) cosγ f 33 = 2β g(2x 3 cosγ x 4 sinγ), f 34 = 2β g(2x 4 cosγ + x 3 sinγ), f 42 = = 60000β g(x x 4 2 ) sinγ 44000β g(x x 4 2 ) sinγ f 43 = 2β g(2x 3 sinγ + x 4 cosγ), f 44 = 2β g(2x 4 sinγ x 3 cosγ) , 0 f 32 f 33 f 34 0 f 42 f 43 f 44 X = X n/n for x 2 <30000 ft, for x ft, for x 2 <30000 ft, for x ft, The Kalman gain K n+l and coariance matrix of X n+1/n and X n+1/n+1, i.e., P n+1/n and P n+1/n+1, respectiely, are K n+1 =P n+1/n H T (HP n+1/n H T +R) 1 (8) P n+1/n =φ n P t n/n φ n +ΓQΓ, (9) P n+1/n+1 =(I K n+l H)P n+l/n (10) and Γ=I t. The unpredictable input term u n in Eq. (6) is the dominant term regarding the accuracy in state estimation. The EKF may conerge well with long time propagation if u n is omitted in Eq. (6). Howeer, a long conergence time is unacceptable for defense considerations, which require a rapid reaction within just a few minutes. Therefore, an algorithm must be deeloped to estimate the input acceleration and achiee a rapid and accurate trajectory estimation. III. Adaptie Filtering with Input Estimation For a situation in which the input term u n in Eq. (6) remains unknown and onset is unpredictable, this section presents a noel adaptie filter scheme which consists of the EKF and a recursie least-squares estimator of input with the test criterion. Adaptie refers to a situation in which the filter sequentially 646
4 Trajectory Estimation for TB Missiles adjusts an input based on measurements. The input estimator attempts to estimate u n using a single radar measurement. Moreoer the estimated input is combined with the EKF if it satisfies the test criterion. 1. Input Estimation Let X n+1/n, X n+1/n+1 denote the predicted and updated states, respectiely, for the EKF with no input at t=(n+l) t. For simplicity, denote X n+1 = X n+1/n+1 and X n+1 = X n+1/n+1, and let M n+1 =(I K n+l H)φ n, N n+1 = (I K n+1 H)ϕ. The updated state can be represented by X n+ l =( n + l M i i = n +1 l 1 Σ l ) X n + ( M n + i ) K n + j Z n + j j =1 + K n + l Z n + l l=1, 2,... (11) Assume that the TBM maneuers with constant lateral acceleration in each axis during kdt t (k+s)dt; then, the abrupt inputs are Then, Allow l 1 Σ Z k + l = H[ ( M k + i ) N k + j u k + j 1 + N k + l u k + l 1 ] j =1 + Z k + l. l l 1 Σ Y k + l = Z k + l H ( M k + i ) N k + j u k + j 1. j =1 Thus, we hae the regression equation where l Y k + l = Φ k + l u k + l 1 + Z k + l, (17) Φ k+l =HN k+l. Therefore, the recursie least-squares estimator of input can be expressed as u = 0 t < k t, t>(k + s) t k, s >0 u k + r k t t (k + s) t r = 0, 1, 2,..., s, (12) where u k + l 1 = u k + l 2 + G k + l (Y k + l Φ k + l u k + l 2 ), (18) where u k+r denotes a constant ector. Similarly, the updated states ector of the EKF formation with input, Eq. (8), can be expressed as X k + l =( k + l M i i = k +1 l 1 Σ l )X k + ( M k + i )(K k + j Z k + j j =1 + N k + j u k + j 1 )+K k + l Z k + l + N k + l u k + l 1. (13) The EKF formations with no input and with an input yield the same results during t k t, so that X k = X k. The difference induced by the abrupt inputs between these two formations during k t t (k+s) t can then be written as X k + l = X k + l X k + l = M k + l X k + l 1 + N k + l u k + l 1. (14) Define the residuals of measurement for the two EKF formations with and without inputs, respectiely, as Z k + l = Z k + l HX k + l, (15) Z k + l = Z k + l HX k + l. (16) Y k + l = Z k + l HM k + l X k + l 1, (19) X k + l 1 = X k + l 1 u = u = M k + l X k + l 2 + N k + l u k + l 2, and the gain G i and ariance of u i and V i are G k+l =V k+l 1 Φ k+l ξ 1, (20) T V k + l 1 = V k + l 2 V k + l 2 Φ k + l T [Φ k + l V k + l 2 Φ k + l + ξ] 1 Φ k + l V k + l 2, (21) where ξ is the ariance of Z. The innoatie measurement Y k + l contains only partial information about Z k+l, leading subsequently to a biased input estimator. Howeer, this biased estimator still functions properly. Note that the estimated input is drien by X k + l and Z k + l. Therefore, it is not only a function of the inputs, but also depends on the estimation errors of the EKF. In fact, X k + l represents the estimated states obtained using the EKF formation with no input, namely, the original filter. Therefore, the actual system must be compared with the simulated system to generate Y k + l during input estimation. 647
5 S.C. Lee et al. In the preious section, we assumed that β was constant, which is not always true. Let u β be the extra acceleration ector induced by iolation of aboe assumption. Then, the state equation can be rewritten as X =F(X,t)+ϕ u+ϕ u β +Iζ =F(X,t)+ϕ U+Iζ, (22) where U=u+u β. Therefore, the augmented input ector U can be estimated in the same manner as described aboe. The partially improper assumption about β may be coered by the proposed filter and is ealuated later by means of simulation. 2. Criterion for Testing In Eq. (12), k and s represent the starting and stopping points of the system input, respectiely. As mentioned earlier, the onset of system input is unpredictable and can not be assumed in adance. Here, we will present a testing method to determine k and s. Assume that ε 1, ε 2,... are independent, identically distributed random ariables, and that ε N(0, R). (23) Since u i, i=3, 4 is biased, according to the central limit theory, we hae u i N(u bi,v ii ) i=3, 4 as t, (24) where u b denotes a certain alue which corresponds to the true input u i and is equal to zero if u i does not exist; in addition, V ii denotes the ariance of u i. Therefore, the hypotheses can be defined to test the existence of an input as H 0 : u bi =0 (absence of input) i=3, 4, (25) H 1 : u bi 0 (existence of input) i=3, 4, (26) and the existence of a normalized statistical ariable as u st = u i V ii i=3, 4. (27) Assume that the randomness of the initial alue of V can be neglected. Then, u st N(0,1) when H 0 is true. Therefore, the test of whether either H 0 or H l holds true becomes Fig. 2. The mechanism of the proposed adaptie filter scheme. u st >t st H 1 is hold, (28) u st t st H 0 is hold, (29) where [ t st,t st ] is the confidence interal, which can be determined by examining the cumulatie normal distribution table for a certain preset confidence leel α. Figure 2 illustrates the mechanism used in the proposed scheme. The estimated input is applied to the EKF when the switch is on, which denotes that H 1 is on hold. Otherwise, the switch is off, indicating that H 0 is on hold. 3. Adaptie Filtering By using the adjusted on-line input, the predicted and the updated states for the adaptie filter at time interal k t t (k+s) t can be obtained as X k+l+1/k+l = φ k + l X k+l/k+l X k+l+1/k+l+1 = X k+l+1/k+l + ϕu k+l 1, (30) + K k+l+1 (Z k+l+1 HX k+l+1/k+l ). (31) The Kalman gain becomes K k+l+1 = P k+l+1/k+l H T (HP k+l+1/k+l H T + R) 1 with the coariance matrix of the adaptie filter at k t t (k+s) t: where P k+l+1/k+l P k+l+1/k+l+1 = P k+l+1/k+l +φ k+l L k+l+1 φ T k+l +ϕv k+l ϕ T = P k+l+1/k+l + P k+l+1/k+l, (32) =(I K k+l+1 H)P k+l+1/k+l, (33) P k+l+1/k+l = φ k+l L k+l+1 φ T k+l + ϕv k+l ϕ T, l Σ l +1 L k+l+1 = ( M k+i 1 ) N k+j V k+j 1 N T T k+j ( M k+i 1 ) j =1 l+1 648
6 Trajectory Estimation for TB Missiles alues. The following simulations are all for a single flight. 1. True Value of β Ealuation T = M k+l L k+l M k+ l. P k+l+1/k+l represents the increment in coariance induced by u i, i=k, k+1,..., k+l. It is easy to proe that the P k+l+1/k+l is positie definite, implying that when u i and i=k, k+1,..., k+l are introduced to reduce the state estimation errors, the coariance of estimated states must increase. For a time beyond the interal t<k t and t>(k+s) t, the estimation states can also be based upon the original EKF. Notably, the initial states and coariance matrices at t>(k+s) t are reinitiated by X k+s/k+s and. P k+s/k+s Fig. 3. Measured position. IV. Computer Simulations This section presents two cases with a confidence leel under α=95% to ealuate the effectieness of the proposed adaptie filter scheme, as shown in Fig. 2. The first case compares the proposed method with the EKF and augmented EKF under the true alue of β to demonstrate the capability of the proposed method. The other case inestigates the robustness of the proposed scheme to β. Consider a scenario corresponding to Eqs. (1) and (2) with an initial altitude of ft, an initial downrange of 1000 ft, an initial elocity of 4500 ft/sec, a re-entry angle of 70, and β=500 1b/ft 2. The lateral accelerations with magnitude 4g and 2g in the X R and Y R axes, respectiely, are applied to the dynamic equations within 20 t 25 sec after re-entry. Figures 3 and 4 display the measured position and elocity generated by Eqs. (4) and (5) with N(0,1) measurement noises, R=0.3I, and Q=0.01I. The first data detected by radar are adopted as the initial Ealuating the adaptie filter with input estimation initially inoles performing 2 dimensional trajectory simulation of the TBM and estimating the trajectory for a gien re-entry course with β= 500 lb/ft 2. This case adopts three methods of trajectory estimation, namely, adaptie filters, the EKF, and the augmented EKF, with acceleration inputs a x and a y taken as two augmented states to be estimated by the EKF. a x and a y can be modeled using a Gauss-Marko process as follows: a x = α x a x + 2α x σ x2 ζ, (34) a y = α y a y + 2α y σ y2 ζ, (35) where 1/α x and 1/α y denote time constants, and σ x 2 and σ y 2 represent the ariance of a x and a y, respectiely. Let the augmented state ector be X =[x 1 x 2 x 3 x 4 x 5 x 6 ] T =[x y x y a x a y ] T ; then, the state equation is gien by where X = F(X,t)+I ζ, (36) Fig. 4. Measured elocity. 649
7 S.C. Lee et al. x 3 x 4 F(X,t)= 2β (x x 4 2 )gcosγ + x 5 2β (x x 4 2 )gsinγ g + x 6 α x α y I I = 2α x σ2 x α y σ2 y. Fig. 6. Detected inputs (upper) and standard deiations (below). The augmented EKF is employed when the state equation is substituted into Eqs. (6) and (7). Figures 5 and 6 show the estimated and detected accelerations and standard deiations (STD) obtained using the augmented EKF and the adaptie filter, respectiely. As mentioned in the preious section, the detected inputs of the adaptie filter are generated from actual inputs and state estimation errors. The larger inputs are required to compensate before the exact inputs are applied. The detected inputs gradually decrease when the lateral accelerations are introduced. The same phenomenon is seen in Fig. 5 for the same reason. Figures 7-9 present comparisons of the resulting estimation errors and STD in position among these Fig. 7. Estimation errors and standard deiation in position without input estimation. Fig. 5. Estimated inputs (upper) and standard deiations (below) obtained using augmented EKF. three methods. As seen in these figures, the proposed filter immediately reduces the estimation errors by more than two orders as compared to the EKF and the augmented EKF. The maximum estimation errors for the EKF are up to ft at t=20 sec and ft at t=21.65 sec in X R and Y R, respectiely. For the augmented EKF, the maximum errors reach 3000 ft and 3000 ft in X R and Y R, respectiely, at t=40 sec. When using the proposed scheme, the errors are bounded by 8.27 ft and 19.5 ft along X R and Y R, respectiely, and rapidly approach a steady state after the inputs are remoed. The same phenomenon appears in Figs , which display the estimation errors and STD in elocity. The maximum errors reach ft/ sec and 5.65 l0 3 ft/sec for the EKF and reach 500 ft/sec and 500 ft/sec for the augmented EKF in X R and 650
8 Trajectory Estimation for TB Missiles Fig. 8. Estimation errors and standard deiation in position obtained using augmented EKF. Fig. 10. Estimation errors and standard deiation in elocity without input estimation. Fig. 9. Estimation errors and standard deiation in position with input estimation. Fig. 11. Estimation errors and standard deiation in elocity obtained using augmented EKF. Y R, respectiely. Howeer, the proposed scheme s errors are limited to 11.2 ft/sec and 17.6 ft/sec and quickly tend toward a steady state. Table 1 compares the estimation errors at impact point, that is, the steady errors, of the adaptie filter with the EKF and the augmented EKF, reealing the much greater accuracy of the proposed scheme. In addition, the augmented EKF requires that α x, α y, σ x, and σ y, which are normally unknown beforehand, be assigned. This difficulty is absent from the proposed method. The contributions and feasibility of the input estimator are, therefore, obious and merit further study. 2. β Variation Ealuation In this simulation study, we present the perfor- Fig. 12. Estimation errors and standard deiation in elocity with input estimation. 651
9 S.C. Lee et al. Table 1. The Steady Errors of All States with β=500 lb/ft 2 States The adaptie The EKF The augmented filter EKF Position in X R 1.57 ft ft 3000 ft Position in Y R 1.03 ft ft 3000 ft Velocity in X R ft/sec ft/sec 325 ft/sec Velocity in Y R ft/sec ft/sec 310 ft/sec Fig. 14. Estimation errors in elocity with β=650 lb/ft 2. Fig. 13. Estimation errors in position with β=650 lb/ft 2. mance of low sensitiity to β of the proposed scheme. As mentioned earlier, β is assumed to be a constant, which is not always true. Two worse scenarios of β ariations of ±30% are inestigated while examining the robustness of the proposed method. Figures show that the estimation errors in position and elocity are under β=650 lb/ft 2 and β=350 lb/ft 2. The proposed method naturally proides better estimation results, roughly more two orders better than those obtained using the original filter. When β increases by 30%, the maximum estimation errors in position of the original filter are ft at t=21.95 sec and ft at t=24.4 sec in X R and Y R, respectiely. In addition, the peak estimation errors of the proposed method in position are 7.53 ft and ft along X R and Y R, respectiely. The errors of the original filter are much larger than those of the proposed method at any time instance, as indicated in Fig. 13. The maximum estimation errors in elocity of the original filter are significantly greater than those of the proposed method in each axis. The original filter s errors in elocity estimation are always larger than those of the proposed method, as depicted in Fig. 14. When β decreases by 30%, the maximum errors of the original filter in position are l ft and ft and in Fig. 15. Estimation errors in elocity with β=350 lb/ft 2. Fig. 16. Estimation errors in elocity with β=350 lb/ft
10 Trajectory Estimation for TB Missiles Table 2. The Steady Errors of All States with β=650 lb/ft 2 Table 3. The Steady Errors of All States with β=350 lb/ft 2 States The proposed scheme The original filter Position in X R 1.71 ft ft Position in Y R 1.59 ft ft Velocity in X R 1.1 ft/sec ft/sec Velocity in Y R 2.12 ft/sec ft/sec States The proposed scheme The original filter Position in X R 1.0 ft ft Position in Y R ft ft Velocity in X R ft/sec ft/sec Velocity in Y R 0.39 ft/sec ft/sec elocity are 5.53 l0 3 ft/sec and 9.84 l0 3 ft/sec along X R and Y R, respectiely. Howeer, the proposed method still produces a satisfactory estimate. The maximum errors are 10 ft and ft in position estimation and are 15 ft/sec and 24.9 ft/sec in elocity estimation along X R and Y R, respectiely. Tables 2 and 3 list the steady alues of errors under ariation of β. From those table, we can infer that the proposed adaptie filter scheme is much lower sensitie to β than is the EKF. V. Conclusions This study has proposed an on-line estimation algorithm composed of the EKF and an input estimator with the test criterion that can achiee accurate and fast trajectory estimation of a maneuering TBM. A recursie least-squares estimator of input estimation has been presented to identify the maneuer forces of the TBM. The method automatically detects the lateral acceleration of the maneuering TBM and also is a more accurate and faster approach than the EKF. Numerical simulation has demonstrated the superior capabilities of this input approach. Moreoer, trajectory estimation for two different alues of the ballistic coefficient has erified the robustness of the proposed scheme in spite of parameter ariation, confirming that this scheme is worth further deelopment. Acknowledgment This work was supported by the National Science Council of the Republic of China under Grant NSC D References Abutaleb, A. S. (1985) Improed trajectory estimation of maneuering reentry ehicles using a nonlinear filter based on the pontryagin minimum principle. IEEE International Radar Conference, pp Arlington, VA, U.S.A. Chang, C. B., R. H. Whiting, and M. Athan (1977) On the state and parameter estimation for maneuering reentry ehicles. IEEE Transactions on Automatic Control, AC-22(1), !"#$"%"#&'!"#$!%&'!"#$%!&'!"#$%&'()*+,-./)*# :;<=>?=@ ABC9+,DE!"#$%&'()*+,-./, :';<=>?@6ABCDBEFGHIJKL!"#$%&'()*+,-./ :;<=>?@ABCDEFGH&456IJK!"#$%&'(#)*+!,-./ *'9:;<3.=>-?@A456BCD!"#$%&'()*+,-./ :; < =>?@A&BCDEFGH9:I!"#$%&'()*+%, 653
Three-dimensional Guidance Law for Formation Flight of UAV
ICCAS25 Three-dimensional Guidance aw for Formation Flight of UAV Byoung-Mun Min * and Min-Jea Tahk ** * Department of Aerospace Engineering KAIST Daejeon Korea (Tel : 82-42-869-3789; E-mail: bmmin@fdcl.kaist.ac.kr)
More informationClarifying Warhead Separation from the Reentry Vehicle Using a Novel Tracking Algorithm
International Journal Clarifying of Warhead Control, Automation, Separation and from Systems, the Reentry ol. Vehicle 4, no. 5, Using pp. 59-538, a Noel racking Octoer Algorithm 006 59 Clarifying Warhead
More informationAstrometric Errors Correlated Strongly Across Multiple SIRTF Images
Astrometric Errors Correlated Strongly Across Multiple SIRTF Images John Fowler 28 March 23 The possibility exists that after pointing transfer has been performed for each BCD (i.e. a calibrated image
More informationProbabilistic Engineering Design
Probabilistic Engineering Design Chapter One Introduction 1 Introduction This chapter discusses the basics of probabilistic engineering design. Its tutorial-style is designed to allow the reader to quickly
More informationDIFFERENTIAL DRAG SPACECRAFT RENDEZVOUS USING AN ADAPTIVE LYAPUNOV CONTROL STRATEGY
(Preprint) IAA-AAS-DyCoSS1-XX-XX DIFFERENIAL DRAG SPACECRAF RENDEZVOUS USING AN ADAPIVE LYAPUNOV CONROL SRAEGY Daid Pérez * and Riccardo Beilacqua I. INRODUCION his paper introduces a noel Lyapuno-based
More informationNotes on Linear Minimum Mean Square Error Estimators
Notes on Linear Minimum Mean Square Error Estimators Ça gatay Candan January, 0 Abstract Some connections between linear minimum mean square error estimators, maximum output SNR filters and the least square
More informationTarget Trajectory Estimation within a Sensor Network
Target Trajectory Estimation within a Sensor Network Adrien Ickowicz IRISA/CNRS, 354, Rennes, J-Pierre Le Cadre, IRISA/CNRS,354, Rennes,France Abstract This paper deals with the estimation of the trajectory
More informationمتلب سایت MatlabSite.com
Vehicle Velocity Estimation based on Data fusion by Kalman Filtering for ABS Melika Amiri Graduate student in electrical engineering, Islamic Azad uniersity, Science and Research ehran, Iran Melika_64@yahoo.com
More informationOptimal Joint Detection and Estimation in Linear Models
Optimal Joint Detection and Estimation in Linear Models Jianshu Chen, Yue Zhao, Andrea Goldsmith, and H. Vincent Poor Abstract The problem of optimal joint detection and estimation in linear models with
More informationNoise constrained least mean absolute third algorithm
Noise constrained least mean absolute third algorithm Sihai GUAN 1 Zhi LI 1 Abstract: he learning speed of an adaptie algorithm can be improed by properly constraining the cost function of the adaptie
More informationERAD THE SEVENTH EUROPEAN CONFERENCE ON RADAR IN METEOROLOGY AND HYDROLOGY
Multi-beam raindrop size distribution retrieals on the oppler spectra Christine Unal Geoscience and Remote Sensing, TU-elft Climate Institute, Steinweg 1, 68 CN elft, Netherlands, c.m.h.unal@tudelft.nl
More informationPurpose of the experiment
Impulse and Momentum PES 116 Adanced Physics Lab I Purpose of the experiment Measure a cart s momentum change and compare to the impulse it receies. Compare aerage and peak forces in impulses. To put the
More informationAn Optimal Split-Plot Design for Performing a Mixture-Process Experiment
Science Journal of Applied Mathematics and Statistics 217; 5(1): 15-23 http://www.sciencepublishinggroup.com/j/sjams doi: 1.11648/j.sjams.21751.13 ISSN: 2376-9491 (Print); ISSN: 2376-9513 (Online) An Optimal
More informationTwo-Dimensional Variational Analysis of Near-Surface Moisture from Simulated Radar Refractivity-Related Phase Change Observations
ADVANCES IN ATMOSPHERIC SCIENCES, VOL. 30, NO. 2, 2013, 291 305 Two-Dimensional Variational Analysis of Near-Surface Moisture from Simulated Radar Refractiity-Related Phase Change Obserations Ken-ichi
More informationParameters Identification of Equivalent Circuit Diagrams for Li-Ion Batteries
Parameters Identification of Equialent Circuit Diagrams for Li-Ion eries Ahmad ahmoun, Helmuth Biechl Uniersity of Applied ciences Kempten Ahmad.ahmoun@stud.fh-empten.de, biechl@fh-empten.de Abstract-eries
More informationNonlinear Trajectory Tracking for Fixed Wing UAVs via Backstepping and Parameter Adaptation. Wei Ren
AIAA Guidance, Naigation, and Control Conference and Exhibit 5-8 August 25, San Francisco, California AIAA 25-696 Nonlinear Trajectory Tracking for Fixed Wing UAVs ia Backstepping and Parameter Adaptation
More informationA study of jacking force for a curved pipejacking
Journal of Rock Mechanics and Geotechnical Engineering. 200, 2 (4): 298 304 A study of jacking force for a cured pipejacking K. J. Shou, J. M. Jiang Department of Ciil Engineering, National Chung Hsing
More informationSEG/San Antonio 2007 Annual Meeting
Characterization and remoal of errors due to magnetic anomalies in directional drilling Nathan Hancock * and Yaoguo Li Center for Graity, Electrical, and Magnetic studies, Department of Geophysics, Colorado
More informationDetection of Critical Driving Situations using Phase Plane Method for Vehicle Lateral Dynamics Control by Rear Wheel Steering
Proceedings of the 7th World Congress The International Federation of Automatic Control Seoul, Korea, July 6-, 28 Detection of Critical Driing Situations using Phase Plane Method for Vehicle Lateral Dynamics
More informationPosition in the xy plane y position x position
Robust Control of an Underactuated Surface Vessel with Thruster Dynamics K. Y. Pettersen and O. Egeland Department of Engineering Cybernetics Norwegian Uniersity of Science and Technology N- Trondheim,
More informationStatement: This paper will also be published during the 2017 AUTOREG conference.
Model Predictie Control for Autonomous Lateral Vehicle Guidance M. Sc. Jochen Prof. Dr.-Ing. Steffen Müller TU Berlin, Institute of Automotie Engineering Germany Statement: This paper will also be published
More informationA Comparative Study of Vision-Based Lateral Control Strategies for Autonomous Highway Driving
A Comparatie Study of Vision-Based ateral Control Strategies for Autonomous Highway Driing J. Košecká, R. Blasi, C. J. Taylor and J. Malik Department of Electrical Engineering and Computer Sciences Uniersity
More informationDesign and Control of Novel Tri-rotor UAV
Design and Control of Noel Tri-rotor UAV Mohamed Kara Mohamed School of Electrical and Electronic Engineering The Uniersity of Manchester Manchester, UK, M 9PL Email: Mohamed.KaraMohamed@postgrad.manchester.ac.uk
More information0 a 3 a 2 a 3 0 a 1 a 2 a 1 0
Chapter Flow kinematics Vector and tensor formulae This introductory section presents a brief account of different definitions of ector and tensor analysis that will be used in the following chapters.
More informationSimulations of Space Probes and their Motions Relative to the Host Orbital Station
Simulations of Space robes and their Motions Relatie to the Host Orbital Station Eugene I. Butiko St. etersburg State Uniersity, St. etersburg, Russia Abstract Launching a space probe from an orbiting
More informationThe Importance of Anisotropy for Prestack Imaging
The Importance of Anisotropy for Prestack Imaging Xiaogui Miao*, Daid Wilkinson and Scott Cheadle Veritas DGC Inc., 715-5th Ae. S.W., Calgary, AB, T2P 5A2 Calgary, AB Xiaogui_Miao@eritasdgc.com ABSTRACT
More informationVISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION
VISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION Predict Obsere Explain Exercise 1 Take an A4 sheet of paper and a heay object (cricket ball, basketball, brick, book, etc). Predict what will
More informationSemi-Empirical Formulas of Drag/Lift Coefficients for High Speed Rigid Body Manoeuvring in Water Column
Semi-Empirical Formulas of Drag/Lift Coefficients for High Speed Rigid Body Manoeuring in Water Column Peter C. Chu and Chenwu Fan Naal Ocean Analysis and Prediction Laboratory Naal Postgraduate School,
More informationCollective circular motion of multi-vehicle systems with sensory limitations
Collectie circular motion of multi-ehicle systems with sensory limitations Nicola Ceccarelli, Mauro Di Marco, Andrea Garulli, Antonello Giannitrapani Dipartimento di Ingegneria dell Informazione Uniersità
More informationMOTION OF FALLING OBJECTS WITH RESISTANCE
DOING PHYSICS WIH MALAB MECHANICS MOION OF FALLING OBJECS WIH RESISANCE Ian Cooper School of Physics, Uniersity of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECORY FOR MALAB SCRIPS mec_fr_mg_b.m Computation
More informationCollective circular motion of multi-vehicle systems with sensory limitations
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seille, Spain, December 12-15, 2005 MoB032 Collectie circular motion of multi-ehicle systems with
More informationGRATING-LOBE PATTERN RETRIEVAL FROM NOISY IRREGULAR BEAM DATA FOR THE PLANCK SPACE TELESCOPE
GRATING-LOBE PATTERN RETRIEVAL FROM NOISY IRREGULAR BEAM DATA FOR THE PLANCK SPACE TELESCOPE Per Heighwood Nielsen (1), Oscar Borries (1), Frank Jensen (1), Jan Tauber (2), Arturo Martín-Polegre (2) (1)
More informationUNDERSTAND MOTION IN ONE AND TWO DIMENSIONS
SUBAREA I. COMPETENCY 1.0 UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS MECHANICS Skill 1.1 Calculating displacement, aerage elocity, instantaneous elocity, and acceleration in a gien frame of reference
More informationRealization of Hull Stability Control System for Continuous Track Vehicle with the Robot Arm
Adanced Science and Technology Letters Vol.86 (Ubiquitous Science and Engineering 05), pp.96-0 http://dx.doi.org/0.457/astl.05.86.0 Realization of Hull Stability Control System for Continuous Track Vehicle
More informationNonlinear Filtering. With Polynomial Chaos. Raktim Bhattacharya. Aerospace Engineering, Texas A&M University uq.tamu.edu
Nonlinear Filtering With Polynomial Chaos Raktim Bhattacharya Aerospace Engineering, Texas A&M University uq.tamu.edu Nonlinear Filtering with PC Problem Setup. Dynamics: ẋ = f(x, ) Sensor Model: ỹ = h(x)
More informationMotion in Two and Three Dimensions
PH 1-A Fall 014 Motion in Two and Three Dimensions Lectures 4,5 Chapter 4 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) 1 Chapter 4 Motion in Two and Three Dimensions In this chapter
More informationarxiv: v1 [physics.comp-ph] 17 Jan 2014
An efficient method for soling a correlated multi-item inentory system Chang-Yong Lee and Dongu Lee The Department of Industrial & Systems Engineering, Kongu National Uniersity, Kongu 34-70 South Korea
More informationSymmetric Functions and Difference Equations with Asymptotically Period-two Solutions
International Journal of Difference Equations ISSN 0973-532, Volume 4, Number, pp. 43 48 (2009) http://campus.mst.edu/ijde Symmetric Functions Difference Equations with Asymptotically Period-two Solutions
More informationDO PHYSICS ONLINE. WEB activity: Use the web to find out more about: Aristotle, Copernicus, Kepler, Galileo and Newton.
DO PHYSICS ONLINE DISPLACEMENT VELOCITY ACCELERATION The objects that make up space are in motion, we moe, soccer balls moe, the Earth moes, electrons moe, - - -. Motion implies change. The study of the
More informationThe hierarchical real-time control of high speed trains for automatic train operation
Computers in Railways XIV 17 The hierarchical real-time control of high speed trains for automatic train operation Q. Y. Wang, P. Wu, Z. C. Liang & X. Y. Feng Southwest Jiaotong Uniersity, China Abstract
More informationUnit 11: Vectors in the Plane
135 Unit 11: Vectors in the Plane Vectors in the Plane The term ector is used to indicate a quantity (such as force or elocity) that has both length and direction. For instance, suppose a particle moes
More informationNAVAL POSTGRADUATE SCHOOL Monterey, California THESIS ANGULAR RATE ESTIMATION FOR MULTI-BODY SPACECRAFT ATTITUDE CONTROL. William J.
NAVAL POSTGRADUATE SCHOOL Monterey, California THESIS ANGULAR RATE ESTIMATION FOR MULTI-BODY SPACECRAFT ATTITUDE CONTROL by William J. Palermo June 2001 Thesis Adisor: Second Reader: Brij N. Agrawal Harold
More informationVariance Reduction for Stochastic Gradient Optimization
Variance Reduction for Stochastic Gradient Optimization Chong Wang Xi Chen Alex Smola Eric P. Xing Carnegie Mellon Uniersity, Uniersity of California, Berkeley {chongw,xichen,epxing}@cs.cmu.edu alex@smola.org
More informationSUPPLEMENTARY MATERIAL. Authors: Alan A. Stocker (1) and Eero P. Simoncelli (2)
SUPPLEMENTARY MATERIAL Authors: Alan A. Stocker () and Eero P. Simoncelli () Affiliations: () Dept. of Psychology, Uniersity of Pennsylania 34 Walnut Street 33C Philadelphia, PA 94-68 U.S.A. () Howard
More informationMotion in Two and Three Dimensions
PH 1-1D Spring 013 Motion in Two and Three Dimensions Lectures 5,6,7 Chapter 4 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) 1 Chapter 4 Motion in Two and Three Dimensions In this chapter
More informationNonlinear Disturbance Decoupling for a Mobile Robotic Manipulator over Uneven Terrain
Nonlinear Disturbance Decoupling for a Mobile Robotic Manipulator oer Uneen Terrain Joel Jimenez-Lozano Bill Goodwine Department of Aerospace and Mechanical Engineering Uniersity of Noe Dame Noe Dame IN
More informationAlternative non-linear predictive control under constraints applied to a two-wheeled nonholonomic mobile robot
Alternatie non-linear predictie control under constraints applied to a to-heeled nonholonomic mobile robot Gaspar Fontineli Dantas Júnior *, João Ricardo Taares Gadelha *, Carlor Eduardo Trabuco Dórea
More informationSection 6: PRISMATIC BEAMS. Beam Theory
Beam Theory There are two types of beam theory aailable to craft beam element formulations from. They are Bernoulli-Euler beam theory Timoshenko beam theory One learns the details of Bernoulli-Euler beam
More informationA Regularization Framework for Learning from Graph Data
A Regularization Framework for Learning from Graph Data Dengyong Zhou Max Planck Institute for Biological Cybernetics Spemannstr. 38, 7076 Tuebingen, Germany Bernhard Schölkopf Max Planck Institute for
More informationEstimation of Efficiency with the Stochastic Frontier Cost. Function and Heteroscedasticity: A Monte Carlo Study
Estimation of Efficiency ith the Stochastic Frontier Cost Function and Heteroscedasticity: A Monte Carlo Study By Taeyoon Kim Graduate Student Oklahoma State Uniersity Department of Agricultural Economics
More information(a) During the first part of the motion, the displacement is x 1 = 40 km and the time interval is t 1 (30 km / h) (80 km) 40 km/h. t. (2.
Chapter 3. Since the trip consists of two parts, let the displacements during first and second parts of the motion be x and x, and the corresponding time interals be t and t, respectiely. Now, because
More informationJournal of Computational and Applied Mathematics. New matrix iterative methods for constraint solutions of the matrix
Journal of Computational and Applied Mathematics 35 (010 76 735 Contents lists aailable at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elseier.com/locate/cam New
More informationDynamic Vehicle Routing with Moving Demands Part II: High speed demands or low arrival rates
ACC 9, Submitted St. Louis, MO Dynamic Vehicle Routing with Moing Demands Part II: High speed demands or low arrial rates Stephen L. Smith Shaunak D. Bopardikar Francesco Bullo João P. Hespanha Abstract
More informationState-space Modelling of Hysteresis-based Control Schemes
European Control Conference (ECC) July 7-9,, Zürich, Switzerland. State-space Modelling of Hysteresis-based Control Schemes Soumya Kundu Ian A. Hiskens Abstract The paper deelops a state-space model for
More informationTools for Investigation of Dynamics of DC-DC Converters within Matlab/Simulink
Tools for Inestigation of Dynamics of DD onerters within Matlab/Simulink Riga Technical Uniersity, Riga, Latia Email: pikulin03@inbox.l Dmitry Pikulin Abstract: In this paper the study of complex phenomenon
More informationarxiv: v1 [stat.ml] 15 Feb 2018
1 : A New Algorithm for Streaming PCA arxi:1802.054471 [stat.ml] 15 Feb 2018 Puyudi Yang, Cho-Jui Hsieh, Jane-Ling Wang Uniersity of California, Dais pydyang, chohsieh, janelwang@ucdais.edu Abstract In
More informationMagnetic Fields Part 3: Electromagnetic Induction
Magnetic Fields Part 3: Electromagnetic Induction Last modified: 15/12/2017 Contents Links Electromagnetic Induction Induced EMF Induced Current Induction & Magnetic Flux Magnetic Flux Change in Flux Faraday
More informationDynamics ( 동역학 ) Ch.2 Motion of Translating Bodies (2.1 & 2.2)
Dynamics ( 동역학 ) Ch. Motion of Translating Bodies (. &.) Motion of Translating Bodies This chapter is usually referred to as Kinematics of Particles. Particles: In dynamics, a particle is a body without
More informationPOSITION CONTROL OF AN INTERIOR PERMANENT MAGNET SYNCHRONOUS MOTOR BY USING ADAPTIVE BACK STEPPING ALGORITHM
Research and Applications (IJERA) ISSN: 2248-9622 www.ijera.com Vol. 2, Issue 3, May-Jun 212, pp.27-276 POSITION CONTROL OF AN INTERIOR PERMANENT MAGNET SYNCHRONOUS MOTOR BY USING ADAPTIVE BACK STEPPING
More informationDoppler shifts in astronomy
7.4 Doppler shift 253 Diide the transformation (3.4) by as follows: = g 1 bck. (Lorentz transformation) (7.43) Eliminate in the right-hand term with (41) and then inoke (42) to yield = g (1 b cos u). (7.44)
More informationHYPOTHESIS TESTING SAMPLING DISTRIBUTION. the sampling distribution for di erences of means is. 2 is known. normal if.
Introduction to Statistics in Psychology PSY Professor Greg Francis Lecture 5 Hypothesis testing for two sample case Why do we let people die? H : µ = a H a : µ 6= a H : = a H a : 6= a always compare one-sample
More informationModeling Highway Traffic Volumes
Modeling Highway Traffic Volumes Tomáš Šingliar1 and Miloš Hauskrecht 1 Computer Science Dept, Uniersity of Pittsburgh, Pittsburgh, PA 15260 {tomas, milos}@cs.pitt.edu Abstract. Most traffic management
More informationUNIVERSITY OF TRENTO ITERATIVE MULTI SCALING-ENHANCED INEXACT NEWTON- METHOD FOR MICROWAVE IMAGING. G. Oliveri, G. Bozza, A. Massa, and M.
UNIVERSITY OF TRENTO DIPARTIMENTO DI INGEGNERIA E SCIENZA DELL INFORMAZIONE 3823 Poo Trento (Italy), Via Sommarie 4 http://www.disi.unitn.it ITERATIVE MULTI SCALING-ENHANCED INEXACT NEWTON- METHOD FOR
More information4. A Physical Model for an Electron with Angular Momentum. An Electron in a Bohr Orbit. The Quantum Magnet Resulting from Orbital Motion.
4. A Physical Model for an Electron with Angular Momentum. An Electron in a Bohr Orbit. The Quantum Magnet Resulting from Orbital Motion. We now hae deeloped a ector model that allows the ready isualization
More information>
The grade length should include 5% of the length of the ertical cures at the start and end of the grade. With two consecutie upgrades, 50% of the length of the ertical cure joining them should be included.
More informationThe challenges associated with Meteosat Third Generation Collocation strategies
SpaceOps 200 ConferenceDeliering on the DreamHosted by NASA Mars 25-30 April 200, Huntsille, Alabama AIAA 200-2255 The challenges associated with Meteosat Third Generation Collocation
More informationChapter 1. The Postulates of the Special Theory of Relativity
Chapter 1 The Postulates of the Special Theory of Relatiity Imagine a railroad station with six tracks (Fig. 1.1): On track 1a a train has stopped, the train on track 1b is going to the east at a elocity
More informationHYPOTHESIS TESTING SAMPLING DISTRIBUTION
Introduction to Statistics in Psychology PSY Professor Greg Francis Lecture 5 Hypothesis testing for two means Why do we let people die? HYPOTHESIS TESTING H : µ = a H a : µ 6= a H : = a H a : 6= a always
More informationInelastic Collapse in One Dimensional Systems with Ordered Collisions
Lake Forest College Lake Forest College Publications Senior Theses Student Publications 4-4-05 Inelastic Collapse in One Dimensional Systems with Ordered Collisions Brandon R. Bauerly Lake Forest College,
More information(a) Taking the derivative of the position vector with respect to time, we have, in SI units (m/s),
Chapter 4 Student Solutions Manual. We apply Eq. 4- and Eq. 4-6. (a) Taking the deriatie of the position ector with respect to time, we hae, in SI units (m/s), d ˆ = (i + 4t ˆj + tk) ˆ = 8tˆj + k ˆ. dt
More informationSemi-implicit Treatment of the Hall Effect in NIMROD Simulations
Semi-implicit Treatment of the Hall Effect in NIMROD Simulations H. Tian and C. R. Soinec Department of Engineering Physics, Uniersity of Wisconsin-Madison Madison, WI 5376 Presented at the 45th Annual
More informationEngineering Notes. CONVENTIONAL streamlined autonomous underwater vehicles. Global Directional Control of a Slender Autonomous Underwater Vehicle
JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS Vol. 3, No., January February 27 Engineering Notes ENGINEERING NOTES are short manuscripts describing new deelopments or important results of a preliminary nature.
More informationResidual migration in VTI media using anisotropy continuation
Stanford Exploration Project, Report SERGEY, Noember 9, 2000, pages 671?? Residual migration in VTI media using anisotropy continuation Tariq Alkhalifah Sergey Fomel 1 ABSTRACT We introduce anisotropy
More informationSmall area estimation under a two-part random effects model with application to estimation of literacy in developing countries
Surey Methodology, December 008 35 Vol. 34, No., pp. 35-49 Small area estimation under a two-part random effects model with application to estimation of literacy in deeloping countries Danny Pfeffermann,
More informationdifferent formulas, depending on whether or not the vector is in two dimensions or three dimensions.
ectors The word ector comes from the Latin word ectus which means carried. It is best to think of a ector as the displacement from an initial point P to a terminal point Q. Such a ector is expressed as
More informationDynamic Distribution State Estimation Using Synchrophasor Data
1 Dynamic Distribution State Estimation Using Synchrophasor Data Jianhan Song, Emiliano Dall Anese, Andrea Simonetto, and Hao Zhu arxi:1901.025541 [math.oc] 8 Jan 2019 Abstract The increasing deployment
More informationOn the influence of horizontal temperature stratification of seawater on the range of underwater sound signals. H. Lichte, Kiel, Germany
On the influence of horizontal temperature stratification of seawater on the range of underwater sound signals Original title: Über den Einfluß horizontaler Temperaturschichtung des Seewassers auf die
More informationAssignment 4 (Solutions) NPTEL MOOC (Bayesian/ MMSE Estimation for MIMO/OFDM Wireless Communications)
Assignment 4 Solutions NPTEL MOOC Bayesian/ MMSE Estimation for MIMO/OFDM Wireless Communications The system model can be written as, y hx + The MSE of the MMSE estimate ĥ of the aboe mentioned system
More informationPeriodic autoregressive stochastic volatility
MPRA Munich Personal RePEc Archie Periodic autoregressie stochastic olatility Abdelhakim Aknouche Uniersity of Science and Technology Houari Boumediene, Qassim Uniersity June 3 Online at https://mpra.ub.uni-muenchen.de/697/
More informationSpace Probe and Relative Motion of Orbiting Bodies
Space robe and Relatie Motion of Orbiting Bodies Eugene I. Butiko Saint etersburg State Uniersity, Saint etersburg, Russia E-mail: e.butiko@phys.spbu.ru bstract. Seeral possibilities to launch a space
More informationRelativity in Classical Mechanics: Momentum, Energy and the Third Law
Relatiity in Classical Mechanics: Momentum, Energy and the Third Law Roberto Assumpção, PUC-Minas, Poços de Caldas- MG 37701-355, Brasil assumpcao@pucpcaldas.br Abstract Most of the logical objections
More informationChapter 1: Kinematics of Particles
Chapter 1: Kinematics of Particles 1.1 INTRODUCTION Mechanics the state of rest of motion of bodies subjected to the action of forces Static equilibrium of a body that is either at rest or moes with constant
More informationAre the two relativistic theories compatible?
Are the two relatiistic theories compatible? F. Selleri Dipartimento di Fisica - Uniersità di Bari INFN - Sezione di Bari In a famous relatiistic argument ("clock paradox") a clock U 1 is at rest in an
More informationInformation Propagation Speed in Bidirectional Vehicular Delay Tolerant Networks
Information Propagation Speed in Bidirectional Vehicular Delay Tolerant Networks Emmanuel Baccelli, Philippe Jacquet, Bernard Mans, Georgios Rodolakis INRIA, France, {firstname.lastname}@inria.fr Macquarie
More informationNAVAL POSTGRADUATE SCHOOL Monterey, California
NPS-MV-- NAVAL POSTGRADUATE SCHOOL Monterey, California Real-Time Tracking and Display of Human Limb Segment Motions Using Sourceless Sensors and a Quaternion-Based Filtering Algorithm Part I: Theory by
More informationFeb 6, 2013 PHYSICS I Lecture 5
95.141 Feb 6, 213 PHYSICS I Lecture 5 Course website: faculty.uml.edu/pchowdhury/95.141/ www.masteringphysics.com Course: UML95141SPRING213 Lecture Capture h"p://echo36.uml.edu/chowdhury213/physics1spring.html
More informationFrames of Reference, Energy and Momentum, with
Frames of Reference, Energy and Momentum, with Interactie Physics Purpose: In this lab we will use the Interactie Physics program to simulate elastic collisions in one and two dimensions, and with a ballistic
More informationOnline Companion to Pricing Services Subject to Congestion: Charge Per-Use Fees or Sell Subscriptions?
Online Companion to Pricing Serices Subject to Congestion: Charge Per-Use Fees or Sell Subscriptions? Gérard P. Cachon Pnina Feldman Operations and Information Management, The Wharton School, Uniersity
More informationUniversity of Babylon College of Engineering Mechanical Engineering Dept. Subject : Mathematics III Class : 2 nd First Semester Year :
Uniersity of Babylon College of Engineering Mechanical Engineering Dept. Subject : Mathematics III Class : nd First Semester Year : 16-17 VECTOR FUNCTIONS SECTION 13. Ideal Projectile Motion Ideal Projectile
More informationReal Gas Thermodynamics. and the isentropic behavior of substances. P. Nederstigt
Real Gas Thermodynamics and the isentropic behaior of substances. Nederstigt ii Real Gas Thermodynamics and the isentropic behaior of substances by. Nederstigt in partial fulfillment of the requirements
More informationA-level Mathematics. MM03 Mark scheme June Version 1.0: Final
-leel Mathematics MM0 Mark scheme 660 June 0 Version.0: Final Mark schemes are prepared by the Lead ssessment Writer and considered, together with the releant questions, by a panel of subject teachers.
More informationFREEWAY WEAVING. Highway Capacity Manual 2000 CHAPTER 24 CONTENTS EXHIBITS
CHAPTER 24 FREEWAY WEAVING CONTENTS I. INTRODUCTION... 24-1 Scope of the Methodology...24-1 Limitations of the Methodology...24-1 II. METHODOLOGY...24-1 LOS...24-2 Weaing Segment Parameters...24-3 Determining
More informationExponential stability of PI control for Saint-Venant equations
Exponential stability of PI control for Saint-Venant equations Georges Bastin and Jean-Michel Coron February 28, 218 Abstract We consider open channels represented by Saint-Venant equations that are monitored
More informationDynamic Vehicle Routing with Heterogeneous Demands
Dynamic Vehicle Routing with Heterogeneous Demands Stephen L. Smith Marco Paone Francesco Bullo Emilio Frazzoli Abstract In this paper we study a ariation of the Dynamic Traeling Repairperson Problem DTRP
More informationOBSERVATIONS ON BAGGING
OBSERVATIONS ON BAGGING Andreas Buja and Werner Stuetzle Uniersity of Pennsylania and Uniersity of Washington Abstract: Bagging is a deice intended for reducing the prediction error of learning algorithms.
More informationRELATIVISTIC DOPPLER EFFECT AND VELOCITY TRANSFORMATIONS
Fundamental Journal of Modern Physics ISSN: 49-9768 Vol. 11, Issue 1, 018, Pages 1-1 This paper is aailable online at http://www.frdint.com/ Published online December 11, 017 RELATIVISTIC DOPPLER EFFECT
More informationNonlinear Lateral Control of Vision Driven Autonomous Vehicles*
Paper Machine Intelligence & Robotic Control, Vol. 5, No. 3, 87 93 (2003) Nonlinear ateral Control of Vision Drien Autonomous Vehicles* Miguel Ángel Sotelo Abstract: This paper presents the results of
More informationA Geometric Review of Linear Algebra
A Geometric Reiew of Linear Algebra The following is a compact reiew of the primary concepts of linear algebra. The order of presentation is unconentional, with emphasis on geometric intuition rather than
More informationarxiv: v2 [hep-ph] 21 Sep 2013
Higher order extensions of the Gaussian effectie potential Fabio Siringo Dipartimento di Fisica e Astronomia dell Uniersità di Catania, INFN Sezione di Catania, Via S.Sofia 64, I-9513 Catania, Italy (Dated:
More informationReversal in time order of interactive events: Collision of inclined rods
Reersal in time order of interactie eents: Collision of inclined rods Published in The European Journal of Physics Eur. J. Phys. 27 819-824 http://www.iop.org/ej/abstract/0143-0807/27/4/013 Chandru Iyer
More information