Robust Fuzzy Control of Electrical Manipulators

Size: px
Start display at page:

Download "Robust Fuzzy Control of Electrical Manipulators"

Transcription

1 J Intell Robot Syst (21) 6: DOI 1.17/s y Robust Fuzzy Control of Electrcal Manpulators Mohammad Mehd Fateh Receved: 1 September 29 / Accepted: 19 Aprl 21 / Publshed onlne: 22 June 21 Sprnger Scence+Busness Meda B.V. 21 Abstract Stablty analyss for fuzzy control of robot manpulators has been a serous challengng problem n lterature. The theoretcal dffcultes are hghly ncreased due to the complexty of both manpulator dynamcs and fuzzy controller structure. Ths paper develops a novel robust fuzzy control approach for electrcal robot manpulators usng the drect method of Lyapunov. We pass analytcal dffcultes by the use of voltage control strategy n replace of torque control strategy. Then, a normalzed and decentralzed Takag Sugeno fuzzy controller s presented n a smple structure. A smple Lyapunov canddate s proposed to apply stablty analyss wthout knowng the explct dynamcs of system. Consequently, fuzzy control s analyzed and desgned as a robust nonlnear control. Roles of scalng factors, gans n consequent parts, and membershp functons n condton parts are consdered n the control desgn. The proposed control approach s appled on a Puma56 robot arm. Keywords Decentralzed control Electrcal robot manpulator Robust control Takag Sugeno fuzzy controller Stablty analyss 1 Introducton Industral robots exhbt hgh accuracy, good resoluton, and sutable repeatablty. Therefore, a smple control strategy such as ndependent jont control [1] usng an ordnary proportonal-dervatve-ntegral (PID) controller shows a good performance on ndustral robots. However, an ndustral robot s constructed wth a hgh expense to have hgh mechancal specfcatons. To save the cost, we may pay more attenton on mprovng the control of robot. In order to use a cheaper robot for performng M. M. Fateh (B) Department of Electrcal and Robotc Engneerng, Shahrood Unversty of Technology, Shahrood, Iran e-mal: mmfateh@shahroodut.ac.r

2 416 J Intell Robot Syst (21) 6: precse applcatons, the vsual feedback of end-effector n a vson control has been proposed to compensate the naccuracy nherent n a cheaply-constructed robot [2]. Snce robust control can guarantee the stablty and show a good performance [3], the robust control of cheaply-constructed robot wth uncertan knematcs and dynamcs wll provde a satsfactory performance. Many valuable robust control approaches have been developed to control robot manpulators n the jont-space [3, 4] and the task-space [5, 6]. However, robust control may nvolve the complexty of manpulator dynamcs. A great attenton has been attracted to remove ths shortage. A proper uncertanty bound parameter has been proposed to smplfy and mprove the robust control of robot manpulators [7]. The voltage control strategy [8] s superor to the torque control strategy n the robust control of electrcal manpulators [9] n terms of smplcty n desgn and performance of control. Alternatvely, fuzzy control as a model-free approach s smply desgned to control complcated systems [1]. To form fuzzy rules, an exact knowledge of model s not requred. Fuzzy controller s an ntellgent controller usng lngustc fuzzy rules to nclude nformaton from experts. Consequently, fuzzy control of robot manpulators has attracted a great deal of researches to overcome uncertanty, nonlnearty and couplng [11 14]. In many applcatons, fuzzy controllers were utlzed and treated as black-box controllers that when constructed properly by tral-and-error method, could produce satsfactory control results [15]. However, an analytcal proof should be gven to guarantee stablty and provde a desred performance. An analytcal well-defned form of fuzzy system wll mprove the analytcal aspect of fuzzy control n nonlnear control theory. The Takag Sugeno (TS) system performs well n an adaptve manner to cope wth nonlnear systems [16]. It performs lke a part-to-part lnear functon that changes smoothly from one part to another. Stablty analyss of fuzzy control to robot manpulators has been a serous challenge n lterature. The theoretcal dffcultes are hghly ncreased due to the complexty of manpulator dynamcs and fuzzy controller together [17]. To pass dffcultes n analyss and desgn, we use the voltage control strategy [8] n replace of the torque control strategy. As a result, the motor dynamcs whch s much smpler than manpulator dynamcs can be employed to decentralze and decouple the robotc system [9]. Ths paper presents stablty analyss wthout knowng the explct dynamcs of system, and smple structures are proposed to the fuzzy controller and the Lyapunov canddate. Then, a novel robust decentralzed TS fuzzy control approach s developed to electrcal robot manpulators based on the nonlnear control theory. Ths paper s organzed as follows. Secton 2 presents the decentralzed system. Secton 3 desgns the controller and presents the stablty analyss. Secton 4 llustrates smulaton results and Secton 5 concludes the paper. 2 Decentralzed System Decentralzed control s a domnant control scheme n the multvarable control because t has many advantages, such as flexblty n operaton, falure tolerance, smplfed desgn and tunng [18]. The decentralzed controller s stll adopted by majorty of modern robots n favor of ts computaton smplcty and low-cost

