ECEN 667 Power System Stability Lecture 21: Modal Analysis
|
|
- Rosemary Hoover
- 5 years ago
- Views:
Transcription
1 ECEN 667 Power System Stablty Lecture 21: Modal Analyss Prof. Tom Overbye Dept. of Electrcal and Computer Engneerng Texas A&M Unversty, 1
2 Announcements Read Chapter 8 Homework 7 s posted; due on Tuesday Nov 28 Fnal s as per TAMU schedule. That s, Frday Dec 8 from 3 to 5pm 2
3 Sngle Machne Infnte Bus A qute useful analyss technque s to consder the small sgnal stablty assocated wth a sngle generator connected to the rest of the system through an equvalent transmsson lne Drvng pont mpedance lookng nto the system s used to calculate the equvalent lne's mpedance The Z value can be calculated qute quckly usng sparse vector methods Rest of the system s assumed to be an nfnte bus wth ts voltage set to match the generator's real and reactve power njecton and voltage 3
4 Small SMIB Example As a small example, consder the 4 bus system shown below, n whch bus 2 really s an nfnte bus GENCLS Bus 4 X=0.1 Bus 1 Bus 2 X=0.2 Infnte Bus slack Deg pu X=0.1 Bus 3 X= Deg pu 4.46 Deg pu 0.00 Deg pu To get the SMIB for bus 4, frst calculate Z Ybus j Z44 j Z 44 s Z th n parallel wth jx' d,4 (whch s j0.3) so Z th s j0.22 4
5 Example: Bus 4 wth GENROU Model The egenvalues can be calculated for any set of generator models Ths example replaces the bus 4 generator classcal machne wth a GENROU model There are now sx egenvalues, wth the domnate response comng from the electro-mechancal mode wth a frequency of 1.83 Hz, and dampng of 6.92% 5
6 Example: Bus 4 wth GENROU Model and Excter Addng an relatvely slow EXST1 excter adds addtonal states (wth K A =200, T A =0.2) As the ntal reactve power output of the generator s decreased, the system becomes unstable Case s saved as B4_GENROU_Sat_SMIB 6
7 Example: Bus 4 wth GENROU Model and Excter Wth Q 4 = 25 Mvar the egenvalues are And wth Q 4 =0 Mvar the egenvalues are 7
8 Example: Bus 4 wth GENROU Model and Excter Graph shows response followng a short fault when Q4 s 0 Mvar gfedcb Rotor Angle_Gen Bus 4 #1 Ths motvates tryng to get addtonal nsght nto how to ncrease system dampng, whch s the goal of modal analyss 8
9 Modal Analyss - Comments Modal analyss (analyss of small sgnal stablty through egenvalue analyss) s at the core of SSA software In Modal Analyss one looks at: Egenvalues Egenvectors (left or rght) Partcpaton factors Mode shape Goal s to determne how the varous parameters affect the response of the system Power System Stablzer (PSS) desgn n a multmachne context s done usng modal analyss method. 9
10 Egenvalues, Rght Egenvectors For an n by n matrx A the egenvalues of A are the roots of the characterstc equaton: det[ A li] A li 0 Assume l 1 l n as dstnct (no repeated egenvalues). For each egenvalue l there exsts an egenvector such that: Av l v v s called a rght egenvector If l s complex, then v has complex entres 10
11 Left Egenvectors For each egenvalue l there exsts a left egenvector w such that: w A w l t t Equvalently, the left egenvector s the rght egenvector of A T ; that s, t A w lw 11
12 Egenvector Propertes The rght and left egenvectors are orthogonal.e. w v t 0, w v 0 ( j) t j We can normalze the egenvectors so that: w v t 1, w v 0 ( j) t j 12
13 Egenvector Example l 4 A, l A I 3 2 l 2 l l l1,2 2 3 (3) 4(10) , Rght Egenvectors l 1 5 v v 4v 5v Av1 5v1 v1 v 3v 2v 5v Smlarly, 4 l2 2 v choose v 1 v 1 v
14 Egenvector Example Left egenvectors t t l w1a w1 w11 w21 w11 w [ ] 5[ ] 3 2 w11 3w21 5w11 Let w21 4, then w11 3 4w 2w 5w w l w v v 3 w 4 w 1 Verfy w v 7, w v 7, w v 0, w v 0 t t t t t wv 1. We would lke to make Ths can be done n many ways. 14
15 1 3 1 Let W Then T W V Egenvector Example I Verfy It can be verfed that W T =V -1. The left and rght egenvectors are used n computng the partcpaton factor matrx. 15
16 Modal Matrces The devaton away from an equlbrum pont can be defned as Δx AΔx From ths equaton t s dffcult to determne how parameters n A affect a partcular x because of the varable couplng To decouple the problem frst defne the matrces of the rght and left egenvectors (the modal matrces) V [ v, v... v ] & W [ w, w,... w ] 1 2 n 1 2 AV VΛ when Λ Dag ( l ) n 16
17 It follows that Modal Matrces To decouple the varables defne z so Then 1 V AV Λ x Vz x Vz AΔx AVz 1 z V AVz WAVz Λz Snce s dagonal, the equatons are now uncoupled wth z l z So x( t) Vz( t) 17
18 Example Assume the prevous system wth 1 4 A V V AV
19 Modal Matrces Thus the response can be wrtten n terms of the ndvdual egenvalues and rght egenvectors as n 1 t x( t) v z ( 0) e l Furthermore wth 1 T Δx= VZ z V x W x Note, we are requrng that the egenvalues be dstnct! So z(t) can be wrtten as usng the left egenvectors as t z( t) W x( t) [ w w... w ] 1 2 n t x1() t xn () t 19
20 Modal Matrces We can then wrte the response x(t) n terms of the modes of the system t z ( t) w x( t) t z ( 0) w x( 0) c so x( t) n 1 v c e l t l1t l2t Expandng x ( t) v c e v c e... v c e n n So c s a scalar that represents the magntude of exctaton of the th mode from the ntal condtons l t n 20
21 Numercal example x1 0 1 x1 1, ( 0) x x x 2 4 Egenvalues are l 4, l Egenvectors are v1, 2 4 v Modal matrx V Normalze so V
22 Numercal example (contd) Left egenvector matrx s: W T V T z = W AVz z1 4 0z1 z 0 2 z
23 Numercal example (contd) 1 z 4z, z( 0) V x( 0) 1 1 z1( 0) z2 2z2, z2( 0) 0 4t 2t T z1( t) z1( 0) e ; z2( t) z2( 0) e, C W x( 0) 0 x = Vz x1( t) 1 1z1( t) x2( t) 4 2 z2( t) c1 z1() t c Because of the ntal condton, the 2 nd mode does not get excted 2 z2( t) c vz ( 0) e l 1 t 23
24 Mode Shape, Senstvty and Partcpaton Factors So we have t x( t) Vz( t), z( t) W x( t) x(t) are the orgnal state varables, z(t) are the transformed varables so that each varable s assocated wth only one mode. From the frst equaton the rght egenvector gves the mode shape.e. relatve actvty of state varables when a partcular mode s excted. For example the degree of actvty of state varable x k n v mode s gven by the element V k of the the rght egenvector matrx V 24
25 Mode Shape, Senstvty and Partcpaton Factors The magntude of elements of v gve the extent of actvtes of n state varables n the th mode and angles of elements (f complex) gve phase dsplacements of the state varables wth regard to the mode. The left egenvector w dentfes whch combnaton of orgnal state varables dsplay only the th mode. 25
26 Egenvalue Parameter Senstvty To derve the senstvty of the egenvalues to the parameters recall Av = l v ; take the partal dervatve wth respect to A kj by usng the chan rule A v l v v A v l A kj Akj Akj Akj Multply by w t A v l v w v w A w v w l t t t t Akj Akj Akj Akj t A t v t l w v w [ A l I] w v A A A kj kj kj 26
27 Egenvalue Parameter Senstvty t Ths s smplfed by notng that w ( A l I) by the defnton of w beng a left egenvector Therefore w t A A kj v Snce all elements of j th column s 1 Thus l WV k j A kj l A kj A A kj 0 are zero, except the k th row, 27
28 Senstvty Example In the prevous example we had A, l1,2 5, 2,, 3 2 V 1 3 W Then the senstvty of l 1 and l 2 to changes n A are l WV l k j, l Akj A A For example wth Aˆ 1 4, ˆ l1,2 5.61, 1.61,
29 Egenvalue Parameter Senstvty t Ths s smplfed by notng that w ( A l I) by the defnton of w beng a left egenvector Therefore w t A A kj v Snce all elements of j th column s 1 Thus l WV k j A kj l A kj A A kj 0 are zero, except the k th row, 29
30 Egenvalue Parameter Senstvty To derve the senstvty of the egenvalues to the parameters recall Av = l v ; take the partal dervatve wth respect to A kj by usng the chan rule A v l v v A v l A kj Akj Akj Akj Multply by w t A v l v w v w A w v w l t t t t Akj Akj Akj Akj t A t v t l w v w [ A l I] w v A A A kj kj kj 30
31 Egenvalue Parameter Senstvty t Ths s smplfed by notng that w ( A l I) by the defnton of w beng a left egenvector Therefore w t A A kj v Snce all elements of j th column s 1 Thus l WV k j A kj l A kj A A kj 0 are zero, except the k th row, 31
32 Senstvty Example In the prevous example we had A, l1,2 5, 2,, 3 2 V 1 3 W Then the senstvty of l 1 and l 2 to changes n A are l WV l k j, l Akj A A For example wth Aˆ 1 4, ˆ l1, ,
33 Partcpaton Factors The partcpaton factors, P k, are used to determne how much the k th state varable partcpates n the th mode P V W k k k The sum of the partcpaton factors for any mode or any varable sum to 1 The partcpaton factors are qute useful n relatng the egenvalues to portons of a model 33
34 Partcpaton Factors For the prevous example wth P k = V k W k and A,, 3 2 V 1 3 W We get P
35 PowerWorld SMIB Partcpaton Factors The magntudes of the partcpaton factors are shown on the PowerWorld SMIB dalog The below values are shown for the four bus example wth Q 4 = 0 Case s saved as B4_GENROU_Sat_SMIB_QZero 35
36 Measurement Based Modal Analyss Wth the advent of large numbers of PMUs, measurement based SSA s ncreasng used The goal s to determne the dampng assocated wth the domnant oscllatory modes n the system Approaches seek to approxmate a sampled sgnal by a seres of exponental functons (usually damped snusodals) Several technques are avalable wth Prony analyss the oldest Method, whch was developed by Gaspard Rche de Prony, dates to 1795; power system applcatons from about 1980's Here we'll consder a newer alternatve, based on the varable projecton method 36
37 Some Useful References J.F. Hauer, C.J. Demeure, and L.L. Scharf, "Intal results n Prony analyss of power system response sgnals," IEEE Trans. Power Systems, vol.5, pp 80-89, Feb 1990 D.J. Trudnowsk, J.M. Johnson, and J.F. Hauer, "Makng Prony analyss more accurate usng multple sgnals," IEEE Trans. Power Systems, vol.14, pp , Feb 1999 A. Borden, B.C. Leseutre, J. Gronqust, "Power System Modal Analyss Tool Developed for Industry Use," Proc North Amercan Power Symposum, Manhattan, KS, Sept
38 Varable Projecton Method (VPM) Idea of all technques s to approxmate a sgnal, y org (t), by the sum of other, smpler sgnals (bass functons) Bass functons are usually exponentals, wth lnear and quadratc functons also added to detrend the sgnal Propertes of the orgnal sgnal can be quantfed from bass functon propertes (such as frequency and dampng) Sgnal s consdered over an nterval wth t=0 at the begnnng Approaches work by samplng the orgnal sgnal y org (t) Vector y conssts of m unformly sampled ponts from y org (t) at a samplng value of T, startng wth t=0, wth values y j for j=1 m Tmes are then t j = (j-1)t 38
39 Varable Projecton Method (VPM) At each tme pont j, where t j = (j-1)t the approxmaton of y j s yˆ ( α) b ( t, α) j j 1 where α n s a vector wth the real and magnary j egenvalue components, wth ( t, α) e for a correspondng to a real egenvalue, and t j j +1 j 1 j ( t, α) e cos( t ) and ( α) e sn( +1t j ) for a complex egenvector value j t t 39
40 Varable Projecton Method (VPM) Error (resdual) value at each pont j s r ( t, α) y yˆ ( t, α) j j j j j s the vector contanng the optmzaton varables Functon beng mnmzed s 1 1 α r α 2 2 m 2 ( y ˆ j y j ( t j, )) ( ) j1 2 2 r() s the resdual vector Method teratvely changes to reduce the mnmzaton functon 40
41 Varable Projecton Method (VPM) A key nsght of the varable projecton method s that yˆ( α) Φ( α) b And then the resdual s mnmzed by selectng b Φ( α) y where Φα ( ) s the m by n matrx wth values j t j ( α) e f corresponds to a real egenvalue, t j j α e 1tj j1 j and ( ) cos( ) and ( α) e sn( t ) for a complex egenvalue; t j 1 T Fnally, Φ( α) s the pseudonverse of Φ( α) j t 1 j 41
ECEN 667 Power System Stability Lecture 23:Measurement Based Modal Analysis, FFT
ECEN 667 Power System Stability Lecture 23:Measurement Based Modal Analysis, FFT Prof. Tom Overbye Dept. of Electrical and Computer Engineering Texas A&M University, overbye@tamu.edu 1 Announcements Read
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More information1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
More informationSalmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
More informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
More informationPARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY
POZNAN UNIVE RSITY OF TE CHNOLOGY ACADE MIC JOURNALS No 86 Electrcal Engneerng 6 Volodymyr KONOVAL* Roman PRYTULA** PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY Ths paper provdes a
More informationU.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017
U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationVector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence.
Vector Norms Chapter 7 Iteratve Technques n Matrx Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematcs Unversty of Calforna, Berkeley Math 128B Numercal Analyss Defnton A vector norm
More informationIntroduction. - The Second Lyapunov Method. - The First Lyapunov Method
Stablty Analyss A. Khak Sedgh Control Systems Group Faculty of Electrcal and Computer Engneerng K. N. Toos Unversty of Technology February 2009 1 Introducton Stablty s the most promnent characterstc of
More informationρ some λ THE INVERSE POWER METHOD (or INVERSE ITERATION) , for , or (more usually) to
THE INVERSE POWER METHOD (or INVERSE ITERATION) -- applcaton of the Power method to A some fxed constant ρ (whch s called a shft), x λ ρ If the egenpars of A are { ( λ, x ) } ( ), or (more usually) to,
More information= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.
Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and
More informationNovember 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project
SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω
More informationHomework Notes Week 7
Homework Notes Week 7 Math 4 Sprng 4 #4 (a Complete the proof n example 5 that s an nner product (the Frobenus nner product on M n n (F In the example propertes (a and (d have already been verfed so we
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationAdvanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)
Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed
More information2.3 Nilpotent endomorphisms
s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms
More information1 Derivation of Point-to-Plane Minimization
1 Dervaton of Pont-to-Plane Mnmzaton Consder the Chen-Medon (pont-to-plane) framework for ICP. Assume we have a collecton of ponts (p, q ) wth normals n. We want to determne the optmal rotaton and translaton
More informationLinear Regression Analysis: Terminology and Notation
ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
More information1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:
68A Solutons to Exercses March 05 (a) Usng a Taylor expanson, and notng that n 0 for all n >, ( + ) ( + ( ) + ) We can t nvert / because there s no Taylor expanson around 0 Lets try to calculate the nverse
More informationMATH Sensitivity of Eigenvalue Problems
MATH 537- Senstvty of Egenvalue Problems Prelmnares Let A be an n n matrx, and let λ be an egenvalue of A, correspondngly there are vectors x and y such that Ax = λx and y H A = λy H Then x s called A
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus
More informationTensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q
For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased
More informationECEN 667 Power System Stability Lecture 20: Oscillations, Small Signal Stability Analysis
ECEN 667 Power System Stability Lecture 20: Oscillations, Small Signal Stability Analysis Prof. Tom Overbye Dept. of Electrical and Computer Engineering Texas A&M University, overbye@tamu.edu 1 Announcements
More informationC/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1
C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned
More informationUsing T.O.M to Estimate Parameter of distributions that have not Single Exponential Family
IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran
More informationLECTURE 9 CANONICAL CORRELATION ANALYSIS
LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationComparison of Regression Lines
STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence
More informationLecture 10: May 6, 2013
TTIC/CMSC 31150 Mathematcal Toolkt Sprng 013 Madhur Tulsan Lecture 10: May 6, 013 Scrbe: Wenje Luo In today s lecture, we manly talked about random walk on graphs and ntroduce the concept of graph expander,
More informationPhysics 1202: Lecture 11 Today s Agenda
Physcs 122: Lecture 11 Today s Agenda Announcements: Team problems start ths Thursday Team 1: Hend Ouda, Mke Glnsk, Stephane Auger Team 2: Analese Bruder, Krsten Dean, Alson Smth Offce hours: Monday 2:3-3:3
More informationConsistency & Convergence
/9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have
More informationQuantum Mechanics I - Session 4
Quantum Mechancs I - Sesson 4 Aprl 3, 05 Contents Operators Change of Bass 4 3 Egenvectors and Egenvalues 5 3. Denton....................................... 5 3. Rotaton n D....................................
More informationAnalysis of inter-area mode oscillations using a PSS for Variable Speed Wind Power Converter and Modal decomposition.
Analyss of nter-area mode oscllatons usng a PSS for Varable Speed Wnd Power Converter and Modal decomposton. C. GALLARDO and P. LEDESMA Electrcal Engneerng department Carlos III Unversty of Madrd Butarque
More informationU.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016
U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and
More informationAGC Introduction
. Introducton AGC 3 The prmary controller response to a load/generaton mbalance results n generaton adjustment so as to mantan load/generaton balance. However, due to droop, t also results n a non-zero
More informationTracking with Kalman Filter
Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationQuantum Mechanics for Scientists and Engineers. David Miller
Quantum Mechancs for Scentsts and Engneers Davd Mller Types of lnear operators Types of lnear operators Blnear expanson of operators Blnear expanson of lnear operators We know that we can expand functons
More informationMTH 263 Practice Test #1 Spring 1999
Pat Ross MTH 6 Practce Test # Sprng 999 Name. Fnd the area of the regon bounded by the graph r =acos (θ). Observe: Ths s a crcle of radus a, for r =acos (θ) r =a ³ x r r =ax x + y =ax x ax + y =0 x ax
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationLecture 8 Modal Analysis
Lecture 8 Modal Analyss 16.0 Release Introducton to ANSYS Mechancal 1 2015 ANSYS, Inc. February 27, 2015 Chapter Overvew In ths chapter free vbraton as well as pre-stressed vbraton analyses n Mechancal
More informationMEM 255 Introduction to Control Systems Review: Basics of Linear Algebra
MEM 255 Introducton to Control Systems Revew: Bascs of Lnear Algebra Harry G. Kwatny Department of Mechancal Engneerng & Mechancs Drexel Unversty Outlne Vectors Matrces MATLAB Advanced Topcs Vectors A
More informationELE B7 Power Systems Engineering. Power Flow- Introduction
ELE B7 Power Systems Engneerng Power Flow- Introducton Introducton to Load Flow Analyss The power flow s the backbone of the power system operaton, analyss and desgn. It s necessary for plannng, operaton,
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationInexact Newton Methods for Inverse Eigenvalue Problems
Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationConvergence of random processes
DS-GA 12 Lecture notes 6 Fall 216 Convergence of random processes 1 Introducton In these notes we study convergence of dscrete random processes. Ths allows to characterze phenomena such as the law of large
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationCSCE 790S Background Results
CSCE 790S Background Results Stephen A. Fenner September 8, 011 Abstract These results are background to the course CSCE 790S/CSCE 790B, Quantum Computaton and Informaton (Sprng 007 and Fall 011). Each
More information1 Convex Optimization
Convex Optmzaton We wll consder convex optmzaton problems. Namely, mnmzaton problems where the objectve s convex (we assume no constrants for now). Such problems often arse n machne learnng. For example,
More informationCALCULUS CLASSROOM CAPSULES
CALCULUS CLASSROOM CAPSULES SESSION S86 Dr. Sham Alfred Rartan Valley Communty College salfred@rartanval.edu 38th AMATYC Annual Conference Jacksonvlle, Florda November 8-, 202 2 Calculus Classroom Capsules
More informationTerm Project - select journal paper and outline. Completed analysis due end
EE 5200 - Lecture 30 Fr ov 4, 2016 Topcs for Today: Announcements Term Project - select journal paper and outlne. Completed analyss due end of Week 12. Submt va e-mal as mn-lecture.ppt wth voce narraton.
More informationDynamic Programming. Preview. Dynamic Programming. Dynamic Programming. Dynamic Programming (Example: Fibonacci Sequence)
/24/27 Prevew Fbonacc Sequence Longest Common Subsequence Dynamc programmng s a method for solvng complex problems by breakng them down nto smpler sub-problems. It s applcable to problems exhbtng the propertes
More informationprinceton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg
prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there
More informationAssortment Optimization under MNL
Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.
More informationFall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede
Fall 0 Analyss of Expermental easurements B. Esensten/rev. S. Errede We now reformulate the lnear Least Squares ethod n more general terms, sutable for (eventually extendng to the non-lnear case, and also
More information8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS
SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationMAE140 - Linear Circuits - Winter 16 Final, March 16, 2016
ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationHongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)
ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of
More informationDeriving the X-Z Identity from Auxiliary Space Method
Dervng the X-Z Identty from Auxlary Space Method Long Chen Department of Mathematcs, Unversty of Calforna at Irvne, Irvne, CA 92697 chenlong@math.uc.edu 1 Iteratve Methods In ths paper we dscuss teratve
More informationNegative Binomial Regression
STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...
More informationSingular Value Decomposition: Theory and Applications
Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real
More informationFor now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.
Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson
More information9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set
More informationFrom Biot-Savart Law to Divergence of B (1)
From Bot-Savart Law to Dvergence of B (1) Let s prove that Bot-Savart gves us B (r ) = 0 for an arbtrary current densty. Frst take the dvergence of both sdes of Bot-Savart. The dervatve s wth respect to
More informationC4B Machine Learning Answers II. = σ(z) (1 σ(z)) 1 1 e z. e z = σ(1 σ) (1 + e z )
C4B Machne Learnng Answers II.(a) Show that for the logstc sgmod functon dσ(z) dz = σ(z) ( σ(z)) A. Zsserman, Hlary Term 20 Start from the defnton of σ(z) Note that Then σ(z) = σ = dσ(z) dz = + e z e z
More informationG4023 Mid-Term Exam #1 Solutions
Exam1Solutons.nb 1 G03 Md-Term Exam #1 Solutons 1-Oct-0, 1:10 p.m to :5 p.m n 1 Pupn Ths exam s open-book, open-notes. You may also use prnt-outs of the homework solutons and a calculator. 1 (30 ponts,
More informationOn a direct solver for linear least squares problems
ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationComposite Hypotheses testing
Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter
More informationBernoulli Numbers and Polynomials
Bernoull Numbers and Polynomals T. Muthukumar tmk@tk.ac.n 17 Jun 2014 The sum of frst n natural numbers 1, 2, 3,..., n s n n(n + 1 S 1 (n := m = = n2 2 2 + n 2. Ths formula can be derved by notng that
More informationThe Dirac Equation for a One-electron atom. In this section we will derive the Dirac equation for a one-electron atom.
The Drac Equaton for a One-electron atom In ths secton we wll derve the Drac equaton for a one-electron atom. Accordng to Ensten the energy of a artcle wth rest mass m movng wth a velocty V s gven by E
More information6.3.7 Example with Runga Kutta 4 th order method
6.3.7 Example wth Runga Kutta 4 th order method Agan, as an example, 3 machne, 9 bus system shown n Fg. 6.4 s agan consdered. Intally, the dampng of the generators are neglected (.e. d = 0 for = 1, 2,
More informationPerron Vectors of an Irreducible Nonnegative Interval Matrix
Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of
More information