State of Charge Estimation of Cells in Series Connection by Using only the Total Voltage Measurement

Size: px
Start display at page:

Download "State of Charge Estimation of Cells in Series Connection by Using only the Total Voltage Measurement"

Transcription

1 213 American Control Conference (ACC) Washington, DC, USA, June 17-19, 213 State of Charge Estimation of Cells in Series Connection by Using only the Total Voltage Measurement Xinfan Lin 1, Anna G. Stefanopoulou 1, Yonghua Li 2 an R. Dyche Anerson 2 Abstract The voltage of lithium ion batteries is usually monitore to prevent overcharge an overischarge. For battery packs consisting of hunres of cells, monitoring the voltage of every single cell as significant cost an complexity to the battery management system (BMS). Reucing voltage sensing by only measuring the total voltage of multiple cells in series connection is esirable if the state of charge () of iniviual cells can be correctly estimate. Such goal cannot be achieve by an extene Kalman filter, because the cell s are not observable in the linearize battery string moel. In this paper, an observer base on solving simultaneously multiple nonlinear equations along the trajectory of evolution is use for the estimation problem. Existence of the solution epens on the nonlinearity of the battery voltage- relationship. The observer is applie to a LiFePO 4 /graphite battery string with 2 cells, where the iniviual cell s are observable in low an high ranges. Experimental results show goo convergence of an voltage estimation, inicating that this new methoology can be applie to, at least, halve the voltage sensing in a battery pack. I. INTRODUCTION Automotive battery systems for hybri electric vehicles (HEV), plug-in hybri electric vehicles (PHEV), an battery electric vehicles (BEV) usually consist of hunres an even thousans of cells. The number of cells in a battery system is a function of the esire total voltage, power an energy, an the characteristics of the selecte cell chemistry. Battery packs in HEVs with Nickel Metal Hyrie (NiMH) batteries typically have cells. Packs using lithium ion cells for the same type of vehicles have about 5-1 cells ue to the higher voltage inherent to the lithium ion chemistry. Prouction HEVs with NiMH batteries typically measure total voltage of moules with 5 to 16 cells in series. Due to the sensitivity of most lithium ion cells to overcharge an overischarge [1], [2], existing battery management systems (BMS) for lithium ion battery packs nee to monitor the state of charge () an voltage of every single cell. This single-cell monitoring strategy as significant cost to the lithium ion battery system incluing sensors, wiring an labor. Therefore, it is highly esirable to reuce the voltage sensing, e.g. by only measuring the total voltage of multiple cells, uner the conition that the an voltage of iniviual cells can still be monitore. *This work is supporte by the For Motor Company. 1 X. Lin an A. Stefanopoulou are with the Department of Mechanical Engineering, University of Michigan, Ann Arbor, MI 4819, USA. xflin@umich.eu an annastef@umich.eu 2 Y. Li an R.D. Anerson are with the Vehicle an Battery Controls Department, Research an Avance Engineering, For Motor Company, Dearborn, MI 48121, USA. yli19@for.com an raner34@for.com In this paper, the possibility of estimating the s of batteries in series connection by only using the total voltage measurement is investigate. First, we try to use the extene Kalman filter (EKF), which have been wiely use to estimate the battery. The EKF works well when the single cell voltage is measure for feeback [3], [4], [5]. However, it has foun that when only the total voltage of a battery string is measure, the EKF can only track the average of the string but not the iniviual cell s. The limitation arises from the lack of observability of the linearize battery string moel iscusse in Section III. In [6], the cell balancing circuit is use to augment the voltage measurement when only the total voltage is measure. In this paper, instea, it is shown that the iniviual cell s are observable in the sense of nonlinear observability for some lithium ion battery chemistry. For example, a 2-cell string with LiFePO 4 /graphite batteries, whose voltage- relationship is highly nonlinear, is shown to have observable iniviual cell s in low an high ranges. Such nonlinear observability provies the possibility of estimating iniviual cell s by using methos other than the EKF. The Newton observer, which is base on solving simultaneously multiple nonlinear equations along the system trajectory, [7], [8], has been propose to aress the state estimation problem with high nonlinearity. In this paper, a similar observer base on the Levenberg-Marquart algorithm [9], [1] is esigne for iniviual cell an voltage estimation uner reuce voltage sensing. Experiments uner a constant charging current have been conucte to valiate the propose algorithm with a 2-cell battery string. It is emonstrate that goo convergence of an voltage estimation can be achieve in the high range, where monitoring overcharge is of critical importance. II. GENERALIZED BATTERY MODEL FORMULATION Some of the most commonly use control-oriente battery moels can be generalize in a state-space representation [11], as ẋ=ax+bi V = g(x,i), where the state matrix A epens on the transient ynamics associate with particle iffusion [3], B on the battery capacity, an g is the nonlinear voltage output function. The state(s) x relates to the energy storage level, the input I is the current an the output V is the voltage. Such moels inclue an are not limite to: Coulomb counting moel (CCM) [12], [13] (1) /$ AACC 74