3 J Intell Robot Syst (21) 6: hardware setup [19]. As a result, mprovng the trackng performance of robot manpulators usng the decentralzed control s stll an attractve research. The decentralzed strategy employs a smple nput output form of fuzzy control for controllng a sngle output varable. Consequently, the number of fuzzy rules n the system s greatly reduced and thus, the computatonal complexty of algorthm s sgnfcantly smplfed. To desgn the decentralzed fuzzy controllers, the robotc system s decomposed. The torque vector s commonly used as the nput of robot manpulator for decentralzng the system [2, 21]. The torque vector of manpulator s gven by D (q) q + C (q, q) q + g (q) = T (1) where q R n s a vector of generalzed jont postons, D(q) R n n s a matrx of manpulator nerta, C (q, q) q R n s a vector of centrfugal and Corols generalzed forces, g(q) R n s a vector of gravtatonal generalzed forces and T R n s the nput generalzed forces. Dynamcs of manpulator as stated by (1) s hghly coupled such that moton of the k-th jont s not only dependent on the k-th nput [1]. Therefore, ths model s not properly suted for decentralzng robotc system. Moreover, the complexty of control desgn and analyss s hghly ncreased snce the robot dynamcs s hghly nonlnear, very large, and uncertan. Consder an electrcal robot drven by permanent magnet dc motors. The electrcal equaton of motor s n the form of u = R a I a + L a İ a + k b r 1 q k (2) where u R s the motor voltage, I a R s the motor current, and q k R s the velocty of k-th jont. Motor parameters, R a, L a k b and r are resstance, nductance, back-emf constant, and the reducton gear coeffcent, respectvely. Equaton 2 s more sutable than (1) for decentralzed control snce velocty of the k-th jont only depends on the voltage and current of the k-th motor that drves the jont. Moreover, as a man advantage, (2) s smple and free of manpulator dynamcs. 3 Control Desgn and Stablty Analyss The control approach s appled on a Puma56 robot arm, whch possesses three degrees of freedom drven by three revolute jonts. Each jont s drven by a dc motor usng the decentralzed fuzzy controller. The motor voltage s consdered as the output of fuzzy controller, and the jont poston error and ts dervatve are the two nputs of fuzzy controller. The consequent part of TS fuzzy rule s a lnear functon of nput varables [22]. The general TS fuzzy systems wth lnear rule consequent are unversal approxmators [23] and nonlnear controllers wth varable gans [24]. They can be nherently nonlnear gan schedulng controllers wth dfferent varable gans n dfferent regons of nput space [25]. Due to the generalty n the confguraton, the expressons of general fuzzy controllers cannot be mathematcally dervable. Thus, stablty analyss should be presented wthout knowng the explct structure of fuzzy controller. To make the dervaton mathematcally feasble, the components of fuzzy controllers must be specfc. All parameters n the rule consequent are adjustable and

4 418 J Intell Robot Syst (21) 6: unknown, and the number of parameters grows exponentally wth ncrease n the number of nputs. The parameters n theory provde a tunng of local control acton, possbly resultng n superor control performance. However, manual tunng of these parameters n practce could be neffectve, neffcent, napproprate, or sometmes even mpossble when the parameters are too many [15]. We desgn the decentralzed TS fuzzy controller by the use of two nputs namely trackng error and ts dervatve, respectvely. If we select three fuzzy sets for each nput, the whole control space wll be covered by nne fuzzy rules. The TS Fuzzy rules are n the form of: Ru l : If x 1 s A l 1 and x 2 s A l 2 Then yl = a l 1 x 1 + a l 2 x 2 + a l (3) where Ru l denotes the fuzzy rule for l = 1,..., 9, and the nputs of rule Ru l are x 1 and x 2. The nput vector s formed as x = [x 1 x 2 ] T U R 2 where U s the unverse of dscourse. In the l-th rule s denoted as Ru l, A l s the fuzzy set, a l s the gan n consequent part for = 1, 2, andy l s the crsp output. The output of TS fuzzy system f (x) V R by usng product nference engne, sngleton fuzzfer and center average defuzzfer s n the form of [1] f (x) = 9 2 μ A l (x ) y l l=1 =1 9 l=1 =1 2 μ A l (x ) (4) where μ A l (x ) [, 1] s the Membershp Functon (MF). An mportant contrbuton of fuzzy systems theory s to provde a systematc procedure for transformng a set of lngustc rules nto a nonlnear mappng. Consequently, the fuzzy controller defned by f (x) s a nonlnear functon. Substtutng y l from (3) nto(4) yelds a smple form for f (x) as where c j (x) for j =, 1, 2 s gven by f (x) = c 1 (x) x 1 + c 2 (x) x 2 + c (x) (5) c j (x) = 9 l=1 =1 9 2 μ A l (x ) a l j l=1 =1 2 μ A l (x ) (6) The obtaned analytcal structure of TS fuzzy controller mproves our study to develop the analyss and desgn. When x 1 s the error and x 2 s the frst tme dervatve of error, (6) presents a decentralzed PD lke TS fuzzy controller. However, we can generalze the proposed structure for other types of varable gan controllers f we select the relevant varable for x. As a result, we smply provde other types of varable gan controller, such as Proportonal (P), Proportonal Integral (PI), or Proportonal Integral Dervatve (PID) fuzzy controller. Moreover, a hgher order decentralzed controller can be obtaned f we use nputs x 1,..., x n. The varable gan c j (x) s a hgh mert of ths controller to make t adaptve n dfferent condtons. To

