Application of Linear Model Predictive Control and Input-Output Linearization to Constrained Control of 3D Cable Robots

Size: px
Start display at page:

Download "Application of Linear Model Predictive Control and Input-Output Linearization to Constrained Control of 3D Cable Robots"

Transcription

1 Moern Mechancal Engneerng, 2011, 1, o: /mme Publshe Onlne November 2011 ( Applcaton of Lnear Moel Prectve Control an Input-Output Lnearzaton to Constrane Control of 3D Cable Robots Abstract Al Ghasem Member of Young Researchers Club, Scence an Research Branch, Islamc Aza Unversty, ehran, Iran E-mal: Receve September 3, 2011; revse October 13, 2011; accepte October 25, 2011 Cable robots are structurally the same as parallel robots but wth the basc fference that cables can only pull the platform an cannot push t. hs feature makes control of cable robots a lot more challengng compare to parallel robots. hs paper ntrouces a controller for cable robots uner force constrant. he controller s base on nput-output lnearzaton an lnear moel prectve control. Performance of nput-output lnearzng (IOL) controllers suffers ue to constrants on nput an output varables. hs problem s successfully tackle by augmentng IOL controllers wth lnear moel prectve controller (LMPC). he effecttveness of the propose metho s llustrate by numercal smulaton. Keywors: Cable Robots, Input-Output Lnearzaton, Lnear Moel Prectve Control 1. Introucton After a moton smulator wth parallel knematc chans was ntrouce n 1965 by D. Stewart [1], parallel manpulators receve more an more attenton because of ther hgh stffness, hgh spee, hgh accuracy, compact an hgh carryng capablty [2]. hey have been use wely n the fels of moton smulators, force/torque sensors, complance evces, mecal evces an machne tools [3,4]. A parallel robot s mae up of an en-effector, wth n egrees of freeom, an a fxe base lnke together by at least two nepenent knematc chans [5]. Actuaton takes place through m smple actuators. Parallel robots rawbacks are ther relatvely small workspace an knematcs complexty. Cable robots are a class of parallel robots n whch the lnks are replace by cables. hey are relatvely smple n form, wth multple cables attache to a moble platform or an en-effector. Cable robots posses a number of esrable characterstcs, nclung: 1) statonary heavy components an few movng parts, resultng n low nertal propertes an hgh payloa-to-weght ratos; 2) ncomparable moton range, much hgher than that obtane by conventonal seral or parallel robots; 3) cables have neglgble nerta an are sutable for hgh acceleraton applcatons; 4) transportablty an ease of sassembly/- reassembly; 5) reconfgurablty by smply relocatng the motors an upatng the control system accorngly; an, 6) economcal constructon an mantenance ue to few movng parts an relatvely smple components [6,7]. Consequently, cable robots are exceptonally well sute for many applcatons such as hanlng of heavy materals, nspecton an repar n shpyars an arplane hangars, hgh-spee manpulaton, raply eployable rescue robots, cleanup of saster areas, an access to remote locatons an nteracton wth hazarous envronments [6-12]. For these applcatons conventonal seral or parallel robots are mpractcal ue to ther lmte workspace. However, cables have the unque property they cannot prove compresson force on an en-effector. Some research has been prevously conucte to guarantee postve tenson n the cables whle the en-effector s movng. he ea of reunancy was utlze n cable system control [13,14]. hs paper ntrouces a controller for cable robots uner force constrant. By conserng, lnear moel prectve covers fferent constrant such as nput constrants. he goal s to apply the lnear moel prectve control to the nput-output lnearze system to account for the constrants. A varety of nonlnear control esgn strategy has been propose n the past two ecaes. Input-output lnearzaton (IOL) an nonlnear moel base control are the most wely stue esgn technques n nonlnear control. he

