Adjoint Transient Sensitivity Analysis in Circuit Simulation
|
|
- Mitchell Riley
- 6 years ago
- Views:
Transcription
1 Ajoint Transient Sensitivity Analysis in Circuit Simulation Z. Ilievski 1, H. Xu 1, A. Verhoeven 1, E.J.W. ter Maten 1,2, W.H.A. Schilers 1,2 an R.M.M. Mattheij 1 1 Technische Universiteit Einhoven; Z.Ilievski@tue.nl 2 NXP Semiconuctors Einhoven Summary. Sensitivity analysis is an important tool that can be use to assess an improve the esign an accuracy of a moel escribing an electronic circuit. Given a moel escription in the form of a set of ifferential-algebraic equations it is possible to observe how a circuit s output reacts to varying input parameters, which are introuce at the requirements stage of esign. In this paper we consier the ajoint metho more closely. This metho is efficient when the number of parameters is large. We exten the transient sensitivity work of Petzol et al., in particular we take into account the parameter epenency of the ynamic term. We also compare the complexity of the irect an ajoint sensitivity an erive some error estimates. Finally we sketch out how Moel Orer Reuction techniques coul be use to improve the efficiency of ajoint sensitivity analysis. Keywors Sensitivity Analysis; Transient Analysis; Ajoint Metho; Moel Orer Reuction 1 Introuction A typical quantity in circuit analysis is the prouct of voltage ifference times the current through an electronic component (power) an, when integrate of time, this reflects the total power that is issipate. Another time omain problem is the etermination of the time moment when a certain unknown, or an expression, crosses a particular value. Such a moment can be the moment at which synchronization is require in co-simulation between a circuit simulator an another simulation tool. In transient analysis, the ajoint metho can be formulate as a convolution of the circuit equations with a carefully constructe function, that, by its nature, requires a backwar integration in time of a relate DAE (an for which a proper initial value has to be etermine). The metho has been popularize in 8, 12] for linear problems. For more general DAEs the metho has been stuie in 7] in a more mathematical way.
2 2 Z. Ilievski et al In 16] the application to the nonlinear DAEs of circuit equations was stuie more closely. Nice applications can be erive for the problem of fining optimal sources in etecting faults in analog circuits 5]. However, in stuying sizing problems (in which for instance the physical area of a capacitor has to be taken into account), it appears that especially parameters of capacitors give rise to terms that require aitional investigation. Here the effect of the inex of the relate DAE shows up. Apart from purposes of optimization, ajoint systems are of interest in etermining optimal reuce orer moels 3, 1], in which case a large number of parameters occurs. Because the ajoint systems are linear the equations themselves can be mae subject to a reuce orer moeling process. We will escribe ways how to calculate sensitivities in a stable an an efficient way. 2 Transient sensitivity analysis Equation (1) is a general Differential-Algebraic Equation (DAE) that can be use to escribe how any circuit behaves over a perio of time. In Moifie Noal Analysis 12], x(t) R N is the state vector an represents the noe voltages an the currents through voltage sources an inuctors, j an q are vector functions that escribe the current an charge (capacitors) or flux (inuctors) behavior. All source values are comprise in s(t) q(x(t))] + j(x(t)) = s(t). (1) The initial solution at t =, the DC-solution x DC, satisfies j(x DC ) = s(). (2) Applying Euler-Backwar time integration between time points t n an t n+1 = t n + t enables to calculate x n+1 as approximation at t n+1 : 1 t q(xn+1 ) q(x n )] + j(x n+1 ) s(t n+1 ) = (3) A Newton-Raphson proceure involves the coefficient matrix Y = 1 t C + G, in which C = q/ x an G = j/ x. Making explicit that the equations an its solution epen on a parameter p R P we will write q(x(t, p), p)] + j(x(t, p), p) = s(t, p). (4) By ajusting these parameters it is possible to optimize the behavior of a require functionality. The sensitivity of x(t, p) with respect to p is enote by
3 Ajoint Transient Sensitivity Analysis 3 ˆx(t, p) x(t, p)/ = ( x i (t, p)/ j ) R N P, an similarly for ˆx DC (p). After solving (3), an saving of the LU-ecomposition of the matrix Y = LU, the sensitivity ˆx n+1 (p) ˆx(t n+1, p) may be calculate by recursion 9, 11] ˆx n+1 (p) = Y 1 f, in which (5) f = 1 q n+1 q n] j n+1 + s n t t Cˆxn (p). (6) The vector f requires O(P N 2 ) operations for the last term in aition to O(P N) evaluations for a term like q etc... For simplicity we assume full matrices. Solving the system requires an aitional O(P N 2 ) operations. A more general basic observation function is enote by F(x(t, p), p) R F from which other observation functions can be obtaine, like From (7) we erive If F x G(x(p), p) = F(x(t, p), p). (7) T ( ) F F G(x(p), p) = ˆx +. (8) x can be etermine rather cheaply in (8), the main emphasis in sensitivity analysis is in the efficient calculation of ˆx, or even in efficiently calculating the inner-prouct F F F x ˆx. Note that, from (5), we erive x ˆx = x Y 1 f = Y T F x ]T ] T f, which can be calculate in O(min(F, P )N 2 +F P N) operations (in aition to those alreay mentione above: the overall leaing N 2 -term still has coefficient P ). This is a irect, forwar, analysis. When, aitionally, some library for evaluating q, j, or s, oes not allow symbolic ifferentiation, here also a symmetric finite ifference will be mae (at the cost of two aitional evaluations for each quantity q q(p+ p) q(p p) 2 p ). This means that at each interior time point of (8) the integran will have an error O( p 2 ) (assuming this iscretization error is ominant). A quarature rule like the Trapezoial Rule as up to O( p 2 / t) leaing to p = o( t) if t an no persistent errors in sensitivities are wante. In the sequel, we now consier an approach base on (backwar) ajoint integration 7]. We ifferentiate (4) w.r.t. p an multiply the result with a function λ (t) R F N (in which the means transpose), yieling = = = λ (t) λ (t) q λ (t) q ] T + ] (9) q + j s ( λ q j + λ s ( λ q x x + q ] T + )] ) + λ ( j x x + j s )]
4 4 Z. Ilievski et al = = λ (t) q λ (t) q ] T + ] T + λ ( λ ( Cˆx + q ) + λ ( Gˆx + j C λ G )ˆx λ q + λ s ( j This result hols for any λ. We now consier some choices. )] s )] (1) 2.1 Backwar, ajoint sensitivity for G(x(p), p) In (8) we encounter the prouct F x ˆx. Equation (1) will enable us to get ri of the ˆx, which oes not nee to be calculate explicitly. We choose λ(t) R N F appropriately an require that λ(t) satisfies the linear ajoint DAE (we assume the inex 1 case, which is not trivial 4]). C λ ( F ). G λ = (11) x This oes not yet make λ(t) unique, because we i not specify the initial value yet. But we are now able to express F λ x ˆx in terms of C G λ] ˆx, after which we can apply (1). T ( ) F F G(x(p), p) = ˆx + x ( ] = λ C + λ G ˆx + F ) = λ (t) q ] ( ( T λ q j t= + λ s ) + F ). The first term involves ˆx() = ˆx DC, which we know. However, if λ (T ) q x (T ), one still nees ˆx(T ), which we i not want to etermine explicitly. Hence goo choices to efine the initial conitions for λ(t) are C(T ) = q (T ) = λ(t ) =, or if (13) x C(t) = q ( F ) x (t) = λ(t ) λ DC, with G λ DC = (14) x t= Note that (13) is only allowe for DAEs of inex up to 1. With this choice (12) simplifies to (12) q G(x(p), p) = λ () x () ˆx DC + q ] () + ( λ q λ ( j s ) + F ). (15)
5 Ajoint Transient Sensitivity Analysis 5 The circuit that escribes the ajoint system (11) is escribe in more etails in 8, 12]. There the above expression is erive by conservation properties base on the Kirchhoff Laws which imply Tellegen s Theorem. Note that (9) is just the time integral of inner-proucts of branch-currents AjointCircuit * branch-voltages OriginalCircuit. The formulation use above follows 7], but reveals more closely the effect of a non-trivial function q in (12) that explicitly epens on p. By this, also the velocity λ will be neee. Alreay a linear function q(x, p) = C(p)x requires this. Because of the DAE-nature of the problem (1), in 16], λ was estimate by symmetric finite ifferences in the interior of the interval an with one-sie approximations at the bounaries. When, as before, a library for evaluating q, j, or s, oes not allow symbolic ifferentiation, here also a symmetric finite ifference will be mae (at the cost of two aitional evaluations for each quantity). This means that at each interior time point the integran will have an error O( t 2 + p 2 ). Any quarature rule (like the Trapezoial Rule) may a these errors up to an error O( t + p 2 / t) which is aitional to the error of the quarature rule, which means that one will require p = O( t) (note that a one-sie ifference will even nee p = O( t 2 )). It also shows that the Trapezoial Rule may not give better results than simple, first orer, Euler integration. System (11) can be etermine by integrating backwars in time after x(t, p) has been etermine in the nominal analysis. This backwar time integration 1 of (11) requires Jacobian matrices tc + G (assuming Euler backwars), similar as in the forwar analysis. When the same step sizes are use an the LU-ecompositions of the Jacobian matrices at the converge values have been save from the nominal analysis, the transpose ecompositions can be reuse for the Jacobians when integrating (11) (one may also save approximative inverse matrices, or preconitioning matrices). However we will assume that one will re-ecompose them uring the backwar integration. Note that in 6] it is remarke that if the same time step is use in the forwar (for x) an backwar analysis (for λ) this may give rise to very inaccurate solutions for λ. This general step size approach introuces effects ue to interpolation (effects which we have not yet stuie). Let W = O(N α ) represent the number of operations for the LU-ecompositions with 1 α 2 for sparse systems an α = 3 for full systems. Note that λ () q x () ˆx DC = C λ()] ˆx DC, which can be hanle similarly as in ajoint sensitivity analysis in the DC-problem, without explicitly calculating ˆx DC, in W + O(min(F, P )N 2 + F P N) operations (assuming that G has been ecompose again at t = ). Each time integration step of (11) requires W + O(F N 2 ) operations, after which the integran in (15) at each time point requires O(P N + F P ) evaluations an O(F P N) aitional operations. In practise F P, which (apart from the W term) makes (15) more efficient as the irect, forwar metho for (5)-(6). Reuction of the W term is iscusse in Section 2.2.
6 6 Z. Ilievski et al 2.2 MOR applie to the global ajoint sensitivity equations The main buren of the backwar ajoint sensitivity equations still is the W = O(N α ) work neee for the LU-ecompositions in the case α = 3 when integrating backwars in time for the ajoint equations for λ(t). In orer to reuce this we observe that in the interior we only nee to know λ(t) for coorinates where q, j s, an are non trivial. However, at t = also for the nontrivial rows of q x () coorinates of λ(t) shoul be known (which may significantly increase the number of neee coorinates). More precisely, in (15) a term like λ () q x () ˆx DC R F P shows that the total number of output is r = F P, an when r N, one may think to apply MOR. We observe that there is a tentative opportunity to apply Proper Orthogonal Decomposition (POD) 1, 13], since the forwar time integration to etermine x(t) elivere a nice series of snapshots {x(t ),..., x(t N )} (an, even cheaply, also of F x (t k)). With POD a matrix V R N r is foun such that x V x, x R r, r N, in which V is time inepenent an similar snapshots coul have been obtaine from a moifie problem VT q(v x(t))] + V T j(v x(t)) = V T s(t). (16) For (16) the matrices are of size r r, which inicates a nice option for MOR. In this paper we apply the V T, V matrices from the POD irectly to the ajoint problem (11). Note that in this approach we can neglect the epenency of V on p. For each p one integrates (1) an calculates the snapshots, resulting in a POD Moel Orer reuction projection matrix V. With V we can reuce the system of DAEs for λ (in which only the error ue to POD matters) resulting in λ P OD (t) λ(t) (please see 13] for error estimates). Note that POD can also be use to reuce (1) itself. In this case the epenency of V on p introuces an aitional error in the proceure that can not be neglecte. This last approach coul further reuce costs an has to be stuie further. In 2] it was shown that (in general) POD not irectly applies to DAEs. Here a Least-Square POD remey was introuce that can be applie to the linear DAE (11). Alternatively to POD, uring the forwar integration one coul aitionally etermine projection matrices for the Trajectory Piecewise-Linear Metho 14, 15]. Next, similarly as to the POD case, a global matrix V is etermine that allows for Piecewise-Linear MOR. In this case the reuce DAE problems can irectly be solve 15].
7 Ajoint Transient Sensitivity Analysis 7 3 Results p=l R (m) G/ G/ withp OD Table 1. Sensitivities, where p is the parameter, in this case the parameter is the length of the secon resistor. Units are (W tm 1 ) Fig. 1. Voltage ifferences for successive sensitivity values, noe 2 Equation (12) was implemente in Matlab for a simple circuit an the sensitivities of the energy G issipate for a resistor R = R(l R ) were observe while changing the length parameter l R G = I(R) V(R). (17) To see the epenency of V R for ifferent values of p, in Fig.1. we plotte, A: (V R (p =.21) V R (p =.2))(t) B: (V R (p =.22) V R (p =.21))(t). In Table 1 the sensitivity G/, calculate by the backwar ajoint metho, are shown for 3 ifferent values of p. To test POD we create a larger, equivalent system by splitting one resistor in to a number of smaller resistors. Doing this enable us to compare the sensitivity of our observe resistor in the reuce large system with the sensitivity in the original system. We generate
8 8 Z. Ilievski et al our projection matrices by applying the singular value ecomposition on the new state snapshot matrix for the enlarge problem. Inee the number of nontrivial singular values for both systems was the same. We applie the projection matrices to reuce only the backwar ajoint calculation steps. The size of the matrices was consierably reuce an the calculate sensitivities from the reuce system were exact to at least 3 significant figures, please see Table 1 for sensitivity results for the POD an for the non POD approach. 4 Conclusions The backwar, ajoint sensitivity methos are immeiately attractive when the original DAE (1) is linear an when the number of parameters P 1. Direct forwar an backwar ajoint approaches impose ifferent accuracy conitions to finite ifference approximations. The irect forwar metho exploits the re-use of LU-ecompositions. The backwar ajoint methos becomes more of interest when MOR can be applie, or when otherwise approximate LU-ecompositions coul have been save uring the forwar time integration. In these cases they can outperform the irect forwar metho when the number of parameters P is large (but still smaller than N, usually F P ). We have shown that applying POD MOR to the backwar ajoint step is possible an works very well. MOR techniques can also be use to reuce the effort in sensitivity calculation in the forwar analysis. In future we want to stuy more closely the application of our sensitivity calculations to larger circuits an with more inustrial relevance. We also want to consier the effect of ifferent time stepping in forwar an backwar analysis. An, finally, we want to stuy the sensitivity of the reuce DAE (16) in which V = V(p) epens on p as well. References 1. P. Astri: Reuction of process simulation moels: a proper orthogonal ecomposition approach, PhD-Thesis, Department of Electrical Engineering, Einhoven University of Technology, P. Astri, A. Verhoeven: Application of least squares MPE technique in the reuce orer moeling of electrical circuits, TU Einhoven, Center for Analysis, Scientific Computing an Applications, CASA Report 11, O. Balima, Y. Favennec, M. Girault, D. Petit: Comparison between the moal ientification metho an the POD-Galerkin metho for moel reuction in nonlinear iffusion systems, Int. J.Numer. Meth. Engng., Vol. 67, pp , K. Balla, R. März : Linear ifferential-algebraic equations of inex-1 an their ajoint equations, Results in Maths, Vol.37, pp.12-35, B. Buriek: Generation of optimal test stimuli for nonlinear analog circuits using nonlinear programming an time-omain sensitivities, Proc. DATE 21, Munchen, pp B. Buriek: Zur Berechnung von Testsignalen für nichtlineare analoge Schaltungen unter Verwenung von Methoen er Optimalsteurungstheorie, Shaker, Achen, 25 (PhD Thesis Univ. of Haannover, 25). 7. Y. Cao, S. Li, L. Petzol, R. Serban: Ajoint sensitivity for ifferential-algebraic equations: the ajoint DAE system an its numerical solution, SIAM J. Sci. Comput., Vol. 24-3, pp , 22.
9 Ajoint Transient Sensitivity Analysis 9 8. A.R. Conn, P.K. Coulman, R.A. Haring, G.L. Morrill, C. Visweswariah, C.W. Wu: JiffyTune: circuit optimization using time-omain sensitivities, IEEE Trans. on CAD of ICs an Systems, Vol , pp , L. Daloss, P. Gubian, M. Quarantelli: Multiparameter time-omain sensitivity computation, IEEE Trans. on Circuits an Systems - I: Fun. Theory an Applics, Vol , pp , Y. Favennec, M. Girault, D. Petit: The ajoint metho couple with the moal ientification metho for nonlinear moel reuction, Inverse Probl. in Science an Engng., Vol. 14, No. 3, , D.E. Hocevar, P. Yang, T.N. Trick, B.D. Epler: Transient sensitivity computation for MOS- FET circuits, IEEE Trans. on CAD of Integr. Circuits an Systems, Vol. CAD-4, Nr. 4, pp , L.T. Pillage, R.A. Rohrer, C. Visweswariah: Electronic circuit an system simulation methos, McGraw-Hill, Inc, New York, USA, ISBN , M. Rathinam, L.R. Petzol: A new look at proper orthogonal ecomposition, SIAM J. Numer. Analysis, Vol. 41-5, pp , M.J. Rewienski: A trajectory piecewise-linear approach to moel orer reuction of nonlinear ynamical systems, PhD-Thesis, Massachutes Institute of Technology, T. Voss: Moel reuction for nonlinear ifferential algebraic equations, MSc. Thesis University of Wuppertal, 25; Unclassifie Report PR-TN-25/919, Philips Research Laboratories, H. Xu: Transient Sensitivity Analysis in Circuit Simulation, MSc-Thesis, Department of Mathematics an Computing Science, Einhoven University of Technology, 24. Online via author search at Library at
Switching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationTransient Sensitivity Analysis CASA Day 13th Nov 2007 Zoran Ilievski. Zoran Ilievski Transient Sensitivity Analysis
CASA Day 13th Nov 2007 Talk Structure Talk Structure Introduction Talk Structure Introduction Recap Sensitivity Talk Structure Introduction Recap Sensitivity Examples and Results Talk Structure Introduction
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationWe G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies
We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationSome vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10
Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the
More informationECE 422 Power System Operations & Planning 7 Transient Stability
ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationDesigning Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations
EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationOptimization of a point-mass walking model using direct collocation and sequential quadratic programming
Optimization of a point-mass walking moel using irect collocation an sequential quaratic programming Chris Dembia June 5, 5 Telescoping actuator y Stance leg Point-mass boy m (x,y) Swing leg x Leg uring
More informationDesigning Information Devices and Systems I Spring 2018 Lecture Notes Note 16
EECS 16A Designing Information Devices an Systems I Spring 218 Lecture Notes Note 16 16.1 Touchscreen Revisite We ve seen how a resistive touchscreen works by using the concept of voltage iviers. Essentially,
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationDomain Decomposition and Model Reduction of Systems with Local Nonlinearities
Domain Decomposition an Moel Reuction of Systems with Local Nonlinearities Kai Sun, Rolan Glowinsi, Matthias Heinenschloss, an Danny C. Sorensen Abstract The goal of this paper is to combine balance truncation
More informationModeling time-varying storage components in PSpice
Moeling time-varying storage components in PSpice Dalibor Biolek, Zenek Kolka, Viera Biolkova Dept. of EE, FMT, University of Defence Brno, Czech Republic Dept. of Microelectronics/Raioelectronics, FEEC,
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationEnergy behaviour of the Boris method for charged-particle dynamics
Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationThermal runaway during blocking
Thermal runaway uring blocking CES_stable CES ICES_stable ICES k 6.5 ma 13 6. 12 5.5 11 5. 1 4.5 9 4. 8 3.5 7 3. 6 2.5 5 2. 4 1.5 3 1. 2.5 1. 6 12 18 24 3 36 s Thermal runaway uring blocking Application
More informationTIME-DELAY ESTIMATION USING FARROW-BASED FRACTIONAL-DELAY FIR FILTERS: FILTER APPROXIMATION VS. ESTIMATION ERRORS
TIME-DEAY ESTIMATION USING FARROW-BASED FRACTIONA-DEAY FIR FITERS: FITER APPROXIMATION VS. ESTIMATION ERRORS Mattias Olsson, Håkan Johansson, an Per öwenborg Div. of Electronic Systems, Dept. of Electrical
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationA. Incorrect! The letter t does not appear in the expression of the given integral
AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationA Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation
A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an
More informationA Path Planning Method Using Cubic Spiral with Curvature Constraint
A Path Planning Metho Using Cubic Spiral with Curvature Constraint Tzu-Chen Liang an Jing-Sin Liu Institute of Information Science 0, Acaemia Sinica, Nankang, Taipei 5, Taiwan, R.O.C., Email: hartree@iis.sinica.eu.tw
More informationMath 1B, lecture 8: Integration by parts
Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationUNIFYING PCA AND MULTISCALE APPROACHES TO FAULT DETECTION AND ISOLATION
UNIFYING AND MULISCALE APPROACHES O FAUL DEECION AND ISOLAION Seongkyu Yoon an John F. MacGregor Dept. Chemical Engineering, McMaster University, Hamilton Ontario Canaa L8S 4L7 yoons@mcmaster.ca macgreg@mcmaster.ca
More information23 Implicit differentiation
23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For
More informationIntroduction to variational calculus: Lecture notes 1
October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationTutorial Test 5 2D welding robot
Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.
More informationAPPPHYS 217 Thursday 8 April 2010
APPPHYS 7 Thursay 8 April A&M example 6: The ouble integrator Consier the motion of a point particle in D with the applie force as a control input This is simply Newton s equation F ma with F u : t q q
More informationAn inductance lookup table application for analysis of reluctance stepper motor model
ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR
More informationExperimental Robustness Study of a Second-Order Sliding Mode Controller
Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationSYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS
SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS GEORGE A HAGEDORN AND CAROLINE LASSER Abstract We investigate the iterate Kronecker prouct of a square matrix with itself an prove an invariance
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationSpace-time Linear Dispersion Using Coordinate Interleaving
Space-time Linear Dispersion Using Coorinate Interleaving Jinsong Wu an Steven D Blostein Department of Electrical an Computer Engineering Queen s University, Kingston, Ontario, Canaa, K7L3N6 Email: wujs@ieeeorg
More informationIPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy
IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationDAE-Index and Convergence Analysis of Lumped Electric Circuits Refined by 3-D Magnetoquasistatic Conductor Models
Bergische Universität Wuppertal Fachbereich Mathematik un Naturwissenschaften Lehrstuhl für Angewante Mathematik un Numerische Mathematik Preprint BUW-AMNA 8/6 Sebastian Schöps DAE-Inex an Convergence
More informationPort-Hamiltonian systems: an introductory survey
Port-Hamiltonian systems: an introuctory survey Arjan van er Schaft Abstract. The theory of port-hamiltonian systems provies a framework for the geometric escription of network moels of physical systems.
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationThe continuity equation
Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal
More informationCOUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND
More informationSystems & Control Letters
Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More informationYear 11 Matrices Semester 2. Yuk
Year 11 Matrices Semester 2 Chapter 5A input/output Yuk 1 Chapter 5B Gaussian Elimination an Systems of Linear Equations This is an extension of solving simultaneous equations. What oes a System of Linear
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationPredictive Control of a Laboratory Time Delay Process Experiment
Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationIERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210
IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc
More informationOrdinary Differential Equations: Homework 1
Orinary Differential Equations: Homework 1 M. Gameiro, J.-P. Lessar, J.D. Mireles James, K. Mischaikow January 12, 2017 2 Chapter 1 Motivation 1.1 Exercises Exercise 1.1.1. (Frictionless spring) Consier
More informationPerturbation Analysis and Optimization of Stochastic Flow Networks
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. Y, MMM 2004 1 Perturbation Analysis an Optimization of Stochastic Flow Networks Gang Sun, Christos G. Cassanras, Yorai Wari, Christos G. Panayiotou,
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationConsider for simplicity a 3rd-order IIR filter with a transfer function. where
Basic IIR Digital Filter The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient ifference equation
More informationTractability results for weighted Banach spaces of smooth functions
Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March
More informationBasic IIR Digital Filter Structures
Basic IIR Digital Filter Structures The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationReal-time economic optimization for a fermentation process using Model Predictive Control
Downloae from orbit.tu.k on: Nov 5, 218 Real-time economic optimization for a fermentation process using Moel Preictive Control Petersen, Lars Norbert; Jørgensen, John Bagterp Publishe in: Proceeings of
More informationTRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS
TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email
More informationLinear Algebra- Review And Beyond. Lecture 3
Linear Algebra- Review An Beyon Lecture 3 This lecture gives a wie range of materials relate to matrix. Matrix is the core of linear algebra, an it s useful in many other fiels. 1 Matrix Matrix is the
More informationDAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1
11 th National Congress on Theoretical an Applie Mechanics, 2-5 Sept. 2009, Borovets, Bulgaria DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1 E. MANOACH Institute of Mechanics,
More informationA new approach to explicit MPC using self-optimizing control
28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June -3, 28 WeA3.2 A new approach to explicit MPC using self-optimizing control Henrik Manum, Sriharakumar Narasimhan an Sigur
More informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
More informationProblem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs
Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable
More informationMonopoly Part III: Multiple Local Equilibria and Adaptive Search
FH-Kiel University of Applie Sciences Prof Dr Anreas Thiemer, 00 e-mail: anreasthiemer@fh-kiele Monopoly Part III: Multiple Local Equilibria an Aaptive Search Summary: The information about market eman
More informationOptimal LQR Control of Structures using Linear Modal Model
Optimal LQR Control of Structures using Linear Moal Moel I. Halperin,2, G. Agranovich an Y. Ribakov 2 Department of Electrical an Electronics Engineering 2 Department of Civil Engineering Faculty of Engineering,
More informationThe Three-dimensional Schödinger Equation
The Three-imensional Schöinger Equation R. L. Herman November 7, 016 Schröinger Equation in Spherical Coorinates We seek to solve the Schröinger equation with spherical symmetry using the metho of separation
More informationNON-SMOOTH DYNAMICS USING DIFFERENTIAL-ALGEBRAIC EQUATIONS PERSPECTIVE: MODELING AND NUMERICAL SOLUTIONS. A Thesis PRIYANKA GOTIKA
NON-SMOOTH DYNAMICS USING DIFFERENTIAL-ALGEBRAIC EQUATIONS PERSPECTIVE: MODELING AND NUMERICAL SOLUTIONS A Thesis by PRIYANKA GOTIKA Submitte to the Office of Grauate Stuies of Texas A&M University in
More informationarxiv:hep-th/ v1 3 Feb 1993
NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,
More informationProof of SPNs as Mixture of Trees
A Proof of SPNs as Mixture of Trees Theorem 1. If T is an inuce SPN from a complete an ecomposable SPN S, then T is a tree that is complete an ecomposable. Proof. Argue by contraiction that T is not a
More information