Real Asymmetric Matrix Eigenvalue Analysis

Size: px
Start display at page:

Download "Real Asymmetric Matrix Eigenvalue Analysis"

Transcription

1 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Rel Asymmetri Mtrix Eigenvlue Anlysis Heewoo Lee omputtionl Menis Lbortory Deprtment of Menil Engineering n Applie Menis University of Miign Ann Arbor, MI Te University of Miign omputtionl Menis Lbortory

2 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis ontents of Presenttion Numeril Algoritms QR QZ Exmples Te University of Miign omputtionl Menis Lbortory

3 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Eigenvlue problem Ax λbx x λx, B A > Hessenberg form O(n ) > O(n ) QR lgoritm Te University of Miign omputtionl Menis Lbortory

4 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Hessenberg Mtrix Elementry similrity trnsformtion Ortogonl trnsformtion Te University of Miign omputtionl Menis Lbortory

5 Te University of Miign omputtionl Menis Lbortory Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Hessenberg Mtrix Mximum i+ m? row & olumn internge i m

6 Te University of Miign omputtionl Menis Lbortory Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Hessenberg Mtrix ,,, mx( ) i y i,,,, + y y ,,, mx( ) i y i,,,, + y y

7 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Hessenberg Mtrix Mximum i+ m? y, i mx( ) i y + y,,,,,,, Te University of Miign omputtionl Menis Lbortory

8 Te University of Miign omputtionl Menis Lbortory Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Single sift QR omplex version of two step QR Double sift QR QR lgoritm λ λ λ λ X λ λ λ λ X X

9 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Single sift QR lgoritm Q + H Q Q is unitry mtrix H Q R R is upper tringulr mtrix H Q ( ζ I) R Te University of Miign omputtionl Menis Lbortory

10 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis omplex version of two step QR Q + H Q H Q ( ζ I) R ζ ζ Q + H + Q + + Q + H + ( + ζ + I) R + Te University of Miign omputtionl Menis Lbortory

11 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Q Q Double sift QR lgoritm Q Q Q Q T T T T + Q ( ζ I)( ζ + I) R R + Q Q, R R + R, ( ζ I)( ζ + ) Γ + I T Q Q +, Q Γ R Te University of Miign omputtionl Menis Lbortory

12 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Double sift QR lgoritm If mtrix H is erive from su tt Qˆ QH ˆ Qˆ T or Qˆ H were Qˆ is ortogonl H is upper Hessenberg If Qˆ s te sme first olumn s Q, ten Te University of Miign Q ˆ Q n + H omputtionl Menis Lbortory

13 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Double sift QR lgoritm P I w r r w T r T T T T P P P P P P P P n n H ( γ γ, γ,,,) T, γ γ γ ( ( ζ + + ζ ) + ζ ζ ζ ζ ) + Te University of Miign omputtionl Menis Lbortory

14 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Double sift QR 5 n n n n nn P T P ' ' ' ' ' ' ' ' ' ' ' ' 5 n n n n nn Te University of Miign omputtionl Menis Lbortory

15 Te University of Miign omputtionl Menis Lbortory Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Double sift QR ( ) T r r r,,,,,,,,, χ β α

16 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Sifting Rule Before e itertion, lulte te roots of te lst prinipl submtrix η n η + Te oie of te origin sifts ζ n ζ + for te itertion epens on; η η η + ρ, ρ η η + η If tey re bot greter tn 5, ζ ζ + if tey re bot less tn 5, ζ η n ζ η + + oterwise we set bot ζ n ζ + to be te rel prt of eiter η or η +, wiever orrespons to te quntity less tn 5 Te University of Miign omputtionl Menis Lbortory

17 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Double sift QR Te University of Miign omputtionl Menis Lbortory

18 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Double sift QR E E Te University of Miign omputtionl Menis Lbortory

19 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Double sift QR lgoritm λ 57, x ± 76, 5 ± 65i i 5 ± 59i 68 ± 8i Te University of Miign omputtionl Menis Lbortory

20 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis MATLAB Eig(A,B) were A is te sme s n B is ientity mtrix Eigenvetors i i i i i i i i i i i i i i i i Eigenvlues E-7i i E-5i i Te University of Miign omputtionl Menis Lbortory

21 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis Nstrn upper Hessenberg meto λ i, x ± 76, i i i 6± 779i ± 65i 7 ± 8i i i i i Te University of Miign omputtionl Menis Lbortory

22 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis omplex Eigenvlue Solver MS/NASTRAN QZ lgoritm Rel Imginry Rel Imginry Te University of Miign omputtionl Menis Lbortory

23 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis QZ lgoritm Ax λbx QAZy λqbzy, x Zy A is reue to upper Hessenberg B is reue to upper tringulr A is reue to qusi-tringulr (generlition of QR) Te University of Miign omputtionl Menis Lbortory

24 Development of Optiml Design Meto for Bre Squel Noise Bse on omplex Eigenvlue Anlysis onlusions For rel symmetri eigenvlue problem, ouble sift QR or QZ lgoritm soul be use Oter metos proue fititious omplex eigenvlues Te University of Miign omputtionl Menis Lbortory

Lecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University

Lecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University Lecture Note 4: Numericl differentition nd integrtion Xioqun Zng Sngi Jio Tong University Lst updted: November, 0 Numericl Anlysis. Numericl differentition.. Introduction Find n pproximtion of f (x 0 ),

More information

4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX. be a real symmetric matrix. ; (where we choose θ π for.

4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX. be a real symmetric matrix. ; (where we choose θ π for. 4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX Some reliminries: Let A be rel symmetric mtrix. Let Cos θ ; (where we choose θ π for Cos θ 4 purposes of convergence of the scheme)

More information

Linear Algebra Introduction

Linear Algebra Introduction Introdution Wht is Liner Alger out? Liner Alger is rnh of mthemtis whih emerged yers k nd ws one of the pioneer rnhes of mthemtis Though, initilly it strted with solving of the simple liner eqution x +

More information

Contents. Outline. Structured Rank Matrices Lecture 2: The theorem Proofs Examples related to structured ranks References. Structure Transport

Contents. Outline. Structured Rank Matrices Lecture 2: The theorem Proofs Examples related to structured ranks References. Structure Transport Contents Structured Rnk Mtrices Lecture 2: Mrc Vn Brel nd Rf Vndebril Dept. of Computer Science, K.U.Leuven, Belgium Chemnitz, Germny, 26-30 September 2011 1 Exmples relted to structured rnks 2 2 / 26

More information

Chapter 2. Determinants

Chapter 2. Determinants Chpter Determinnts The Determinnt Function Recll tht the X mtrix A c b d is invertible if d-bc0. The expression d-bc occurs so frequently tht it hs nme; it is clled the determinnt of the mtrix A nd is

More information

"Add"-operator "Mul"-operator "Pow"-operator. def. h b. def

Add-operator Mul-operator Pow-operator. def. h b. def Opertors A sort review of opertors. Te isussions out tetrtion le me to two impressions. ) It my e etter to see opertors using prmeters, inste of two, s it is ommon use upte 4 ) Sering for noter onsistent

More information

Tangent Lines-1. Tangent Lines

Tangent Lines-1. Tangent Lines Tangent Lines- Tangent Lines In geometry, te tangent line to a circle wit centre O at a point A on te circle is defined to be te perpendicular line at A to te line OA. Te tangent lines ave te special property

More information

Lecture Note 9: Orthogonal Reduction

Lecture Note 9: Orthogonal Reduction MATH : Computtionl Methods of Liner Algebr 1 The Row Echelon Form Lecture Note 9: Orthogonl Reduction Our trget is to solve the norml eution: Xinyi Zeng Deprtment of Mthemticl Sciences, UTEP A t Ax = A

More information

18.06 Problem Set 4 Due Wednesday, Oct. 11, 2006 at 4:00 p.m. in 2-106

18.06 Problem Set 4 Due Wednesday, Oct. 11, 2006 at 4:00 p.m. in 2-106 8. Problem Set Due Wenesy, Ot., t : p.m. in - Problem Mony / Consier the eight vetors 5, 5, 5,..., () List ll of the one-element, linerly epenent sets forme from these. (b) Wht re the two-element, linerly

More information

Section 2.1 Special Right Triangles

Section 2.1 Special Right Triangles Se..1 Speil Rigt Tringles 49 Te --90 Tringle Setion.1 Speil Rigt Tringles Te --90 tringle (or just 0-60-90) is so nme euse of its ngle mesures. Te lengts of te sies, toug, ve very speifi pttern to tem

More information

arxiv: v2 [math.nt] 2 Feb 2015

arxiv: v2 [math.nt] 2 Feb 2015 rxiv:407666v [mthnt] Fe 05 Integer Powers of Complex Tridigonl Anti-Tridigonl Mtrices Htice Kür Duru &Durmuş Bozkurt Deprtment of Mthemtics, Science Fculty of Selçuk University Jnury, 08 Astrct In this

More information

Probabilistic Investigation of Sensitivities of Advanced Test- Analysis Model Correlation Methods

Probabilistic Investigation of Sensitivities of Advanced Test- Analysis Model Correlation Methods Probbilistic Investigtion of Sensitivities of Advnced Test- Anlysis Model Correltion Methods Liz Bergmn, Mtthew S. Allen, nd Dniel C. Kmmer Dept. of Engineering Physics University of Wisconsin-Mdison Rndll

More information

Topic 6b Finite Difference Approximations

Topic 6b Finite Difference Approximations /8/8 Course Instructor Dr. Rymond C. Rump Oice: A 7 Pone: (95) 747 6958 E Mil: rcrump@utep.edu Topic 6b Finite Dierence Approximtions EE 486/5 Computtionl Metods in EE Outline Wt re inite dierence pproximtions?

More information

Introduction to Arnoldi method

Introduction to Arnoldi method Introduction to Arnoldi method SF2524 - Matrix Computations for Large-scale Systems KTH Royal Institute of Technology (Elias Jarlebring) 2014-11-07 KTH Royal Institute of Technology (Elias Jarlebring)Introduction

More information

Numerical Linear Algebra Assignment 008

Numerical Linear Algebra Assignment 008 Numericl Liner Algebr Assignment 008 Nguyen Qun B Hong Students t Fculty of Mth nd Computer Science, Ho Chi Minh University of Science, Vietnm emil. nguyenqunbhong@gmil.com blog. http://hongnguyenqunb.wordpress.com