2 70 A. GHASEMI central ea of the nput output lnearzaton approach s to algebracally transform the nonlnear system nto lnear one an apply a sutable lnear control esgn technque [15,16]. LMPC s prmarly evelope for process control. herefore ts applcaton n robot control has less been reporte. he ncpent nterest n the applcatons of MPC ates back to the late 1970s. In 1978, Rchalet et al. [17], presente the Moel Prectve Heurstc Control (MPHC) metho n whch an mpulse response moel was use to prect the effect at the output of the future control actons. Lnear moel prectve control refers to a class of control algorthms that compute a manpulate varable profle by utlzng a lnear process moel to optmze a lnear or quaratc open-loop performance objectve subject to lnear constrants over a future tme horzon. he frst move of ths open loop optmal manpulate varable profle s then mplemente. hs proceure s repeate at each control nterval wth the process measurements use to upate the optmzaton problem. Durng 1980s, MPC quckly became popular partcularly n chemcal process nustry ue to the smplcty of the algorthm an to the use of the mpulse or step response moel, whch s preferre, as beng more ntutve an requrng no prevous nformaton for ts entfcaton [18]. A cable-suspene robot s actuate by servo motors that control the tensons n the cables. A major savantage of cable robots s that each cable can only exert tenson. hs constrant leas to performance eteroraton an even nstablty, f not properly accounte for n the control esgn proceure. Due to ths feature, well known results n robotcs for trajectory plannng an control are not rectly applcable to them. Several approaches nclung a lyapunov base controller wth varable gans an a feeback lnearzng controller wth varable gans [19], feeback lnearzaton together wth metho of reference sgnal management [20], lyapunov base slng controller wth metho of sgnal management [21] have been suggeste to satsfy the postve tenson n the cables whle the platform s movng. In ths paper a lnear moel prectve control s apple to lnearze moel. Moel prectve control, a computer control algorthm that utlzes an explct moel to prect the future response of a system s an effectve tool for hanlng constrane control problems. cable connecton n the platform frame. herefore, = a R b c s the vector representng the length of each cable an l s the recton of tenson force along each cable, where c s the poston vector of mass center of platform parameterze. R s the rotaton matrx between the two frames, the base an the movng, Fgure System Dynamcs he nerta of each lnk of cable robots s neglgble compare to that of the platform because the so-calle lnk s just a cable or wre. herefore, the ynamcs of the lnks can be gnore whch wll sgnfcantly smplfy the ynamc moel of the manpulator. One can erve the Newton-Euler equatons of moton of the manpulator wth respect to the center of mass on the platform as follows [10] F mamc f mg ( n s number of cables) M I I n ext 1 ext n 1 Rb f l I I where m an I are the mass an nerta tensor of the platform nclung any attache payloa; g s the gravty acceleraton vector; an f ext an τ ext are external force an moment vectors apple to the platform. c an α are the lnear an angular acceleraton vectors of the platform; an f are the force vector an force value of the th cable. Equaton (1) can be rewrtten nto a compact form as: M x Cx D Ju (2) where mi M,, f D C 033 I 033 I ext mg ext (1) 2. Knematcs Moelng of the Cable Robots he knematc notaton of a spatal cable-rven manpulator s presente n Fgure 1, where P an B are two attachng ponts of the th cable to the platform an the base, respectvely. a represents the poston vector of B n the base frame an b shows the poston vector of the Fgure 1. General knematcs of a cable robot.

3 A. GHASEMI 71 l1 ln J, J jacopan matrx Rb1l1 Rbnln f1 u f n I 3*3 s a 3*3 entty matrx an 0 I I z 0 z I I 0 y I I I x y x Equaton (2) can be wrtten nto a steay state form as: where y I x X F Gu y E X h X x x 06* n X, F, G 1 1 x M Cx D M J wth constrants 0 u u max 4. Input-Output Lnearzaton he ffculty of the trackng control esgn can be reuce f we can fn a rect an smple relaton between the system output y an the control nput u. Inee, ths ea consttutes the ntutve bass for the so-calle nput -output lnearzaton approach to nonlnear control esgn. hs artcle s ame to use the LMPC whch covers fferent constrants such as nput constrants. Because of LMPC s usually use for lnear screte systems. In the begnnng we wll lnearze the ynamcs equatons base on n Input-Output Lnearzaton. o generate a rect relatonshp between the output y an the nput u, Dfferentate the output y wth respect to tme t n Equaton (3), we have y L h X L h u t F gj j j1 where L F h (X) an L gj h (X) are the Le ervatves of h (X) wth respect to F(X) an g j (X), respectvely. If L gj h (X) = 0, y = L F h (X) then the r th ervatve of y can be represent n the followng form. n r r r1 F gj F j j1 In ths way, we can wrte all plant s nput-output equatons as n y L h X L L h u t (3) r1 r r 1 11 r11 y 1 LF h L 1 g L 1 F h1 Lg L h n F 1 u r q r q rq1 rq1 yq LF hq Lg L h 1 F qlg L h n F q where q an r are the number of egree of freeom of the robot an the relatve orer of the plant, respectvely. Equaton (4) can be represente n the followng compact form: v P Wu Now u can be obtane as: u W W W vp where 1 W s the pseuo nverse of W. 1 W W WW For a cable robotc system, t can be easly shown that the system has no zero ynamcs. herefore, the ecouplng matrx W s full rank an the control law s well X efne an a sutable change of coornates X y y y y yq y q where yel s a close- n the normal form [15]. loop system y xr 2 1,,q r1 rq n A Bv y H 5. Lnear Moel Prectve Control he basc structure of MPC to mplement s shown n Fgure 2 A moel s use to prect the future plant outputs, base on past an current values an on propose future control actons. hese actons are calculate by optmzer takng nto account the cost functon as well as constrants. he optmzer s another funamental part of the strategy as t proves the control acton. each component of ths structure s escrbe n more etal In the followng of ths artcle. Fgure 2. Basc structure of MPC. (4)

