Derivative-Free Kalman Filtering-Based Control of Prosthetic Legs*

Size: px
Start display at page:

Download "Derivative-Free Kalman Filtering-Based Control of Prosthetic Legs*"

Transcription

1 7 Amercan Control Conerence Sheraton Seattle Hotel May 4 6, 7, Seattle, USA Dervatve-Free Kalman Flterng-Based Control o Prosthetc Legs* S Mahmoud Moosav, Seyed Abolazl Fakooran, Vahd Azm, Hanz Rchter and Dan Smon Abstract A dervatve-ree method or state estmatonbased control o a robot/prosthess system s presented The system s the combnaton o a test robot that emulates human hp and thgh moton, and a powered transemoral prosthetc leg The robot/prosthess combnaton s modeled as a three degreeo-reedom (DOF) robot: vertcal hp dsplacement, thgh angle, and knee angle We develop a dervatve-ree Kalman lter (DKF) or state estmaton-based control or an n-dof robotc system We then propose a method to make the DKF robust when the robot dynamcs nclude dsturbances In the robust DKF, we use two derent methods or dsturbance rejecton: PD and PI These dsturbance compensators are used or supervsory control to make the DKF robust n the presence o dsturbances The smulaton results show the advantages o the DKF and the robust DKF or the three-dof robot/prosthess system or state estmaton-based control I INTRODUCTION The desgn and development o transemoral leg prostheses has receved consderable attenton due to the ncreasng number o above-knee amputees Recent advances n mcroelectroncs and robotc technologes have led to the development o new powered prosthetc legs Transemoral amputees who use powered (actve) prosthetc legs expend less energy compared wth those who use passve prostheses snce powered prostheses are able to generate net power at the jonts [] Some recent research on the development o lower lmb prostheses wth powered jonts ncludes [], [3], [4], [] Robotc testng n the sagttal plane o transemoral prostheses s presented n [6], [7], where nonlnear control technques are proposed on a test robot combned wth a prosthetc leg Passvty-based robust control, robust adaptve mpedance control, hybrd control, and robust composte adaptve mpedance control are some recent attempts or the control o a three degree-o-reedom robot/prosthess system [7], [8], [9], [] The proposed methods have demonstrated good trackng perormance n the presence o parameter uncertantes and ground reacton orce (GRF) as an external dsturbance However, GRF s treated n the aorementoned research as an nput wth known dynamc propertes and known bounds Thus, trackng perormance may break down the system s aected by dsturbances In ths paper we propose observer-based prosthess controllers n order to address ths shortcomng Note that unknown *Ths work was supported by NSF Grant Department o Electrcal Engneerng and Computer Scence, Cleveland State Unversty, Cleveland, OH, USA Correspondng author e-mal address: smoosav7@csuohoedu Department o Mechancal Engneerng, Cleveland State Unversty, Cleveland, OH, USA nputs can also be consdered as dsturbances n the context o ths paper A dervatve-ree Kalman lter (DKF) s a new method or estmaton-based control loop whch uses exact eedback lnearzaton to control nonlnear systems [] The DKF s adapted and mplemented n ths paper or estmatonbased control o a class o nput-output lnearzable nonlnear systems The nonlnear system s lnearzed and controlled by exact eedback lnearzaton, and then the standard Kalman lter s appled to the lnearzed model to estmate the state vector However, unlke the extended Kalman lter (EKF), the DKF provdes state estmaton o the nonlnear system wthout on-lne dervatve and Jacoban calculatons The eedback lnearzaton step s nstead part o the problem ormulaton, and so on-lne dervatve calculatons are not necessary Thereore, the problems o the EKF whch result rom the local lnearzaton o the nonlnear system by rst-order Taylor seres expansons can be avoded n the DKF These EKF problems can aect the accuracy o the state estmaton and consequently the stablty o the state estmaton-based controller, and are allevated wth the DKF The man contrbuton o ths paper s the desgn o a DKF or a leg prosthess test robot that s used or transemoral prosthess desgn and test Ths research s the rst attempt or closed-loop estmaton-based control or the robot/prosthess system usng a DKF We rst consder an n-dof robot and desgn a DKF when the robot s model s known n parameters, nputs, and GRF dynamcs Then we examne the case n whch the robot model s perturbed wth dsturbances as parametrc uncertantes or GRF In both cases nose terms are modeled n an ane state space model or the robot The dsturbances are assumed to be bounded but ther bounds and dynamcs are unknown System perormance may degrade when the system s perturbed by dsturbances snce the DKF s based on eedback lnearzaton and the accuracy o exact eedback lnearzaton deterorates under dsturbances Thereore, the new contrbutons n ths paper nclude the use o two methods, namely PD and PI dsturbance compensators, as supervsory control terms or the rejecton o dsturbances The use o these methods or dsturbance rejecton n the DKF s proposed or the rst tme n ths paper These methods combne wth the DKF to ncrease robustness The combned DKF and dsturbance compensator s called the robust DKF snce t can tolerate dsturbances or estmaton-based control We consder a smulated three-dof prosthess test robot to show the results o the DKF; the DKF s used or onlne estmaton-based control o hp dsplacement, thgh angle, knee angle, and correspondng veloctes /$3 7 AACC

