Orthogonal Function Solution of Differential Equations

Similar documents
Orthogonal Functions

Lecture #3. Math tools covered today

Perturbation Theory, Zeeman Effect, Stark Effect

AP Calculus BC Review Chapter 12 (Sequences and Series), Part Two. n n th derivative of f at x = 5 is given by = x = approximates ( 6)

Summer MA Lesson 13 Section 1.6, Section 1.7 (part 1)

Indefinite Integral. Lecture 21 discussed antiderivatives. In this section, we introduce new notation and vocabulary. The notation f x dx

X. Perturbation Theory

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES

8.3 Perturbation theory

Null Spaces, Column. Transformations. Remarks. Definition ( ) ( ) Null A is the set of all x's in R that go to 0 in R 3

Chapter 4. Fourier Series

Apply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.

Uncertainty Principle of Mathematics

THE CLOSED FORMS OF CONVERGENT INFINITE SERIES ESTIMATION OF THE SERIES SUM OF NON-CLOSED FORM ALTERNATING SERIES TO A HIGH DEGREE OF PRECISION.

AVERAGE MARKS SCALING

Binomial transform of products

Application of Homotopy Analysis Method for Solving various types of Problems of Ordinary Differential Equations

k=1 s k (x) (3) and that the corresponding infinite series may also converge; moreover, if it converges, then it defines a function S through its sum

September 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1

19.1 The dictionary problem

Lecture 10: Bounded Linear Operators and Orthogonality in Hilbert Spaces

and then substitute this into the second equation to get 5(11 4 y) 3y

) is a square matrix with the property that for any m n matrix A, the product AI equals A. The identity matrix has a ii

6.4 Binomial Coefficients

A Generalization of Ince s Equation

The Differential Transform Method for Solving Volterra s Population Model

3. Calculus with distributions

Lecture 19. Curve fitting I. 1 Introduction. 2 Fitting a constant to measured data

A New Type of q-szász-mirakjan Operators

THE LEGENDRE POLYNOMIALS AND THEIR PROPERTIES. r If one now thinks of obtaining the potential of a distributed mass, the solution becomes-

Chapter 2. Asymptotic Notation

Lommel Polynomials. Dr. Klaus Braun taken from [GWa] source. , defined by ( [GWa] 9-6)

Name Period ALGEBRA II Chapter 1B and 2A Notes Solving Inequalities and Absolute Value / Numbers and Functions

SphericalHarmonicY. Notations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation

Chapter 4 Postulates & General Principles of Quantum Mechanics

The structure of Fourier series