4 72 A. GHASEMI he goal s to apply the lnear moel prectve control to the nput-output lnearze system to account for the constrants. Snce the lnear moel prectve s more naturally formulate n screte tme, the lnear subsystem n (5) s scretze wth a samplng pero to yel k1 A kb vk (6) yk H k where A,B an H are obtane rectly from the contnuous-tme matrces [22]. Also, the state-epenent relaton between u(k) an v( k) s obtane as r vk L h r F X k LGLFhX k uk hs mappng can be rewrtten n the followng form v k PX k W X k u k 5.1. Constrant Mappng When lnear moel prectve control s apple to the system, t s necessary to map the constrants from the orgnal nput space to lnearze system. By conserng v s a new nput to be etermne, o obtan constrants on the new nput, he nput constrant mappng s performe usng nput output lnearzaton law an the current state measurement x(k). he transforme constrants can be etermne on each samplng pero by solvng the followng optmzaton problem: mn mn k vmn k j k P X k j max ukjk W X k j k u k j k 0 j N 1 uk ( jk ) u u k j j u max vmax k j k P X k j k W X k j k u k j k 0 j N 1 (7) where Z (k + j k) s the precte value of the system state (Z) at tme k + j base on the nformaton avalable at tme k. Note that the varable X (k + j k) cannot be calculate prove that the nput sequence s calculate, whch s not possble untl the constrants are specfe. herefore, at the begnnng (k = 0), the nput constrant over the entre control horzon can be presente by: v k j k v k k mn mn j1,, N1 max v k j k v k k max j1,, N1 hen, we use nputs calculate at last samplng tme to etermne future constrants at the current samplng tme. (8) herefore, Equaton (7) wll be change nto Equaton (9). mn max mn uk ( jk ) 1 max 1 uk ( jk ) max v k j k P X k j k W X k j k1 ] u k j k 0 jn1 v k j k P X k j k W X k j k1 u k j k 0 jn1 umn u k j j u (9) Now Equaton (9) can be solve to obtan v mn (k+ j k) an v max (k + j k). If W(,j) s postve, the control u(j) must be the smallest value for v mn an the largest value for v max an f W(,j) s negatve, then t must be the largest value for v mn an the smallest value for v max Lnear Moel Prectve Control Desgn he goal s to apply LMPC to the lnearze system to account for these constrants. Now the moel (6) s use n the nfnte horzon lnear moel prectve strategy propose by Muske an Rawlngs [23]. herefore, the openloop optmal control problem that the nput control foun by mnmzng the nfnte horzon crteron, can be expresse as V( kk ) j1 k j k H QH k j k mn ( ) ( ) vk j k vk j 1 k R vk j k vk j 1 k v k j k v S v k j k v (10) where ξ an v are target values for ξ an v, respectvely, an Q,S, an R, are postve sem efnte matrces. In orer to obtan value v(k + j k) t s necessary to mnmze the functonal of Equaton (10) to o ths value of the precte output are calculate as functon of pas values of nputs an outputs an future control sgnals obtan an expresson whose mnmzaton leas to the looke for values. he ecson vector s efne as V(k k) = [v(k k) v( k+1 k) v(k+n-1 k)], where N s the control horzon. All future moves beyon the control horzon are set equal to the target value v. As scusse n [23], the matrx A s unstable an n orer for the lnear moel prectve problem to have a feasble soluton t s necessary to mpose the equalty constrant ξ(k + N k) = ξ. o obtan a fnte set of ecson varables, nputs beyon the control horzon are set equal to the esre value: v(k + j k) = v, j N. herefore, the nfnte horzon lnear