2 Equivalent circuit moel (ECM) [13], [14] Simplifie electrochemical moel featuring lithium ion transport an iffusion uring battery charge an ischarge (SECM) [3], [15]. Physical meaning of the variables an parameters in (1) for the above three moels are summarize in Table I. More etails on the moel generalization can be foun in [11]. It is note that there are more complicate partial ifferential equation (PDE) base battery moels, which capture etaile electrochemical processes uring the battery operation [16], [17]. These moels may not be written in the form of (1) with linear state equations an a nonlinear output equation. For a battery string connecte in series, the current input is the same for each cell, an the measure output, the total voltage, is the summation of all the single cell voltages. In this paper, we consier the situation that the cells may only iffer in their s, an all the cells have the same moel parameters. Such imbalance can be cause by factors such as ifference among cells in self-ischarge rate, which is typically relate to the temperature variability along the string. Inee, since the battery egraation is also temperatureepenent, this variability may also result in imbalance in capacity an resistance among cells. Such imbalance is more likely to be evelope over longer time-span an will be aresse in the future work. Base on the single cell moel in (1), a string with r cells in series can be moele as, where ẋ str =F(x str,i )=A str x str + B str I V str =G(x str,i)= r i=1 g(x i,i), x str = [ x 1 x r ] T, A str = iag(a,,a), B str = [ B T B T B T] T. (3) This general form will be use in the analysis in the subsequent sections. The state estimation of (1), where the single cell voltage is measure, is referre to as the estimation uner full voltage sensing, an the state estimation of (2), where only the total voltage is measure, is calle the estimation uner reuce voltage sensing. III. EXTENDED KALMAN FILTER AND OBSERVABILITY OF THE LINEARIZED BATTERY STRING MODEL The extene Kalman filter is one of the most commonly use methos for nonlinear state estimation problems, an has been wiely aopte to estimate the battery uner full voltage sensing [3], [4], [5]. However, it is shown here that the EKF cannot be use uner reuce sensing because the linearize battery string moel is not observable. The EKF for (1) takes the form ˆx=A ˆx+Bu+K(V ˆV) ˆV = g( ˆx,I), where ˆx an ˆV are the estimate state(s), e.g., an voltage, V is the measure voltage an K is the observer (2) (4) gain. The estimate state of the EKF will converge to the real state asymptotically, if the linearize system moel is observable along the system trajectory [18]. Furthermore, the linear observability practically guarantees that the errors in estimation by the EKF will be boune by some constant value even if the moel parameters are not precise [11]. Such observability conition is satisfie if the single cell voltage is measure, leaing to the popularity of EKF uner full voltage sensing. When only the total voltage of the string is measure, the linearize battery string moel can be obtaine by linearizing the output voltage function g of each single cell, as g(x i,i)= C i x+d i I. (5) As liste in Table I, C is mainly the slope of the open circuit voltage (OCV) curve an D is the ohmic resistance. Since g is a nonlinear function of the state, the values of C an D can be ifferent for cells with ifferent s. Hence, the linearize battery string moel will be ẋ str = A str x+b str I V str = C str x+d str I C str = [ ] C 1 C 2 C r (6) D str = r i=1 D i, with C i C j an D i D j if i j. The moel is observable if an only if its observability matrix C str C 1 C 2 C r U = C str A str = C 1 A C 2 A C r A C str A n 1 str C 1 A n 1 C 2 A n 1 C r A n 1 (7) }{{} U 1 }{{} U 2 }{{} U r is of full rank [18], [19], where n is the imension of the A str matrix. Accoring to Tab.I, for the coulomb counting an the simplifie electrochemical moels, C i s are proportional to each other since they only have one non-zero element, C j = η i, j C i, (8) where η i, j = α j α i for CCM, an η i, j = α j+β j α i +β j for SECM. Hence the U i s in (7) is linearly epenent on each other, C j η i, j C i C i U j = C j A = η i, j C i A = η C i A i, j (9) C j A n 1 η i, j C i A n 1 C i A n 1 = η i, j U i. As a result, the observability matrix will be short of full rank, inicating the the linearize battery string moel is unobservable. The same conclusion can be rawn for the equivalent circuit moel upon similar analysis. Simulation has been conucte to show the inaequacy of the EKF for the estimation problem uner reuce voltage sensing. A coulomb counting moel of the LiFePO 4 /graphite 75

3 TABLE I MODEL GENERALIZATION Coulomb counting moel (CCM) Equivalent circuit moel (ECM) Simplifie electrochemical moel (SECM) State x, spatially istribute lithium concentrations RC voltages: V c,1v c,m c s,1c s,q Input u current I current I current I Output y voltage V voltage V voltage V A R c,1 C c,1 D s 2 r... 1 R q 1 c,mc c,m q q 1 q R c,i : equivalent resistance D s : iffusion coefficient C c,i : equivalent capacitance r : spatial iscretization parameter B Q 1 Q 1 (A e δfa s r ) 1 [,,, q+1 q ]T Q: battery capacity C 1 c,i A e :electroe area, δ: electroe thickness, F: Faraay constant, a s : active particle surface g OCV()+IR, OCV(c s,q )+OP(c s,q,i)+ir f OCV : open circuit voltage OCV() Σ m i V c,i+ IR OP: over potential R: internal resistance [ ] [ R f : ohmic resistance ] C α α α+ β α: local slope of OCV w.r.t. β: local slope of OP w.r.t. D R R R f + γ γ: local slope of OP w.r.t. I battery is use to emulate a battery string with 2 cells. The output voltage function of the LiFePO 4 /graphite battery uner a constant current is shown in Fig. 1. In simulation, SÔC SÔC ˆ ˆ t (s) Fig. 1. Voltage- relationship of LiFePO 4 batteries uner a constant charging current Fig. 2. Estimation of single cell an string voltage by the EKF one cell is initialize to 5%, an the other to 4%. The two cells are charge with the constant current until the cell with higher reaches 1% mark. The total voltage of the two cells is fe back to the extene Kalman filter. Simulation results are shown in Fig. 2. It can be seen that as the two batteries are charge up, their s increase linearly with time an the ifference between them stays at 1%. The two estimate s, however, instea of tracing 1 an 2 respectively, both converge to the intermeiate values that woul make the estimate voltage match the average voltage of the string. IV. OBSERVABILITY ANALYSIS OF THE NONLINEAR BATTERY STRING MODEL Although the unobservable linearize battery string moel prevents the EKF from estimating the iniviual cell s uner reuce voltage sensing, it oes not rule out the possibility of solving the estimation problem with a nonlinear estimation technique. The existence of the solution epens on the observability of the nonlinear battery moel, which is to be investigate in this section. A nonlinear system is observable at a certain state if the graient (matrix) of the Lie erivative vector of the output function is of full rank [2]. For the string moel in (2), the ith orer Lie erivative of the output function G is L i F(G)= V (i) str = G str F(x str,i), (1) where the current I is fixe at a constant value [2]. The Lie erivative vector is efine as LF (G) l(x str )=, (11) (G) L n 1 F 76

4 an the graient of l(x str ) takes the form l(x str )= l = str L F (G) str LF n 1 (G) str. (12) Consier the coulomb counting moel as an example to analyze the moel observability, whose l(x str ) matrix is g g x1 xr 2 g x1 I l(x str )= 2 Q 2 g xr I 2 Q r g x1 r ( Q I )n 1 r, (13) g xr r ( Q I )n 1 with x i being the of the ith cell in the string. It can be seen that the rank of (13) epens on the graients of g with respect to the. Clearly, if g is linear, all its graients higher than 2n orer will be zero, an the rank of (13) will be 1. The system will thus be unobservable uner reuce voltage sensing. Nevertheless, it is possible to have a full-rank matrix in (13) if the voltage- relationship g is nonlinear. For example, the g(x) function of the LiFePO 4 /graphite battery, is shown in Fig. 1 an its graients in Fig. 3. In the range between 1% an 9%, g(x) is 2 V 2 V Fig. 3. Graients of g(x) for LiFePO 4 /graphite batteries uner a constant charging current close to linear, making the battery string moel unobservable. However, in the low ( 1%) an high (9 1%) ranges, as shown in Fig. 1 an Fig. 3, g(x) is highly nonlinear with non-zero 1st an 2n (an even higher) orer graients. It can be checke numerically that for battery strings with two cells, the l matrix is of full rank for any combination of iniviual cell s along the evolution trajectory. Therefore, it is possible to istinguish the iniviual cell s when they are in those ranges. In general cases, for strings with r cells, up to (r 1)th orer erivatives nee to be calculate to etermine the rank of (13) an hence the observability. Base on the above conclusion, a nonlinear observer will be esigne to estimate iniviual cell s in a 2-cell battery string, while the methoology is applicable to strings with more cells. V. NONLINEAR OBSERVER BASED ON THE LEVENBERG-MARQUARDT ALGORITHM Apart from extening the linear estimation theory, e.g. EKF, nonlinear observers can also be constructe base on solving multiple nonlinear equations along the system evolution trajectory, such as the newton observer [7], [8]. Such algorithms, instea of hanling the output ata one point at a time, process the output trajectory over a certain time span simultaneously. As a result, compare with methos extene from the linear estimation theory base on moel linearization, they retain the nonlinearity of the system an have a wier range of application. Since the nonlinear observer is intene for onboar BMS, it will be esigne in iscrete-time omain to accommoate sample measurements. In iscrete-time omain, the moel in (2) can be written as x str,k+1 = F (x str,k,i k ) V str,k = G(x str,k,i k ), where k is the time step. The function F equations in iscrete-time omain, as F (x str,k,i k )=A str, x str,k + B str, I k A str, = e A strt, B str, = T e A strτ B str τ, (14) is the state (15) where T is the sampling perio. At each estimation step, a set of N consecutive measurements an inputs, V str,[k N+1,k] = V str,k N+1,I [k N+1,k] = V str,k are use for estimation. Let us efine I k N+1 I k, (16) F I k+1 F I k (x ( ) str,k) := F F (x str,k,i k ),I k+1, (17) as stans for function composition. It can be erive from (14) that G I k N+1(x str,k N+1 ) V str,[k N+1,k] = G I k N+2 F I k N+1 (x str,k N+1 ) (18) G I k F I k 1 F I k N+1 (x str,k N+1 ) = H(x str,k N+1,I [k N+1,k] ). The problem is then reuce to solving (18) for x str,k N+1 given V str,[k N+1,k] an I str,[k N+1,k]. Evolution of the state can then be etermine base on (14). The Newton observer has been propose in [7] an [8] to solve (18). At each estimation step, x str,k N+1 is calculate through multiple iterations base on the Newton Raphson algorithm, [ x j+1 str,k N+1 = x j str,k N+1 + H ( x j str,k N+1,I [k N+1,k]) str,k N+1 ( Vstr,[k N+1,k] G( x j str,k N+1,I [k N+1,k]) ). (19) ] 1 77

5 The resulting x str,k N+1 at the final iteration step, is taken as the estimation for x str,k N+1, as ˆx str,k N+1 = x str,k N+1. (2) Estimates of the subsequent states, ˆx str,k N+2,..., ˆx str,k, can then be etermine by ˆx str,k i = F I k i... F I k N+1 ( ˆx str,k N+1 ). (21) At the next estimation step, V str,[k N+1,k] an I [k N+1,k] will be upate by the newly acquire measurements, an the initial guess of the state estimation is etermine base on the previous estimation. Although the Newton s metho usually converges fast, it might not be robust uner some circumstances. For example, when the graient of g(x) is small, H in (19) becomes illconitione, making it ifficult to calculate its inverse. One convenient remey is to replace the Newton iteration in (19) with the Levenberg-Marquart iteration [9], [1], as [ x j+1 str,k N+1 = x j str,k N+1 + H ( x j str,k N+1,I [k N+1,k])+ str,k N+1 λi] 1( V str,[k N+1,k] G( x j str,k N+1,I [k N+1,k]) ), (22) where λ is a scalar use to lower the conition number of the inverte matrix when g is small, an I is the ientity matrix. VI. EXPERIMENTAL VALIDATION Experiments have been conucte to valiate the esigne nonlinear observer. Two A LiFePO 4 /graphite batteries are use for testing. Before the experiment, the two cells were first initialize with ifferent s, aroun 5% an % respectively. The two cells were then connecte in series an charge up by a single current source with a constant current of 2 A. The voltage of each single cell was monitore to prevent overcharge. The actual cell s were calculate base on current integration for valiation. The current was cut off when the voltage of the cell with higher hit V. The actual s an voltages are shown in Fig. 4. The nonlinear observer in (22) is use to estimate the iniviual cell s an voltages by only using the total voltage measurement. The coulomb counting moel is use for the battery string moel, which is aequate uner the constant current charging scenario. Estimation starts at the 2th secon, an the initial guess of the is taken as the same for both cells by inverting the average measure voltage. At each estimation step, 15 secons length of ata are use, corresponing to an span of 4.15%. The time interval between ata points is 1 s an a total of 15 sample ata are processe at each step. The estimation results are plotte in Fig. 5 an Fig. 6. The plotte value correspons to the last point at each estimation step. The final values are shown in Tab.(II). It can be seen from Fig. 5 that the iniviual cell s are not istinguishable before the 12th estimation step, when t (s) t (s) Fig. 4. Actual s an voltages of iniviual cells SÔC 1 1 SÔC Estimation step Fig. 5. estimation of iniviual cells 1 2 SÔC 2 SÔC Estimation step Fig. 6. ˆ ˆ ˆ Voltage estimation of iniviual cells ˆ 78

6 TABLE II FINAL ESTIMATION RESULTS Estimation Measurement error (%) % 94.36% % 1%.19 6 V 4 V V V.28 both of them are below 85%. This has been preicte by the nonlinear observability analysis, since the moel is not observable in that range ue to almost linear voltage- relationship, as shown in Fig. 1. When the cell s get above 9%, the estimates start to converge graually to respective actual values as the moel becomes observable. It can be observe that the an voltage estimation of cell 1, which is the cell with the higher, is more accurate than that of cell 2. The reason can be attribute to the graients of g(x) at the s of the two cells. The change in total voltage can be linearize locally as δ( + )= g δ 1 + g 1 δ 2. (23) 2 Since cell 1 has higher, g 1 is larger than g 2 accoring to Fig. 3. Consequently, the change in total voltage will be more sensitive to the change in 1. This will lea to more accurate estimation of 1, which is important when the cell is nearly fully charge. VII. CONCLUSIONS This paper investigates the estimation of the iniviual cell s with only the total voltage measurement for cells in series connection. It is pointe out that the existence of the solution relies on the observability of the nonlinear battery string moel. For battery chemistry with linear voltage- relationship, the iniviual cell s are not observable uner reuce voltage sensing. Nevertheless, for some battery chemistry, such as LiFePO 4, the nonlinearity of the voltage- relationship reners moel observability in certain ranges. A nonlinear observer base on the Levenberg-Marquart algorithm is then esigne to estimate the iniviual cell s an voltages. The algorithm has been implemente to a LiFePO 4 /graphite battery string with 2 cells. As inicate by the observability analysis, the estimate s converge faster an are much more accurate at the observable high ens (than in the unobservable mile range), where the estimation is more critical. In principle, the methoology can be extene to cell strings with more cells an of other chemistry given proper voltage- relationship. REFERENCES [1] S. Hossain, Y. Saleh, an R. Loutfy, Carbon-carbon composite as anoes for lithium-ion battery systems, Journal of Power Sources, vol. 96, pp. 5 13, 21. [2] Y.-S. Lee an M.-W. Cheng, Intelligent control battery equalization for series connecte lithium-ion battery strings, IEEE Transactions on Inustrial Electronics, vol. 52, pp , 25. [3] D. D. Domenico, A. Stefanopoulou, an G. Fiengo, Lithium-ion battery state of charge an critical surface charge estimation using an electrochemical moel-base extene kalman filter, Journal of Dynamic Systems, Measurement, an Control, vol. 132, pp , 21. [4] S. Santhanagopalan an R. E. White, Online estimation of the state of charge of a lithium ion cell, Journal of Power Sources, vol. 16, no. 2, pp , 26. [5] G. L. Plett, Extene kalman filtering for battery management systems of lipb-base hev battery packs part 3. state an parameter estimation, Journal of Power Sources, vol. 134, pp , 24. [6] L. Y. Wang, M. Polis, G. Yin, W. Chen, Y. Fu, an C. Mi, Battery cell ientification an soc estimation using string terminal voltage measurements, IEEE Transactions on Vehicular Technology, vol. 61, pp , 212. [7] P. E. Moraal an J. W. Grizzle, Observer esign for nonlinear systems with iscrete-time measurements, IEEE Transactions on Automatic Control, vol. 4, pp , [8], Asymptotic observers for etectable an poorly observable systems, in IEEE Proceeings of the 34th Conference on Decision & Control, [9] K. Levenberg, A metho for the solution of certain non-linear problems in least squares, Quarterly of Applie Mathematics, vol. 2, pp , [1] D. W. Marquart, An algorithm for least-squares estimation of nonlinear parameters, Journal of the Society for Inustrial an Applie Mathematics, vol. 11, pp , [11] X. Lin, A. G. Stefanopoulou, P. Laskowsky, J. S. Freuenberg, Y. Li, an R. D. Anerson, State of charge estimation error ue to parameter mismatch in a generalize explicit lithium ion battery moel, in Proceeings of ASME Dynamic Systems an Control Conference, no. DSCC , 211, pp [12] Y. Cairci an Y. Ozkazanc, Microcontroller-base on-line state-ofcharge estimator for seale leacaci batteries, Journal of Power Sources, vol. 129, pp , 24. [13] V. Johnson, A. Pesaran, an T. Sack, Temperature-epenent battery moels for high-power lithium-lon batteries, in Proceeings of the 17th Electric Vehicle Symposium, 2. [14] H. E. Perez, J. B. Siegel, X. Lin, A. G. Stefanopoulou, Y. Ding, an M. P. Castanier, Parameterization an valiation of an integrate electro-thermal lfp battery moel, in ASME Dynamic Systems an Control Conference (DSCC), 212. [15] K. A. Smith, C. D. Rahn, an C.-Y. Wang, Moel-base electrochemical estimation an constraint management for pulse operation of lithium ion batteries, IEEE Transactions on Control Systems Technology, vol. 18, pp , 21. [16] T. F. Fuller, M. Doyle, an J. Newman, Simulation an optimization of the ual lithium ion insertion cell, Journal of the Electrochemical Society, vol. 141, no. 1, pp. 1 1, January [17] C. Y. Wang, W. B. Gu, an B. Y. Liaw, Micro-macroscopic couple moeling of batteries an fuel cells i. moel evelopment, Journal of the Electrochemical Society, vol. 145, no. 1, pp , October [18] Y. Song an J. W. Grizzle, The extene kalman filter as a local asymptotic observable for iscrete-time nonlinear systems, Journal of Mathematical Systems, Estimations an Control, vol. 5, pp , [19] R. Williams an D. Lawrence, Linear state-space control systems. Wiley, 27. [2] R. Hermann an A. J. Krener, Nonlinear controllability an observability, IEEE Transactions on Automatic Control, vol. 5, pp , ACKNOWLEDGMENT This work has been supporte by the For Motor Company (For/UMICH Alliance Project). 79

BATTERY pack, as a common energy storage system,

BATTERY pack, as a common energy storage system, 1052 IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 23, NO. 3, MAY 2015 State of Charge Imbalance Estimation for Battery Strings Under Reduced Voltage Sensing Xinfan Lin, Anna G. Stefanopoulou,

More information

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain

More information

Efficient Macro-Micro Scale Coupled Modeling of Batteries

Efficient Macro-Micro Scale Coupled Modeling of Batteries A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

Neural Network Training By Gradient Descent Algorithms: Application on the Solar Cell

Neural Network Training By Gradient Descent Algorithms: Application on the Solar Cell ISSN: 319-8753 Neural Networ Training By Graient Descent Algorithms: Application on the Solar Cell Fayrouz Dhichi*, Benyounes Ouarfi Department of Electrical Engineering, EEA&TI laboratory, Faculty of

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

Predictive Control of a Laboratory Time Delay Process Experiment

Predictive Control of a Laboratory Time Delay Process Experiment Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

PCCP PAPER. 1 Introduction. A. Nenning,* A. K. Opitz, T. M. Huber and J. Fleig. View Article Online View Journal View Issue

PCCP PAPER. 1 Introduction. A. Nenning,* A. K. Opitz, T. M. Huber and J. Fleig. View Article Online View Journal View Issue PAPER View Article Online View Journal View Issue Cite this: Phys. Chem. Chem. Phys., 2014, 16, 22321 Receive 4th June 2014, Accepte 3r September 2014 DOI: 10.1039/c4cp02467b www.rsc.org/pccp 1 Introuction

More information

A simplified macroscopic urban traffic network model for model-based predictive control

A simplified macroscopic urban traffic network model for model-based predictive control Delft University of Technology Delft Center for Systems an Control Technical report 9-28 A simplifie macroscopic urban traffic network moel for moel-base preictive control S. Lin, B. De Schutter, Y. Xi,

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

A Novel Unknown-Input Estimator for Disturbance Estimation and Compensation

A Novel Unknown-Input Estimator for Disturbance Estimation and Compensation A Novel Unknown-Input Estimator for Disturbance Estimation an Compensation Difan ang Lei Chen Eric Hu School of Mechanical Engineering he University of Aelaie Aelaie South Australia 5005 Australia leichen@aelaieeuau

More information

Nonlinear Observer Designs for State-of-Charge Estimation of Lithium-ion Batteries

Nonlinear Observer Designs for State-of-Charge Estimation of Lithium-ion Batteries 2014 American Control Conference (ACC) June 4-6, 2014. Portland, Oregon, USA Nonlinear Observer Designs for State-of-Charge Estimation of Lithium-ion Batteries Satadru Dey and Beshah Ayalew Abstract State-of-Charge

More information

PD Controller for Car-Following Models Based on Real Data

PD Controller for Car-Following Models Based on Real Data PD Controller for Car-Following Moels Base on Real Data Xiaopeng Fang, Hung A. Pham an Minoru Kobayashi Department of Mechanical Engineering Iowa State University, Ames, IA 5 Hona R&D The car following

More information

Two Dimensional Numerical Simulator for Modeling NDC Region in SNDC Devices

Two Dimensional Numerical Simulator for Modeling NDC Region in SNDC Devices Journal of Physics: Conference Series PAPER OPEN ACCESS Two Dimensional Numerical Simulator for Moeling NDC Region in SNDC Devices To cite this article: Dheeraj Kumar Sinha et al 2016 J. Phys.: Conf. Ser.

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

Optimal operating strategies for semi-batch reactor used for chromium sludge regeneration process

Optimal operating strategies for semi-batch reactor used for chromium sludge regeneration process Latest Trens in Circuits, Automatic Control an Signal Processing Optimal operating strategies for semi-batch reactor use for chromium sluge regeneration process NOOSAD DAID, MACKŮ LUBOMÍR Tomas Bata University

More information

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Yin Li, Stuent Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract This paper proposes a robust power system stabilizer (PSS)

More information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental Robustness Study of a Second-Order Sliding Mode Controller Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans

More information

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Preprints of the 9th Worl Congress The International Feeration of Automatic Control Cape Town, South Africa. August -9, Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Zhengqiang

More information

Solutions to Math 41 Second Exam November 4, 2010

Solutions to Math 41 Second Exam November 4, 2010 Solutions to Math 41 Secon Exam November 4, 2010 1. (13 points) Differentiate, using the metho of your choice. (a) p(t) = ln(sec t + tan t) + log 2 (2 + t) (4 points) Using the rule for the erivative of

More information

Text S1: Simulation models and detailed method for early warning signal calculation

Text S1: Simulation models and detailed method for early warning signal calculation 1 Text S1: Simulation moels an etaile metho for early warning signal calculation Steven J. Lae, Thilo Gross Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, 01187 Dresen, Germany

More information

Left-invariant extended Kalman filter and attitude estimation

Left-invariant extended Kalman filter and attitude estimation Left-invariant extene Kalman filter an attitue estimation Silvere Bonnabel Abstract We consier a left-invariant ynamics on a Lie group. One way to efine riving an observation noises is to make them preserve

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Experimental Determination of Mechanical Parameters in Sensorless Vector-Controlled Induction Motor Drive

Experimental Determination of Mechanical Parameters in Sensorless Vector-Controlled Induction Motor Drive Experimental Determination of Mechanical Parameters in Sensorless Vector-Controlle Inuction Motor Drive V. S. S. Pavan Kumar Hari, Avanish Tripathi 2 an G.Narayanan 3 Department of Electrical Engineering,

More information

Lithium-ion Battery Electro-thermal Model and Its Application in the Numerical Simulation of Short Circuit Experiment

Lithium-ion Battery Electro-thermal Model and Its Application in the Numerical Simulation of Short Circuit Experiment EVS27 Barcelona, Spain, November 17-20, 2013 Lithium-ion Battery Electro-thermal Moel an Its Application in the Numerical Simulation of Short Circuit Experiment Chengtao Lin 1, Can Cui 2, Xiaotian Xu 3

More information

Perturbation Analysis and Optimization of Stochastic Flow Networks

Perturbation Analysis and Optimization of Stochastic Flow Networks IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. Y, MMM 2004 1 Perturbation Analysis an Optimization of Stochastic Flow Networks Gang Sun, Christos G. Cassanras, Yorai Wari, Christos G. Panayiotou,

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

IEEE TRANSACTIONS ON COMPUTER-AIDED DESIGN OF INTEGRATED CIRCUITS AND SYSTEMS, VOL. 33, NO. 8, AUGUST

IEEE TRANSACTIONS ON COMPUTER-AIDED DESIGN OF INTEGRATED CIRCUITS AND SYSTEMS, VOL. 33, NO. 8, AUGUST IEEE TRANSACTIONS ON COMPUTER-AIDED DESIGN OF INTEGRATED CIRCUITS AND SYSTEMS, VOL. 33, NO. 8, AUGUST 214 1145 A Semi-Analytical Thermal Moeling Framework for Liqui-Coole ICs Arvin Srihar, Member, IEEE,

More information

Adjoint Transient Sensitivity Analysis in Circuit Simulation

Adjoint Transient Sensitivity Analysis in Circuit Simulation Ajoint Transient Sensitivity Analysis in Circuit Simulation Z. Ilievski 1, H. Xu 1, A. Verhoeven 1, E.J.W. ter Maten 1,2, W.H.A. Schilers 1,2 an R.M.M. Mattheij 1 1 Technische Universiteit Einhoven; e-mail:

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Control of a PEM Fuel Cell Based on a Distributed Model

Control of a PEM Fuel Cell Based on a Distributed Model 21 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July 2, 21 FrC1.6 Control of a PEM Fuel Cell Base on a Distribute Moel Michael Mangol Abstract To perform loa changes in proton

More information

How the potentials in different gauges yield the same retarded electric and magnetic fields

How the potentials in different gauges yield the same retarded electric and magnetic fields How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Thermal runaway during blocking

Thermal runaway during blocking Thermal runaway uring blocking CES_stable CES ICES_stable ICES k 6.5 ma 13 6. 12 5.5 11 5. 1 4.5 9 4. 8 3.5 7 3. 6 2.5 5 2. 4 1.5 3 1. 2.5 1. 6 12 18 24 3 36 s Thermal runaway uring blocking Application

More information

ECE 422 Power System Operations & Planning 7 Transient Stability

ECE 422 Power System Operations & Planning 7 Transient Stability ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of

More information

A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL

A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL 6th European Signal Processing Conference EUSIPCO 28, Lausanne, Switzerlan, August 25-29, 28, copyright by EURASIP A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL Leonaro Tomazeli

More information

Event based Kalman filter observer for rotary high speed on/off valve

Event based Kalman filter observer for rotary high speed on/off valve 28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June 11-13, 28 WeC9.6 Event base Kalman filter observer for rotary high spee on/off valve Meng Wang, Perry Y. Li ERC for Compact

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

Stable and compact finite difference schemes

Stable and compact finite difference schemes Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long

More information

State observers and recursive filters in classical feedback control theory

State observers and recursive filters in classical feedback control theory State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent

More information

Nested Saturation with Guaranteed Real Poles 1

Nested Saturation with Guaranteed Real Poles 1 Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Conservation laws a simple application to the telegraph equation

Conservation laws a simple application to the telegraph equation J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness

More information

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E.

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E. AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 14 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering an Applie Science Gonzaga University SDC

More information

Non-Linear Bayesian CBRN Source Term Estimation

Non-Linear Bayesian CBRN Source Term Estimation Non-Linear Bayesian CBRN Source Term Estimation Peter Robins Hazar Assessment, Simulation an Preiction Group Dstl Porton Down, UK. probins@stl.gov.uk Paul Thomas Hazar Assessment, Simulation an Preiction

More information

Dot trajectories in the superposition of random screens: analysis and synthesis

Dot trajectories in the superposition of random screens: analysis and synthesis 1472 J. Opt. Soc. Am. A/ Vol. 21, No. 8/ August 2004 Isaac Amiror Dot trajectories in the superposition of ranom screens: analysis an synthesis Isaac Amiror Laboratoire e Systèmes Périphériques, Ecole

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

UNIFYING PCA AND MULTISCALE APPROACHES TO FAULT DETECTION AND ISOLATION

UNIFYING PCA AND MULTISCALE APPROACHES TO FAULT DETECTION AND ISOLATION UNIFYING AND MULISCALE APPROACHES O FAUL DEECION AND ISOLAION Seongkyu Yoon an John F. MacGregor Dept. Chemical Engineering, McMaster University, Hamilton Ontario Canaa L8S 4L7 yoons@mcmaster.ca macgreg@mcmaster.ca

More information

An Optimal Algorithm for Bandit and Zero-Order Convex Optimization with Two-Point Feedback

An Optimal Algorithm for Bandit and Zero-Order Convex Optimization with Two-Point Feedback Journal of Machine Learning Research 8 07) - Submitte /6; Publishe 5/7 An Optimal Algorithm for Banit an Zero-Orer Convex Optimization with wo-point Feeback Oha Shamir Department of Computer Science an