2 The paper s organzed as ollows In Secton II the DKF s desgned or estmaton-based control or general robotc systems, under some assumptons In Secton III methods or dsturbance rejecton are presented to make the DKF more robust n the presence o dsturbances Secton IV provdes smulaton results o the DKF on a three-dof robot/prosthess system and analyzes the perormance o the DKF n the presence o dsturbances Secton V concludes the paper and suggests uture work II DERIVATIVE-FREE KALMAN FILTERING FOR ROBOTS The Kalman lter apples drectly only to lnear systems However, we can lnearze a nonlnear system and then use lnear estmaton technques Applyng the standard Kalman lter to a nonlnear system through the transormaton o the nonlnear system to the observer canoncal orm s called dervatve-ree Kalman lterng (DKF) The lnearzaton transormaton o the nonlnear system s based on a deomorphsm and, unlke the EKF, does not nvolve the computaton o Jacobans The DKF can be used or estmaton-based control o robots the robot model s subjected to a lnearzaton transormaton or exact eedback lnearzaton control and then state estmaton s perormed on the lnearzed model wth a standard Kalman lter A dynamc model or a mult-nput mult-output (MIMO) n-dof robot s gven as M (q) q +C (q, q) q + J e T F e + g(q) = u () where q s the vector o jont dsplacements, M (q) s the mass matrx, C (q, q) s a matrx accountng or centrpetal and Corols eects, J e s the knematc Jacoban relatve to the pont o the applcaton o external orces F e, g(q) s the gravty vector, and u s the vector o control sgnals [6], [7], [], [] The dynamc model () can be wrtten n ane state space orm whle also ncludng nose terms n the process and measurement equatons as ollows: ẋ = (x,t) + g(x,t)u + ω(t) y = h(x,t) + ϖ(t) () where state vector x(t), nput u(t), and output y(t) are n R n, R p and R m respectvely; and the smooth vector elds () and g() are n R n and R n p respectvely The nose terms ω(t) and ϖ(t) are assumed to be bandlmted, zero-mean, and uncorrelated Note that the process nose s partally due to parameter uncertantes, unmodeled dynamcs, and nput uncertantes n the robot model, and the measurement nose corresponds to errors n the measurement equpment The nonlnear system () s subjected to a lnearzaton transormaton to obtan an observer canoncal orm n order to perorm dervatve-ree Kalman lterng (DKF) [] In ths secton the nput-output eedback lnearzaton s presented to lnearze the model n () Snce t s desred to obtan an nvertble decouplng matrx n the lnearzaton transormaton, we assume that the number o nputs s equal to the number o the outputs; that s, p = m The robot model s nput-output lnearzable under the condtons gven n [3, Chapter 6] The classcal nput-output lnearzaton method can be appled by derentatng each output uncton y ( =,,,m) untl the nputs appear Remark: We assume that the lat output o the system [4, Chapter ] s a lnear uncton o the state vector elements In ths case the appled transormaton does not have nonlnear eects on the nose sgnals We also assume that we know the mean and covarance o the nose so that we can calculate the mean and covarance o the derentated nose [, Chapter 9] Assume that r s the smallest nteger such that at least one o the nputs appears n y (r ) ; then y (r ) (t) = L r h (x) + m j= L g j L r h (x)u j (t) + γ (t) (3) where L h s the Le dervatve L h = ( h), wth L g j L r h It can be shown that γ (t) s nose wth zero mean: γ (t) = ω (t) + d(r ) ϖ (t) (4) dt Note that one o the man advantages o transormng the robot dynamcs () nto the ane state-space orm () s to allow treatment o the nose terms n the deomorphsm Based on (3) and exact eedback lnearzaton, the system model () (gnorng nose or now) can be transormed nto the ollowng orm [3, Chapter 6]: v (t) L r h (x) = + Γ(x)u(t) () v m (t) L r m h m (x) where Γ(x) s the m m nvertble decouplng matrx and the scalar r = r + r + + r m s called the total relatve degree o the system The vrtual control term v (t) wll be desgned later n ths paper The partal relatve degrees r are called well-dened L g j L r h In ths case, the nput transormaton v (t) L r u(t) = Γ (x) h (x) v m (t) L r m h m (x) (6) yelds m smple lnear systems rom (3) n new coordnates: Ż (t) = Z (t) Ż r (t) = y (r ) (t) = v (t) + γ (t) (7) or =,,,m Now we can wrte each decoupled subsystem n state space orm: Ż (t) = A Z (t) + B v (t) + M γ (t) y (t) = C Z (t) + ϖ (t) (8) 6