5 A. GHASEMI 73 moel prectve problem Equaton (10) can be wrtten as a fnte horzon problem. V( kk ) vk N kv 1 mn 1 N 1 k j k H QH k j k j1 vk j k vk j 1 k 1 S v kn k v v k jk v S v k jk v R v k j k v k j k (11) hs optmzaton problem must be solve subjecte to the followng constrants: k N k (12) v k j k v k j k v k j k mn mn N 1 N 2 EE A B A B B N K N A HQH A 0 he soluton of Equaton (13) belongs to the regular system. o fn the solutons of the trackng one, that of the regular system shoul be shfte nto the orgn of the system to the steay state escrbe by ξ, an v. he esre values must le wthn the feasble regon efne by nput constrants for lnear moel prectve control an mnmze the control effort (Equaton (15)). Straghtforwar algebrac manpulaton of quaratc subject to : objectve functon of the corresponng regular system presente n Equaton (11) results n the followng stan- I A B 0 U ar quaratc program form: H 0 y mn V HHV 2V GG kffvk 1 (13) 0 DD U C Vkk where S B KN 1 A 0 GG, FF B K0 A 0 Equaton (13) shoul be solve subjecte to the followng constrants: DDV CC N (14) EE V A k where by conserng vmax k k I vmax k N 1 k 6* N DD, CC, I 6* N vmn k k vmn kn 1 k 2 mn ( 0 U 0m2m Imm U U U v m m Im m U U Qs 2 m*2m (15) U an Qs are esre value of the nput an a postve efnte matrx, respectvely. herefore control law s the summaton of the answers of Equaton (13) an Equaton (15) whch can be shown n the followng form. ( W X k k v k kp( X k k u k W X k k v k k P X k k (16) Fgure 3 shows schematc of the propose structure control. Fgure 3. Schematc of control structure. N 1 B K N1B R2S B AKN2 B S B A K0 B N 2 B KN 2 AB S B KN 2 B R 2S B A K0 B HH N1 N2 B K0 A B B K0 A B B K0 B R2S

6 74 A. GHASEMI 6. Smulaton In ths secton, smulaton results of applyng moel prectve control on a 3D cable robot wll be presente. able 1. shows the mensons of the cable robot. We conser a efnte movement of the platform from X 0 = [ 0.1, 0.1,1.5,0.002,.001, 0.001] to X = [0.2sn(t), 0.4 sn(t), t t 3,0,0,0], on a esre trajectory shown n otte lnes by Fgures 4(a) to 4(f), for poston/orentaton of the platform. Also, we conser the parameters of the moel precttve controller as: the control horzon N = 25, S, Q, an R of Equaqon (10), R = I, Q = I, S = 0.01I, an Q s of Equaton (15), Q s = I. Snce the controller nees platform s poston/orenttaton, at frst, we must solve forwar knematcs of the robot. hs has been carre out by the authors usng neural network algorthm [24]. Fgure 4 shows the moel prectve controller wth nput-out lnearzng worke qute well an a trajectory trackng are one. Fgure 5 shows the sx tensons n the cables vs. tme. As t can be seen, all of them reman postve urng the moton. able 1. Dmensons of the cable robot. Poston vector X(m) Y(m) Z(m) Poston vector x(m) y(m) z(m) a b a b a b a b a b a b 6 (a) (b) (c) () (e) Fgure 4. Plots of esre an actual poston an orentaton of the platform, (a)-(c) poston n X-Y-Z rectons, respectvely, ()-(f) orentaton aroun X-Y-Z recton, respectvely. (f)

7 A. GHASEMI 75 (a) (b) (c) () (e) (f) Fgure 5. Plots of tenson trajectores. 7. Conclusons In ths paper, a lnear moel prectve controller together wth an nput-output lnearzng control strategy for a constrane robotc system, a 3D cable robot, was evelope an evaluate. he control system s comprse of: 1) an nput-output lnearzng controller that accounts for cable robot nonlneartes; 2) a constrant mappng scheme that transforms the actual nput constrants nto nput constrants on the feeback lnearze system; an 3) a lnear moel prectve controller that proves explct compensaton for nput constrants. he smulaton re- of the propose metho. It sults showe the effectveness s worth nothng that ths approach can be extene for the reunant cable robots. 8. References [1] D. Stewart, A Platform wth Sx Degrees of Freeom, Proceengs of the Insttuton of Mechancal Engneer, Vol. 180, No. 15, 1965, pp o: /pime_proc_1965_180_029_02 [2] J. -P. Merlet, Parallel Robots, Sol Mechancs an Its Applcatons, Kluwer, Norwell, [3] J. -P. Merlet, Stll a Long Way to Go on the Roa for Parallel Mechansms, A Keynote Speech at Desgn Engneerng echncal Conferences, Montreal, 29 September-2 October [4] J. -P. Merlet, Parallel Robots, Open Problems, INRIA Sopha-Antpols, France. [5] L. -W. sa, Robot Analyss, he Mechancs of Seral an Parallel Manpulators, Wley, New York, [6] P. Bosscher, A.. Rechel an I. Ebert-Uphoff, Wrench- Feasble Workspace Generaton for Cablerven Robots, Journal of Intellgent an Robotc Systems, Vol. 22, No. 4, 2006, pp [7] A.. Rechel an I. Ebert-Uphoff, Force-Feasble Workspace Analyss for Unerconstrane Pontmass Cable Robots, Proceengs of IEEE Internatonal Conference on Robotcs an Automaton, New Orleans, 26 Aprl-1 May 2004, pp [8] S. Kawamura, W. Choe, S. anaka an S. R. Panan, Development of an Ultrahgh Spee FALCON Usng Wre Drve System, In Proceengs of the 1995 IEEE Internatonal Conference on Robotcs an Automaton, May 2003, pp [9] P. Lafourcae, M. Llbre an C. Reboulet, Desgn of a Par-

8 76 A. GHASEMI [10] X. Dao, O. Ma an R. Paz, Stuy of 6-DOF Cable Robots for Potental Applcaton of HIL Mcrogravty Con- tact-dynamcs Smulaton, In Proceengs of the AIAA Moelng an Smulaton echnologes Conference an Exhbt, Keystone, August 2006, pp [11] P. Gallna, G. Rosat an A. Ross, 3-D.O.F. Wre Drven Planar Haptc Interface, Journal of Intellgent an Robotc Systems,Vol. 32, No. 1, 2001, pp o: /a: [12] J. Albus, R. Bostelman an N. Dagalaks, he NIS Robocrane, Journal of Robotc Systems, Vol. 10, No. 5, 1993, pp o: /rob [13] Y. Q. Zheng, Workspace Analyss of a Sx DOF Wre- Drven Parallel Manpulator, Proceengs of the WORK- SHOP on Funamental Issues an Future Research Derecton for Parallel Mechansms an Manpulators, Quebec, 3-4 October 2002, pp [14] W. J. Shang, D. Cannon an J. Gorman, Dynamc Analyss of the Cable Array Robotc Crane, Proceengs of the IEEE Internatonal Conference on Robotcs an Automaton, Detrot, May 1999, pp [15] J. J. Slotne an W. Lepng, Apple Nonlnear Control, Prentce Hall, Englewoo Clffs, [16] H. Khall, Nonlnear Systems, 3r Eton, Prentce-Hall, Upper Sale Rver, [17] J. Rchalet, A. Raault, J. L. estu an J. Papon, Moel allel Wre-Drven Manpulator for Wn unnels, In Proceengs of the Workshop on Funa mental Issues an Future Research Drectons for Parallel Mechansms an Manpulators, Quebec, 3-4 October 2002, pp Prectve Heurstc Control: A pplcaton to Inustry Processes, Automatca, Vol. 14. No. 2, 1978, pp o: / (78) [18] C. E. Garca, D. M. Prett an Morar, Moel Prectve Control: heory an Practce a Survey, Automatca, Vol. 25, No , pp o: / (89) [19] A. B. Alp an A. K. Agrawal, Cable Suspene Robots: Desgn, Plannng an Control, Internatonal Conference on Robotcs Robotcs & Automaton, Washngton, DC, 9-13 May 2002, pp [20] S. R. Oh an A. K. Agrawal, Controller Desgn for a Non-reunant Cable Robot Uner Input Constrant, ASME Internatonal Mechancal Engneerng Congress & Exposton, November 2003, Washngton, DC. [21] S. R. Oh an A. K. Agrawal, A Control Lyapunov Approach for Feeback Control of Cable-Suspene Robots, IEEE Internatonal Conference on Robotcs an Automaton, Aprl 2007, pp [22] F. Frankln, J. Powell an L. Workman Dgtal Control of Dynamc Systems 2n Eton, Ason Wesley, Boston, 1994, pp [23] K. R. Muske an J. B. Rawlngs, Moel prectve control wth lnear moels, AIChE Journal, Vol. 39, No. 2, 1993, pp o: /ac [24] A. Ghasem, M. Eghtesa an M. Far Neural Network Soluton for Forwar Knematcs Problem of Cable Robot, Journal of Intellgent an Robotc Systems, Vol. 60, No. 2, 2010, pp o: /s z

