Chapter 11: Feedback and PID Control Theory

Size: px
Start display at page:

Download "Chapter 11: Feedback and PID Control Theory"

Transcription

1 Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by measuring the current state of a physical system, etermining how far the current state is from the esire state, an then automatically applying a control signal to bring the system closer to the esire state. This process is repeate iteratively to bring the system to the esire state an keep it there. Feeback can be use very effectively to stabilize the state of a system, while also improving its performance: Engineers use feeback to control otherwise unstable esigns; op-amps use feeback to stabilize an linearize their gain; an physicists use feeback to stabilize an improve the performance of their instruments. A. Feeback in engineering Feeback is ubiquitous in engineering. ts application has le to evice features an machines which woul not otherwise function. Here are few examples: Climate control: A sensor measures the temperature an humiity in a room an then heats or cools an humiifies or ehumiifies accoringly. Automobile cruise control: The car measures its spee an then applies the accelerator or not epening on whether the spee must be increase or ecrease to maintain the target spee. Highly maneuverable fighter jets: The F-16 Falcon fighter jet is an inherently unstable aircraft (i.e. the airframe will not glie on its own). The F-16 oes fly because 5 onboar computers constantly measure the aircraft s flight characteristics an then apply corrections to the control surfaces (i.e. ruer, flaps, ailerons, etc ) to keep it from tumbling out of control. The avantage of this technique is that the aircraft has the very rapi response an maneuverability of a naturally unstable airframe, while also being able to fly. B. Feeback in electronics: Op-amps use feeback to achieve very high linearity an preictability for their close-loop gain by sacrificing some of their extremely high open-loop gain. Another common application of feeback in electronics is in precision, fast- response power supplies. Constant current an constant voltage power supplies which have a high egree of stability use feeback to regulate their current or their voltage, by measuring the current an voltage across a precision shunt resistors an then using feeback to automatically correct for any eviations from the esire output. Feeback also allows the power supply to ajust its voltage or current very quickly an controllably in response to a change in loa

2 Chapter 11: Feeback an D Control Theory C. Feeback in physics Feeback has become a familiar tool for experimental physicists to improve the stability of their instruments. n particular, physicists use feeback for precise control of temperature, for stabilizing an cooling particle beams in accelerators, for improving the performance of atomic force microscopes, for locking the optical frequency of lasers to atomic transitions, an referencing quartz oscillators to groun state atomic hyperfine microwave transitions in atomic clocks, to name just a few example. Temperature control: Many elicate physics evices, such as crystals, lasers, RF oscillators, an amplifiers, require their temperature to be very stable in orer to guarantee their performance. For example, the wavelength of ioe lasers generally has a temperature epenence on the orer of. nm/ C, but requires a stability of 1-6 nm for experiments. tochastic cooling: n a particle accelerator, the transverse momentum sprea of particles must be reuce to a minimum. The reuce momentum sprea increases the particle ensity, or beam luminosity, an consequently the probability of collisions with a similar counter-propagating particle beam in the etector area. tochastic cooling works by measuring the transverse positions an momenta of the particles as they pass through a section of the accelerator, an then applying appropriate momentum kicks to some of the particles at other points in the accelerator ring to reuce the overall transverse momentum sprea. The process is repeate until the momentum sprea is sufficiently reuce. The 1984 Nobel rize in hysics was aware in part to imon van er Meer for his invention of stochastic cooling which contribute to the iscovery of the W an Z bosons (weak force meiators) at CERN. Atomic force microscope: An atomic force microscope uses a very sharp tip (just a few nanometers in size at the very tip) which is scanne back an forth just a few nanometers above the surface to be image. nstea of scanning the tip at a constant height above the surface, which coul lea to the tip actually running into a bump on the surface, the microscope uses feeback to ajust the tip height such that the force (from the surface atoms) on the tip is constant. Laser locking: Many experiments in atomic an optical physics require lasers which have a very stable optical frequency. The optical frequency of the laser is locke by measuring the optical frequency ifference between the laser an an atomic transition an using feeback to set this ifference to a constant value. Lasers can be routinely stabilize with feeback to better than 1 MHz out of 3x1 14 Hz (about 1 part per billion), though stabilities close to 1 Hz have been reporte after heroic efforts. Atomic clocks: n an atomic clock, the frequency of an RF oscillator (a quartz crystal for example) is compare to that of a groun state atomic hyperfine microwave transition (6.8 or 9. GHz). The frequency ifference is measure an the frequency of the RF oscillator is correcte by feeback. The process is constantly repeate to eliminate any rift in the frequency of the RF oscillator. Atomic fountain clocks can achieve accuracies in the range of 1 part in 1 15, an plans are unerway to construct optical atomic clocks with accuracies an stabilities of about 1 part in