More information

Engineering Analysis ENG 3420 Fall Dan C. Marinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00

Engineering Analysis ENG 3420 Fall Dan C. Marinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00 Engineering Anlysis ENG 3420 Fll 2009 Dn C. Mrinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00 Lecture 13 Lst time: Problem solving in preprtion for the quiz Liner Algebr Concepts Vector Spces,

More information

Math 61CM - Solutions to homework 9

Math 61CM - Solutions to homework 9 Mth 61CM - Solutions to homework 9 Cédric De Groote November 30 th, 2018 Problem 1: Recll tht the left limit of function f t point c is defined s follows: lim f(x) = l x c if for ny > 0 there exists δ

More information

Eigenvectors and Eigenvalues

Eigenvectors and Eigenvalues MTB 050 1 ORIGIN 1 Eigenvets n Eigenvlues This wksheet esries the lger use to lulte "prinipl" "hrteristi" iretions lle Eigenvets n the "prinipl" "hrteristi" vlues lle Eigenvlues ssoite with these iretions.

More information

Position Analysis: Review (Chapter 2) Objective: Given the geometry of a mechanism and the input motion, find the output motion

Position Analysis: Review (Chapter 2) Objective: Given the geometry of a mechanism and the input motion, find the output motion Position Anlysis: Review (Chpter Ojetive: Given the geometry of mehnism n the input motion, fin the output motion Grphil pproh Algeri position nlysis Exmple of grphil nlysis of linges, four r linge. Given

More information

MATRICES AND VECTORS SPACE

MATRICES AND VECTORS SPACE MATRICES AND VECTORS SPACE MATRICES AND MATRIX OPERATIONS SYSTEM OF LINEAR EQUATIONS DETERMINANTS VECTORS IN -SPACE AND -SPACE GENERAL VECTOR SPACES INNER PRODUCT SPACES EIGENVALUES, EIGENVECTORS LINEAR

More information

MTH 5102 Linear Algebra Practice Exam 1 - Solutions Feb. 9, 2016

MTH 5102 Linear Algebra Practice Exam 1 - Solutions Feb. 9, 2016 Nme (Lst nme, First nme): MTH 502 Liner Algebr Prctice Exm - Solutions Feb 9, 206 Exm Instructions: You hve hour & 0 minutes to complete the exm There re totl of 6 problems You must show your work Prtil

More information

Pythagorean Theorem and Trigonometry

Pythagorean Theorem and Trigonometry Ptgoren Teorem nd Trigonometr Te Ptgoren Teorem is nient, well-known, nd importnt. It s lrge numer of different proofs, inluding one disovered merin President Jmes. Grfield. Te we site ttp://www.ut-te-knot.org/ptgors/inde.stml

More information

E E I M (E, I) E I 2 E M I I X I Y X Y I X, Y I X > Y x X \ Y Y {x} I B E B M E C E C C M r E X E r (X) X X r (X) = X E B M X E Y E X Y X B E F E F F E E E M M M M M M E B M E \ B M M 0 M M M 0 0 M x M

More information

1 ode.mcd. Find solution to ODE dy/dx=f(x,y). Instructor: Nam Sun Wang

1 ode.mcd. Find solution to ODE dy/dx=f(x,y). Instructor: Nam Sun Wang Fin solution to ODE /=f(). Instructor: Nam Sun Wang oe.mc Backgroun. Wen a sstem canges wit time or wit location, a set of ifferential equations tat contains erivative terms "/" escribe suc a namic sstem.

More information

K e sub x e sub n s i sub o o f K.. w ich i sub s.. u ra to the power of m i sub fi ed.. a sub t to the power of a

K e sub x e sub n s i sub o o f K.. w ich i sub s.. u ra to the power of m i sub fi ed.. a sub t to the power of a - ; ; ˆ ; q x ; j [ ; ; ˆ ˆ [ ˆ ˆ ˆ - x - - ; x j - - - - - ˆ x j ˆ ˆ ; x ; j κ ˆ - - - ; - - - ; ˆ σ x j ; ˆ [ ; ] q x σ; x - ˆ - ; J -- F - - ; x - -x - - x - - ; ; 9 S j P R S 3 q 47 q F x j x ; [ ]

More information

Differential Calculus: Differentiation (First Principles, Rules) and Sketching Graphs (Grade 12) *

Differential Calculus: Differentiation (First Principles, Rules) and Sketching Graphs (Grade 12) * OpenStax-CNX moule: m39313 1 Differential Calculus: Differentiation (First Principles, Rules) an Sketcing Graps (Grae 12) * Free Hig Scool Science Texts Project Tis work is prouce by OpenStax-CNX an license

More information

Precalculus Notes: Unit 6 Law of Sines & Cosines, Vectors, & Complex Numbers. A can be rewritten as

Precalculus Notes: Unit 6 Law of Sines & Cosines, Vectors, & Complex Numbers. A can be rewritten as Dte: 6.1 Lw of Sines Syllus Ojetie: 3.5 Te student will sole pplition prolems inoling tringles (Lw of Sines). Deriing te Lw of Sines: Consider te two tringles. C C In te ute tringle, sin In te otuse tringle,

More information

Matrices. Elementary Matrix Theory. Definition of a Matrix. Matrix Elements:

Matrices. Elementary Matrix Theory. Definition of a Matrix. Matrix Elements: Mtrices Elementry Mtrix Theory It is often desirble to use mtrix nottion to simplify complex mthemticl expressions. The simplifying mtrix nottion usully mkes the equtions much esier to hndle nd mnipulte.

More information

REGULARITY OF NONLOCAL MINIMAL CONES IN DIMENSION 2

REGULARITY OF NONLOCAL MINIMAL CONES IN DIMENSION 2 EGULAITY OF NONLOCAL MINIMAL CONES IN DIMENSION 2 OVIDIU SAVIN AND ENICO VALDINOCI Abstrct. We show tht the only nonlocl s-miniml cones in 2 re the trivil ones for ll s 0, 1). As consequence we obtin tht

More information

University of Houston, Department of Mathematics Numerical Analysis II

University of Houston, Department of Mathematics Numerical Analysis II University of Houston, Deprtment of Mtemtics Numericl Anlysis II 6 Glerkin metod, finite differences nd colloction 6.1 Glerkin metod Consider sclr 2nd order ordinry differentil eqution in selfdjoint form

More information

Matrix Eigenvalues and Eigenvectors September 13, 2017

Matrix Eigenvalues and Eigenvectors September 13, 2017 Mtri Eigenvlues nd Eigenvectors September, 7 Mtri Eigenvlues nd Eigenvectors Lrry Cretto Mechnicl Engineering 5A Seminr in Engineering Anlysis September, 7 Outline Review lst lecture Definition of eigenvlues

More information

Edexcel Level 3 Advanced GCE in Mathematics (9MA0) Two-year Scheme of Work

Edexcel Level 3 Advanced GCE in Mathematics (9MA0) Two-year Scheme of Work Eexel Level 3 Avne GCE in Mthemtis (9MA0) Two-yer Sheme of Work Stuents stuying A Level Mthemtis will tke 3 ppers t the en of Yer 13 s inite elow. All stuents will stuy Pure, Sttistis n Mehnis. A level

More information

Department of Physical Pharmacy and Pharmacokinetics Poznań University of Medical Sciences Pharmacokinetics laboratory

Department of Physical Pharmacy and Pharmacokinetics Poznań University of Medical Sciences Pharmacokinetics laboratory Deprtment of Physicl Phrmcy nd Phrmcoinetics Poznń University of Medicl Sciences Phrmcoinetics lbortory Experiment 1 Phrmcoinetics of ibuprofen s n exmple of the first-order inetics in n open one-comprtment

More information

Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Quantum Optical Communication

Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Quantum Optical Communication Msschusetts Institute of Technology Deprtment of Electricl Engineering nd Computer Science 6.453 Quntum Opticl Communiction Problem Set 6 Fll 2004 Issued: Wednesdy, October 13, 2004 Due: Wednesdy, October

More information

Matrix Solution to Linear Equations and Markov Chains

Matrix Solution to Linear Equations and Markov Chains Trding Systems nd Methods, Fifth Edition By Perry J. Kufmn Copyright 2005, 2013 by Perry J. Kufmn APPENDIX 2 Mtrix Solution to Liner Equtions nd Mrkov Chins DIRECT SOLUTION AND CONVERGENCE METHOD Before

More information

Matrix & Vector Basic Linear Algebra & Calculus

Matrix & Vector Basic Linear Algebra & Calculus Mtrix & Vector Bsic Liner lgebr & lculus Wht is mtrix? rectngulr rry of numbers (we will concentrte on rel numbers). nxm mtrix hs n rows n m columns M x4 M M M M M M M M M M M M 4 4 4 First row Secon row

More information

Fractals on non-euclidean metric

Fractals on non-euclidean metric Frctls on non-eucliden metric Yery Cchón Sntn April, 8 As fr s I know, there is no study on frctls on non eucliden metrics.this pper proposes rst pproch method bout generting frctls on non-eucliden metric.

More information

Math 270A: Numerical Linear Algebra

Math 270A: Numerical Linear Algebra Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner

More information

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these. Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra

More information

POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY

POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 1 Sturm-Liouville Theory In the two preceing lectures I emonstrte the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series re just the tip of the iceerg of the theory

More information

Polynomial Functions. Linear Functions. Precalculus: Linear and Quadratic Functions

Polynomial Functions. Linear Functions. Precalculus: Linear and Quadratic Functions Concepts: definition of polynomial functions, linear functions tree representations), transformation of y = x to get y = mx + b, quadratic functions axis of symmetry, vertex, x-intercepts), transformations