ENGI9496 Lecture Notes Multiport Models in Mechanics

ENGI9496 Lecture Notes Multiport Models in Mechanics ENGI9496 Moellng an Smulaton of Dynamc Systems Mechancs an Mechansms ENGI9496 Lecture Notes Multport Moels n Mechancs (New text Secton 4..3; Secton 9.1 generalzes to 3D moton) Defntons Generalze coornates

More information

New Liu Estimators for the Poisson Regression Model: Method and Application

New Liu Estimators for the Poisson Regression Model: Method and Application New Lu Estmators for the Posson Regresson Moel: Metho an Applcaton By Krstofer Månsson B. M. Golam Kbra, Pär Sölaner an Ghaz Shukur,3 Department of Economcs, Fnance an Statstcs, Jönköpng Unversty Jönköpng,

More information

Advanced Mechanical Elements

Advanced Mechanical Elements May 3, 08 Advanced Mechancal Elements (Lecture 7) Knematc analyss and moton control of underactuated mechansms wth elastc elements - Moton control of underactuated mechansms constraned by elastc elements

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

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

High-Order Hamilton s Principle and the Hamilton s Principle of High-Order Lagrangian Function

High-Order Hamilton s Principle and the Hamilton s Principle of High-Order Lagrangian Function Commun. Theor. Phys. Bejng, Chna 49 008 pp. 97 30 c Chnese Physcal Socety Vol. 49, No., February 15, 008 Hgh-Orer Hamlton s Prncple an the Hamlton s Prncple of Hgh-Orer Lagrangan Functon ZHAO Hong-Xa an

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

Analytical classical dynamics

Analytical classical dynamics Analytcal classcal ynamcs by Youun Hu Insttute of plasma physcs, Chnese Acaemy of Scences Emal: yhu@pp.cas.cn Abstract These notes were ntally wrtten when I rea tzpatrck s book[] an were later revse to

More information

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,

More information

Dynamic Modeling of a Synchronous Generator Using T-S Fuzzy Approach

Dynamic Modeling of a Synchronous Generator Using T-S Fuzzy Approach e-issn : 0975-0 Hee-Jn Lee / Internatonal Journal of Engneerng an echnology (IJE) Dynamc oelng of a Synchronous Generator Usng -S Fuzzy Approach Hee-Jn Lee Department of Electronc Engneerng Kumoh Natonal

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

ENTROPIC QUESTIONING

ENTROPIC QUESTIONING ENTROPIC QUESTIONING NACHUM. Introucton Goal. Pck the queston that contrbutes most to fnng a sutable prouct. Iea. Use an nformaton-theoretc measure. Bascs. Entropy (a non-negatve real number) measures

More information

CHAPTER 4 HYDROTHERMAL COORDINATION OF UNITS CONSIDERING PROHIBITED OPERATING ZONES A HYBRID PSO(C)-SA-EP-TPSO APPROACH

CHAPTER 4 HYDROTHERMAL COORDINATION OF UNITS CONSIDERING PROHIBITED OPERATING ZONES A HYBRID PSO(C)-SA-EP-TPSO APPROACH 77 CHAPTER 4 HYDROTHERMAL COORDINATION OF UNITS CONSIDERING PROHIBITED OPERATING ZONES A HYBRID PSO(C)-SA-EP-TPSO APPROACH 4.1 INTRODUCTION HTC consttutes the complete formulaton of the hyrothermal electrc

More information

Spin-rotation coupling of the angularly accelerated rigid body

Spin-rotation coupling of the angularly accelerated rigid body Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s

More information

On the Error Modeling of a Novel Mobile Hybrid Parallel Robot

On the Error Modeling of a Novel Mobile Hybrid Parallel Robot On the Error Moelng of a Novel Moble Hbr Parallel Robot Yongbo Wang (,2), Huapeng Wu (), Hekk Hanroos (), Bngku Chen (2) () Department of Mechancal Engneerng IMVE, Lappeenranta Unverst of echnolog P.O.

More information

Chapter 7: Conservation of Energy

Chapter 7: Conservation of Energy Lecture 7: Conservaton o nergy Chapter 7: Conservaton o nergy Introucton I the quantty o a subject oes not change wth tme, t means that the quantty s conserve. The quantty o that subject remans constant

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Yukawa Potential and the Propagator Term

Yukawa Potential and the Propagator Term PHY304 Partcle Physcs 4 Dr C N Booth Yukawa Potental an the Propagator Term Conser the electrostatc potental about a charge pont partcle Ths s gven by φ = 0, e whch has the soluton φ = Ths escrbes the