More information

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device Close an Open Loop Optimal Control of Buffer an Energy of a Wireless Device V. S. Borkar School of Technology an Computer Science TIFR, umbai, Inia. borkar@tifr.res.in A. A. Kherani B. J. Prabhu INRIA

More information

Modeling time-varying storage components in PSpice

Modeling time-varying storage components in PSpice Moeling time-varying storage components in PSpice Dalibor Biolek, Zenek Kolka, Viera Biolkova Dept. of EE, FMT, University of Defence Brno, Czech Republic Dept. of Microelectronics/Raioelectronics, FEEC,

More information

Separation Principle for a Class of Nonlinear Feedback Systems Augmented with Observers

Separation Principle for a Class of Nonlinear Feedback Systems Augmented with Observers Proceeings of the 17th Worl Congress The International Feeration of Automatic Control Separation Principle for a Class of Nonlinear Feeback Systems Augmente with Observers A. Shiriaev, R. Johansson A.

More information

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,

More information

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS Yannick DEVILLE Université Paul Sabatier Laboratoire Acoustique, Métrologie, Instrumentation Bât. 3RB2, 8 Route e Narbonne,

More information

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign

More information

An algebraic expression of stable inversion for nonminimum phase systems and its applications

An algebraic expression of stable inversion for nonminimum phase systems and its applications Proceeings of the 17th Worl Congress The International Feeration of Automatic Control An algebraic expression of stable inversion for nonminimum phase systems an its applications Takuya Sogo Chubu University,

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