5 J Intell Robot Syst (21) 6: normalze the controller, we use the nput scalng factors k 1 > and k 2 >, andthe output scalng factor k o >. The nput vector s formed as where for the k-th jont z 1 and z 2 are defned as x = [ k 1 z 1 k 2 z 2 ] T (7) z 1 = q dk q k (8) z 2 = q dk q k (9) where q dk and q k are the desred and actual trajectores, respectvely. Substtutng (8)and(9) nto(7)forx 1 = k 1 z 1, x 2 = k 2 z 2,andż 1 = z 2 yelds where α = k 1 /k 2 >. Fuzzy controller by the use of scalng factors s formed as ẋ 1 = αx 2 (1) u (x) = k o (c 1 (x) x 1 + c 2 (x) x 2 + c (x)) (11) Ths general structure shows a nonlnear varable gan controller that fnds many applcatons n control. The nonlnear gan c j (x) covers the nonlnearty of controller by parameters n hand. The control purposes are smply descrbed by lngustc rules n fuzzy controller transformed to a nonlnear functon as stated by (11). c j (x) s comprsed by MFs μ A l (x ) n the rule condton parts and the gans a l j n the rule consequent parts as descrbed by (6). The MFs play a sgnfcant role n the performance of control system. The center, range and shape of MF are the most sgnfcant parameters. The operatng range of varables should be covered by the selected fuzzy sets characterzed by MFs. The MFs may be classfed accordng to have ether a lmted range or an unlmted range. For example, a trangular type covers a lmted range whle a sgmod type covers an unlmted range. In the trackng control, to cover a large unwanted error caused by dsturbances, ntal trackng error, or large scaled error, we should cover a wde range by MFs of nputs. For ths purpose, the MFs located n the left and rght sdes of nput can be sgmod type whle the MF n the center can be a Gaussan or trangular shape. When a fuzzy controller s desgned as a normalzed one, the scalng factors are used to scale the actual range of varables n the range of MFs. The archtecture of TS controller n (11) ncludes the scalng factors. Substtutng (11)nto(1) forms the closed loop system as follows: (c 1 (x) x 1 + c 2 (x) x 2 + c (x)) k o = R a I a + L a İ a + k b r 1 q k (12) Assume that the motor voltage u expressed by (2) s lmted such that R a I a + L a İ a + k b r 1 q k u max (13) where < u max s a maxmum permtted voltage for the motor. Ths assumpton s a techncal regard to protect motor aganst over voltages. The complexty of desgn and analyss has been changed to smplcty for usng the model of motor n place of model of manpulator. Here, we should know only the upper lmts for the motor voltages as nputs of robotc system. Because electrcal motors drve the electrcal manpulator, the motor voltages are the system nputs. The desred trajectory should

6 42 J Intell Robot Syst (21) 6: be planned wth regardng the maxmum permtted voltages for motors somehow each motor s so strong such that can track the desred trajectory under the permtted voltage. Moreover, the desred trajectory should be smooth such that ts dervatves up to the requred order are avalable and lmted. The drect method of Lyapunov s utlzed for asymptotc stablty analyss of TS fuzzy controller. To fnd a control law for a guaranteed stablty, a Lyapunov canddate s proposed as V (x) = x1 c 2 (x) x 1 dx 1, V () = (14) where V(x) s a postve defnte functon f c 2 (x) s postve. c 2 (x) s stated by (6)as c 2 (x) = 9 l=1 =1 9 2 μ A l (x ) a l 2 l=1 =1 2 μ A l (x ) To satsfy < c 2 (x) t s suffcent that a l 2 and l such that < 2 μ A l (x ) a l 2. Proof Assume that C 2 c 2 (x) where C 2 s a postve constant. Thus x1 C 2 x 1 dx 1 x1 =1 (15) c 2 (x) x 1 dx 1 (16) x1 We have C 2 x 1 dx 1 =.5C 2 x 2 1 and <.5C 2x 2 1 for x 1 =. Thus, (16) mples that < V(x) for x =. Snce c 2 (x) x 1 dx 1 = and c 2 (x)x 1 s lmted, V() =. Thus, V(x) s a postve defnte functon. For smplcty, the feature < c 2 (x) s used n replace of the explct form of c 2 (x). The tme dervatve of V(x) by the use of (1)scalculatedas V = c 2 (x) x 1 ẋ 1 = αc 2 (x) x 1 x 2 (17) From (12)wecanwrte c 2 (x) x 2 = c 1 (x) x 1 c (x) + ( R a I a + L a İ a + k b r 1 q k ) /ko (18) Substtutng (18)nto(17) yelds V = α ( c 1 (x) x 2 1 x ( 1c (x) + x 1 Ra I a + L a İ a + k b r 1 ) ) q k /ko (19) Snce c 1 (x) x 2 1 for c 1(x), tosatsfy V for stablty, t s requred that ( x 1 c (x) + x 1 Ra I a + L a İ a + k b r 1 ) q k /ko (2) It can be wrtten as (( x 1 Ra I a + L a İ a + k b r 1 ) q k ko c (x) ) (21) Assume that c (x) =. To establsh stablty, k o s selected as x 1 u max / (x 1 c (x)) k o (22)

7 J Intell Robot Syst (21) 6: Proof Equaton 21 can be wrtten as The use of (13) yelds x 1 ( Ra I a + L a İ a + k b r 1 q k ) ko x 1 c (x) (23) x 1 ( Ra I a + L a İ a + k b r 1 q k ) x1. R a I a + L a İ a + k b r 1 q k = x 1. u (24) To satsfy (23), t s suffcent that x 1. u k o x 1 c (x) (25) Snce < k o, to guarantee stablty < x 1 c (x). Thsmeansthatc (x) must be desgned wth the same sgn as x 1. Ths condton s smply satsfed f a l s selected wth the same sgn as x 1. Ths fact s understood by studyng the behavor of dc motor as stated by (2) where q k has a drect relaton wth u. Ths means that a postve u ncreases q k. Thus, when x 1 = k 1 (q dk q k ) >,.e. the actual poston q k s under the desred poston q kd then u should be postve to reduce x 1.Wth the same reasonng when x 1 < thenu should be negatve to reduce the sze of x 1. Therefore, the sgn u s selected the same as sgn of x 1 to reduce the sze of x 1. From (25)andx 1 c (x) >, we obtan That s x 1 u / (x 1 c (x)) k o (26) u / c (x) k o (27) To know c (x), we pay attenton on c (x) as gven by (6)as c (x) = 9 l=1 =1 9 2 μ A l (x ) a l l=1 =1 2 μ A l (x ) Snce c (x) s an average functon due to 2 μ A l (x ) 1, thus where a,mn = mn 9 { } a l and a,max = 9 l=1 =1 (28) a,mn c (x) a,max (29) max l=1 { a l }. To select a constant value, we should select a value for k o that satsfes (27) n all cases. On the other hand, a larger value for k o may provde an over voltage. Thus, to protect the motor from over voltages we should regard u u max. The maxmum permtted value for k o s then selected as k o = u max / a,max (3) where u max and a l are known n advance. If the selected value s not suffcently large to satsfy (27), we should select an stronger motor wth a larger u max that n the operatng range u / c (x) u max / a,max = k o (31)