3 We can then wrte each state space equaton n canoncal orm as ollows: A Ż (t) {}}{ Z Ż (t) (t) {}}{{}}{ = Z (t) + v (t) + γ (t) Żr (t) Zr (t) y (t) = [ ] Z (t) + ϖ (t) (9) }{{} C Ater nput-output lnearzaton and model transormaton or each subsystem, a lnear control law (vrtual control) v can be desgned to control each transormed subsystem o (9) The controller s am s to make the system s output ollow a gven desred trajectory Z d The vrtual control s obtaned or each subsystem as v (t) = Z (r ) d (t) θ T e (t) () where the gan matrx θ T = [ θr,,θ] s desgned by pole placement and the error dynamcs are dened on the bass o the states e (t) = Ẑ Z d wth e (t) = [ e,ė,,e (r ) ] T such that the polynomal e (r ) B M +θ e(r ) + + θ r e s Hurwtz To obtan the estmate Ẑ o the state vector, we apply a standard Kalman lter to the lnearzed model o the robot (9) as ollows [6, Chapter 8]: Ẑ () = E[Z ()] [ (Z P () = E () Ẑ () )( Z () Ẑ () ) ] T K (t) = P (t)c T R (t) Ẑ (t) = A Ẑ (t) + B v (t) + K (t)(y (t) C Ẑ (t)) Ṗ (t) = A P (t) + P (t)a T (t) + Q (t) P (t)c T R (t)c P (t) () where K (t) and P (t) are the Kalman gan and estmaton error covarance matrx; and the covarances o the process and measurement nose are Q (t) = E [ γ (t)γ T (t) ] and R (t) = E [ ϖ (t)ϖ T (t) ], respectvely We see that a Kalman lter can be separately desgned or each decoupled subsystem, =,,,m These two steps transormaton by eedback lnearzaton, and state estmaton by standard Kalman lters comprse the DKF Remark: Measured values o all states are requred n eedback lnearzaton However, the Kalman lter can be used to estmate unmeasured states n the transormed system Thus, n the DKF, we can use the states to obtan the control nput u(t) III DERIVATIVE-FREE KALMAN FILTERING IN THE PRESENCE OF DISTURBANCES Parameter uncertantes and dsturbances aect the robot Parameter uncertantes are oten modeled as process nose (as dscussed n the prevous secton), whle determnstc dsturbances are modeled as external nputs In robotc systems, dsturbances can be an order o magntude or more greater than the process nose The DKF typcally handles process nose well, but s not desgned to handle dsturbances In ths secton we extend the DKF to handle dsturbances The dsturbances can arse due to the ollowng [4]: () unknown torques on the robot jonts, () unknown orces on the robot (eg, rcton) In the DKF we use eedback lnearzaton to lnearze the nonlnear model o the robot to apply the standard Kalman lter However, the perormance o the DKF may break down snce eedback lnearzaton s not robust to dsturbances The robot s dynamc model n ane state space orm () wth dsturbances can be wrtten as ẋ = (x,t) + g(x,t)u + ω(t) + E(x,t)d(x,t) y = h(x,t) + ϖ(t) () where d(x,t) R l s the dsturbance vector and E(x,t)d(x,t) comprses the eect o the dsturbances on each jont, where E(x,t) R n l can be calculated rom the dynamc model () It should be noted that the dsturbance and ts dervatves must be bounded or stablty reasons, but we do not assume any pror bounds and we do not assume any normaton or the dynamcs o the dsturbances We use the same approach here or the transormaton lnearzaton o the nonlnear model () as we dd n Secton II Exact eedback lnearzaton s used and the nput transormaton s obtaned as n (6), whch results n the transormaton o () nto the ollowng m decoupled subsystems: Ż (t) = Z (t) Ż r (t) = y (r ) (t) = v (t) + γ (t) + d (x,t) (3) or =,,,m Note that we stll have the assumpton about the number o nputs p beng equal to the number o outputs m, as n Secton II, or the exstence o u The components ( d, d,) whch represent the transormed eects o the dsturbances can be obtaned as ollows [7]: d (x,t) = d m (x,t) L r E d h + + dt d (r ) (LE d h ) L r m E d h m + + dt d (4) (r m ) (LE d h m ) The state space canoncal orm can be wrtten or every subsystem rom the transormed model (3) and the bass o (9): Ż (t) = A Z (t) + B v (t) + M γ (t) + N d (x,t) y (t) = C Z (t) + ϖ (t) () or =,,,m, where N = M In order to acheve state estmaton-based control o the robot wth the DKF, we apply the standard Kalman lter () to the lnearzed model () However, the DKF s not able to estmate the states o the robot n the presence o dsturbances So we need to desgn a vrtual control v as 7

4 n () wth dsturbance rejecton to stablze the transormed system dynamcs v (t) = Z (r ) d (t) θ T e (t) ˆ d (x,t) (6) where ˆ d (x,t) s the estmate o the dsturbance whch s used or the rejecton/compensaton o the dsturbance We propose two derent methods to obtan ˆ d (x,t) ) Use a PD dsturbance compensator or dsturbance rejecton; ˆ d (x,t) = L T E, wth the gan L T = [L P,L D ] and E = [e,ė ] T ) Use a PI dsturbance compensator or dsturbance rejecton; ˆ d (x,t) = J T Ξ, wth the gan J T = [J P,J I ] and Ξ = [ e, t e(τ)dτ] T The vrtual control v (6) ams to make the DKF robust aganst dsturbances The DKF or state estmaton and the proposed methods or dsturbance rejecton results n the augmented system [ ] [ ][ ] [ ] [ ] Ẑ (t) A N = Ẑ (t) B K (t)(y(t) C + v ˆ d (x,t) ˆ d (x,t) + Ẑ (t)) L T E or J T Ξ (7) We dene the state estmaton errors as ξ (t) = Ẑ (t) Z (t) and the dsturbance estmaton errors as ϕ(t) = ˆ d (x,t) d (x,t) I L P L D or J P J I, the estmaton errors ξ (s) F or ϕ(s) F converge to a postve small value ε s d (s) F s bounded, where F represents the Frobenus norm and s s the Laplace transorm varable The proposed method or dsturbance rejecton and the DKF or estmatonbased control o the robot s depcted n Fg IV SIMULATION RESULTS In ths secton, we show the eectveness o the DKF and the robust DKF or state estmaton-based control usng smulaton studes o the three-dof prosthess test robot Fg shows a schematc o the hp robot and prosthess combnaton, whch s modeled as a three-lnk robot The dynamc model s gven by (), and the vector o generalzed jont dsplacements s q = [q,q,q 3 ] T, where q s vertcal hp dsplacement, q s thgh angle, and q 3 s knee angle The state vector and outputs are gven as x = [q, q,q, q,q 3, q 3 ] T, y = [q,q,q 3 ] T We thus have three states or the generalzed dsplacements and three states or the generalzed veloctes Moreover, the number o nputs u = [u,u,u 3 ] T s the same as the number o outputs In the dynamcs o the prosthetc leg robot, the nput u wll appear ater two dervatves o the correspondng output y, =,,3 So the total relatve degree s r = r + r + r 3 = 6 as gven n (3), whch s equal to the system order n It should be noted that there are no nternal dynamcs assocated wth the nput-output lnearzaton o the robot/prosthess system snce r = n We assume that the robot model s known except or process nose In Subsecton A we assume that the GRF s modeled as n [8], whch results n the external orce vector F e n (), and that the GRF model s known to the controller In Subsecton B we make no assumptons about the controller s knowledge o the GRF or ts model, and we treat GRF as dsturbances We consder the perormance o the robot/prosthess system durng our steps o walkng at a normal speed, whch s approxmately our seconds The reerence data were provded by research partcpants at the Moton Studes Laboratory o the Cleveland Department o Veterans Aars Medcal Center [8] A The robot/prosthess system wthout dsturbances: As mentoned earler, the robot dynamcs o () s nput-output lnearzable Based on the assumpton o the equalty o the number o nputs and outputs, the robot dynamcs can be transormed nto three decoupled subsystems as [ ] [ ] [ ] Ż (t) = Z (t) + v + γ (t) y (t) = [ ] Z (t) + ϖ (t) (8) or =,,3, where Z (t) = [q, q ] T The unts o Z are m and m/s, the unts o Z are rad and rad/sec, and the unts o Z 3 are also rad and rad/sec We apply the DKF to each lnearzed subsystem (8) to obtan state estmaton-based control The covarances o the process and measurement nose are consdered as Q = 3 and R = 3 or =,, 3 These covarances are based on pror experence wth the accuracy o system dynamc modelng and measurement uncertanty Nose Dsturbances Exact Feedback lnearzaton Inputs u = [u,, u m] Robot model Outputs y (r) = [y (r),, y (rm) ] m decoupled subsystems through nput-output lnearzed model v = [v,, v m] Vrtual control term Dsturbances rejecton d Desred trajectores PD or PI dsturbance compensator State estmaton DKF Fg Block dagram o the proposed methods or dsturbance rejecton and or state estmaton-based control o the robotc system, (robust DKF) Fg Dagram o the robot/prosthess system wth rgd ankle jont 8