On Characterizing the Delay-Performance of Wireless Scheduling Algorithms

On Characterizing the Delay-Performance of Wireless Scheduling Algorithms On Characterizing the Delay-Performance of Wireless Scheuling Algorithms Xiaojun Lin Center for Wireless Systems an Applications School of Electrical an Computer Engineering, Purue University West Lafayette,

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

YVES F. ATCHADÉ, GARETH O. ROBERTS, AND JEFFREY S. ROSENTHAL

YVES F. ATCHADÉ, GARETH O. ROBERTS, AND JEFFREY S. ROSENTHAL TOWARDS OPTIMAL SCALING OF METROPOLIS-COUPLED MARKOV CHAIN MONTE CARLO YVES F. ATCHADÉ, GARETH O. ROBERTS, AND JEFFREY S. ROSENTHAL Abstract. We consier optimal temperature spacings for Metropolis-couple

More information

New Statistical Test for Quality Control in High Dimension Data Set

New Statistical Test for Quality Control in High Dimension Data Set International Journal of Applie Engineering Research ISSN 973-456 Volume, Number 6 (7) pp. 64-649 New Statistical Test for Quality Control in High Dimension Data Set Shamshuritawati Sharif, Suzilah Ismail

More information

Designing Information Devices and Systems I Spring 2018 Lecture Notes Note 16