8 422 J Intell Robot Syst (21) 6: Therefore, stablty s guaranteed under assumptons < c 1 (x),< c 2 (x),< x 1 c (x), and k o = u max / a,max. The output scalng factor plays an effectve role n stablty of closed loop system. We may nterest n the role of nput scalng factors. It can be wrtten that c (x) x = c 1 (x) x 1 + c 2 (x) x 2 (32) where c(x) = [ c 1 (x) c 2 (x) ]. Then, (5) srewrttenas c (x) x = ( ) R a I a + L a İ a + k b q k /ko c (x) (33) By the use of (7), we can wrte [ ][ ] k1 z1 x = = kz (34) k 2 z 2 [ ] k1 where k = and z = [ ] T.Hence(33)sformedas z k 1 z 2 2 c (x) kz = ( ) R a I a + L a İ a + k b q k /ko c (x) (35) Thus c (x). k. z R a I a + L a İ a + k b q k /k o + c (x) (36) Therefore, by the use of (13), the norm of error z s bounded as z < (u max /k o + c (x) ) / ( k. c (x) ) (37) Membershp functons of x 1 1 N Z P.8 Degree of membershp x 1 Fg. 1 The membershp functons for the nputs

9 J Intell Robot Syst (21) 6: Table 1 Fuzzy rules f (x) x 2 N Z P x 1 P Z.5 15x 1 + 1x N Ths mples that selectng larger values for k o, k 1 and k 2 provdes a smaller norm of error for z. Fuzzy rules n the 9 subsectons for l = 1,...,9 are desgned where the followng cases occur: Case 1 Assume that x 1 x 2 < resultng n asymptotc stablty by V < n (17). Thus, a small control effort s gven to u. Case 2 Assume that x 1 = or x 2 = resultng n Lyapunov stablty by V = n (17). Thus, a medum control effort s gven to u. Case 3 Assume that x 1 and x 2 both are postve or negatve resultng n nstablty by V > n (17). Thus, a very hgh effort s gven to u. Assumptons < c 1 (x), < c 2 (x), < x 1 c (x), andk o = u max / a,max must be satsfed, as well. In the l-th subsecton, a l s selected wth the same sgn as x 1 to satsfy < x 1 c (x). Wecanselecta l 1 and al 2 n all subsectons but l that a l 1 > and al 2 > to satsfy < c 2 (x), < x 1 c (x), respectvely. Fuzzy controller as a nonlnear functon 1.5 f(x) x x Fg. 2 The normalzed TS fuzzy controller as a nonlnear functon

10 424 J Intell Robot Syst (21) 6: The robust normalzed TS fuzzy controller s then smply desgned as follows: Choose the two nputs x 1 and x 2, for error and ts frst tme dervatve, and y or the output. a) Gve three MFs for each nput as shown n a range of [ 1 1]nFg.1. b) Defne nne fuzzy rules n the form of (3)asfollows: Rule 1 If x 1 s P and x 2 s P Then y = 1 Rule 2 If x 1 s P and x 2 s Z Then y =.75 Rule 3 If x 1 s P and x 2 s N Then y =.25 Rule 4 If x 1 s Z and x 2 s P Then y =.5 Rule 5 If x 1 s Z and x 2 s Z Then y = a 1 x 1 + a 2 x 2 Rule 6 If x 1 s Z and x 2 s N Then y =.5 Rule 7 If x 1 s N and x 2 s P Then y =.25 Rule 8 If x 1 s N and x 2 s Z Then y =.75 Rule 9 If x 1 s N and x 2 s N Then y = 1 Table 1 summarzed the fuzzy rules. In Rule 5, a 1 > and a 2 >. Rule 5 s specal snce t governs the equlbrum subsecton. The fuzzy system as a nonlnear functon of mappng from nput vector to output n the case of a 1 = 1 and a 2 = 1 s shown n Fg Smulaton The control system s smulated for trackng control of PUMA 56 Arm drven by brushed permanent magnet dc motors wth specfcatons [26] gven n Table 2. The nductances of motors are ncluded to consder a more complcated model n smulatons. Motors that drve the frst, second and thrd jonts are 4 V and 16 W [27]. Thus, the maxmum permtted voltage s gven 4 V. The smulaton model of PUMA 56 n MATLAB [28] s used n the smulatons. The desred trajectory s suffcently smooth such that all requred dervatves up to requred order s avalable as shown n Fg. 3. It starts from zero to reach a jont angle of 2rad n a perod of 4 S. The fuzzy controller desgned n Secton 3 s appled on the system whle the same motors, controllers and desred trajectores are used for the jonts. The dsturbance s nserted to the nput of each motor as a perodc pulse functon wth a perod of 2 S, ampltude of 2 V, tme delay of.7 S, and pulse wdth 3% of perod. Ths form of dsturbance s an example of any form that can be appled but t ncludes jumps to cover the complex cases. The output scalng factor k o s gven 4 as calculated from (31) to protect motors from over voltages Table 2 Parameters of permanent magnet dc motors Jont R L K b J B 1/r

