Time dependent weight functions for the Trajectory Piecewise-Linear approach?
|
|
- Gary Lee
- 5 years ago
- Views:
Transcription
1 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 ynamcal systems (.e an more equatons). One technque for non-lnear MOR s the trajectory pecewse-lnear approach (TPWL). TPWL approxmates a non-lnear fferental system by a weghte sum of lnear systems whch have a sgnfcantly reuce number of equatons. One open queston s whch weghtng schemes are approprate. We scuss whether a tme epenent weghtng scheme whch s computatonally faster as well as computatonally cheaper than the orgnally propose one s approprate. 1 The Trajectory Pecewse-Lnear approach (TPWL) Moel orer reucton (MOR) spees up the smulaton of hgh mensonal ynamcal systems by reucng ts menson whle keepng the output error small. MOR for lnear tme-nvarant systems s well unerstoo an use n nustry, however for non-lnear systems there s stll much research to o. The trajectory pecewse-lnear approach [4] s a numercal metho for the moel orer reucton of non-lnear systems escrbe by frst-orer, non-lnear ornary fferental equatons an non-lnear fferental-algebrac equatons. Here we conser the ornary fferental equaton: Juan Pablo Amorocho D. Insttut Computatonal Mathematcs, Technsche Unverstät Braunschweg, D Braunschweg, e-mal: jamoroch@tu-braunschweg.e Heke Faßbener Insttut Computatonal Mathematcs, Technsche Unverstät Braunschweg, D Braunschweg, e-mal: h.fassbener@tu-braunschweg.e 1
2 2 Juan Pablo Amorocho an Heke Faßbener t x(t) = f(x(t))+dv(t), 0 t T, x(0) = x 0, (1a) y(x(t)) = L T x(t), (1b) where x(t) R N s the state vector, f : R N R N s a non-lnear functon, D R N M s the nput matrx, v(t) R M s the nput, L T R P N s the matrx mappng the state vector nto the output, an y(x(t)) R P s the output of the system. The TPWL reuces the hgh computatonal cost of the evaluaton of the non-lnearty by approxmatng equaton (1a) through a non-lnear convex combnaton of s fferent lnear moels aroun (t,x ), = 0,...,s 1 where x = x(t ). Let ˆf(x(t)) be the frst orer Taylor approxmaton of f(x(t)) aroun x : ˆf(x(t)) = f + F (x(t) x ), where f = f(x ) an F R N N s the Jacoban matrx of f(x(t)) evaluate at x. Wth ths, the pecewse approxmaton of system (1) s gven by s 1 t ˆx(t) = ŵ ( ˆx(t))(F ˆx(t)+B u(t)) =0 (2a) ŷ( ˆx(t)) = L T ˆx(t) (2b) wth the weghtng functons ŵ : R N [0,1], an s 1 =0 ŵ(x(t)) = 1 for all t [0,T], B = [D, f F x ] R N (M+1), an u(t) = [v(t),1] T R M+1. The tuples (t,x ) are foun va a smulaton of system (1a), see [4] for more etals. The resultng trajectory s calle the tranng trajectory. Next, lnear moel orer reucton s apple to the system (2) as follows: Bul the projector π q = V q W T q, such that q N an V q,w q R N q,w T q V q = I q (I q : entty matrx of menson q). Next, use π q to apply a Petrov-Galerkn projecton to system (2). The resultng ynamcal system s calle the TPWL moel of system (1) gven by: s 1 t x r(t) = =0 y r (x r (t)) = (L r ) T x r (t), w (x r (t))(f r x r(t)+b r u(t)), (3a) (3b) where x r (t) R q, an w : R q [0,1], s 1 =0 w (x r (t)) = 1 for all t [0,T]. Matrces F r,br an L r are gven by F r = Wq T F V q,b r = W q T B,(L r ) T = L T V q. Note that x r (t) = Wq T ˆx(t). The use of the weghts w (x r (t)) n (3a) nstea of ŵ ( ˆx(t)) ams to reuce the cost of the computaton of the weght functon from O(Ns) to O(qs). In case the non-lnearty of (1a) s not too pronounce, the TPWL moel (3a) can prect the evoluton of the soluton trajectory even for ntal contons an/or nput sgnal fferent from the ones use to set up the TPWL moel. In [4] two algorthms to bul the projector π q are propose: One that nclues all lnear moels an another that only uses the lnear moel at the ntal conton. The latter s faster than the former, but the former yels a better approxmaton than the latter. Furthermore, [4] uses Krylov subspace methos to bul π q. However ths s not the only opton, [7] employs balance truncaton methos, an [1] combnes the proper orthogonal ecomposton metho (POD) together wth the TPWL approach.
3 Tme epenent weght functons for the Trajectory Pecewse-Lnear approach? 3 In our mplementaton we use the algorthm gven n [3] usng only the lnear moel at the ntal conton an settng W q = V q whch yels an orthogonal projector. In bulng the pecewse approxmaton (2) two questons arse: () how to select the lnear moels along the tranng trajectory? () how to choose sutable weght functons? In [4] the extracton starts wth the lnear moel foun at the ntal conton x 0 an contnues to smulate an select lnear moels keepng the tranng trajectory covere by balls B α x, = 0,...,s 1 of constant raus α each centere at x. Startng at the lnear moel at x 0, a new lnear moel s selecte at (t,x ) when x(t ) oes not longer belong to B α x 1. Ths process s contnue untl x(t) s reache. Here we aopt the proceure from [4]. Concernng the secon queston, t s far to say that sutable weght functons shoul aequately approxmate the system s non-lnearty by the resultng convex combnaton of the lnearze moels. In [4] the weghts are gven by w ( x r (t)) = e β /m s 1 =0 e β /m, = 0,...,s 1 (4) where = x r (t) x 2, an m = mn an β s a postve parameter usually =0,...,s 1 set to 25. Ths ensures that the weghts change raply as x r (t) evolves n the state space. Here, x are gven by [ x 0, x 1,..., x s 1 ] = [W T x 0,W T x 1,...,W T x s 1 ], where ether W = W q as n (3a) or W results from the orthonormalzaton of the columns of W = {W q,x 0,..., x s 1 } an x r (t) s gven by x r (t) = W T ˆx(t). In the latter case, n (3a) W s use nstea of W q. However, ths ncreases the menson of the reuce space whch s unesrable. An alternatve to compute wthout havng to bul W s to project x r (t) back to the full state space,.e. V q x r (t) an compute the as V q x r (t) x 2 [5]. We use ths approach n our experments. We refer to the reaer to [2, 6] for other weght functons. 2 Tme epenent weght functons State-epenent weght functon are qute expensve to evaluate. Therefore, one mght want to conser weght functons whch are cheaper to evaluate, e.g., weght functons that are only a functon of tme. The ea behn these weght functons s smple: At each t we efne a pecewse weght functon w (t) as w (l) (t) f t [t 1,t ] w (t) = w (r) (t) f t [t,t +1 ] for = 0,...,s 1, 0 otherwse where w (l) 0 (t) = w(r) s 1 (t) an w (r) (t) are cubc poly- (t 1 ) = t w(l) (t ) = 0 nomals such that w (l) (t) = 0. Here, the functons w(l) (t ) = w (r) (t ) = 1, w (l) (5) (t 1 ) = t w(l)
4 4 Juan Pablo Amorocho an Heke Faßbener Fg. 1 Top: State-epenent weghts (4) from the tranng nput. Mle: Tme-epenent weghts (5) from tranng an test nput. Bottom: state-epenent weghts (4) from test nput. an w (r) (t +1 ) = t w(r) (t ) = t w(r) (t +1 ) = 0. If t s 1 < T, then we nee an atonal conton for w s 1 : w s 1 (t) = 1,t t s 1. Wth these cubc splnes, (3a) approxmates (1a) for each t [t,t +1 ] by at most two lnear fferental equatons, t x r(t) = w (r) (x r (t))(f r x r (t)+b r u(t))+w (l) +1 (x r(t))(f+1x r r (t)+b r +1u(t)). (6) Fgure 1 shows an example of tme-epenent an state-epenent weghts for the example consere n Secton 3. The TPWL moel constructe s smulate once for a tranng nput (whch was also use for smulatng the orgnal problem when settng up the TPWL moel) an once for a test nput. The tme-epenent weght functon oes not change for fferent nput functons, whle the stateepenent one oes change. As can be seen, consere as functons of tme, n case of the tranng nput, both weghtng functons o have a smlar behavor, the state-epenent weghts are almost pecewse constant (ether 0 or 1), whle the tme-epenent weght functons are smoother an allow (n the sense of the above equaton) for a better overlap of the two lnear moels whch lve on the nterval [t,t +1 ]. But for the test nput, the state-epenent weghts allow for a much stronger swtch between the fferent lnear moels than the tme-epenent weghts o. The tme-epenent weghts ffer from those n [2, 4, 6] n several ways. Frst, the computaton of the tme-epenent weghts s O(s) nstea of O(qs), as the computaton of the tme epenent weghts requres an evaluaton of a scalar polynomal of egree 3 per tme step. In contrast, n orer to compute the weght functons (4) sums of vectors of menson q have to be etermne n every weght w per tme step. Tme-epenent weght functons not only spee up the TPWL, but also reuce the storage, as there s no nee to bul the bass W as escrbe at the en of the prevous secton. Atonally, the use of these new weghts makes the pecewse system a lnear tme epenent one as the weghts are no longer state epene. Ths woul enable the use of the theory of lnear tme epenent systems to analyze the TPWL moel. But (6) has a serous problem when the TPWL moel s smulate wth ntal contons an/or nput sgnal fferent from the ones use to set up the TPWL
5 Tme epenent weght functons for the Trajectory Pecewse-Lnear approach? 5 1 N Fg. 2 Non-lnear transmsson lne. [4] moel. In that case, there s no guarantee that, for nstance, n the frst tme nterval [0,t 1 ] the lnear moels F 0 an F 1 are the best lnear moels to approxmate the ynamcs. It mght be better to use, say, F j an F k. State-epenent weght functons mght be able to allow for ths choce, a tme-epenent weght functon can not acheve ths, as t wll always be local wth respect to tme. Therefore, the use of tme-epenent weght functons s very lmte, unless the TPWL moel s bul from several smulaton runs wth fferent parameter settng for the ntal conton an nput. In that case one mght be able to mofy the tme-epenent weghts approach above to be comparatve to state-epenent weghts. 3 Numercal Experments We have teste our weght representatons on a non-lnear transmsson lne [4] shown n Fg. 2. Applyng the mofe noal analyss to ths crcut yels the ornary fferental equaton x 1 (t) = 2x 1(t)+x 2 (t) 1,0 1,2 + v(t) x j(t) = x j 1 (t) 2x j (t)+x j+1 (t)+ j 1, j j, j+1, j = 2,...,N 1 x N (t) = x N 1 x N + N 1,N, where x k s the kth noe voltage an k,l (t) = exp(40(x k (t) x l (t))) 1 s the oe s current between noes k an l. l = 0 means the groun noe. We have smulate the above crcut usng the TPWL wth two fferent nputs, one for the tranng nput 1 f 0 t < 0.5 v(t) = 1 t otherwse an one for the test nput v(t) = cos( π 2 t), an wth two fferent weght functons: the state epenent weght functon (4) an the tme epenent weght functon (5). In the fgures an tables below we refer to each case as Rewensk an splnes, respectvely. The smulatons run from 0 to 10 tme unts, the ntal conton s zero, N = 100, an the output y(t) s the voltage at noe 1. The lnear moels of the TPWL moel are extracte usng the algorthm escrbe n [4], the number of extracte lnear moels s 10 an the matrx V q s bul usng the algorthm of [3] an has rank 10. The smulaton was one usng MATLAB R2007b mplct ODE solver
6 6 Juan Pablo Amorocho an Heke Faßbener TPWL smulaton of non lnear transmsson lne usng fferent weghts 0.04 Voltage noe 1 (V) Splnes tme Exact soluton Rewensk state Tme Fg. 3 Smulaton of a non-lnear transmsson lne usng TPWL. oe15s on a Intel Core 2 DUO 1.80Ghz wth 2.0Gb of RAM memory runnng a Lnux kernel Fgure 3 shows the output of the smulaton, where the exact soluton s obtane by a smulaton of the full system. Ths s a worse case scenaro where the shape of the tranng nput s fferent from the test nput. The ba performance of the tmeepenent weghts can be unerstoo by lookng at the weght functons n Fgure 1. The test nput requres the TPWL to strongly swtch between the lnear moels, somethng moel (6) oes not o. One mght be tempte to reverse the orer of the nputs so that the TPWL s set up wth a more oscllatng nput. However our experments show that moel (6) wll stll perform poorly, but more nterestngly moel (1) usng weghts (4) ntally follows the exact soluton, but eventually verges. Acknowlegements Ths work, part of the SyreNe network, s supporte by the German Feeral Mnstry of Eucaton an Research (BMBF) grant no. 03FAPAE2. References 1. Bechtol, T., Strebel, M., Mohaghegh, K., ter Maten, E.J.W.: Nonlnear moel orer reucton n nanoelectroncs: Combnaton of POD an TPWL. In: PAMM, vol. 8, pp. 10,057 10,060. WILEY-VCH Verlag GmbH & Co. (2008). 79th Annual Meetng of the Internatonal Assocaton of Apple Mathematcs an Mechancs (GAMM) 2. Dong, N., Roychowhury, J.: Automate nonlnear macromoelng of output buffers for hghspee gtal applcatons. In: ACM/IEEE Desgn Automaton Conference, pp (2005). 3. Oabasoglu, A., Celk, M., Plegg, L.T.: PRIMA: passve reuce-orer nterconnect macromoelng algorthm. IEEE Transactons on Computer-Ae Desgn of Integrate Crcuts an Systems 17(8), (1998). DOI / Rewensky, M.: A trajectory pecewse-lnear approach to moel orer reucton of nonlnear ynamcal systems. Ph.D. thess, Massachusetts Insttute of Technology (2003) 5. Strebel, M., Rommes, J.: Moel orer reucton of non-lnear systems: status, open ssues, an applcatons. Tech. Rep. CSC/08-07, Technsche Unverstät Chemntz (2008) 6. Twary, S., Rutenbar, R.A.: Scalable trajectory methos for on-eman analog macromoel extracton. In: ACM/IEEE Desgn Automaton Conference, pp (2005). 7. Vaslev, D., Rewensky, M., Whte, J.: A TBR-base trajectory pecewse-lnear algorthm for generatng accurate low-orer moels for analog crcuts an mems. DAC pp (2003)
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 informationTransfer 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 informationConvexity 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 informationNumerical 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 informationVisualization of 2D Data By Rational Quadratic Functions
7659 Englan UK Journal of Informaton an Computng cence Vol. No. 007 pp. 7-6 Vsualzaton of D Data By Ratonal Quaratc Functons Malk Zawwar Hussan + Nausheen Ayub Msbah Irsha Department of Mathematcs Unversty
More informationSTORM: A Nonlinear Model Order Reduction Method via Symmetric Tensor Decomposition
STORM: A Nonlnear Moel Orer Reucton Metho va Symmetrc Tensor Decomposton Jan Deng, Haotan Lu, Km Batseler, Yu-Kwong Kwok an Nga Wong Department of Electrcal an Electronc Engneerng The Unversty of Hong
More informationModule 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 informationChapter 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 informationNew 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 informationSOLVING NON-LINEAR SYSTEMS BY NEWTON s METHOD USING SPREADSHEET EXCEL Tay Kim Gaik Universiti Tun Hussein Onn Malaysia
SOLVING NON-LINEAR SYSTEMS BY NEWTON s METHOD USING SPREADSHEET EXCEL Tay Km Gak Unverst Tun Hussen Onn Malaysa Kek Se Long Unverst Tun Hussen Onn Malaysa Rosmla Abdul-Kahar
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More informationOn 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 informationPolynomial 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 informationCSci 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 informationOn 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 information4DVAR, 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 informationLeast squares cubic splines without B-splines S.K. Lucas
Least squares cubc splnes wthout B-splnes S.K. Lucas School of Mathematcs and Statstcs, Unversty of South Australa, Mawson Lakes SA 595 e-mal: stephen.lucas@unsa.edu.au Submtted to the Gazette of the Australan
More informationHongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)
ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of
More informationProblem 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 informationGrover s Algorithm + Quantum Zeno Effect + Vaidman
Grover s Algorthm + Quantum Zeno Effect + Vadman CS 294-2 Bomb 10/12/04 Fall 2004 Lecture 11 Grover s algorthm Recall that Grover s algorthm for searchng over a space of sze wors as follows: consder the
More informationNON-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 informationENGI9496 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 informationNON-LINEAR CONVOLUTION: A NEW APPROACH FOR THE AURALIZATION OF DISTORTING SYSTEMS
NON-LINEAR CONVOLUTION: A NEW APPROAC FOR TE AURALIZATION OF DISTORTING SYSTEMS Angelo Farna, Alberto Belln and Enrco Armellon Industral Engneerng Dept., Unversty of Parma, Va delle Scenze 8/A Parma, 00
More informationOne-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 information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationSupporting Information
Supportng Informaton The neural network f n Eq. 1 s gven by: f x l = ReLU W atom x l + b atom, 2 where ReLU s the element-wse rectfed lnear unt, 21.e., ReLUx = max0, x, W atom R d d s the weght matrx to
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationReview of Taylor Series. Read Section 1.2
Revew of Taylor Seres Read Secton 1.2 1 Power Seres A power seres about c s an nfnte seres of the form k = 0 k a ( x c) = a + a ( x c) + a ( x c) + a ( x c) k 2 3 0 1 2 3 + In many cases, c = 0, and the
More informationNote 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 informationAnnexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances
ec Annexes Ths Annex frst llustrates a cycle-based move n the dynamc-block generaton tabu search. It then dsplays the characterstcs of the nstance sets, followed by detaled results of the parametercalbraton
More informationAppendix 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 informationOn 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 informationMarkov 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 informationExplicit bounds for the return probability of simple random walk
Explct bouns for the return probablty of smple ranom walk The runnng hea shoul be the same as the ttle.) Karen Ball Jacob Sterbenz Contact nformaton: Karen Ball IMA Unversty of Mnnesota 4 Ln Hall, 7 Church
More informationSome 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 informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationStructure 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 informationCubic 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 informationParallel MAC Based On Radix-4 & Radix-8 Booth Encodings
Shankey Goel. et al. / Internatonal Journal of Engneerng Scence an Technology (IJEST) Parallel MAC Base On Rax-4 & Rax-8 Booth Encongs SHAKE GOEL Stuent, Department of Electroncs & Communcaton Engneerng
More informationC/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1
C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned
More informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More informationCHAPTER 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 informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationParameter Estimation for Dynamic System using Unscented Kalman filter
Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,
More informationSuppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl
RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com
More informationLarge-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 information10-701/ Machine Learning, Fall 2005 Homework 3
10-701/15-781 Machne Learnng, Fall 2005 Homework 3 Out: 10/20/05 Due: begnnng of the class 11/01/05 Instructons Contact questons-10701@autonlaborg for queston Problem 1 Regresson and Cross-valdaton [40
More informationComparison of Wiener Filter solution by SVD with decompositions QR and QLP
Proceedngs of the 6th WSEAS Int Conf on Artfcal Intellgence, Knowledge Engneerng and Data Bases, Corfu Island, Greece, February 6-9, 007 7 Comparson of Wener Flter soluton by SVD wth decompostons QR and
More informationStabilizing Schemes for Piecewise-Linear Reduced Order Models via Projection and Weighting Functions
Stablzng Schemes for Pecewse-Lnear Reduced Order Models va Projecton and Weghtng Functons Bradley N Bond Research Laboratory n Electroncs, Massachusetts Insttute of Technology, 77 Massachusetts Ave, Cambrdge,
More informationModeling curves. Graphs: y = ax+b, y = sin(x) Implicit ax + by + c = 0, x 2 +y 2 =r 2 Parametric:
Modelng curves Types of Curves Graphs: y = ax+b, y = sn(x) Implct ax + by + c = 0, x 2 +y 2 =r 2 Parametrc: x = ax + bxt x = cos t y = ay + byt y = snt Parametrc are the most common mplct are also used,
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus
More informationNice plotting of proteins II
Nce plottng of protens II Fnal remark regardng effcency: It s possble to wrte the Newton representaton n a way that can be computed effcently, usng smlar bracketng that we made for the frst representaton
More informationChaper 2: Stress in beams
Chaper : Stress n eams FLEURE Beams suject to enng wll fle COPRESSON TENSON On the lower surface the eam s stretche lengthwse. Ths sujects t to tensle stress. N.A. N.A. s the neutral as On the upper surface
More informationMultiline Overview. Karl Bois, PhD Lead SI Engineer Servers & Blades. December 5th, 2016
Multlne Overvew Karl Bos, PhD Lead SI Engneer Servers & Blades December 5th, 216 Why are we here? Orgnal Equpment Manufacturers (OEMs) and Orgnal Desgn Manufacturers (ODMs) have desgns that requre that
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationLecture 9 Sept 29, 2017
Sketchng Algorthms for Bg Data Fall 2017 Prof. Jelan Nelson Lecture 9 Sept 29, 2017 Scrbe: Mtal Bafna 1 Fast JL transform Typcally we have some hgh-mensonal computatonal geometry problem, an we use JL
More informationOn a direct solver for linear least squares problems
ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationNeuro-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 informationMMA 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 informationMechanics Physics 151
Mechancs Physcs 5 Lecture 3 Contnuous Systems an Fels (Chapter 3) Where Are We Now? We ve fnshe all the essentals Fnal wll cover Lectures through Last two lectures: Classcal Fel Theory Start wth wave equatons
More informationPop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing
Advanced Scence and Technology Letters, pp.164-168 http://dx.do.org/10.14257/astl.2013 Pop-Clc Nose Detecton Usng Inter-Frame Correlaton for Improved Portable Audtory Sensng Dong Yun Lee, Kwang Myung Jeon,
More informationNUMERICAL 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 informationDistance-Based Approaches to Inferring Phylogenetic Trees
Dstance-Base Approaches to Inferrng Phylogenetc Trees BMI/CS 576 www.bostat.wsc.eu/bm576.html Mark Craven craven@bostat.wsc.eu Fall 0 Representng stances n roote an unroote trees st(a,c) = 8 st(a,d) =
More informationREAL 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 informationLarge-Scale Kernel Methods - I
Large-Scale Kernel Methos - I Sanjv Kumar, Google Research, NY EECS-6898, Columba Unversty - Fall, 2010 EECS6898 Large Scale Machne Learnng 1 Lnear Moels Popular n machne learnng / Statstcs ue to ther
More informationMatrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD
Matrx Approxmaton va Samplng, Subspace Embeddng Lecturer: Anup Rao Scrbe: Rashth Sharma, Peng Zhang 0/01/016 1 Solvng Lnear Systems Usng SVD Two applcatons of SVD have been covered so far. Today we loo
More informationFinding Dense Subgraphs in G(n, 1/2)
Fndng Dense Subgraphs n Gn, 1/ Atsh Das Sarma 1, Amt Deshpande, and Rav Kannan 1 Georga Insttute of Technology,atsh@cc.gatech.edu Mcrosoft Research-Bangalore,amtdesh,annan@mcrosoft.com Abstract. Fndng
More informationU.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 informationUsing the estimated penetrances to determine the range of the underlying genetic model in casecontrol
Georgetown Unversty From the SelectedWorks of Mark J Meyer 8 Usng the estmated penetrances to determne the range of the underlyng genetc model n casecontrol desgn Mark J Meyer Neal Jeffres Gang Zheng Avalable
More informationOn the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros
Appled Mathematcal Scences, Vol. 5, 2011, no. 75, 3693-3706 On the Interval Zoro Symmetrc Sngle-step Procedure for Smultaneous Fndng of Polynomal Zeros S. F. M. Rusl, M. Mons, M. A. Hassan and W. J. Leong
More informationExample: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,
The greatest common dvsor of two ntegers a and b (not both zero) s the largest nteger whch s a common factor of both a and b. We denote ths number by gcd(a, b), or smply (a, b) when there s no confuson
More informationCurve Fitting with the Least Square Method
WIKI Document Number 5 Interpolaton wth Least Squares Curve Fttng wth the Least Square Method Mattheu Bultelle Department of Bo-Engneerng Imperal College, London Context We wsh to model the postve feedback
More informationDesign, Control and Application of Modular Multilevel Converters for HVDC Transmission Systems by Kamran Sharifabadi, Lennart Harnefors, Hans Peter
Desgn, Control an Applcaton o Moular Multleel Conerters or HVDC Transmsson Systems by Kamran Sharaba, Lennart Harneors, Hans Peter Nee, Staan Norrga, Remus Teoorescu ISBN 0: 885560 Copyrght Wley 06 Chapter
More informationSIMPLIFIED 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 informationBoise State University Department of Electrical and Computer Engineering ECE 212L Circuit Analysis and Design Lab
Bose State Unersty Department of Electrcal and omputer Engneerng EE 1L rcut Analyss and Desgn Lab Experment #8: The Integratng and Dfferentatng Op-Amp rcuts 1 Objectes The objectes of ths laboratory experment
More informationParametric 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 informationSimulated 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 informationGeneralized 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 informationReport 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 informationErrors 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 informationBoise State University Department of Electrical and Computer Engineering ECE 212L Circuit Analysis and Design Lab
Bose State Unersty Department of Electrcal and omputer Engneerng EE 1L rcut Analyss and Desgn Lab Experment #8: The Integratng and Dfferentatng Op-Amp rcuts 1 Objectes The objectes of ths laboratory experment
More informationCHAPTER 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 informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationDO 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 informationHigh-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 informationAn efficient method for computing single parameter partial expected value of perfect information
An effcent metho for computng sngle parameter partal expecte value of perfect nformaton Mark Strong,, Jeremy E. Oakley 2. School of Health an Relate Research ScHARR, Unversty of Sheffel, UK. 2. School
More informationChapter 2 Transformations and Expectations. , and define f
Revew for the prevous lecture Defnton: support set of a ranom varable, the monotone functon; Theorem: How to obtan a cf, pf (or pmf) of functons of a ranom varable; Eamples: several eamples Chapter Transformatons
More informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
More informationHigh 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 informationA Note on the Numerical Solution for Fredholm Integral Equation of the Second Kind with Cauchy kernel
Journal of Mathematcs an Statstcs 7 (): 68-7, ISS 49-3644 Scence Publcatons ote on the umercal Soluton for Freholm Integral Equaton of the Secon Kn wth Cauchy kernel M. bulkaw,.m.. k Long an Z.K. Eshkuvatov
More informationNon-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 informationENTROPIC 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 informationMAE140 - Linear Circuits - Winter 16 Final, March 16, 2016
ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have
More informationCalculation of time complexity (3%)
Problem 1. (30%) Calculaton of tme complexty (3%) Gven n ctes, usng exhaust search to see every result takes O(n!). Calculaton of tme needed to solve the problem (2%) 40 ctes:40! dfferent tours 40 add
More informationNumerical modeling of a non-linear viscous flow in order to determine how parameters in constitutive relations influence the entropy production
Technsche Unverstät Berln Fakultät für Verkehrs- un Maschnensysteme, Insttut für Mechank Lehrstuhl für Kontnuumsmechank un Materaltheore, Prof. W.H. Müller Numercal moelng of a non-lnear vscous flow n
More informationC4B Machine Learning Answers II. = σ(z) (1 σ(z)) 1 1 e z. e z = σ(1 σ) (1 + e z )
C4B Machne Learnng Answers II.(a) Show that for the logstc sgmod functon dσ(z) dz = σ(z) ( σ(z)) A. Zsserman, Hlary Term 20 Start from the defnton of σ(z) Note that Then σ(z) = σ = dσ(z) dz = + e z e z
More informationChapter 9: Statistical Inference and the Relationship between Two Variables
Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,
More informationStudy 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