5 The ntal value o the state vector s gven as Z () = [9,93] T, Z () = [3,779] T, Z 3 () = [9,4] T We set the ntal value o the state vector to provde an arbtrary but nonzero ntal estmaton error: Ẑ () = [39,6] T, Ẑ () = [9,3] T, Ẑ 3 () = [3,884] T The gans θ n () are desgned by pole placement to provde good estmaton-based control perormance and are obtaned as θ T = [,] or =,,3 The results or the state estmaton-based control o the robot/prosthess system are shown n Fg 3 Although sgncant ntal estmaton errors are mposed on the jont dsplacements and veloctes, the DKF quckly converges to the states B The robot/prosthess system wth dsturbances: In the robot/prosthess system, parameters may change sgncantly rom ther nomnal values Also, the external orces F e nclude unmodeled GRF whch can be consdered as a dsturbance sensors are not avalable We consder these two sources o uncertanty as dsturbances d(x,t) n the robot dynamcs We assume that the ndependent dsturbances d,d,d 3 are ncluded n the robot dynamcs The dsturbances that we use n our smulaton are d = sn(t) N, d = N, and d 3 = sn(t + π) N These nputs aect each jont o the robot as gven by the transormaton lnearzaton (4); the transormed nputs are denoted as d, d, d 3 n the lnearzed model () The robot dynamcs are transormed nto three decoupled subsystems as shown n () The DKF can be desgned usng dsturbance rejecton as shown n Secton III The ntal state vector and estmate, and the covarances o the nose terms, are dentcal wth those n Subsecton A above In the ollowng, we rst use the PD compensator to estmate d or dsturbance rejecton, and then we use the PI compensator Fg 4 shows the perormance o the robust DKF usng PD compensaton or state estmaton-based control o the robot n the presence o these large dsturbances The dsturbances are as a sde benet o ths approach (see Fg (a)) We next use the PI compensator or dsturbance rejecton Estmaton results are shown n Fg (b) Table I compares the PD and PI compensators n terms o root mean square error (RMSE) o the state estmates, and RMS values o the control sgnals Table I shows that the PD compensator outperorms the PI compensator n the ollowng ways: () The PD compensator has aster dynamcs than the PI compensator due to the error dervatve term ė; ths provdes predcton whch mproves the compensator s bandwdth or dsturbance estmaton () The ntegral term n the PI compensator nvolves a pole at the orgn whch results not only n overshoot durng transents, but also slowness n dsturbance estmaton Snce the system has been lnearzed based on eedback lnearzaton, a pole at the orgn may cause oscllaton or nstablty n some cases V CONCLUSION AND FUTURE WORK We developed and appled a DKF to robotc systems or estmaton-based control We proposed a PD or PI dsturbance compensator as a supervsory control term or hp dsplasment (m) thgh angle (rad) knee angle (rad) 4 3 tme tme tme hp velocty (m/s) 4 4 angular velocty o thgh (rad/s) angular velocty o knee (rad/s) 4 4 Fg 3 State estmaton-based control o the robot/prosthess system usng dervatve-ree Kalman lterng n the absence o dsturbances hp dsplasment (m) thgh angle (rad) knee angle (rad) 4 3 tme tme tme hp velocty (m/s) 4 4 angular velocty o thgh (rad/s) angular velocty o knee (rad/s) 4 4 Fg 4 State estmaton-based control o the robot/prosthess system n the presence o dsturbances usng robust dervatve-ree Kalman lterng wth PD dsturbance compensaton 9