11 J Intell Robot Syst (21) 6: Desred jont angle Jont angle (rad) tme (s) Fg. 3 Desred trajectory and overcome uncertantes n the operatng range of motors under the permtted voltages of motors. Smulaton 1 The control system s smulated for trackng control wth zero ntal errors. The nput scalng factors are set to k 1 = 1 and k 2 = 1. The gans of consequent part of Rule 5 are gven a 1 = 15 and a 2 = 1. The control system performs well as shownnfg.4 whle the maxmum values of errors for jonts 1 3 are , and rad, respectvely. The effects of dsturbances are represented n the left and rght sdes of Fg. 4 as small jumps on the curves of trackng errors. Snce the ntal errors are zero, the control system starts from equlbrum subsecton governed by Rule 5 and stays there as confrmed by the trackng errors. Therefore, the gans of a 1 and a 2 n ths rule, play a sgnfcant role n the performance of controller. They exhbt a behavor smlar as the proportonal and dervatve gans n an ordnary PD controller. The control system overcomes dsturbances very well as shownn Fg.4. The motor voltages behave well under the maxmum permtted value of 4V as shown n Fg. 5. The control efforts compensate the dsturbances such that the jumps on the curves respond to the jumps nserted by dsturbances otherwse the

12 426 J Intell Robot Syst (21) 6: x Performance of control system jont 1 jont 2 jont 3 5 Trackng erorr (rad) tme (s) Fg. 4 Performance of robust fuzzy control system curves are smooth. There are not any chatterng, oscllatons, and over voltages n the control efforts. Smulaton 2 Three dfferent values of 15, 3 and 6 are gven to a 1 n Rule 5 and the value of 1 s gven to a 2 for consderng the role of gan a 1 on trackng error as shown n Fg.6. Other condtons are the same as Smulaton 1. The control system behaves well n these cases, as well. The maxmum value of trackng error for jont 2 n the smulatons s reached to the value of rad, rad and rad n respondng to the values of 15, 3 and 6 gven to a 1, respectvely. Snce the control system s around the equlbrum pont, a larger a 1 provdes smaller trackng error. Smulaton 3 A very large dsturbance of 1 tmes as the value n Smulaton 1 s nserted n the control system whle other condtons are the same as Smulaton 1. The control system s robust aganst the large dsturbances as verfed by performance of trackng error and the voltages of motors shown n Fgs. 7 and 8, respectvely. It s very nterestng that the system does not go out of the equlbrum even n the case of a very large dsturbance. The control efforts behave well to compensate the

13 J Intell Robot Syst (21) 6: Control efforts n trackng control u 1 u 2 u 3 15 motor voltage (V) tme (S) Fg. 5 Voltages of motors to perform trackng control dsturbances. The maxmum value of trackng error for jont 2 s ncreased from rad n Smulaton 1 to a value of rad n Smulaton 3. Smulaton 4 The effect of ntal error s studed whle a value of 2 rad s gven as ntal errors to jonts. The nserted dsturbance s removed to pay attenton on the effect of ntal error whereas other condtons are the same as smulaton 1. Control system starts from a pont that s far away from the equlbrum; however t goes to the equlbrum well. The trackng errors approach the value of about rad whle they start from 2 rad as shown n Fg. 9. The responses behave smoothly to reduce the sze of trackng errors. The robustness of system n terms of stablty and trackng performance s presented by ths smulaton n drectng the trackng system to the equlbrum. The control efforts show jumps when startng to control the effect of nonzero ntal errors. They behave well under the maxmum permtted voltages as shown n Fg. 1. We can see the role of MFs of nputs namely P and N n the performance of control system. If we select the trangular shape n replace of the sgmod shape for the P and N MFs of nputs, then the startng pont wll not n the range of nputs. As a result, the control system wll not operate snce the control efforts wll be zero. Smulaton 5 The control system s smulated for set pont control where a desred value of 1 rad s gven to the jonts. Other condtons are the same as Smulaton 1. Set pont applcaton s used for pont-to-pont moton control of robot manpulators. In ndustry, the set pont control s a domnant approach n the process control, as well. The control system responses very well wthout over shoot and relatvely fast wth an gnorable steady state error and robust aganst dsturbances. The actual jont

14 428 J Intell Robot Syst (21) 6: x 1-3 Role of nput scalng factors a1 a2.8.7 Norm of trackng error (rad) tme (s) Fg. 6 Increasng nput scalng factors reduces the trackng error. a1 k 1 = 2, k 2 =.2, k o = 16, a2 k 1 = 1, k 2 = 1, k o = 16 poston reaches to a value of about 1 rad at tme of 4S as shown n Fg. 11. The effect of dsturbances s compensated such that the responses are not affected by them. The control efforts show jumps when startng to control the effect of hgh value of errors. They behave well to compensate the effects of dsturbances and reduce the errors under the maxmum permtted voltages as shown n Fg. 12. Smulaton 6 The role of nput scalng factors becomes more mportant for usng the trangular MFs for nputs. The nput scalng factor s employed to take the nput nto the operatng range covered by MFs otherwse the controller wll not respond to the nput. When the ntal error s very large ether n trackng control or n set pont control, the nput scalng factor s selected small. If the error s very small, then a large nput scalng factor can be selected to operate effcently. In ths smulaton, the trangular shape MFs of nputs are arranged n [ 1 1] as shown n Fg. 13. To control the system for a desred set pont of 2 rad, the gans of nput scalng factors are gven k 1 =.5 and k 2 = 1, and the gans of Rule 5 are gven a 1 = 1 and a 2 = 1. Other condtons are the same as Smulaton 5. The control system responses very