More information

Underactuated Control and Distribution of Multi-Agent Systems in Stratified Flow Environments

Underactuated Control and Distribution of Multi-Agent Systems in Stratified Flow Environments 49th IEEE Conference on Decson an Control December 15-17, 1 Hlton Atlanta Hotel, Atlanta, GA, USA Uneractuate Control an Dstrbuton of Mult-Agent Systems n Stratfe Flow Envronments Robert H Krohn an Thomas

More information

Active Vibration Control Based on a 3-DOF Dual Compliant Parallel Robot Using LQR Algorithm

Active Vibration Control Based on a 3-DOF Dual Compliant Parallel Robot Using LQR Algorithm The 009 IEEE/RSJ Internatonal Conference on Intellgent Robots an Systems October -5, 009 St. Lous, USA Actve Vbraton Control Base on a 3-DOF Dual Complant Parallel Robot Usng LQR Algorthm Yuan Yun an Yangmn

More information

WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION?

WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION? ISAHP 001, Berne, Swtzerlan, August -4, 001 WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION? Masaak SHINOHARA, Chkako MIYAKE an Kekch Ohsawa Department of Mathematcal Informaton Engneerng College

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

REAL TIME OPTIMIZATION OF a FCC REACTOR USING LSM DYNAMIC IDENTIFIED MODELS IN LLT PREDICTIVE CONTROL ALGORITHM

REAL TIME OPTIMIZATION OF a FCC REACTOR USING LSM DYNAMIC IDENTIFIED MODELS IN LLT PREDICTIVE CONTROL ALGORITHM REAL TIME OTIMIZATION OF a FCC REACTOR USING LSM DYNAMIC IDENTIFIED MODELS IN LLT REDICTIVE CONTROL ALGORITHM Durask, R. G.; Fernandes,. R. B.; Trerweler, J. O. Secch; A. R. federal unversty of Ro Grande

More information

Large-Scale Data-Dependent Kernel Approximation Appendix

Large-Scale Data-Dependent Kernel Approximation Appendix Large-Scale Data-Depenent Kernel Approxmaton Appenx Ths appenx presents the atonal etal an proofs assocate wth the man paper [1]. 1 Introucton Let k : R p R p R be a postve efnte translaton nvarant functon

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

The Study of Teaching-learning-based Optimization Algorithm

The Study of Teaching-learning-based Optimization Algorithm Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute

More information

SIMPLIFIED MODEL-BASED OPTIMAL CONTROL OF VAV AIR- CONDITIONING SYSTEM

SIMPLIFIED MODEL-BASED OPTIMAL CONTROL OF VAV AIR- CONDITIONING SYSTEM Nnth Internatonal IBPSA Conference Montréal, Canaa August 5-8, 2005 SIMPLIFIED MODEL-BASED OPTIMAL CONTROL OF VAV AIR- CONDITIONING SYSTEM Nabl Nassf, Stanslaw Kajl, an Robert Sabourn École e technologe

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

On Liu Estimators for the Logit Regression Model

On Liu Estimators for the Logit Regression Model CESIS Electronc Workng Paper Seres Paper No. 59 On Lu Estmators for the Logt Regresson Moel Krstofer Månsson B. M. Golam Kbra October 011 The Royal Insttute of technology Centre of Excellence for Scence

More information

Design and Analysis of Landing Gear Mechanic Structure for the Mine Rescue Carrier Robot

Design and Analysis of Landing Gear Mechanic Structure for the Mine Rescue Carrier Robot Sensors & Transducers 214 by IFSA Publshng, S. L. http://www.sensorsportal.com Desgn and Analyss of Landng Gear Mechanc Structure for the Mne Rescue Carrer Robot We Juan, Wu Ja-Long X an Unversty of Scence

More information

Kinematics of Fluid Motion

Kinematics of Fluid Motion Knematcs of Flu Moton R. Shankar Subramanan Department of Chemcal an Bomolecular Engneerng Clarkson Unversty Knematcs s the stuy of moton wthout ealng wth the forces that affect moton. The scusson here

More information

PHZ 6607 Lecture Notes