3 Chapter 11: Feeback an D Control Theory. Feeback n this section we introuce the main elements of a generic feeback moel. A. ystem Consier a simple system characterize by a single variable. Uner normal conitions the system has a steay state value of = which may vary an rift somewhat over time ue to the variation of environmental variables v which we cannot measure or are unaware of. We possess a mechanism for measuring the state of the system as well as a control input u with which we can use to moify the state of the system. n summary, the system has the following functional form (u; v;. We will make the final assumption that is monotonic with u in the vicinity of (i.e. that the plot of vs. u oes not have any maxima or minima, an that /u is either always positive or always negative). Figure 1 shows a conceptual schematic of the relationship between the system, the variables u an v, an the measurement of the system state. unknowns v moify ystem tate: = Measurement of Control u Moifies Figure 1: Conceptual schematic of system B. Objective Our objective is to set or lock the state of the system to a esire value = an keep it there without letting it rift or vary over time, regarless of variations in the environmental variables v. C. Feeback moel We will set or lock the state of the system to = with the following proceure (see also figure ): 1. Measure the state of the system

4 Chapter 11: Feeback an D Control Theory. Determine how far the system is from its esire set point by efining an error variable, e=-. 3. Calculate a trial control value u=u(e). 4. Fee the calculate control value, u(e), back into the control input of the system. 5. The state of the system changes in response to the change in the control value. 6. Return to step. f we repeat this feeback cycle inefinitely with an appropriately calculate control value u(e), then the system will converge to the state = an remain there even uner the influence of small changes to other variables (i.e. v) which influence the value of the state. This feeback moel can be aapte to inclue several state variables an several feeback variables. unknowns v moify ystem tate: Measurement of Calculate error e=- Control u Moifies Calculate u=u(e) Figure : Conceptual schematic of system with the feeback loop. n section, we iscuss a frequently use expression for calculating the feeback control variable u(e).. D Feeback Control The most popular type of feeback stabilization control, u(e), is roportional-ntegral- Derivative (D) gain feeback. D is very effective an easy to implement. The expression for u(e) epens only on the error signal e=- an is given by u( e; = g e( + g t e( t + g D t e( (1) - 4 -

5 Chapter 11: Feeback an D Control Theory where g, g, an g D are respectively the proportional, integral, an erivative gains. We also note that g, g, an g D o not have the same units. We will assume for simplicity that g is imensionless in which case u(e) has the same units as. A. Time evolution of the system with D feeback control We are now in a position to calculate the time evolution of the system uner the influence of feeback. Without feeback, the system woul remain in the state : ( no feeback = () may vary in time, but we will ignore this effect until part.c. n the presence of feeback, the state of the system at time t+δt (step 5) epens on the state of the system without feeback,, which has been moifie by the control input variable u(e). We now make the following simplifying assumption that the control input variable, u(e), controls or moifies the state of the system through the process of aition. n this case, the system state variable evolves accoring to the following equation: ( t Δ = + u( e; + (3) We can convert this equation to an integro-ifferential equation, if we assume that the system has a characteristic reponse time τ (small). n this case, equation 3 becomes ( + τ ( = + ge( + g e( t + gd e( (4) t t t B. pecial case: pure proportional gain feeback As a limiting case we consier pure proportional gain feeback (g = an g D =). We stuy this special case, because it is the basis for op-amp feeback an is also the simplest form of feeback. For g = an g D =, equation 4 becomes ( ( = + g e( t +τ (5) We can solve this 1 st orer ifferential equation, for the initial conition (t=)=, with the same technique we use in chapter 3 (equations 17-1). After a little bit of integration an algebra, which is left as an exercise to the reaer, we fin the following solution: 1 g g t τ ( = e + g 1 g 1 g (6) - 5 -

6 Chapter 11: Feeback an D Control Theory Equation 6 shows that the system will converge to the state =( -g )/(1-g ) when feeback control is applie, so long as the exponential exponent is negative (i.e. g <1), otherwise will iverge. We note that g < correspons to negative feeback. Figure 3 shows the response of a system for a imensionless gain of g =-1 an state values =.5 an =1, with time measure in units of τ (the system characteristic response time). ( Figure 3: ystem response with pure proportional gain feeback control an parameters g =-1, =.5, an =1. Time is measure in units of τ (the system characteristic response time). As figure 3 makes clear, the system oes not converge to the esire state =, though it oes reach its final steay state value relatively quickly. f we restrict ourselves to negative feeback, then accoring to equation 1, the system will converge to the steay state value ss of ss g = 1 g Equation 7 inicates that the system can be mae to converge to a steay state value ss which is arbitrarily close to = just by increasing the gain. n fact, for infinite proportional gain (i.e g - ) the system oes converge to ss = : This is the limit in which op-amps feeback operates. A note of caution: On its own, equation 7 is a little misleaing since it woul seem to imply that large positive feeback, g -, woul also prouce ss =. Of course, this is (7) - 6 -

7 Chapter 11: Feeback an D Control Theory not true since accoring to equation 6, the system will never achieve a steay state, but instea will iverge forever. C. olution for feeback control A large majority of D feeback controllers are actually just controllers (i.e. proportional an integral gain, but no erivative gain), an so for simplicity we solve equation 4 without the erivative gain term (g D =). The inclusion of the erivative gain term is conceptually simple an follows the same treatment as feeback an is left as an exercise to reaer. Derivative gain is use to improve the time response of the feeback, so that the system converges more quickly to its steay state value. With the erivative gain term omitte, equation (4) becomes ( + τ ( = + g e( + g e( t (8) t We can convert this integro-ifferential equation to a n orer linear ifferential equation with constant coefficients by taking the time erivative of equation (8) to obtain t τ = g + g g, (9) t t + where we employe the substitution e(=-. After combining similar terms, equation (9) becomes τ + (1 g ) g = g (1) t t Equation 1 is an inhomogeneous n orer ifferential equation with constant coefficients. The full solution to equation (7) is given by + λ+ t λ t ( = A+ e + A e (11) ( g 1) ± ( g 1) + 4g τ λ = τ ± (11a) The first two terms of equation 11 represent the homogeneous solution to equation 1, while the 3 r term is the inhomogeneous solution to the equation (it oes not epen on the initial contions). A + an A - are constants to be etermine from the initial conitions. Equation 11 shows that the system will converge to the state = when feeback control is applie, so long as λ + an λ - are negative (i.e. negative feeback), otherwise will iverge (exactly the opposite of what we want to accomplish with feeback). t - 7 -

8 Chapter 11: Feeback an D Control Theory f we choose (t=)= an (t=)/t= as our initial conitions, we can calculate the constants A + an A -. After a little bit of algebra, we fin that A + λ = λ λ+ ( ) (1a) A λ + = λ λ+ ( ) (1b) n figure 4, the behavior of the system uner feeback control is plotte for several ifferent parameters configurations. ( ( Figure 4: Time-evolution of a generic system with control feeback for =.5 an =1. For the left han plot the gain parameters are g =-1, g =-3, while for the right han plot the parameters are g =-1, g =-4. The small overshoot in the left han plot is ue to a small imaginary part in the exponential exponent of equation 11(a). The primary purpose of integral gain is to provie essentially infinite gain at DC ( Hz), which guarantees that ss =, as can be seen in figure 4. Figure 4 also shows that the larger the gain, the faster the correction time of the feeback control loop. D. Fourier space analysis of noise suppression One of the primary objectives of feeback is to make the system insensitive to noise on the system state, so that the system state stays locke to = regarless of external influences. n the absence of corrective feeback, external noise will cause the system state to eviate from =. External noise at a frequency ω will cause the system state to oscillate aroun its natural steay state such that = + N cos(ω, where N is the amplitue of the oscillations. Following the stanar Fourier space recipe of chapter 3, we replace cos(ω with exp(iω, an then take the real part at the en of our calculations. n essence, we must re-solve equations 8 with the following moification: - 8 -

9 Chapter 11: Feeback an D Control Theory e iωt + N (13) Using the substitution of equation 13, equation 1 becomes τ + (1 g ) t t g = iω N e iωt g Equation 14 has the same homogeneous solution as equation 1, but the inhomogeneous solution, ih (, iffers an is given by the following expression iω t = e + (14) iωt ih ( ) N (15) (1 g ) iω τω g We see that the noise term is present in the inhomogeneous solution, but with an aitional factor moifying the amplitue of the noise. n the case of negative feeback the moulus of this suppression factor, which we will call A N, is always less than unity an is given by the following expression A N = ω ( τω + g ) + (1 g ) ω (16) The plot in figure 5 shows the epenence of the suppression factor, A N, on frequency for ifferent feeback schemes. The plot shows that a combination of proportional an integral control gives the best suppression of noise, except in the vicinity of the resonant frequency ω = g / τ. The high frequency rop-off of the suppression factor is not ue to feeback but simply the natural response time τ of the system which also suppresses noise. ω, frequency (ras/τ) A N, suppression factor Figure 5: Comparison of the suppression factor, A N, for ifferent feeback schemes. The feeback control loop parameters are g =-1 an g =

10 Chapter 11: Feeback an D Control Theory V. Reality n practice, feeback is not quite as straightforwar as presente in the previous section. A. Gain vs, Frequency n the theoretical treatment of part, we assume that the proportional gain was inepenent of frequency. n practice, gain will generally fall off at higher frequencies ue to natural low-pass RC filtering in an amplifier an the larger circuit. As an example, figure 6 shows a plot of the open-loop gain of an op-amp as function of frequency, which has a clear rop-off in gain at higher frequencies. Figure 6: Open-loop gain of the O7 op-amp (Analog Devices O7 atasheet revision F, p. 1 (6)). B. hase shifts an positive feeback The natural or stray RC filtering of an amplifier not only rolls off the gain at high frequencies, but also introuces a -π/ phase shift. f the feeback loop has a secon stray unintentional RC filter present (for example, the natural time response of the system), then a secon -π/ phase shift is introuce. f the feeback gain is larger than 1 at the frequency at which the total accumulate phase is -π, then the feeback loop goes into positive feeback which causes the state of the system to iverge or sometimes oscillate out of control

11 Chapter 11: Feeback an D Control Theory C. tray RC positive feeback compensation One way to avoi having the system go into positive feeback is to purposely introuce an aitional RC low-pass filter into the feeback loop. f this RC filter has a f 3B frequency which is sufficiently smaller than the frequency at which the positive feeback occurs then the attenuation of the filter can bring the gain below 1 when the -π phase shift occurs. This way the feeback loop will no longer go into positive feeback above a certain frequency (of course there will not be any noise suppression or feeback action above this frequency either). Design Exercises: Design exercise 11-1: Consier an LED facing a photoioe in a manner similar to what you i last week in lab 1. Design a circuit which will maintain a constant optical power incient on the photoioe, even in the presence of external fluctuations in the room lighting. You shoul use the op-amp circuit of chapter 1, figure 11 an use feeback control to stabilize the intensity of the LED. Your circuit shoul be able to provie at least 1 ma at ~ V to the LED. Design exercise 11-: The non-inverting op-amp amplifier with finite open-loop gain. n this exercise you will NOT use the op-amp golen rules to solve the problem, unless explicitly inicate. Consier the non-inverting op-amp amplifier in the circuit below. V N + V OUT R 1 R n the following parts, V + an V - refer respectively to the voltages at the non-inverting an inverting terminals of the op-amp. The open-loop gain of the op-amp is A. a. Write own the funamental op-amp relation between V +, V -, V out, an A when no feeback is present (i.e. when R 1 an R are not presen. b. Assuming that the V - input raws no current, erive an expression for V - in terms of V out, R 1, an R

12 Chapter 11: Feeback an D Control Theory c. Obtain an expression for V out in terms of V in, A, R 1, an R, an etermine the gain G of the amplifier. Calculate V out in the limit of A +. Calculate V out in the limit A - an comment on its physical meaning.. uppose that the open-loop gain, A, is not very constant with frequency an changes by ΔA between frequencies f 1 an f. Derive an expression for the resulting relative variation in the amplifier gain ΔG/G in terms of the relative variation in the openloop gain ΔA/A. Calculate ΔG/G an ΔA/A for A=1 6, ΔA=1 5, R =1 kω, an R 1 =1 kω. e. Most op-amps feature a significant rop-off in their open-loop gain at frequencies above ~1 Hz. The rop-off follows a well establishe curve given by Af=constant, where the constant is calle the gain-banwith prouct (f is frequency in Hz). The gainbanwith prouct of the O-7 is 8 MHz. On a same log-log graph, plot the open-loop gain A vs frequency an the close loop gain G vs. frequency (for R =1 kω an R 1 =1 kω) for an O-7-base non-inverting amplifier

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by

More information

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an ID Control Theory Chapter 11: Feeback an ID Control Theory I. Introuction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works

More information

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

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

EE 370L Controls Laboratory. Laboratory Exercise #7 Root Locus. Department of Electrical and Computer Engineering University of Nevada, at Las Vegas

EE 370L Controls Laboratory. Laboratory Exercise #7 Root Locus. Department of Electrical and Computer Engineering University of Nevada, at Las Vegas EE 370L Controls Laboratory Laboratory Exercise #7 Root Locus Department of Electrical an Computer Engineering University of Nevaa, at Las Vegas 1. Learning Objectives To emonstrate the concept of error

More information

6.003 Homework #7 Solutions

6.003 Homework #7 Solutions 6.003 Homework #7 Solutions Problems. Secon-orer systems The impulse response of a secon-orer CT system has the form h(t) = e σt cos(ω t + φ)u(t) where the parameters σ, ω, an φ are relate to the parameters

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

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

Lecture 6: Control of Three-Phase Inverters

Lecture 6: Control of Three-Phase Inverters Yoash Levron The Anrew an Erna Viterbi Faculty of Electrical Engineering, Technion Israel Institute of Technology, Haifa 323, Israel yoashl@ee.technion.ac.il Juri Belikov Department of Computer Systems,

More information

Module 2. DC Circuit. Version 2 EE IIT, Kharagpur

Module 2. DC Circuit. Version 2 EE IIT, Kharagpur Moule 2 DC Circuit Lesson 9 Analysis of c resistive network in presence of one non-linear element Objectives To unerstan the volt (V ) ampere ( A ) characteristics of linear an nonlinear elements. Concept

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

APPPHYS 217 Thursday 8 April 2010

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

SYNCHRONOUS SEQUENTIAL CIRCUITS

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

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

Linear First-Order Equations

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

State observers and recursive filters in classical feedback control theory

State observers and recursive filters in classical feedback control theory State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent

More information

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables* 51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*

More information

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

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

Table of Common Derivatives By David Abraham

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

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

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

State-Space Model for a Multi-Machine System

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

ECE 422 Power System Operations & Planning 7 Transient Stability

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

CURRENT ELECTRICITY Q.1

CURRENT ELECTRICITY Q.1 CUENT EECTCTY Q. Define Electric current an its unit.. Electric Current t can be efine as the time rate of flow of charge in a conuctor is calle Electric Current. The amount of flow of charge Q per unit

More information

PD Controller for Car-Following Models Based on Real Data

PD Controller for Car-Following Models Based on Real Data PD Controller for Car-Following Moels Base on Real Data Xiaopeng Fang, Hung A. Pham an Minoru Kobayashi Department of Mechanical Engineering Iowa State University, Ames, IA 5 Hona R&D The car following

More information

PES 1120 Spring 2014, Spendier Lecture 36/Page 1

PES 1120 Spring 2014, Spendier Lecture 36/Page 1 PES 0 Spring 04, Spenier ecture 36/Page Toay: chapter 3 - R circuits: Dampe Oscillation - Driven series R circuit - HW 9 ue Wenesay - FQs Wenesay ast time you stuie the circuit (no resistance) The total

More information

On Using Unstable Electrohydraulic Valves for Control

On Using Unstable Electrohydraulic Valves for Control Kailash Krishnaswamy Perry Y. Li Department of Mechanical Engineering, University of Minnesota, 111 Church St. SE, Minneapolis, MN 55455 e-mail: kk,pli @me.umn.eu On Using Unstable Electrohyraulic Valves

More information

Determine Power Transfer Limits of An SMIB System through Linear System Analysis with Nonlinear Simulation Validation

Determine Power Transfer Limits of An SMIB System through Linear System Analysis with Nonlinear Simulation Validation Determine Power Transfer Limits of An SMIB System through Linear System Analysis with Nonlinear Simulation Valiation Yin Li, Stuent Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract This paper extens

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Yin Li, Stuent Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract This paper proposes a robust power system stabilizer (PSS)

More information

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique J. Basic. Appl. Sci. Res., ()9-99,, TextRoa Publication ISSN 9-434 Journal of Basic an Applie Scientific Research www.textroa.com A Comparison between a Conventional Power System Stabilizer (PSS) an Novel

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

ELEC3114 Control Systems 1

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

Chapter 6. Electromagnetic Oscillations and Alternating Current

Chapter 6. Electromagnetic Oscillations and Alternating Current hapter 6 Electromagnetic Oscillations an Alternating urrent hapter 6: Electromagnetic Oscillations an Alternating urrent (hapter 31, 3 in textbook) 6.1. Oscillations 6.. The Electrical Mechanical Analogy

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

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

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: 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 information

Applications of First Order Equations

Applications of First Order Equations Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

Alpha Particle scattering

Alpha Particle scattering Introuction Alpha Particle scattering Revise Jan. 11, 014 In this lab you will stuy the interaction of α-particles ( 4 He) with matter, in particular energy loss an elastic scattering from a gol target

More information

Separation of Variables

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

Experimental Robustness Study of a Second-Order Sliding Mode Controller

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

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR 5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR Tatsunori Kaneko* an Yasuto Asano* * Department of Mechanical Engineering,

More information

Schrödinger s equation.

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

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations

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

Differentiability, Computing Derivatives, Trig Review. Goals:

Differentiability, Computing Derivatives, Trig Review. Goals: Secants vs. Derivatives - Unit #3 : Goals: Differentiability, Computing Derivatives, Trig Review Determine when a function is ifferentiable at a point Relate the erivative graph to the the graph of an

More information

ELECTRONIC SYSTEMS. Real operational amplifiers. Electronic Systems - C2 28/04/ DDC Storey 1

ELECTRONIC SYSTEMS. Real operational amplifiers. Electronic Systems - C2 28/04/ DDC Storey 1 Politecnico i Torino ICT school Lesson C2 Real Op mp parameters an moel ELECTRONIC SYSTEMS C OPERTIONL MPLIFIERS C.2 Real Op mp an examples» Real Op mp parameters an moel» Static an ynamic parameters»

More information

involve: 1. Treatment of a decaying particle. 2. Superposition of states with different masses.

involve: 1. Treatment of a decaying particle. 2. Superposition of states with different masses. Physics 195a Course Notes The K 0 : An Interesting Example of a Two-State System 021029 F. Porter 1 Introuction An example of a two-state system is consiere. involve: 1. Treatment of a ecaying particle.

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

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

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask 5.4 FUNDAMENTAL THEOREM OF CALCULUS Do you remember the Funamental Theorem of Algebra? Just thought I' ask The Funamental Theorem of Calculus has two parts. These two parts tie together the concept of

More information

CONTROL CHARTS FOR VARIABLES

CONTROL CHARTS FOR VARIABLES UNIT CONTOL CHATS FO VAIABLES Structure.1 Introuction Objectives. Control Chart Technique.3 Control Charts for Variables.4 Control Chart for Mean(-Chart).5 ange Chart (-Chart).6 Stanar Deviation Chart

More information

Quantum optics of a Bose-Einstein condensate coupled to a quantized light field

Quantum optics of a Bose-Einstein condensate coupled to a quantized light field PHYSICAL REVIEW A VOLUME 60, NUMBER 2 AUGUST 1999 Quantum optics of a Bose-Einstein conensate couple to a quantize light fiel M. G. Moore, O. Zobay, an P. Meystre Optical Sciences Center an Department

More information

PARALLEL-PLATE CAPACITATOR

PARALLEL-PLATE CAPACITATOR Physics Department Electric an Magnetism Laboratory PARALLEL-PLATE CAPACITATOR 1. Goal. The goal of this practice is the stuy of the electric fiel an electric potential insie a parallelplate capacitor.

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

Calculus in the AP Physics C Course The Derivative

Calculus in the AP Physics C Course The Derivative Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

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

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

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

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

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

Chapter 6: Energy-Momentum Tensors

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

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,

More information

New Simple Controller Tuning Rules for Integrating and Stable or Unstable First Order plus Dead-Time Processes

New Simple Controller Tuning Rules for Integrating and Stable or Unstable First Order plus Dead-Time Processes Proceeings of the 3th WSEAS nternational Conference on SYSTEMS New Simple Controller Tuning Rules for ntegrating an Stable or Unstable First Orer plus Dea-Time Processes.G.ARVANTS Department of Natural

More information

In the usual geometric derivation of Bragg s Law one assumes that crystalline

In the usual geometric derivation of Bragg s Law one assumes that crystalline Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

05 The Continuum Limit and the Wave Equation

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

Physics 115C Homework 4

Physics 115C Homework 4 Physics 115C Homework 4 Problem 1 a In the Heisenberg picture, the ynamical equation is the Heisenberg equation of motion: for any operator Q H, we have Q H = 1 t i [Q H,H]+ Q H t where the partial erivative

More information

State Space Analysis of Power System Stability Enhancement with Used the STATCOM

State Space Analysis of Power System Stability Enhancement with Used the STATCOM tate pace Analysis of Power ystem tability Enhancement with Use the ACOM M. Mahavian () - G. hahgholian () () Department of Electrical Engineering, Islamic Aza University, Naein Branch, Esfahan, Iran ()

More information

FET Inrush Protection

FET Inrush Protection FET Inrush Protection Chris Pavlina https://semianalog.com 2015-11-23 CC0 1.0 Universal Abstract It is possible to use a simple one-transistor FET circuit to provie inrush protection for low voltage DC

More information

A Model of Electron-Positron Pair Formation

A Model of Electron-Positron Pair Formation Volume PROGRESS IN PHYSICS January, 8 A Moel of Electron-Positron Pair Formation Bo Lehnert Alfvén Laboratory, Royal Institute of Technology, S-44 Stockholm, Sween E-mail: Bo.Lehnert@ee.kth.se The elementary

More information

Differentiability, Computing Derivatives, Trig Review

Differentiability, Computing Derivatives, Trig Review Unit #3 : Differentiability, Computing Derivatives, Trig Review Goals: Determine when a function is ifferentiable at a point Relate the erivative graph to the the graph of an original function Compute

More information

Thermal runaway during blocking

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

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

4. CONTROL OF ZERO-SEQUENCE CURRENT IN PARALLEL THREE-PHASE CURRENT-BIDIRECTIONAL CONVERTERS

4. CONTROL OF ZERO-SEQUENCE CURRENT IN PARALLEL THREE-PHASE CURRENT-BIDIRECTIONAL CONVERTERS 4. CONRO OF ZERO-SEQUENCE CURREN IN PARAE HREE-PHASE CURREN-BIDIRECIONA CONVERERS 4. A NOVE ZERO-SEQUENCE CURREN CONRO 4.. Zero-Sequence Dynamics he parallel boost rectifier moel in Figure.4 an the parallel

More information

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.

More information

water adding dye partial mixing homogenization time

water adding dye partial mixing homogenization time iffusion iffusion is a process of mass transport that involves the movement of one atomic species into another. It occurs by ranom atomic jumps from one position to another an takes place in the gaseous,

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

Diagonalization of Matrices Dr. E. Jacobs

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

Nested Saturation with Guaranteed Real Poles 1

Nested Saturation with Guaranteed Real Poles 1 Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically

More information

A Review of Feedforward Control Approaches in Nanopositioning for High-Speed SPM

A Review of Feedforward Control Approaches in Nanopositioning for High-Speed SPM Garrett M. Clayton Department of Mechanical Engineering, Villanova University, Villanova, PA 1985 e-mail: garrett.clayton@villanova.eu Szuchi Tien Department of Mechanical Engineering, National Cheng Kung

More information

PHYS 414 Problem Set 2: Turtles all the way down

PHYS 414 Problem Set 2: Turtles all the way down PHYS 414 Problem Set 2: Turtles all the way own This problem set explores the common structure of ynamical theories in statistical physics as you pass from one length an time scale to another. Brownian

More information

UNIT 4:Capacitors and Dielectric

UNIT 4:Capacitors and Dielectric UNIT 4:apacitors an Dielectric SF7 4. apacitor A capacitor is a evice that is capable of storing electric charges or electric potential energy. It is consist of two conucting plates separate by a small

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Linear and quadratic approximation

Linear and quadratic approximation Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function

More information

Experiment 2, Physics 2BL

Experiment 2, Physics 2BL Experiment 2, Physics 2BL Deuction of Mass Distributions. Last Upate: 2009-05-03 Preparation Before this experiment, we recommen you review or familiarize yourself with the following: Chapters 4-6 in Taylor

More information

3.0 DETECTION THEORY. Figure 3-1 IF Receiver/Signal Processor Representation

3.0 DETECTION THEORY. Figure 3-1 IF Receiver/Signal Processor Representation 3.0 DEECON HEORY 3. NRODUCON n some of our raar range equation problems we looke at fining the etection range base on NRs of 3 an 0 B. We now want to evelop some of the theory that explains the use of

More information

Tutorial Test 5 2D welding robot

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

Numerical Integrator. Graphics

Numerical Integrator. Graphics 1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,

More information

The Press-Schechter mass function

The Press-Schechter mass function The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for

More information

Situation awareness of power system based on static voltage security region

Situation awareness of power system based on static voltage security region The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran

More information

Chapter 31: RLC Circuits. PHY2049: Chapter 31 1

Chapter 31: RLC Circuits. PHY2049: Chapter 31 1 Chapter 31: RLC Circuits PHY049: Chapter 31 1 LC Oscillations Conservation of energy Topics Dampe oscillations in RLC circuits Energy loss AC current RMS quantities Force oscillations Resistance, reactance,

More information

Conservation laws a simple application to the telegraph equation

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

Dynamics of the synchronous machine

Dynamics of the synchronous machine ELEC0047 - Power system ynamics, control an stability Dynamics of the synchronous machine Thierry Van Cutsem t.vancutsem@ulg.ac.be www.montefiore.ulg.ac.be/~vct These slies follow those presente in course

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

A Parametric Device Study for SiC Power Electronics

A Parametric Device Study for SiC Power Electronics A Parametric evice Stuy for SiC Power Electronics Burak Ozpineci urak@ieee.org epartment of Electrical an Computer Engineering The University of Tennessee Knoxville TN 7996- Leon M. Tolert tolert@utk.eu

More information

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010 NTNU Page of 6 Institutt for fysikk Fakultet for fysikk, informatikk og matematikk This solution consists of 6 pages. Solution to the exam in TFY423 STATISTICAL PHYSICS Wenesay ecember, 2 Problem. Particles

More information

Mechanical Vibrations Misc Topics Base Excitation & Seismic

Mechanical Vibrations Misc Topics Base Excitation & Seismic Mechanical Vibrations Misc Topics Base Excitation & Seismic Peter Avitabile Mechanical Engineering Department University of Massachusetts Lowell 1 Dr. Peter Avitabile Seismic events are often moele with

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

TIME-DELAY ESTIMATION USING FARROW-BASED FRACTIONAL-DELAY FIR FILTERS: FILTER APPROXIMATION VS. ESTIMATION ERRORS

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