More information

a a a a a a a a a a a a a a a a a a a a a a a a In this section, we introduce a general formula for computing determinants.

a a a a a a a a a a a a a a a a a a a a a a a a In this section, we introduce a general formula for computing determinants. Section 9 The Lplce Expnsion In the lst section, we defined the determinnt of (3 3) mtrix A 12 to be 22 12 21 22 2231 22 12 21. In this section, we introduce generl formul for computing determinnts. Rewriting

More information

( ) 2. ( ) is the Fourier transform of! ( x). ( ) ( ) ( ) = Ae i kx"#t ( ) = 1 2" ( )"( x,t) PC 3101 Quantum Mechanics Section 1

( ) 2. ( ) is the Fourier transform of! ( x). ( ) ( ) ( ) = Ae i kx#t ( ) = 1 2 ( )( x,t) PC 3101 Quantum Mechanics Section 1 1. 1D Schrödinger Eqution G chpters 3-4. 1.1 the Free Prticle V 0 "( x,t) i = 2 t 2m x,t = Ae i kxt "( x,t) x 2 where = k 2 2m. Normliztion must hppen: 2 x,t = 1 Here, however: " A 2 dx " " As this integrl

More information

Discrete Least-squares Approximations

Discrete Least-squares Approximations Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve

More information