6 TABLE I RMSE AND RMS VALUES OF CONTROL SIGNALS FOR THE ROBUST DKF: COMPARISON BETWEEN PD DISTURBANCE COMPENSATOR AND PI x (m) x (m/s) DISTURBANCE COMPENSATOR x 3 (rad) x 4 (rad/s) x (rad) x 6 (rad/s) u (N) u (Nm) PD compensator PI compensator u 3 (Nm) dstrubance/unkown nput estmaton dstrubance/unkown nput estmaton dstrubance/unkown nput estmaton (a) dstrubance/unkown nput estmaton dstrubance/unkown nput estmaton dstrubance/unkown nput estmaton Fg (a) Estmates o dsturbances usng PD compensaton and (b) estmates o dsturbances usng PI compensaton dsturbance rejecton, whch makes the DKF robust A three- DOF robot/prosthess system or transemoral amputees was smulated to demonstrate the DKF or estmaton-based control o generalzed jont dsplacements and veloctes Results show that the DKF estmates quckly converge to the states, even wth sgncant nose and ntal condton errors In addton, we presented robust state estmatonbased control n the presence o dsturbances We showed that the robust DKF acheves smaller estmaton errors wth PD compensaton than wth PI compensaton For uture work, we wll analytcally prove the stablty o the DKF or the robot/prosthess system The DKF can be appled to an underactuated robot as well, whch wll requre an nvestgaton o the stablty o the nternal dynamcs Future work wll nclude expermental hardware mplementaton and human trals (b) [] Z Harvey, Z Benjamn, J Vandersea, and E Wol, Prosthetc advances, Journal o Surgcal Orthopedc Advances, vol, no, pp 8 64, [3] H Zhao, S Kolathaya, and A D Ames, Quadratc programmng and mpedance control or transemoral prosthess, n Internatonal Conerence on Robotcs and Automaton (ICRA), (Hong Kong), pp , 4 [4] R Gregg, T Lenz, L Hargrove, and J Sensnger, Vrtual constrant control o a powered prosthetc leg: rom smulaton to experments wth transemoral amputees, IEEE Transactons on Robotcs, vol 3, no 6, pp 4 47, 4 [] A J Young, A M Smon, and L J Hargrove, A tranng method or locomoton mode predcton usng powered lower lmb prostheses, IEEE Transactons on Neural Systems and Rehabltaton Engneerng, vol, no 3, pp , 4 [6] H Rchter, D Smon, W Smth, and S Samorezov, Dynamc modelng, parameter estmaton and control o a leg prosthess test robot, Appled Mathematcal Modellng, vol 39, no, pp 9 73, [7] H Rchter and D Smon, Robust trackng control o a prosthess test robot, Journal o Dynamc Systems, Measurement, and Control, vol 36, no 3, p 3, 4 [8] V Azm, D Smon, and H Rchter, Stable robust adaptve mpedance control o a prosthetc leg, n Proc o the ASME Dynamc Systems and Control Conerence, (Columbus, Oho), pp VT9A3 VT9A3, [9] D Ebegbe, D Smon, and H Rchter, Hybrd uncton approxmaton based control wth applcaton to prosthetc legs, n IEEE Int Systems Conerence, (Orlando, Florda), pp 6, 6 [] V Azm, D Smon, H Rchter, and S A Fakooran, Robust composte adaptve transemoral prosthess control wth non-scalar boundary layer trajectores, n Proc Amercan Control Conerence, (Boston, Massachusetts), 6 [] G G Rgatos, A dervatve-ree Kalman lterng approach to state estmaton-based control o nonlnear systems, IEEE Trans Ind Electron, vol 9, no, pp , [] S A Fakooran, D Smon, H Rchter, and V Azm, Ground reacton orce estmaton n prosthetc legs wth an extended Kalman lter, n IEEE Int Systems Conerence, (Orlando, Florda), pp , 6 [3] J-J E Slotne and W L, Appled Nonlnear Control Prentce-Hall Englewood Cls, NJ, 99 [4] G G Rgatos, Nonlnear Control and Flterng Usng Derental Flatness Approaches: Applcatons to Electromechancal Systems Sprnger, [] E C Rasmussen and C K I Wllams, Gaussan Processes or Machne Learnng The MIT press, 6 [6] D Smon, Optmal State Estmaton: Kalman, H-Innty, and Nonlnear Approaches John Wley & Sons, 6 [7] Y Lu and D Söker, Robust control approach or nputoutput lnearzable nonlnear systems usng hgh-gan dsturbance observer, Internatonal Journal o Robust and Nonlnear Control, vol 4, no, pp , 4 [8] G Khadem, H Mohammad, E C Hardn, and D Smon, Evolutonary optmzaton o user ntent recognton or transemoral amputees, n Proc Bomedcal Crcuts and Systems Conerence, (Atlanta, Georga), pp 4, REFERENCES [] F Sub and H A Varol, Upslope walkng wth a powered knee and ankle prosthess: Intal results wth an amputee subject, IEEE Trans Neural Systems and Rehabltaton Engneerng, vol 9, no, pp 7 78,

Chapter 3 Differentiation and Integration

Chapter 3 Differentiation and Integration MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton

More information

, rst we solve te PDE's L ad L ad n g g (x) = ; = ; ; ; n () (x) = () Ten, we nd te uncton (x), te lnearzng eedbac and coordnates transormaton are gve

, rst we solve te PDE's L ad L ad n g g (x) = ; = ; ; ; n () (x) = () Ten, we nd te uncton (x), te lnearzng eedbac and coordnates transormaton are gve Freedom n Coordnates Transormaton or Exact Lnearzaton and ts Applcaton to Transent Beavor Improvement Kenj Fujmoto and Tosaru Suge Dvson o Appled Systems Scence, Kyoto Unversty, Uj, Kyoto, Japan suge@robotuassyoto-uacjp

More information

Parameter Estimation for Dynamic System using Unscented Kalman filter

Parameter Estimation for Dynamic System using Unscented Kalman filter Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,

More information

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI 2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Prof. Paolo Colantonio a.a

Prof. Paolo Colantonio a.a Pro. Paolo olantono a.a. 3 4 Let s consder a two ports network o Two ports Network o L For passve network (.e. wthout nternal sources or actve devces), a general representaton can be made by a sutable

More information

Trajectory Optimization of Flexible Mobile Manipulators Using Open-Loop Optimal Control method

Trajectory Optimization of Flexible Mobile Manipulators Using Open-Loop Optimal Control method rajectory Optmzaton o Flexble Moble Manpulators Usng Open-Loop Optmal Control method M.H. Korayem and H. Rahm Nohooj Robotc Lab, Mechancal Department, Iran Unversty o Scence and echnology hkorayem@ust.ac.r,

More information

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming OPTIMIATION Introducton ngle Varable Unconstraned Optmsaton Multvarable Unconstraned Optmsaton Lnear Programmng Chapter Optmsaton /. Introducton In an engneerng analss, sometmes etremtes, ether mnmum or

More information

EE 330 Lecture 24. Small Signal Analysis Small Signal Analysis of BJT Amplifier

EE 330 Lecture 24. Small Signal Analysis Small Signal Analysis of BJT Amplifier EE 0 Lecture 4 Small Sgnal Analss Small Sgnal Analss o BJT Ampler Eam Frda March 9 Eam Frda Aprl Revew Sesson or Eam : 6:00 p.m. on Thursda March 8 n Room Sweene 6 Revew rom Last Lecture Comparson o Gans

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to ChE Lecture Notes - D. Keer, 5/9/98 Lecture 6,7,8 - Rootndng n systems o equatons (A) Theory (B) Problems (C) MATLAB Applcatons Tet: Supplementary notes rom Instructor 6. Why s t mportant to be able to

More information

Work is the change in energy of a system (neglecting heat transfer). To examine what could

Work is the change in energy of a system (neglecting heat transfer). To examine what could Work Work s the change n energy o a system (neglectng heat transer). To eamne what could cause work, let s look at the dmensons o energy: L ML E M L F L so T T dmensonally energy s equal to a orce tmes

More information

Force = F Piston area = A

Force = F Piston area = A CHAPTER III Ths chapter s an mportant transton between the propertes o pure substances and the most mportant chapter whch s: the rst law o thermodynamcs In ths chapter, we wll ntroduce the notons o heat,

More information

Spring Force and Power