Designing Information Devices and Systems I Spring 2018 Lecture Notes Note 16 EECS 16A Designing Information Devices an Systems I Spring 218 Lecture Notes Note 16 16.1 Touchscreen Revisite We ve seen how a resistive touchscreen works by using the concept of voltage iviers. Essentially,

More information

An inductance lookup table application for analysis of reluctance stepper motor model

An inductance lookup table application for analysis of reluctance stepper motor model ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR

More information

4. Important theorems in quantum mechanics

4. Important theorems in quantum mechanics TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this

More information

APPPHYS 217 Thursday 8 April 2010

APPPHYS 217 Thursday 8 April 2010 APPPHYS 7 Thursay 8 April A&M example 6: The ouble integrator Consier the motion of a point particle in D with the applie force as a control input This is simply Newton s equation F ma with F u : t q q

More information

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1 Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close

More information

Damage identification based on incomplete modal data and constrained nonlinear multivariable function

Damage identification based on incomplete modal data and constrained nonlinear multivariable function Journal of Physics: Conference Series PAPER OPEN ACCESS Damage ientification base on incomplete moal ata an constraine nonlinear multivariable function To cite this article: S S Kourehli 215 J. Phys.:

More information

Estimation of District Level Poor Households in the State of. Uttar Pradesh in India by Combining NSSO Survey and