PHZ 6607 Lecture Notes NOTE PHZ 6607 Lecture Notes 1. Lecture 2 1.1. Defntons Books: ( Tensor Analyss on Manfols ( The mathematcal theory of black holes ( Carroll (v Schutz Vector: ( In an N-Dmensonal space, a vector s efne

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

The Noether theorem. Elisabet Edvardsson. Analytical mechanics - FYGB08 January, 2016

The Noether theorem. Elisabet Edvardsson. Analytical mechanics - FYGB08 January, 2016 The Noether theorem Elsabet Evarsson Analytcal mechancs - FYGB08 January, 2016 1 1 Introucton The Noether theorem concerns the connecton between a certan kn of symmetres an conservaton laws n physcs. It

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

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

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15 NGN 40 ynamcs and Vbratons Homework # 7 ue: Frday, Aprl 15 1. Consder a concal hostng drum used n the mnng ndustry to host a mass up/down. A cable of dameter d has the mass connected at one end and s wound/unwound

More information

4DVAR, according to the name, is a four-dimensional variational method.

4DVAR, according to the name, is a four-dimensional variational method. 4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The

More information

On a one-parameter family of Riordan arrays and the weight distribution of MDS codes

On a one-parameter family of Riordan arrays and the weight distribution of MDS codes On a one-parameter famly of Roran arrays an the weght strbuton of MDS coes Paul Barry School of Scence Waterfor Insttute of Technology Irelan pbarry@wte Patrck Ftzpatrck Department of Mathematcs Unversty

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

Robust Dynamic Programming for Discounted Infinite-Horizon Markov Decision Processes with Uncertain Stationary Transition Matrice

Robust Dynamic Programming for Discounted Infinite-Horizon Markov Decision Processes with Uncertain Stationary Transition Matrice roceengs of the 2007 IEEE Symposum on Approxmate Dynamc rogrammng an Renforcement Learnng (ADRL 2007) Robust Dynamc rogrammng for Dscounte Infnte-Horzon Markov Decson rocesses wth Uncertan Statonary Transton

More information

On the Coordinated Control of Multiple HVDC Links: Modal Analysis Approach

On the Coordinated Control of Multiple HVDC Links: Modal Analysis Approach R. Erksson an V. Knazkns / GMSARN nternatonal Journal 2 (2008) 15-20 On the Coornate Control of Multple HVDC Lnks: Moal Analyss Approach Robert Erksson an Valerjs Knazkns Abstract There are several possbltes

More information

Field and Wave Electromagnetic. Chapter.4

Field and Wave Electromagnetic. Chapter.4 Fel an Wave Electromagnetc Chapter.4 Soluton of electrostatc Problems Posson s s an Laplace s Equatons D = ρ E = E = V D = ε E : Two funamental equatons for electrostatc problem Where, V s scalar electrc

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 3: Element Properties Lecture 1: Natural Coordinates Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force. Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act

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

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

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

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

Application of Gravitational Search Algorithm for Optimal Reactive Power Dispatch Problem

Application of Gravitational Search Algorithm for Optimal Reactive Power Dispatch Problem Applcaton of Gravtatonal Search Algorthm for Optmal Reactve Power Dspatch Problem Serhat Duman Department of Electrcal Eucaton, Techncal Eucaton Faculty, Duzce Unversty, Duzce, 8620 TURKEY serhatuman@uzce.eu.tr

More information

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Solutions to exam in SF1811 Optimization, Jan 14, 2015 Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable

More information

MEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut

MEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut MEV442 Introducton to Robotcs Module 2 Dr. Santhakumar Mohan Assstant Professor Mechancal Engneerng Natonal Insttute of Technology Calcut Jacobans: Veloctes and statc forces Introducton Notaton for tme-varyng

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 Multi-Axis Force Measurement System for a Space Docking Mechanism

A Multi-Axis Force Measurement System for a Space Docking Mechanism 3rd Internatonal Conference on Materal, Mechancal and Manufacturng Engneerng (IC3ME 215) A Mult-Axs orce Measurement System for a Space Dockng Mechansm Gangfeng Lu a*, Changle L b and Zenghu Xe c Buldng

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

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

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

1 Convex Optimization

1 Convex Optimization Convex Optmzaton We wll consder convex optmzaton problems. Namely, mnmzaton problems where the objectve s convex (we assume no constrants for now). Such problems often arse n machne learnng. For example,

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 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

Time dependent weight functions for the Trajectory Piecewise-Linear approach?

Time dependent weight functions for the Trajectory Piecewise-Linear approach? Tme epenent weght functons for the Trajectory Pecewse-Lnear approach? Juan Pablo Amorocho an Heke Faßbener Abstract Moel orer reucton (MOR) has become an ubqutous technque n the smulaton of large-scale

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

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

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

Passive Bilateral Teleoperation with Constant Time Delays

Passive Bilateral Teleoperation with Constant Time Delays Proceengs of the 6 IEEE Internatonal Conference on Robotcs an Automaton Orlano, Flora - May 6 Passve Blateral eleoperaton wth Constant me Delays Dongjun Lee an Mark W Spong Coornate Scence Laboratory Unversty

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

Robust Model Predictive Control

Robust Model Predictive Control Robust Model Predctve Control Formulatons of robust control [1] he robust control problem concerns to the control of plants that are only approxmately nown. Usually, t s assumed that the plant les n a

More information

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two World Academy of Scence Engneerng and echnology Internatonal Journal of Mathematcal and omputatonal Scences Vol: No:0 00 ubc rgonometrc B-Splne Appled to Lnear wo-pont Boundary Value Problems of Order

More information

Optimization Techniques for Natural Resources

Optimization Techniques for Natural Resources Optmzaton Technques for Natural Resources SEFS 540 / ESRM 490 B Lecture 1 (3/27/2017) About the Instructor Teachng: nspre stuents to be curous but crtcal learners who can thnk for themselves an nurture

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

( ) = : a torque vector composed of shoulder torque and elbow torque, corresponding to

( ) = : a torque vector composed of shoulder torque and elbow torque, corresponding to Supplementary Materal for Hwan EJ, Donchn O, Smth MA, Shamehr R (3 A Gan-Fel Encon of Lmb Poston an Velocty n the Internal Moel of Arm Dynamcs. PLOS Boloy, :9-. Learnn of ynamcs usn bass elements he nternal

More information

Lossy Compression. Compromise accuracy of reconstruction for increased compression.

Lossy Compression. Compromise accuracy of reconstruction for increased compression. Lossy Compresson Compromse accuracy of reconstructon for ncreased compresson. The reconstructon s usually vsbly ndstngushable from the orgnal mage. Typcally, one can get up to 0:1 compresson wth almost

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

Non-negative Matrices and Distributed Control

Non-negative Matrices and Distributed Control Non-negatve Matrces an Dstrbute Control Yln Mo July 2, 2015 We moel a network compose of m agents as a graph G = {V, E}. V = {1, 2,..., m} s the set of vertces representng the agents. E V V s the set of

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

The Chaotic Robot Prediction by Neuro Fuzzy Algorithm (2) = θ (3) = ω. Asin. A v. Mana Tarjoman, Shaghayegh Zarei

The Chaotic Robot Prediction by Neuro Fuzzy Algorithm (2) = θ (3) = ω. Asin. A v. Mana Tarjoman, Shaghayegh Zarei The Chaotc Robot Predcton by Neuro Fuzzy Algorthm Mana Tarjoman, Shaghayegh Zare Abstract In ths paper an applcaton of the adaptve neurofuzzy nference system has been ntroduced to predct the behavor of

More information

Phys 331: Ch 7,.2 Unconstrained Lagrange s Equations 1

Phys 331: Ch 7,.2 Unconstrained Lagrange s Equations 1 Phys 33: Ch 7 Unconstrane agrange s Equatons Fr0/9 Mon / We /3 hurs /4 7-3 agrange s wth Constrane 74-5 Proof an Eaples 76-8 Generalze Varables & Classcal Haltonan (ecoen 79 f you ve ha Phys 33) HW7 ast

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

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

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

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

Appendix B: Resampling Algorithms

Appendix B: Resampling Algorithms 407 Appendx B: Resamplng Algorthms A common problem of all partcle flters s the degeneracy of weghts, whch conssts of the unbounded ncrease of the varance of the mportance weghts ω [ ] of the partcles

More information

Physics 111: Mechanics Lecture 11

Physics 111: Mechanics Lecture 11 Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton

More information

Modeling and Simulation of a Hexapod Machine Tool for the Dynamic Stability Analysis of Milling Processes. C. Henninger, P.

Modeling and Simulation of a Hexapod Machine Tool for the Dynamic Stability Analysis of Milling Processes. C. Henninger, P. Smpack User Meetng 27 Modelng and Smulaton of a Heapod Machne Tool for the Dynamc Stablty Analyss of Mllng Processes C. Hennnger, P. Eberhard Insttute of Engneerng project funded by the DFG wthn the framework

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

PLANAR TRANSLATIONAL CABLE-DIRECT-DRIVEN ROBOTS

PLANAR TRANSLATIONAL CABLE-DIRECT-DRIVEN ROBOTS PLANA TANSLATIONAL CABLE-DIECT-DIVEN OBOTS obert L. Wllams II Oho Unversty Athens, Oho Paolo Gallna 2 Unversty of Treste Treste, Italy Jgar Vada 3 Oho Unversty Athens, Oho Journal of obotc Systems Vol.

More information

High resolution entropy stable scheme for shallow water equations

High resolution entropy stable scheme for shallow water equations Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal

More information

Amiri s Supply Chain Model. System Engineering b Department of Mathematics and Statistics c Odette School of Business

Amiri s Supply Chain Model. System Engineering b Department of Mathematics and Statistics c Odette School of Business Amr s Supply Chan Model by S. Ashtab a,, R.J. Caron b E. Selvarajah c a Department of Industral Manufacturng System Engneerng b Department of Mathematcs Statstcs c Odette School of Busness Unversty of

More information

Distributed Exponential Formation Control of Multiple Wheeled Mobile Robots

Distributed Exponential Formation Control of Multiple Wheeled Mobile Robots Proceedngs of the Internatonal Conference of Control, Dynamc Systems, and Robotcs Ottawa, Ontaro, Canada, May 15-16 214 Paper No. 46 Dstrbuted Exponental Formaton Control of Multple Wheeled Moble Robots

More information

Design of Optimum Controllers for Gas Turbine Engines

Design of Optimum Controllers for Gas Turbine Engines Desgn of Optmum Controllers for Gas Turbne Engnes Junxa Mu 1, Dav Rees 1, Cer Evans 1 an Neophytos Chras 1 School of Electroncs, Unversty of Glamorgan Pontypr, CF37 1DL, Wales, UK Phone: +44() 1443 4859

More information

Mathematical Preparations

Mathematical Preparations 1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity 1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum

More information

A Particle Swarm approach for the Design of Variable Structure Stabilizer for a Nonlinear Model of SMIB System

A Particle Swarm approach for the Design of Variable Structure Stabilizer for a Nonlinear Model of SMIB System A Partcle Swarm approach for the Desgn of Varable Structure Stablzer for a Nonlnear Moel of SMIB System NAI A AL-MUSABI**, ZAKARIYA M AL-HAMOUZ*, HUSSAIN N AL-DUWAISH* ** The Petroleum Insttute, Electrcal

More information

Lecture 23: Newton-Euler Formulation. Vaibhav Srivastava

Lecture 23: Newton-Euler Formulation. Vaibhav Srivastava Lecture 23: Newton-Euler Formulaton Based on Chapter 7, Spong, Hutchnson, and Vdyasagar Vabhav Srvastava Department of Electrcal & Computer Engneerng Mchgan State Unversty Aprl 10, 2017 ECE 818: Robotcs

More information

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010 Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton

More information

Report on Image warping

Report on Image warping Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.

More information