Spring Force and Power Lecture 13 Chapter 9 Sprng Force and Power Yeah, energy s better than orces. What s net? Course webste: http://aculty.uml.edu/andry_danylov/teachng/physcsi IN THIS CHAPTER, you wll learn how to solve problems

More information

Chapter 6. Operational Amplifier. inputs can be defined as the average of the sum of the two signals.

Chapter 6. Operational Amplifier.  inputs can be defined as the average of the sum of the two signals. 6 Operatonal mpler Chapter 6 Operatonal mpler CC Symbol: nput nput Output EE () Non-nvertng termnal, () nvertng termnal nput mpedance : Few mega (ery hgh), Output mpedance : Less than (ery low) Derental

More information

Modeling of Dynamic Systems

Modeling of Dynamic Systems Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how

More information

Department of Electrical and Computer Engineering FEEDBACK AMPLIFIERS

Department of Electrical and Computer Engineering FEEDBACK AMPLIFIERS Department o Electrcal and Computer Engneerng UNIT I EII FEEDBCK MPLIFIES porton the output sgnal s ed back to the nput o the ampler s called Feedback mpler. Feedback Concept: block dagram o an ampler

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

Lecture 2 Solution of Nonlinear Equations ( Root Finding Problems )

Lecture 2 Solution of Nonlinear Equations ( Root Finding Problems ) Lecture Soluton o Nonlnear Equatons Root Fndng Problems Dentons Classcaton o Methods Analytcal Solutons Graphcal Methods Numercal Methods Bracketng Methods Open Methods Convergence Notatons Root Fndng

More information

Note 10. Modeling and Simulation of Dynamic Systems

Note 10. Modeling and Simulation of Dynamic Systems Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada

More information

Physics 2A Chapter 3 HW Solutions

Physics 2A Chapter 3 HW Solutions Phscs A Chapter 3 HW Solutons Chapter 3 Conceptual Queston: 4, 6, 8, Problems: 5,, 8, 7, 3, 44, 46, 69, 70, 73 Q3.4. Reason: (a) C = A+ B onl A and B are n the same drecton. Sze does not matter. (b) C

More information

Neuro-Adaptive Design II:

Neuro-Adaptive Design II: Lecture 37 Neuro-Adaptve Desgn II: A Robustfyng Tool for Any Desgn Dr. Radhakant Padh Asst. Professor Dept. of Aerospace Engneerng Indan Insttute of Scence - Bangalore Motvaton Perfect system modelng s

More information

Tracking with Kalman Filter

Tracking with Kalman Filter Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,

More information

FEEDBACK AMPLIFIERS. v i or v s v 0

FEEDBACK AMPLIFIERS. v i or v s v 0 FEEDBCK MPLIFIERS Feedback n mplers FEEDBCK IS THE PROCESS OF FEEDING FRCTION OF OUTPUT ENERGY (VOLTGE OR CURRENT) BCK TO THE INPUT CIRCUIT. THE CIRCUIT EMPLOYED FOR THIS PURPOSE IS CLLED FEEDBCK NETWORK.