15 J Intell Robot Syst (21) 6: x Control performance under a very large dsturbance jont1 jont2 jont3 5 trackng error (rad) tme (S) Fg. 7 Control performance under a very large dsturbance 5 4 Control efforts under a very large dsturbance u 1 u 2 u 3 3 motor voltge (V) tme (S) Fg. 8 Control efforts under a very large dsturbance

16 43 J Intell Robot Syst (21) 6: jont1 jont2 jont3 trackng error (rad) tme (S) Fg. 9 Trackng performance under a large ntal error 3 u 1 2 u 2 u 3 1 voltage of motor (V) tme (S) Fg. 1 Control effort under a large ntal error

17 J Intell Robot Syst (21) 6: q d 1.8 q 1 q 2 q 3 jont angle (rad) tme (S) Fg. 11 Performance of set pont control u 1 u 2 u 3 25 motor voltge (V) tme (S) Fg. 12 Control efforts n set pont control

18 432 J Intell Robot Syst (21) 6: Trangular MFs for x 1 n Smulaton 6 N Z P 1.8 Degree of membershp x 1 Fg. 13 Trangular MFs for nputs n Smulaton Set pont control usng trangular MFs and scalng factors q d q 1 q 2 q jont angle (rad) tme (S) Fg. 14 Performance of set pont control usng trangular MFs and scalng factors

19 J Intell Robot Syst (21) 6: Control efforts n set pont control usng trangular MFs u 1 u 2 u 3 25 motor voltage (V) tme (S) Fg. 15 Control efforts n set pont control usng trangular MFs and scalng factors well wthout over shoot and relatvely fast wth an gnorable steady state error and robust aganst dsturbances. The actual jont poston reaches to a value of about 2 rad at tme of 4S as shown n Fg. 14. The control efforts behave well to compensate the effects of dsturbances and reduce the errors under the maxmum permtted voltages as shown n Fg Concluson A robust normalzed TS fuzzy controller was desgned for both trackng and set pont control of electrcal manpulators. The analyss and desgn of fuzzy control as a robust and nonlnear control approach was presented by nonlnear control theory usng the drect method of Lyapunov. The stablty was analyzed wthout knowng the explct dynamcs of system. Maxmum permtted voltages of motors are the only requred knowledge about the robotc system. The complextes n analyss and control desgn have been removed well usng the voltage-based control. Fnally, a smple algorthm was provded to desgn a robust normalzed TS fuzzy controller. References 1. Spong, M.W., Hutchnson, S., Vdyasagar, M.: Robot Modelng and Control. Wley, Hoboken (26)

20 434 J Intell Robot Syst (21) 6: Hodges, S.E.: Lookng for a cheaper robot: vsual feedback for automated PCB manufacture. Ph.D. thess n Unversty of Cambrdge (1996) 3. Qu, Z., Dawson, D.M.: Robust Trackng Control of Robot Manpulators. IEEE, New York (1996) 4. Abdallah, C., Dawson, D., Dorato, P., Jamshd, M.: Survey of robust control for rgd roots. Control Syst. Mag. 11, 24 3 (1991) 5. Cheah, C.C., Hrano, M., Kawamura, S., Armoto, S.: Approxmate Jacoban control for robots wth uncertan knematcs and dynamcs. IEEE Trans. Robot. Autom. 19(4), (23) 6. Fateh, M.M., Soltanpour, M.R.: Robust task-space control of robot manpulators under mperfect transformaton of control space. Int. J. Innov. Comput. Info. Control. 5(11A), (29) 7. Fateh, M.M.: Proper uncertanty bound parameter to robust control of electrcal manpulators usng nomnal model. Nonlnear Dyn. (21). do:1.17/s Fateh, M.M.: On the voltage-based control of robot manpulators. Int. J. Control. Autom. Syst. 6(5), (28) 9. Fateh, M.M.: Robust voltage control of electrcal manpulators n task-space. Int. J. Innov. Comput. Info. Control. 6(6), (21) 1. Wang, L.X.: A Course n Fuzzy Systems and Control. Prentce Hall, New York (1996) 11. Lm, C.M., Hyama, T.: Applcaton of fuzzy logc control to a manpulator. IEEE Trans. Robot. Autom. 1(5), (1991) 12. Ham, C., Qu, Z., Johnson, R.: Robust fuzzy control for robot manpulators. IEE Proc., Control Theory Appl. 147(2), (2) 13. Km, E.: Output feedback trackng control of robot manpulator wth model uncertanty va adaptve fuzzy logc. IEEE Trans. Fuzzy Syst. 12(3), (24) 14. Hwang, J.P., Km, E.: Robust trackng control of an electrcally drven robot: adaptve fuzzy logc approach. IEEE Trans. Fuzzy Syst. 14(2), (26) 15. Yng, H.: The Takag Sugeno fuzzy controllers usng the smplfed lnear control rules are nonlnear varable gan controllers. Automatca 34(2), (1998) 16. Tsay, D.L., Chung, H.Y., Lcc, C.J.: The adaptve control of nonlnear system usng the Sugenotype of fuzzy logc. IEEE Trans. Fuzzy Syst. 7(2), (1999) 17. Tsa, C.H., Wang, C.H., Ln, W.S.: Robust fuzzy model-followng control of robot manpulators. IEEE Trans. Fuzzy Syst. 8(4), (2) 18. Shen, Y., Ca, W.J., L, S.: Multvarable process control: decentralzed, decouplng, or sparse. Ind. Eng. Chem. Res. 49, (21) 19. Hsua, S.H., Fua, L.C.: A fully adaptve decentralzed control of robot manpulators. Automatca 42, (26) 2. Jn, Y.: Decentralzed adaptve fuzzy control of robot manpulators. IEEE Trans. Syst. Man Cybern., Part B, Cybern. 28(1), (1998) 21. Km, V.T.: Independent jont adaptve fuzzy control of robot manpulator. In: The 5th Bannual World Automaton Congress, vol. 14, pp (22) 22. Takag, T., Sugeno, M.: Fuzzy dentfcaton of systems and ts applcatons to modelng and control. IEEE Trans. Syst. Man. Cybern. 15, (1985) 23. Yng, H.: Suffcent condtons on unform approxmaton of multvarate functons by general Takag Sugeno fuzzy systems wth lnear rule consequent. IEEE Trans. Syst., Man. Cybernet. 28, (1998) 24. Yng, H.: An analytcal study on structure, stablty and desgn of general Takag Sugeno fuzzy control systems. Automatca 34, (1998) 25. Dng, Y., Yng, H., Shao, S.: Typcal Takag Sugeno PI and PD fuzzy controllers: analytcal structures and stablty analyss. Inf. Sc. 151, (23) 26. Corke, P.I., Armstrong-Hlouvry, B.: A search for consensus among model parameters reported for the PUMA 56 robot. Proc. IEEE Int. Conf. Robot. Autom. 1, (1994) 27. Wyeth, G.F., Kennedy, J., Lllywhte, J.: Dstrbuted dgtal control of a robot arm. In: Proceedngs of the Australan Conference on Robotcs and Automaton (ACRA 2), pp (2) 28. Corke, P.I.: Robotcs toolbox for {MATLAB}. IEEE Robot. Autom. Mag. 3(1), (1996)

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

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

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

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

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm Desgn and Optmzaton of Fuzzy Controller for Inverse Pendulum System Usng Genetc Algorthm H. Mehraban A. Ashoor Unversty of Tehran Unversty of Tehran h.mehraban@ece.ut.ac.r a.ashoor@ece.ut.ac.r Abstract:

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

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

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

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

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

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

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

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. EE 539 Homeworks Sprng 08 Updated: Tuesday, Aprl 7, 08 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. For full credt, show all work. Some problems requre hand calculatons.