KENDRIYA VIDYALAYA IIT KANPUR HOME ASSIGNMENTS FOR SUMMER VACATIONS CLASS - XII MATHEMATICS (Relations and Functions & Binary Operations)

KENDRIYA VIDYALAYA IIT KANPUR HOME ASSIGNMENTS FOR SUMMER VACATIONS CLASS - XII MATHEMATICS (Relations and Functions & Binary Operations) KENDRIYA VIDYALAYA IIT KANPUR HOME ASSIGNMENTS FOR SUMMER VACATIONS 6-7 CLASS - XII MATHEMATICS (Reltions nd Funtions & Binry Opertions) For Slow Lerners: - A Reltion is sid to e Reflexive if.. every A

More information

Now we must transform the original model so we can use the new parameters. = S max. Recruits

Now we must transform the original model so we can use the new parameters. = S max. Recruits MODEL FOR VARIABLE RECRUITMENT (ontinue) Alterntive Prmeteriztions of the pwner-reruit Moels We n write ny moel in numerous ifferent ut equivlent forms. Uner ertin irumstnes it is onvenient to work with

More information

TANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION

TANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION Indin Journl of Mthemtics nd Mthemticl Sciences Vol. 7, No., (June ) : 9-38 TANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION

More information

L 2 STABILITY ANALYSIS OF THE CENTRAL DISCONTINUOUS GALERKIN METHOD AND A COMPARISON BETWEEN THE CENTRAL AND REGULAR DISCONTINUOUS GALERKIN METHODS

L 2 STABILITY ANALYSIS OF THE CENTRAL DISCONTINUOUS GALERKIN METHOD AND A COMPARISON BETWEEN THE CENTRAL AND REGULAR DISCONTINUOUS GALERKIN METHODS ESAIM: MAN 008 593 607 DOI: 0.05/mn:00808 ESAIM: Mtemticl Modelling nd Numericl Anlysis www.esim-mn.org L STABILITY ANALYSIS OF THE CENTRAL DISCONTINUOUS GALERKIN METHOD AND A COMPARISON BETWEEN THE CENTRAL

More information

REPRESENTATION THEORY OF PSL 2 (q)

REPRESENTATION THEORY OF PSL 2 (q) REPRESENTATION THEORY OF PSL (q) YAQIAO LI Following re notes from book [1]. The im is to show the qusirndomness of PSL (q), i.e., the group hs no low dimensionl representtion. 1. Representtion Theory

More information

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40 Mth B Prof. Audrey Terrs HW # Solutions by Alex Eustis Due Tuesdy, Oct. 9 Section 5. #7,, 6,, 5; Section 5. #8, 9, 5,, 7, 3; Section 5.3 #4, 6, 9, 3, 6, 8, 3; Section 5.4 #7, 8,, 3, 5, 9, 4 5..7 Since

More information

SOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS

SOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS ELM Numericl Anlysis Dr Muhrrem Mercimek SOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS ELM Numericl Anlysis Some of the contents re dopted from Lurene V. Fusett, Applied Numericl Anlysis using MATLAB.

More information

A P P E N D I X POWERS OF TEN AND SCIENTIFIC NOTATION A P P E N D I X SIGNIFICANT FIGURES

A P P E N D I X POWERS OF TEN AND SCIENTIFIC NOTATION A P P E N D I X SIGNIFICANT FIGURES A POWERS OF TEN AND SCIENTIFIC NOTATION In science, very lrge nd very smll deciml numbers re conveniently expressed in terms of powers of ten, some of wic re listed below: 0 3 0 0 0 000 0 3 0 0 0 0.00

More information

TRANSVERSE VIBRATION OF A BEAM VIA THE FINITE ELEMENT METHOD Revision E

TRANSVERSE VIBRATION OF A BEAM VIA THE FINITE ELEMENT METHOD Revision E RANSVERSE VIBRAION OF A BEAM VIA HE FINIE ELEMEN MEHOD Revision E B om Irvine Emil: tomirvine@ol.com November 8 8 Introuction Mn structures re too complex for nlsis vi clssicl meto. Close-form solutions

More information

Stochastic Programming Project Konrad Borys. Model for Optical Fiber Manufacturing

Stochastic Programming Project Konrad Borys. Model for Optical Fiber Manufacturing Stochstic Progrmming Project Konrd Borys Model for Opticl Fiber Mnufcturing. Introduction Opticl fibers re mde of solid rods of glss clled preforms. he s of the preforms re heted nd fibers re drwn from

More information

Symmetry Labeling of Molecular Energies