More information

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 NMT EE 589 & UNM ME 48/58 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 48/58 7. Robot Dynamcs 7.5 The Equatons of Moton Gven that we wsh to fnd the path q(t (n jont space) whch mnmzes the energy

More information

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN Int. J. Chem. Sc.: (4), 04, 645654 ISSN 097768X www.sadgurupublcatons.com COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN R. GOVINDARASU a, R. PARTHIBAN a and P. K. BHABA b* a Department

More information

11. Dynamics in Rotating Frames of Reference

11. Dynamics in Rotating Frames of Reference Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

A NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT

A NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT Internatonal Conference Mathematcal and Computatonal ology 0 Internatonal Journal of Modern Physcs: Conference Seres Vol. 9 0 68 75 World Scentfc Publshng Company DOI: 0.4/S009450059 A NUMERICAL COMPARISON

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

A Comparison Study of Advanced State Observer Design Techniques

A Comparison Study of Advanced State Observer Design Techniques A Comparson Study of Advanced State Observer Desgn Technques Wewen Wang and Zhqang Gao The Appled Control Research Laboratory Department of Electrcal and Computer Engneerng Cleveland State Unversty, Cleveland

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

ECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals

ECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals ECEN 5005 Crystals, Nanocrystals and Devce Applcatons Class 9 Group Theory For Crystals Dee Dagram Radatve Transton Probablty Wgner-Ecart Theorem Selecton Rule Dee Dagram Expermentally determned energy

More information

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t 8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes

More information

Chapter 3 and Chapter 4

Chapter 3 and Chapter 4 Chapter 3 and Chapter 4 Chapter 3 Energy 3. Introducton:Work Work W s energy transerred to or rom an object by means o a orce actng on the object. Energy transerred to the object s postve work, and energy

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method Stablty Analyss A. Khak Sedgh Control Systems Group Faculty of Electrcal and Computer Engneerng K. N. Toos Unversty of Technology February 2009 1 Introducton Stablty s the most promnent characterstc of

More information

Physics 181. Particle Systems

Physics 181. Particle Systems Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system

More information

: Numerical Analysis Topic 2: Solution of Nonlinear Equations Lectures 5-11:

: Numerical Analysis Topic 2: Solution of Nonlinear Equations Lectures 5-11: 764: Numercal Analyss Topc : Soluton o Nonlnear Equatons Lectures 5-: UIN Malang Read Chapters 5 and 6 o the tetbook 764_Topc Lecture 5 Soluton o Nonlnear Equatons Root Fndng Problems Dentons Classcaton

More information

Controller Design of High Order Nonholonomic System with Nonlinear Drifts

Controller Design of High Order Nonholonomic System with Nonlinear Drifts Internatonal Journal of Automaton and Computng 6(3, August 9, 4-44 DOI:.7/s633-9-4- Controller Desgn of Hgh Order Nonholonomc System wth Nonlnear Drfts Xu-Yun Zheng Yu-Qang Wu Research Insttute of Automaton,

More information

CHAPTER 4d. ROOTS OF EQUATIONS

CHAPTER 4d. ROOTS OF EQUATIONS CHAPTER 4d. ROOTS OF EQUATIONS A. J. Clark School o Engneerng Department o Cvl and Envronmental Engneerng by Dr. Ibrahm A. Assakka Sprng 00 ENCE 03 - Computaton Methods n Cvl Engneerng II Department o

More information

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites 7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT

More information

Unknown input extended Kalman filter-based fault diagnosis for satellite actuator

Unknown input extended Kalman filter-based fault diagnosis for satellite actuator Internatonal Conference on Computer and Automaton Engneerng (ICCAE ) IPCSI vol 44 () () IACSI Press, Sngapore DOI: 776/IPCSIV44 Unnown nput extended Kalman flter-based fault dagnoss for satellte actuator

More information

Physics 2A Chapters 6 - Work & Energy Fall 2017

Physics 2A Chapters 6 - Work & Energy Fall 2017 Physcs A Chapters 6 - Work & Energy Fall 017 These notes are eght pages. A quck summary: The work-energy theorem s a combnaton o Chap and Chap 4 equatons. Work s dened as the product o the orce actng on

More information

A NOVEL DESIGN APPROACH FOR MULTIVARIABLE QUANTITATIVE FEEDBACK DESIGN WITH TRACKING ERROR SPECIFICATIONS

A NOVEL DESIGN APPROACH FOR MULTIVARIABLE QUANTITATIVE FEEDBACK DESIGN WITH TRACKING ERROR SPECIFICATIONS A OVEL DESIG APPROACH FOR MULTIVARIABLE QUATITATIVE FEEDBACK DESIG WITH TRACKIG ERROR SPECIFICATIOS Seyyed Mohammad Mahd Alav, Al Khak-Sedgh, Batool Labb Department of Electronc and Computer Engneerng,

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k) ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

Army Ants Tunneling for Classical Simulations

Army Ants Tunneling for Classical Simulations Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Supplemental Material: Causal Entropic Forces

Supplemental Material: Causal Entropic Forces Supplemental Materal: Causal Entropc Forces A. D. Wssner-Gross 1, 2, and C. E. Freer 3 1 Insttute for Appled Computatonal Scence, Harvard Unversty, Cambrdge, Massachusetts 02138, USA 2 The Meda Laboratory,

More information

FE REVIEW OPERATIONAL AMPLIFIERS (OP-AMPS)

FE REVIEW OPERATIONAL AMPLIFIERS (OP-AMPS) FE EIEW OPEATIONAL AMPLIFIES (OPAMPS) 1 The Opamp An opamp has two nputs and one output. Note the opamp below. The termnal labeled wth the () sgn s the nvertng nput and the nput labeled wth the () sgn

More information

TREATMENT OF THE TURNING POINT IN ADK-THEORY INCLUDING NON-ZERO INITIAL MOMENTA

TREATMENT OF THE TURNING POINT IN ADK-THEORY INCLUDING NON-ZERO INITIAL MOMENTA 41 Kragujevac J. Sc. 5 (00) 41-46. TREATMENT OF THE TURNING POINT IN ADK-THEORY INCLUDING NON-ZERO INITIAL MOMENTA Vladmr M. Rstć a and Tjana Premovć b a Faculty o Scence, Department o Physcs, Kragujevac

More information

A Tale of Friction Basic Rollercoaster Physics. Fahrenheit Rollercoaster, Hershey, PA max height = 121 ft max speed = 58 mph

A Tale of Friction Basic Rollercoaster Physics. Fahrenheit Rollercoaster, Hershey, PA max height = 121 ft max speed = 58 mph A Tale o Frcton Basc Rollercoaster Physcs Fahrenhet Rollercoaster, Hershey, PA max heght = 11 t max speed = 58 mph PLAY PLAY PLAY PLAY Rotatonal Movement Knematcs Smlar to how lnear velocty s dened, angular

More information

6.3.7 Example with Runga Kutta 4 th order method

6.3.7 Example with Runga Kutta 4 th order method 6.3.7 Example wth Runga Kutta 4 th order method Agan, as an example, 3 machne, 9 bus system shown n Fg. 6.4 s agan consdered. Intally, the dampng of the generators are neglected (.e. d = 0 for = 1, 2,

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

Neuro-Adaptive Design - I:

Neuro-Adaptive Design - I: Lecture 36 Neuro-Adaptve Desgn - I: A Robustfyng ool for Dynamc Inverson Desgn Dr. Radhakant Padh Asst. Professor Dept. of Aerospace Engneerng Indan Insttute of Scence - Bangalore Motvaton Perfect system

More information

2 STATISTICALLY OPTIMAL TRAINING DATA 2.1 A CRITERION OF OPTIMALITY We revew the crteron of statstcally optmal tranng data (Fukumzu et al., 1994). We

2 STATISTICALLY OPTIMAL TRAINING DATA 2.1 A CRITERION OF OPTIMALITY We revew the crteron of statstcally optmal tranng data (Fukumzu et al., 1994). We Advances n Neural Informaton Processng Systems 8 Actve Learnng n Multlayer Perceptrons Kenj Fukumzu Informaton and Communcaton R&D Center, Rcoh Co., Ltd. 3-2-3, Shn-yokohama, Yokohama, 222 Japan E-mal:

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

Output Feedback Adaptive Robust Control of Uncertain Linear Systems With Disturbances

Output Feedback Adaptive Robust Control of Uncertain Linear Systems With Disturbances Bn Yao 1 School of Mechancal Engneerng, Purdue Unversty, West Lafayette, IN 47907 e-mal: byao@purdue.edu L Xu Research and Advanced Engneerng, Ford Motor Company, Dearborn, MI 48124 Output Feedback Adaptve

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

SOC Estimation of Lithium-ion Battery Packs Based on Thevenin Model Yuanqi Fang 1,a, Ximing Cheng 1,b, and Yilin Yin 1,c. Corresponding author

SOC Estimation of Lithium-ion Battery Packs Based on Thevenin Model Yuanqi Fang 1,a, Ximing Cheng 1,b, and Yilin Yin 1,c. Corresponding author Appled Mechancs and Materals Onlne: 2013-02-13 ISSN: 1662-7482, Vol. 299, pp 211-215 do:10.4028/www.scentfc.net/amm.299.211 2013 Trans Tech Publcatons, Swtzerland SOC Estmaton of Lthum-on Battery Pacs

More information

PHYS 1441 Section 002 Lecture #15

PHYS 1441 Section 002 Lecture #15 PHYS 1441 Secton 00 Lecture #15 Monday, March 18, 013 Work wth rcton Potental Energy Gravtatonal Potental Energy Elastc Potental Energy Mechancal Energy Conservaton Announcements Mdterm comprehensve exam

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Why working at higher frequencies?

Why working at higher frequencies? Advanced course on ELECTRICAL CHARACTERISATION OF NANOSCALE SAMPLES & BIOCHEMICAL INTERFACES: methods and electronc nstrumentaton. MEASURING SMALL CURRENTS When speed comes nto play Why workng at hgher

More information

THE Kalman filter (KF) rooted in the state-space formulation

THE Kalman filter (KF) rooted in the state-space formulation Proceedngs of Internatonal Jont Conference on Neural Networks, San Jose, Calforna, USA, July 31 August 5, 211 Extended Kalman Flter Usng a Kernel Recursve Least Squares Observer Pngpng Zhu, Badong Chen,

More information

Lecture 5: Operational Amplifiers and Op Amp Circuits

Lecture 5: Operational Amplifiers and Op Amp Circuits Lecture 5: peratonal mplers and p mp Crcuts Gu-Yeon We Dson o Engneerng and ppled Scences Harard Unersty guyeon@eecs.harard.edu We erew eadng S&S: Chapter Supplemental eadng Background rmed wth our crcut

More information

Testing for seasonal unit roots in heterogeneous panels

Testing for seasonal unit roots in heterogeneous panels Testng for seasonal unt roots n heterogeneous panels Jesus Otero * Facultad de Economía Unversdad del Rosaro, Colomba Jeremy Smth Department of Economcs Unversty of arwck Monca Gulett Aston Busness School

More information

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω

More information

Translational Equations of Motion for A Body Translational equations of motion (centroidal) for a body are m r = f.

Translational Equations of Motion for A Body Translational equations of motion (centroidal) for a body are m r = f. Lesson 12: Equatons o Moton Newton s Laws Frst Law: A artcle remans at rest or contnues to move n a straght lne wth constant seed there s no orce actng on t Second Law: The acceleraton o a artcle s roortonal

More information

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient 58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),

More information

Iterative General Dynamic Model for Serial-Link Manipulators

Iterative General Dynamic Model for Serial-Link Manipulators EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for Seral-Lnk Manpulators In ths set of notes, we are gong to develop a method for computng a general

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

P R. Lecture 4. Theory and Applications of Pattern Recognition. Dept. of Electrical and Computer Engineering /

P R. Lecture 4. Theory and Applications of Pattern Recognition. Dept. of Electrical and Computer Engineering / Theory and Applcatons of Pattern Recognton 003, Rob Polkar, Rowan Unversty, Glassboro, NJ Lecture 4 Bayes Classfcaton Rule Dept. of Electrcal and Computer Engneerng 0909.40.0 / 0909.504.04 Theory & Applcatons

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

Lesson 5: Kinematics and Dynamics of Particles

Lesson 5: Kinematics and Dynamics of Particles Lesson 5: Knematcs and Dynamcs of Partcles hs set of notes descrbes the basc methodology for formulatng the knematc and knetc equatons for multbody dynamcs. In order to concentrate on the methodology and

More information

Lesson 16: Basic Control Modes

Lesson 16: Basic Control Modes 0/8/05 Lesson 6: Basc Control Modes ET 438a Automatc Control Systems Technology lesson6et438a.tx Learnng Objectves Ater ths resentaton you wll be able to: Descrbe the common control modes used n analog

More information

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 4. Moton Knematcs 4.2 Angular Velocty Knematcs Summary From the last lecture we concluded that: If the jonts

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Robust Fuzzy Control of Electrical Manipulators

Robust Fuzzy Control of Electrical Manipulators J Intell Robot Syst (21) 6:415 434 DOI 1.17/s1846-1-943-y Robust Fuzzy Control of Electrcal Manpulators Mohammad Mehd Fateh Receved: 1 September 29 / Accepted: 19 Aprl 21 / Publshed onlne: 22 June 21 Sprnger

More information

Designing of Combined Continuous Lot By Lot Acceptance Sampling Plan

Designing of Combined Continuous Lot By Lot Acceptance Sampling Plan Internatonal Journal o Scentc Research Engneerng & Technology (IJSRET), ISSN 78 02 709 Desgnng o Combned Contnuous Lot By Lot Acceptance Samplng Plan S. Subhalakshm 1 Dr. S. Muthulakshm 2 1 Research Scholar,

More information

PHYS 1443 Section 004 Lecture #12 Thursday, Oct. 2, 2014

PHYS 1443 Section 004 Lecture #12 Thursday, Oct. 2, 2014 PHYS 1443 Secton 004 Lecture #1 Thursday, Oct., 014 Work-Knetc Energy Theorem Work under rcton Potental Energy and the Conservatve Force Gravtatonal Potental Energy Elastc Potental Energy Conservaton o

More information

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud Resource Allocaton wth a Budget Constrant for Computng Independent Tasks n the Cloud Wemng Sh and Bo Hong School of Electrcal and Computer Engneerng Georga Insttute of Technology, USA 2nd IEEE Internatonal

More information

A revised adaptive fuzzy sliding mode controller for robotic manipulators

A revised adaptive fuzzy sliding mode controller for robotic manipulators A revsed adaptve fuzzy sldng mode controller for robotc manpulators Xaosong Lu* Department of Systems and Computer Engneerng, Carleton Unversty, 5 Colonel By Drve, Ottawa, Ontaro, Canada E-mal: luxaos@sce.carleton.ca

More information

Numerical Solution of Ordinary Differential Equations

Numerical Solution of Ordinary Differential Equations Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples

More information

Chapter 07: Kinetic Energy and Work

Chapter 07: Kinetic Energy and Work Chapter 07: Knetc Energy and Work Conservaton o Energy s one o Nature s undamental laws that s not volated. Energy can take on derent orms n a gven system. Ths chapter we wll dscuss work and knetc energy.

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information