More information

Adaptive Consensus Control of Multi-Agent Systems with Large Uncertainty and Time Delays *

Adaptive Consensus Control of Multi-Agent Systems with Large Uncertainty and Time Delays * Journal of Robotcs, etworkng and Artfcal Lfe, Vol., o. (September 04), 5-9 Adaptve Consensus Control of Mult-Agent Systems wth Large Uncertanty and me Delays * L Lu School of Mechancal Engneerng Unversty

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

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

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

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

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

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

The Analysis of Coriolis Effect on a Robot Manipulator

The Analysis of Coriolis Effect on a Robot Manipulator Internatonal Journal of Innovatons n Engneerng and echnology (IJIE) he Analyss of Corols Effect on a Robot Manpulator Pratap P homas Assstant Professor Department of Mechancal Engneerng K G Reddy college

More information

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several

More information

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &

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

Robust observed-state feedback design. for discrete-time systems rational in the uncertainties

Robust observed-state feedback design. for discrete-time systems rational in the uncertainties Robust observed-state feedback desgn for dscrete-tme systems ratonal n the uncertantes Dmtr Peaucelle Yosho Ebhara & Yohe Hosoe Semnar at Kolloquum Technsche Kybernetk, May 10, 016 Unversty of Stuttgart

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

Power law and dimension of the maximum value for belief distribution with the max Deng entropy

Power law and dimension of the maximum value for belief distribution with the max Deng entropy Power law and dmenson of the maxmum value for belef dstrbuton wth the max Deng entropy Bngy Kang a, a College of Informaton Engneerng, Northwest A&F Unversty, Yanglng, Shaanx, 712100, Chna. Abstract Deng

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

Fuzzy Boundaries of Sample Selection Model

Fuzzy Boundaries of Sample Selection Model Proceedngs of the 9th WSES Internatonal Conference on ppled Mathematcs, Istanbul, Turkey, May 7-9, 006 (pp309-34) Fuzzy Boundares of Sample Selecton Model L. MUHMD SFIIH, NTON BDULBSH KMIL, M. T. BU OSMN

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

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

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

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

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

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

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

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

Control of Uncertain Bilinear Systems using Linear Controllers: Stability Region Estimation and Controller Design

Control of Uncertain Bilinear Systems using Linear Controllers: Stability Region Estimation and Controller Design Control of Uncertan Blnear Systems usng Lnear Controllers: Stablty Regon Estmaton Controller Desgn Shoudong Huang Department of Engneerng Australan Natonal Unversty Canberra, ACT 2, Australa shoudong.huang@anu.edu.au

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

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

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

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,

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

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

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set

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

Some modelling aspects for the Matlab implementation of MMA

Some modelling aspects for the Matlab implementation of MMA Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton

More information

Automatic PID Controller Tuning for Robots with Nonlinear Friction at the Joints

Automatic PID Controller Tuning for Robots with Nonlinear Friction at the Joints Automatc PID Controller unng for Robots wth Nonlnear Frcton at the Jonts Abílo Azenha, Ph.D. Abstract hs paper descrbes an approach to automatc tunng of PID poston controllers for manpulators wth nonlnear

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

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v10157-01-00-y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION

More information

Complement of Type-2 Fuzzy Shortest Path Using Possibility Measure

Complement of Type-2 Fuzzy Shortest Path Using Possibility Measure Intern. J. Fuzzy Mathematcal rchve Vol. 5, No., 04, 9-7 ISSN: 30 34 (P, 30 350 (onlne Publshed on 5 November 04 www.researchmathsc.org Internatonal Journal of Complement of Type- Fuzzy Shortest Path Usng

More information

Adaptive sliding mode reliable excitation control design for power systems

Adaptive sliding mode reliable excitation control design for power systems Acta Technca 6, No. 3B/17, 593 6 c 17 Insttute of Thermomechancs CAS, v.v.. Adaptve sldng mode relable exctaton control desgn for power systems Xuetng Lu 1, 3, Yanchao Yan Abstract. In ths paper, the problem

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

x i1 =1 for all i (the constant ).

x i1 =1 for all i (the constant ). Chapter 5 The Multple Regresson Model Consder an economc model where the dependent varable s a functon of K explanatory varables. The economc model has the form: y = f ( x,x,..., ) xk Approxmate ths by

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

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations Physcs 178/278 - Davd Klenfeld - Wnter 2015 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons We consder a network of many neurons, each of whch obeys a set of conductancebased,

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

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

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

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

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

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

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

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

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

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

THE GUARANTEED COST CONTROL FOR UNCERTAIN LARGE SCALE INTERCONNECTED SYSTEMS

THE GUARANTEED COST CONTROL FOR UNCERTAIN LARGE SCALE INTERCONNECTED SYSTEMS Copyrght 22 IFAC 5th rennal World Congress, Barcelona, Span HE GUARANEED COS CONROL FOR UNCERAIN LARGE SCALE INERCONNECED SYSEMS Hroak Mukadan Yasuyuk akato Yoshyuk anaka Koch Mzukam Faculty of Informaton

More information

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse

More information

Modelling of the precise movement of a ship at slow speed to minimize the trajectory deviation risk

Modelling of the precise movement of a ship at slow speed to minimize the trajectory deviation risk Computatonal Methods and Expermental Measurements XIV 29 Modellng of the precse movement of a shp at slow speed to mnmze the trajectory devaton rsk J. Maleck Polsh Naval Academy, Poland Faculty of Mechancs

More information

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning Journal of Machne Learnng Research 00-9 Submtted /0; Publshed 7/ Erratum: A Generalzed Path Integral Control Approach to Renforcement Learnng Evangelos ATheodorou Jonas Buchl Stefan Schaal Department of

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

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

Polynomial Regression Models

Polynomial Regression Models LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance

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

Research Article Green s Theorem for Sign Data

Research Article Green s Theorem for Sign Data Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of

More information

Markov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement

Markov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement Markov Chan Monte Carlo MCMC, Gbbs Samplng, Metropols Algorthms, and Smulated Annealng 2001 Bonformatcs Course Supplement SNU Bontellgence Lab http://bsnuackr/ Outlne! Markov Chan Monte Carlo MCMC! Metropols-Hastngs

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

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results. Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

arxiv:cs.cv/ Jun 2000

arxiv:cs.cv/ Jun 2000 Correlaton over Decomposed Sgnals: A Non-Lnear Approach to Fast and Effectve Sequences Comparson Lucano da Fontoura Costa arxv:cs.cv/0006040 28 Jun 2000 Cybernetc Vson Research Group IFSC Unversty of São

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

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem H.K. Pathak et. al. / (IJCSE) Internatonal Journal on Computer Scence and Engneerng Speedng up Computaton of Scalar Multplcaton n Ellptc Curve Cryptosystem H. K. Pathak Manju Sangh S.o.S n Computer scence

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

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1]

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1] DYNAMIC SHORTEST PATH SEARCH AND SYNCHRONIZED TASK SWITCHING Jay Wagenpfel, Adran Trachte 2 Outlne Shortest Communcaton Path Searchng Bellmann Ford algorthm Algorthm for dynamc case Modfcatons to our algorthm

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

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

Chapter 2 A Class of Robust Solution for Linear Bilevel Programming

Chapter 2 A Class of Robust Solution for Linear Bilevel Programming Chapter 2 A Class of Robust Soluton for Lnear Blevel Programmng Bo Lu, Bo L and Yan L Abstract Under the way of the centralzed decson-makng, the lnear b-level programmng (BLP) whose coeffcents are supposed

More information

Group Analysis of Ordinary Differential Equations of the Order n>2

Group Analysis of Ordinary Differential Equations of the Order n>2 Symmetry n Nonlnear Mathematcal Physcs 997, V., 64 7. Group Analyss of Ordnary Dfferental Equatons of the Order n> L.M. BERKOVICH and S.Y. POPOV Samara State Unversty, 4430, Samara, Russa E-mal: berk@nfo.ssu.samara.ru

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

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

Modelli Clamfim Equazione del Calore Lezione ottobre 2014 CLAMFIM Bologna Modell 1 @ Clamfm Equazone del Calore Lezone 17 15 ottobre 2014 professor Danele Rtell danele.rtell@unbo.t 1/24? Convoluton The convoluton of two functons g(t) and f(t) s the functon (g

More information

Markov Chain Monte Carlo Lecture 6

Markov Chain Monte Carlo Lecture 6 where (x 1,..., x N ) X N, N s called the populaton sze, f(x) f (x) for at least one {1, 2,..., N}, and those dfferent from f(x) are called the tral dstrbutons n terms of mportance samplng. Dfferent ways

More information