Estimation of District Level Poor Households in the State of. Uttar Pradesh in India by Combining NSSO Survey and Int. Statistical Inst.: Proc. 58th Worl Statistical Congress, 2011, Dublin (Session CPS039) p.6567 Estimation of District Level Poor Househols in the State of Uttar Praesh in Inia by Combining NSSO Survey

More information

A Course in Machine Learning

A Course in Machine Learning A Course in Machine Learning Hal Daumé III 12 EFFICIENT LEARNING So far, our focus has been on moels of learning an basic algorithms for those moels. We have not place much emphasis on how to learn quickly.

More information

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables* 51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*

More information

Invariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem

Invariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem Invariant Extene Kalman Filter: Theory an application to a velocity-aie estimation problem S. Bonnabel (Mines ParisTech) Joint work with P. Martin (Mines ParisTech) E. Salaun (Georgia Institute of Technology)

More information

State estimation for predictive maintenance using Kalman filter

State estimation for predictive maintenance using Kalman filter Reliability Engineering an System Safety 66 (1999) 29 39 www.elsevier.com/locate/ress State estimation for preictive maintenance using Kalman filter S.K. Yang, T.S. Liu* Department of Mechanical Engineering,

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

CURRENT ELECTRICITY Q.1

CURRENT ELECTRICITY Q.1 CUENT EECTCTY Q. Define Electric current an its unit.. Electric Current t can be efine as the time rate of flow of charge in a conuctor is calle Electric Current. The amount of flow of charge Q per unit

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