Symmetry Labeling of Molecular Energies Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry

More information

THE QUADRATIC RECIPROCITY LAW OF DUKE-HOPKINS. Circa 1870, G. Zolotarev observed that the Legendre symbol ( a p

THE QUADRATIC RECIPROCITY LAW OF DUKE-HOPKINS. Circa 1870, G. Zolotarev observed that the Legendre symbol ( a p THE QUADRATIC RECIPROCITY LAW OF DUKE-HOPKINS PETE L CLARK Circ 1870, Zolotrev observed tht the Legendre symbol ( p ) cn be interpreted s the sign of multipliction by viewed s permuttion of the set Z/pZ

More information

Digital Filter Structures

Digital Filter Structures Digital Filter Structures Te convolution sum description of an LTI discrete-time system can, in principle, be used to implement te system For an IIR finite-dimensional system tis approac is not practical

More information

Chapter 3 Polynomials

Chapter 3 Polynomials Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling

More information

Chapter 18 Two-Port Circuits

Chapter 18 Two-Port Circuits Cpter 8 Two-Port Circuits 8. Te Terminl Equtions 8. Te Two-Port Prmeters 8.3 Anlysis of te Terminted Two-Port Circuit 8.4 nterconnected Two-Port Circuits Motivtion Tévenin nd Norton equivlent circuits

More information

(a 1 m. a n m = < a 1/N n

(a 1 m. a n m = < a 1/N n Notes on a an log a Mat 9 Fall 2004 Here is an approac to te eponential an logaritmic functions wic avois any use of integral calculus We use witout proof te eistence of certain limits an assume tat certain

More information

1 Orthogonalisation in finite precision arithmetic

1 Orthogonalisation in finite precision arithmetic 1 Orthogonlistion in finite precision rithmetic We investigte the differences nd similrities between the following four wys to compute the QR-decomposition of given rectngulr mtrix A C m n in Mtlb: (CGS)

More information

Introduction to Determinants. Remarks. Remarks. The determinant applies in the case of square matrices

Introduction to Determinants. Remarks. Remarks. The determinant applies in the case of square matrices Introduction to Determinnts Remrks The determinnt pplies in the cse of squre mtrices squre mtrix is nonsingulr if nd only if its determinnt not zero, hence the term determinnt Nonsingulr mtrices re sometimes

More information

Path product and inverse M-matrices

Path product and inverse M-matrices Electronic Journl of Liner Algebr Volume 22 Volume 22 (2011) Article 42 2011 Pth product nd inverse M-mtrices Yn Zhu Cheng-Yi Zhng Jun Liu Follow this nd dditionl works t: http://repository.uwyo.edu/el

More information

Fuzzy transform to approximate solution of boundary value problems via optimal coefficients

Fuzzy transform to approximate solution of boundary value problems via optimal coefficients 217 Interntionl Conference on Hig Performnce Computing & Simultion Fuzzy trnsform to pproximte solution of boundry vlue problems vi optiml coefficients Zr Alijni Mtemtics nd Sttistics Dept University of

More information

Matrices SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics (c) 1. Definition of a Matrix

Matrices SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics (c) 1. Definition of a Matrix tries Definition of tri mtri is regulr rry of numers enlosed inside rkets SCHOOL OF ENGINEERING & UIL ENVIRONEN Emple he following re ll mtries: ), ) 9, themtis ), d) tries Definition of tri Size of tri

More information

232 Calculus and Structures

232 Calculus and Structures 3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE

More information

SOLUTIONS FOR ANALYSIS QUALIFYING EXAM, FALL (1 + µ(f n )) f(x) =. But we don t need the exact bound.) Set