x+ 2 + c p () x c p () x is an arbitrary function. ( sin x ) dx p f() d d f() dx = x dx p cosx = cos x+ 2 d p () x + x-a r (1.

COMP 2804 Solutions Assignment 1

x !1! + 1!2!

Matrix Operators and Functions Thereof

A string of not-so-obvious statements about correlation in the data. (This refers to the mechanical calculation of correlation in the data.

Topic 1 2: Sequences and Series. A sequence is an ordered list of numbers, e.g. 1, 2, 4, 8, 16, or

Inhomogeneous Poisson process

Representing Functions as Power Series. 3 n ...

Jacobi symbols. p 1. Note: The Jacobi symbol does not necessarily distinguish between quadratic residues and nonresidues. That is, we could have ( a

The Non-homogeneous Diffusion Equation

CALCULUS BASIC SUMMER REVIEW

New method for evaluating integrals involving orthogonal polynomials: Laguerre polynomial and Bessel function example

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t

Bernoulli Polynomials Talks given at LSBU, October and November 2015 Tony Forbes

Research Article Sums of Products of Cauchy Numbers, Including Poly-Cauchy Numbers

Definition 2.1 (The Derivative) (page 54) is a function. The derivative of a function f with respect to x, represented by. f ', is defined by

Lecture 11. Solution of Nonlinear Equations - III

Integrals of Functions of Several Variables

5.6 Binomial Multi-section Matching Transformer

Name Date Class. Think: Use the Quotient Property. Rationalize the denominator. Use the Product Property.

Sigma notation. 2.1 Introduction

Machine Learning for Data Science (CS 4786)

Numerical Solution of Nonlinear Integro-Differential Equations with Initial Conditions by Bernstein Operational Matrix of Derivative

! " * (x,t) " (x,t) dx =! #(x,t) dx = 1 all space

Mathematical Series (You Should Know)

ENGINEERING COMPUTATION Lecture 5. Approximate representation of data and functions

( a) ( ) 1 ( ) 2 ( ) ( ) 3 3 ( ) =!

Automated Proofs for Some Stirling Number Identities

MATHEMATICS. 61. The differential equation representing the family of curves where c is a positive parameter, is of

Fourier Series and the Wave Equation

SOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.

INTEGRATION BY PARTS (TABLE METHOD)

An Insight into Differentiation and Integration

Complete Solutions to Supplementary Exercises on Infinite Series

ECE 901 Lecture 4: Estimation of Lipschitz smooth functions

A generalization of the Leibniz rule for derivatives

Question 1: The magnetic case

Generating Functions and Their Applications

Data Analysis and Statistical Methods Statistics 651

Math 216A Notes, Week 3

Contents Two Sample t Tests Two Sample t Tests

Properties and Tests of Zeros of Polynomial Functions

42 Dependence and Bases

Math 4707 Spring 2018 (Darij Grinberg): midterm 2 page 1. Math 4707 Spring 2018 (Darij Grinberg): midterm 2 with solutions [preliminary version]

(3) If you replace row i of A by its sum with a multiple of another row, then the determinant is unchanged! Expand across the i th row:

d dx where k is a spring constant

Seunghee Ye Ma 8: Week 5 Oct 28

q-fibonacci polynomials and q-catalan numbers Johann Cigler [ ] (4) I don t know who has observed this well-known fact for the first time.

Acoustic Field inside a Rigid Cylinder with a Point Source

5.6 Binomial Multi-section Matching Transformer

CHAPTER 10 INFINITE SEQUENCES AND SERIES

GROUPOID CARDINALITY AND EGYPTIAN FRACTIONS

Answer Key, Problem Set 1, Written

MA131 - Analysis 1. Workbook 3 Sequences II

A PROBABILITY PROBLEM

Ma/CS 6a Class 22: Power Series

FUZZY RELIABILITY ANALYSIS OF COMPOUND SYSTEM BASED ON WEIBULL DISTRIBUTION

We have also learned that, thanks to the Central Limit Theorem and the Law of Large Numbers,

Linear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy

Proof of two divisibility properties of binomial coefficients conjectured by Z.-W. Sun

Algebra II Notes Unit Seven: Powers, Roots, and Radicals

Taylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH

6.451 Principles of Digital Communication II Wednesday, March 9, 2005 MIT, Spring 2005 Handout #12. Problem Set 5 Solutions

A widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α

Transcription:

Royal Holloway Uiversity of Loo Departet of Physics Orthogoal Fuctio Solutio of Differetial Equatios trouctio A give oriary ifferetial equatio will have solutios i ters of its ow fuctios Thus, for eaple, the solutio of the SHO Schröiger equatio is epresse i ters of Herite polyoials ultiplyig a gausia: SHO Schröiger equatio H 16 e / We are talkig here about the eigefuctios of the equatio a their eigevalues ( i this case) Now we are goig to look ito usig the wrog fuctios to solve a ifferetial equatio Why ight we wat to o this? There are a uber of reasos 16 / We coul use the set of H e like a quatu SHO to see how uch the behaviour of the syste was Siilarly, we coul use the set e i0k ωt5 to see how uch the behaviour was like that of free travelig waves These proceures coul be iportat i their ow right, or as part of a perturbatio or approiatio schee The practical cosequece of this proceure is that ifferetial equatios are trasfore ito atri equatios (A of course this is at the heart of the equivalece of the Schröiger a the Heiseberg approaches to quatu echaics The fuaetal ieas The fuctios which arise fro the solutio of a ifferetial equatio ofte have soe iportat properties The equatios which arise i physical probles ca be assue to coprise a coplete orthogoal set of fuctios (Forally, this result follows fro the Stur-Liouville theore) Orthogoal We kow what orthogoal eas For the set of fuctios orthogoality itegral 16 16 (it ight also have a weight fuctio) < 16A there is a PH130 Matheatical Methos 1

Royal Holloway Uiversity of Loo Departet of Physics The itegral ust be evaluate over the correct iterval a it is zero uless = f the fuctios are oralise the the itegral is uity whe = Coplete f the set of fuctios is coplete, this eas that ay fuctio f16ca be epresse as a liear su of the fuctios f = f This epressio for f 16 16 1 6 will be vali over the orthogoality iterval You shoul reeber how to fi the = coefficiets f : f16 f 16 16 The eoiator is eee whe the fuctios are ot oralise Differetiatio A iportat a highly relevat eaple of the copleteess property is the ifferetiatio of a fuctio Startig fro a fuctio f16, if we ifferetiate it we obtai a ew fuctio f / Now copleteess tells us that ay fuctio f16 ca be epresse as a liear su of the basis fuctios 16 But sice f / is just aother fuctio, this too ca be epresse as a liear su of the basis fuctios We have 16= 16 f f a we siilarly epect that f = 16 g The questio, the, is how to fi the coefficiets g for the erivative Give the basis set < 16Athe fuctio f16 is specifie by the coefficiets f We wat to fi the coefficiets g which correspoigly specify the erivative f / Let s procee i the followig way Start fro f16: a ifferetiate it 16= 16 f f PH130 Matheatical Methos

Royal Holloway Uiversity of Loo Departet of Physics f = 16 f Copleteess assures us that each of the 16/ will be a liear fuctio of the basis fuctios For eaple, the erivative of Herite polyoials is epresse as H H 16= 1 16 geeral we will have 16 = 16 D The coefficiets D ee two iices, oe to iicate which basis fuctio we are ifferetiatig (), a the other to iicate how uch of each of the basis fuctios we ee i the erivative () For the Herite polyoial case we the have D δ, Check that you uersta this = + 1 Usig the above result we ca retur to the task of ifferetiatig our fuctio f f = = f D 16 16 16 Df = f We frae our questio by sayig that if we write the erivative as the epasio f = 16 g the what are the coefficiets g? We ow have the aswer: g = D f 16 You shoul observe that this has the for of a atri prouct f we represet g a f as colu atrices the g D D D : f g g : 1 3 = 11 1 13 D1 D : : f D31 : : : f3 : : : : : 1 PH130 Matheatical Methos 3

Royal Holloway Uiversity of Loo Departet of Physics Notatio a 3 vector aalogy With 3 vectors, oce oe has a give basis set = i, j, k B the ay vector ca be epresse as a liear su of these v = vi+ v i jj+ vkk Oce the basis set is ecie upo the ay physical vector v is specifie by givig its coefficiets v i, v j, a v k a these coefficiets ca be asseble ito a colu atri vi v j vk Matheaticias ight call such a colu of ubers a vector We ay eve occasioally succub to this ieactitue Strictly speakig, this colu of ubers is a particular represetatio of the physical vector v if we chage the basis set i, j, k = B the we ee a ifferet colu of coefficiets to represet the sae ol vector v Now cosier fuctios With a give basis set < 16 Aay fuctio f16 ca be epresse as a liear su of these f16 = 16 f Oce the basis set is ecie upo the ay fuctio f16 is specifie by givig its coefficiets;@a f these ubers ca be asseble ito a colu atri f1 f f This colu of coefficiets (a atheaticia s vector) is a particular represetatio of the fuctio f16 f we chage the basis set < 16 Awe were usig the we woul ee a ifferet colu of coefficiets to represet the sae ol fuctio Takig the otatioal aalogy further we coul eote the colu of coefficiets by a vector sybol f That is, we coul ake the associatio f 16~ f As a eaple of usig this otatioal aalogy cosier the ifferetiatio result of before g = D f 3 : Here the f are the coefficiets of the fuctio f16a the g are the coefficiets of its erivative f we eote these colus of coefficiets by the vectors f a f the the ifferetiatio result becoes the vector-atri equatio f = Df PH130 Matheatical Methos 4

Royal Holloway Uiversity of Loo Departet of Physics where D represets the atri whose eleets are D This shows us that ultiplyig with the atri D has the effect of ifferetiatio Matri represetatio of ifferetial equatios The above aalogies suggest that ifferetial equatios ay be represete by atri equatios This ay be oe forally i the followig way Startig fro the siple eaple 16 16 g = f, i ters of the basis set < 16 Athis ay be epresse as 16 g 16 = Df (Reeber that the suatio iices are copletely arbitrary) The trick to tur this ito a atri equatio is to ultiply both sies by oe of the basis fuctios a the perfor a orthogoality itegral Let s ultiply by p 16a itegrate O the left ha sie as varies through its rage of values the orthogoality itegral of 16with p 16will always be zero ecept whe = p Thus we pick out just the p ter fro the su O the right ha sie, as varies through its rage we will siilarly pick out the p ter We the e up with the equatio g = D f, p p or g= Df which is the atri represetatio of the origial equatio Differetial equatios with costat coefficiets Hoogeeous equatios Let s cosier a seco orer (hoogeeous) ifferetial equatio with costat coefficiets f16 A B f 1 6 + + Cf 0 16= This is trasfore ito a atri equatio i four steps: 1 Epress f as a liear su of the basis fuctios: f = f 16 16 PH130 Matheatical Methos 5

Royal Holloway Uiversity of Loo Departet of Physics Substitute this ito the ifferetial equatio: 16 16 A f + B f C f + 16 = 0 3 Multiply by p A & ' p 1 6 a o the orthogoality itegral: 16 16 ( f B p16 16 ) f C > + p16 16 Cf * &' ( )* + = 0 4 etify the atri for of each ter By copariso with the erivatio of the ifferetiatio atri D above, the ifferetial equatio ay be writte i atri for as ADf+ BDf+ 1f=0 or AD + BD+ 1 f =0 where 1 is the uit atri < A Now AD + BD+ 1 is just aother atri, which we ca eote by M Thus our ifferetial equatio has bee trasfore, with the basis fuctio set < 16A, ito the atri equatio Mf = 0 The hoogeeous ifferetial equatio has bee trasfore ito a hoogeeous atri equatio This woul the be solve by the usual proceures for solvig hoogeeous atri equatios hoogeeous equatios Let us ow cosier a ihoogeeous equatio such as f A + B f + Cf = s where s ust be epresse as liear sus of the basis fuctios: f = f 16 16 1 6 1 6 1 6 is the source ter, akig the equatio ihoogeeous Now both f a s 16 16 s 16= 16 s We use the vector otatio aalogy to eote the colu of coefficiets of f16 by the vector sybol f a siilarly the coefficiets of s 16by s 16 f s 16 ~ f ~ s The, usig the sae proceure as i the previous case, we arrive at the ihoogeeous atri equatio PH130 Matheatical Methos 6

Royal Holloway Uiversity of Loo Departet of Physics Mf = s The ihoogeeous ifferetial equatio has bee trasfore ito a ihoogeeous atri equatio This woul the be solve by the usual proceures for solvig ihoogeeous atri equatios The foral solutio is particularly straightforwar: we ultiply both equatios, fro the left, by the iverse of M The we have 1 1 M Mf = M s or f = M 1 s This solves the proble; f16is give i ters of s Differetial equatios with variable coefficiets Whe the coefficiets of the ifferetial equatio are variable, that is whe we have a equatio of the for f f A 16 1 6 B 16 1 6 + + C 1616 f = 0 thigs are slightly ore coplicate We procee to trasfor this ito a atri equatio usig the sae four steps as before: 1 Epress f as a liear su of the basis fuctios: f = f 16 16 Substitute this ito the ifferetial equatio: A 16 1 6 f B 16 1 6 + f C 16 16 f + = 0 16 3 Multiply by p 16a o the orthogoality itegral: p1a 6 1 6 1 6 f p16b16 1 6 & ( f ' ) + * &' ( )* + + > p1616 C16 Cf = 0 4 etify the atri for of each ter this case the we o t have a irect coectio with the ifferetiatio atri D we ha before Nevertheless this still has the for of a atri equatio Mf = 0 where the eleets of the atri M are ow Mp = p1616 1 6 A + p1616 1 6 B + p1616 C 16 PH130 Matheatical Methos 7