Constraint Reformulation and a Lagrangian Relaxation based Solution Algorithm for a Least Expected Time Path Problem Abstract 1.

Constraint Reformulation and a Lagrangian Relaxation based Solution Algorithm for a Least Expected Time Path Problem Abstract 1. Constraint Reformulation an a Lagrangian Relaation base Solution Algorithm for a Least Epecte Time Path Problem Liing Yang State Key Laboratory of Rail Traffic Control an Safety, Being Jiaotong University,

More information

Time-Optimal Motion Control of Piezoelectric Actuator: STM Application

Time-Optimal Motion Control of Piezoelectric Actuator: STM Application Time-Optimal Motion Control of Piezoelectric Actuator: STM Application Yongai Xu, Peter H. Mecl Abstract This paper exaes the problem of time-optimal motion control in the context of Scanning Tunneling

More information

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation

More information

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 13. Khyruddin Akbar Ansari, Ph.D., P.E.

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 13. Khyruddin Akbar Ansari, Ph.D., P.E. AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 13 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering Gonzaga University SDC PUBLICATIONS Schroff

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

Distributed coordination control for multi-robot networks using Lyapunov-like barrier functions

Distributed coordination control for multi-robot networks using Lyapunov-like barrier functions IEEE TRANSACTIONS ON 1 Distribute coorination control for multi-robot networks using Lyapunov-like barrier functions Dimitra Panagou, Dušan M. Stipanović an Petros G. Voulgaris Abstract This paper aresses