SOLUTIONS FOR ANALYSIS QUALIFYING EXAM, FALL (1 + µ(f n )) f(x) =. But we don t need the exact bound.) Set SOLUTIONS FOR ANALYSIS QUALIFYING EXAM, FALL 28 Nottion: N {, 2, 3,...}. (Tht is, N.. Let (X, M be mesurble spce with σ-finite positive mesure µ. Prove tht there is finite positive mesure ν on (X, M such

More information

Songklanakarin Journal of Science and Technology SJST R1 Akram. N-Fuzzy BiΓ -Ternary Semigroups

Songklanakarin Journal of Science and Technology SJST R1 Akram. N-Fuzzy BiΓ -Ternary Semigroups ongklnkrin Journl of cience nd Technology JT-0-0.R Akrm N-Fuzzy Bi -Ternry emigroups Journl: ongklnkrin Journl of cience nd Technology Mnuscript ID JT-0-0.R Mnuscript Type: Originl Article Dte ubmitted

More information

Matrices 13: determinant properties and rules continued

Matrices 13: determinant properties and rules continued Mtrices : determinnt properties nd rules continued nthony Rossiter http://controleduction.group.shef.c.uk/indexwebbook.html http://www.shef.c.uk/cse Deprtment of utomtic Control nd Systems Engineering

More information

The Perron-Frobenius operators, invariant measures and representations of the Cuntz-Krieger algebras

The Perron-Frobenius operators, invariant measures and representations of the Cuntz-Krieger algebras The Perron-Frobenius opertors, invrint mesures nd representtions of the Cuntz-Krieger lgebrs Ktsunori Kwmur Reserch Institute for Mthemticl Sciences Kyoto University, Kyoto 606-8502, Jpn For trnsformtion

More information

Lecture 11 Binary Decision Diagrams (BDDs)

Lecture 11 Binary Decision Diagrams (BDDs) C 474A/57A Computer-Aie Logi Design Leture Binry Deision Digrms (BDDs) C 474/575 Susn Lyseky o 3 Boolen Logi untions Representtions untion n e represente in ierent wys ruth tle, eqution, K-mp, iruit, et

More information

Qubit and Quantum Gates

Qubit and Quantum Gates Quit nd Quntum Gtes Shool on Quntum omputing @Ygmi Dy, Lesson 9:-:, Mrh, 5 Eisuke Ae Deprtment of Applied Physis nd Physio-Informtis, nd REST-JST, Keio University From lssil to quntum Informtion is physil

More information

Design of Members. Shear Force. Example : Shear resistance of webs without and with stiffeners

Design of Members. Shear Force. Example : Shear resistance of webs without and with stiffeners TALAT Lecture 0 Design o Members Ser Force Exmple 6. 6.6 : Ser resistnce o ebs itout nd it stieners pges Advnced Level prepred by Torsten Höglund, Royl Institute o Tecnology, Stockolm Dte o Issue: 999

More information

MORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)

MORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.) MORE FUNCTION GRAPHING; OPTIMIZATION FRI, OCT 25, 203 (Lst edited October 28, 203 t :09pm.) Exercise. Let n be n rbitrry positive integer. Give n exmple of function with exctly n verticl symptotes. Give

More information

4.8 Arnoldi Iteration, Krylov Subspaces and GMRES

4.8 Arnoldi Iteration, Krylov Subspaces and GMRES 48 Arnoldi Iteration, Krylov Subspaces and GMRES We start with the problem of using a similarity transformation to convert an n n matrix A to upper Hessenberg form H, ie, A = QHQ, (30) with an appropriate

More information

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008 Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of

More information

Modified midpoint method for solving system of linear Fredholm integral equations of the second kind

Modified midpoint method for solving system of linear Fredholm integral equations of the second kind Americn Journl of Applied Mtemtics 04; (5: 55-6 Publised online eptember 30, 04 (ttp://www.sciencepublisinggroup.com/j/jm doi: 0.648/j.jm.04005. IN: 330-0043 (Print; IN: 330-006X (Online Modified midpoint

More information

0.1 Differentiation Rules

0.1 Differentiation Rules 0.1 Differentiation Rules From our previous work we ve seen tat it can be quite a task to calculate te erivative of an arbitrary function. Just working wit a secon-orer polynomial tings get pretty complicate

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

Joint Distribution of any Record Value and an Order Statistics

Joint Distribution of any Record Value and an Order Statistics Interntionl Mthemticl Forum, 4, 2009, no. 22, 09-03 Joint Distribution of ny Record Vlue nd n Order Sttistics Cihn Aksop Gzi University, Deprtment of Sttistics 06500 Teknikokullr, Ankr, Turkey entelpi@yhoo.com

More information

Power System Representation and Equations. A one-line diagram of a simple power system

Power System Representation and Equations. A one-line diagram of a simple power system Power ystem epresenttion nd Equtions Lod B Lod A Bus Bus A oneline digrm of simple power system Oil or liquid iruit reker otting mhine Twowinding power trnsformer Wye onnetion, neutrl ground PerPhse, Per

More information

Prerna Tower, Road No 2, Contractors Area, Bistupur, Jamshedpur , Tel (0657) ,

Prerna Tower, Road No 2, Contractors Area, Bistupur, Jamshedpur , Tel (0657) , R rern Tower, Rod No, Contrctors Are, Bistupur, Jmshedpur 800, Tel 065789, www.prernclsses.com IIT JEE 0 Mthemtics per I ART III SECTION I Single Correct Answer Type This section contins 0 multiple choice

More information

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems 5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we

More information

2011 Fermat Contest (Grade 11)

2011 Fermat Contest (Grade 11) Te CENTRE for EDUCATION in MATHEMATICS and COMPUTING 011 Fermat Contest (Grade 11) Tursday, February 4, 011 Solutions 010 Centre for Education in Matematics and Computing 011 Fermat Contest Solutions Page

More information

Chapter 5 FINITE DIFFERENCE METHOD (FDM)

Chapter 5 FINITE DIFFERENCE METHOD (FDM) MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential

More information

Key words. Numerical quadrature, piecewise polynomial, convergence rate, trapezoidal rule, midpoint rule, Simpson s rule, spectral accuracy.

Key words. Numerical quadrature, piecewise polynomial, convergence rate, trapezoidal rule, midpoint rule, Simpson s rule, spectral accuracy. O SPECTRA ACCURACY OF QUADRATURE FORMUAE BASED O PIECEWISE POYOMIA ITERPOATIO A KURGAOV AD S TSYKOV Abstrct It is well-known tt te trpezoidl rule, wile being only second-order ccurte in generl, improves

More information

Lecture 19: Continuous Least Squares Approximation

Lecture 19: Continuous Least Squares Approximation Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for

More information

Solving large scale eigenvalue problems

Solving large scale eigenvalue problems arge scale eigenvalue problems, Lecture 5, March 23, 2016 1/30 Lecture 5, March 23, 2016: The QR algorithm II http://people.inf.ethz.ch/arbenz/ewp/ Peter Arbenz Computer Science Department, ETH Zürich

More information

Online Appendix for Lerner Symmetry: A Modern Treatment

Online Appendix for Lerner Symmetry: A Modern Treatment Online Appendix or Lerner Symmetry: A Modern Treatment Arnaud Costinot MIT Iván Werning MIT May 2018 Abstract Tis Appendix provides te proos o Teorem 1, Teorem 2, and Proposition 1. 1 Perect Competition

More information

Designing Information Devices and Systems I Fall 2016 Babak Ayazifar, Vladimir Stojanovic Homework 6. This homework is due October 11, 2016, at Noon.

Designing Information Devices and Systems I Fall 2016 Babak Ayazifar, Vladimir Stojanovic Homework 6. This homework is due October 11, 2016, at Noon. EECS 16A Designing Informtion Devices nd Systems I Fll 2016 Bk Ayzifr, Vldimir Stojnovic Homework 6 This homework is due Octoer 11, 2016, t Noon. 1. Homework process nd study group Who else did you work

More information

Problem Set 3

Problem Set 3 14.102 Problem Set 3 Due Tuesdy, October 18, in clss 1. Lecture Notes Exercise 208: Find R b log(t)dt,where0

More information

Lecture Notes Di erentiating Trigonometric Functions page 1

Lecture Notes Di erentiating Trigonometric Functions page 1 Lecture Notes Di erentiating Trigonometric Functions age (sin ) 7 sin () sin 8 cos 3 (tan ) sec tan + 9 tan + 4 (cot ) csc cot 0 cot + 5 sin (sec ) cos sec tan sec jj 6 (csc ) sin csc cot csc jj c Hiegkuti,

More information

CONNECTEDNESS OF SELF-AFFINE SETS WITH PRODUCT DIGIT SETS. 1. Introduction. T = d DA 1 (T +d). (1.1)

CONNECTEDNESS OF SELF-AFFINE SETS WITH PRODUCT DIGIT SETS. 1. Introduction. T = d DA 1 (T +d). (1.1) CONNECTEDNESS OF SELF-AFFINE SETS WITH PRODUCT DIGIT SETS JING-CHENG LIU 1, JUN JASON LUO 2 AND KE TANG 1 rxiv:1603.06087v1 [mth.gn] 19 Mr 2016 Abstrct. Let T(A, D) be self-ffine set generted by n expnding

More information

Tutorial 4. b a. h(f) = a b a ln 1. b a dx = ln(b a) nats = log(b a) bits. = ln λ + 1 nats. = log e λ bits. = ln 1 2 ln λ + 1. nats. = ln 2e. bits.

Tutorial 4. b a. h(f) = a b a ln 1. b a dx = ln(b a) nats = log(b a) bits. = ln λ + 1 nats. = log e λ bits. = ln 1 2 ln λ + 1. nats. = ln 2e. bits. Tutoril 4 Exercises on Differentil Entropy. Evlute the differentil entropy h(x) f ln f for the following: () The uniform distribution, f(x) b. (b) The exponentil density, f(x) λe λx, x 0. (c) The Lplce

More information

Math Week 5 concepts and homework, due Friday February 10

Math Week 5 concepts and homework, due Friday February 10 Mt 2280-00 Week 5 concepts nd omework, due Fridy Februry 0 Recll tt ll problems re good for seeing if you cn work wit te underlying concepts; tt te underlined problems re to be nded in; nd tt te Fridy

More information

INTERNATIONAL CENTRE FOR THEORETICAL PHYSICS THE ALGEBRAIC APPROACH TO THE SCATTERING PROBLEM ABSTRACT

INTERNATIONAL CENTRE FOR THEORETICAL PHYSICS THE ALGEBRAIC APPROACH TO THE SCATTERING PROBLEM ABSTRACT IC/69/7 INTERNAL REPORT (Limited distribution) INTERNATIONAL ATOMIC ENERGY AGENCY INTERNATIONAL CENTRE FOR THEORETICAL PHYSICS THE ALGEBRAIC APPROACH TO THE SCATTERING PROBLEM Lot. IXARQ * Institute of

More information

Differentiation Rules and Formulas

Differentiation Rules and Formulas Differentiation Rules an Formulas Professor D. Olles December 1, 01 1 Te Definition of te Derivative Consier a function y = f(x) tat is continuous on te interval a, b]. Ten, te slope of te secant line

More information

308K. 1 Section 3.2. Zelaya Eufemia. 1. Example 1: Multiplication of Matrices: X Y Z R S R S X Y Z. By associativity we have to choices:

308K. 1 Section 3.2. Zelaya Eufemia. 1. Example 1: Multiplication of Matrices: X Y Z R S R S X Y Z. By associativity we have to choices: 8K Zely Eufemi Section 2 Exmple : Multipliction of Mtrices: X Y Z T c e d f 2 R S X Y Z 2 c e d f 2 R S 2 By ssocitivity we hve to choices: OR: X Y Z R S cr ds er fs X cy ez X dy fz 2 R S 2 Suggestion

More information

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x

More information