More information

An Introduction to Event-triggered and Self-triggered Control

An Introduction to Event-triggered and Self-triggered Control An Introuction to Event-triggere an Self-triggere Control W.P.M.H. Heemels K.H. Johansson P. Tabuaa Abstract Recent evelopments in computer an communication technologies have le to a new type of large-scale

More information

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique J. Basic. Appl. Sci. Res., ()9-99,, TextRoa Publication ISSN 9-434 Journal of Basic an Applie Scientific Research www.textroa.com A Comparison between a Conventional Power System Stabilizer (PSS) an Novel

More information

CONTROL CHARTS FOR VARIABLES

CONTROL CHARTS FOR VARIABLES UNIT CONTOL CHATS FO VAIABLES Structure.1 Introuction Objectives. Control Chart Technique.3 Control Charts for Variables.4 Control Chart for Mean(-Chart).5 ange Chart (-Chart).6 Stanar Deviation Chart

More information

Linear and quadratic approximation

Linear and quadratic approximation Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Relative Position Sensing by Fusing Monocular Vision and Inertial Rate Sensors

Relative Position Sensing by Fusing Monocular Vision and Inertial Rate Sensors Proceeings of ICAR 2003 The 11th International Conference on Avance Robotics Coimbra, Portugal, June 30 - July 3, 2003 Relative Position Sensing by Fusing Monocular Vision an Inertial Rate Sensors Anreas

More information

On Using Unstable Electrohydraulic Valves for Control

On Using Unstable Electrohydraulic Valves for Control Kailash Krishnaswamy Perry Y. Li Department of Mechanical Engineering, University of Minnesota, 111 Church St. SE, Minneapolis, MN 55455 e-mail: kk,pli @me.umn.eu On Using Unstable Electrohyraulic Valves

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.

More information