ECE Spring Prof. David R. Jackson ECE Dept. Notes 20

Similar documents
ECE Spring Prof. David R. Jackson ECE Dept. Notes 8

Notes 19 Bessel Functions

Chapter 4. Fourier Series

BESSEL EQUATION and BESSEL FUNCTIONS

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0,

A second look at separation of variables

Notes 8 Singularities

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES

De Moivre s Theorem - ALL

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

Different kinds of Mathematical Induction

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

x !1! + 1!2!

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

CHAPTER 10 INFINITE SEQUENCES AND SERIES

1. Hydrogen Atom: 3p State

Chapter 9: Numerical Differentiation

Eigenvalues and Eigenvectors

Ma 530 Infinite Series I

Chapter Vectors

UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014

PHYC - 505: Statistical Mechanics Homework Assignment 4 Solutions

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

Integrals of Functions of Several Variables

Problem 4: Evaluate ( k ) by negating (actually un-negating) its upper index. Binomial coefficient

The Binomial Multi- Section Transformer

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

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

Orthogonal Functions

6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.

CHAPTER 5. Theory and Solution Using Matrix Techniques

The Binomial Multi-Section Transformer

PHYSICS 116A Homework 2 Solutions

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

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

5.6 Binomial Multi-section Matching Transformer

5.6 Binomial Multi-section Matching Transformer

Appendix F: Complex Numbers

MIDTERM 3 CALCULUS 2. Monday, December 3, :15 PM to 6:45 PM. Name PRACTICE EXAM SOLUTIONS

Notes 12 Asymptotic Series

42 Dependence and Bases

RADICAL EXPRESSION. If a and x are real numbers and n is a positive integer, then x is an. n th root theorems: Example 1 Simplify

For use only in [the name of your school] 2014 FP2 Note. FP2 Notes (Edexcel)

Chapter 1. Complex Numbers. Dr. Pulak Sahoo

6.4 Binomial Coefficients

JEE ADVANCED 2013 PAPER 1 MATHEMATICS

Presentation of complex number in Cartesian and polar coordinate system

Lecture 2 Appendix B: Some sample problems from Boas, Chapter 1. Solution: We want to use the general expression for the form of a geometric series

Inverse Matrix. A meaning that matrix B is an inverse of matrix A.

Created by T. Madas SERIES. Created by T. Madas

3sin A 1 2sin B. 3π x is a solution. 1. If A and B are acute positive angles satisfying the equation 3sin A 2sin B 1 and 3sin 2A 2sin 2B 0, then A 2B

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

In algebra one spends much time finding common denominators and thus simplifying rational expressions. For example:

Physics 219 Summary of linear response theory

AVERAGE MARKS SCALING

Application 10.5D Spherical Harmonic Waves

X. Perturbation Theory

Math 12 Final Exam, May 11, 2011 ANSWER KEY. 2sinh(2x) = lim. 1 x. lim e. x ln. = e. (x+1)(1) x(1) (x+1) 2. (2secθ) 5 2sec2 θ dθ.

SEQUENCE AND SERIES NCERT

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

ECE 901 Lecture 4: Estimation of Lipschitz smooth functions

Solutions to Problem Set 8

Chapter 8. Euler s Gamma function

(A sequence also can be thought of as the list of function values attained for a function f :ℵ X, where f (n) = x n for n 1.) x 1 x N +k x N +4 x 3

8.3 Perturbation theory

Math 155 (Lecture 3)

Recursive Algorithms. Recurrences. Recursive Algorithms Analysis

Evaluation of Some Non-trivial Integrals from Finite Products and Sums

Introduction to Extreme Value Theory Laurens de Haan, ISM Japan, Erasmus University Rotterdam, NL University of Lisbon, PT

The Non-homogeneous Diffusion Equation

Subject: Differential Equations & Mathematical Modeling-III

[ 11 ] z of degree 2 as both degree 2 each. The degree of a polynomial in n variables is the maximum of the degrees of its terms.

18.01 Calculus Jason Starr Fall 2005

Sequences and Limits

ON SUPERSINGULAR ELLIPTIC CURVES AND HYPERGEOMETRIC FUNCTIONS

MATH 6101 Fall 2008 Newton and Differential Equations

Sequences of Definite Integrals, Factorials and Double Factorials

Numerical Methods in Fourier Series Applications

Introduction to Astrophysics Tutorial 2: Polytropic Models


Random Models. Tusheng Zhang. February 14, 2013

Linearly Independent Sets, Bases. Review. Remarks. A set of vectors,,, in a vector space is said to be linearly independent if the vector equation

( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.

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

Solutions to Final Exam Review Problems

Ma/CS 6a Class 22: Power Series

SEQUENCES AND SERIES

Linear Regression Demystified

Bertrand s postulate Chapter 2

CSE 1400 Applied Discrete Mathematics Number Theory and Proofs

THE KENNESAW STATE UNIVERSITY HIGH SCHOOL MATHEMATICS COMPETITION PART II Calculators are NOT permitted Time allowed: 2 hours

Math 113 Exam 3 Practice

Complex Numbers Solutions

Expected Number of Level Crossings of Legendre Polynomials

The natural exponential function

Enumerative & Asymptotic Combinatorics

Math 475, Problem Set #12: Answers

Chapter 2 The Solution of Numerical Algebraic and Transcendental Equations

Transcription:

ECE 6341 Sprig 016 Prof. David R. Jackso ECE Dept. Notes 0 1

Spherical Wave Fuctios Cosider solvig ψ + k ψ = 0 i spherical coordiates z φ θ r y x

Spherical Wave Fuctios (cot.) I spherical coordiates we have ψ = ψ ψ ψ ψ r siθ 1 1 1 = + + r r r r siθ θ θ r si θ φ Hece we have 1 1 1 r ψ ψ ψ siθ k 0 (1) + ψ + + = r r r r siθ θ θ r si θ φ Usig separatio of variables, let ( ) ( ) ψ = Rr ( ) H θ Φ φ () 3

Spherical Wave Fuctios (cot.) After substitutig Eq. () ito Eq. (1), divide by ψ : 1 R 1 H r siθ + rr r r rsiθh θ θ 1 1 Φ + + k = r si θ Φ φ 0 At this poit, we caot yet say that all of the depedece o ay give variable is oly withi oe ter. 4

Spherical Wave Fuctios (cot.) Next, ultiply by r si θ : si θ R siθ H r + siθ R r r H θ θ 1 Φ + + kr θ = Φ φ si 0 (3) Sice the uderlied ter is the oly oe which depeds o, It ust be equal to a costat, φ Hece, set 1 Φ φ Φ = (4) 5

Spherical Wave Fuctios (cot.) Hece, cos φ Φ ( φ ) = (5) si φ I geeral, w (ot a iteger). Now divide Eq. (3) by si θ ad use Eq. (4), to obtai 1 R 1 H r + siθ R r r Hsiθ θ θ si + kr = θ 0 (6) 6

Spherical Wave Fuctios (cot.) The uderlied ters are the oly oes that ivolve θ ow. This tie, the separatio costat is custoarily chose as ( 1) + 1 d siθ dh = ( + 1 ) (7) Hsiθ dθ dθ si θ I geeral, (ot a iteger) To siplify this, let x = cosθ dx = siθ dθ d d = siθ dθ dx ad deote yx ( ) = H( θ ) 7

Spherical Wave Fuctios (cot.) Usig x = cosθ d d siθ dθ = dx ( ) 1/ siθ = 1 x 1 d siθ dh = + 1 Hsiθ dθ dθ si θ ( ) y 1 ( 1 x ) 1/ ( ) 1/ d ( ) 1/ ( )( ) 1/ 1 x 1 x 1 1 x y dx + ( + 1) = 0 ( 1 x ) 8

Spherical Wave Fuctios (cot.) Cacelig ters, y 1 ( 1 x ) 1/ ( ) 1/ d ( ) 1/ ( )( ) 1/ 1 x 1 x 1 1 x y dx + ( + 1) = 0 ( 1 x ) Multiplyig by y, we have d dx ( ) 1 x y + ( + 1) y= 0 (8) ( 1 x ) 9

Spherical Wave Fuctios (cot.) Eq. (8) is the associated Legedre equatio. The solutios are represeted as yx ( ) P = Q ( x) ( x) Associated Legedre fuctio of the first kid. Associated Legedre fuctio of the secod kid. = order, = degree If = 0, Eq. (8) is called the Legedre equatio, i which case yx ( ) 0 P ( x) P( x) = = 0 Q ( ) Q ( x) x Legedre fuctio of the first kid. Legedre fuctio of the secod kid. 10

Spherical Wave Fuctios (cot.) Hece: H ( θ ) P = Q (cos θ ) (cos θ ) To be as geeral as possible: w H ( θ ) P = Q w w (cos θ ) (cos θ ) 11

Spherical Wave Fuctios (cot.) Substitutig Eq. (7) ito Eq. (6) ow yields 1 d R dr r + + kr = dr dr ( ) 1 0 Next, let x = kr dx d dr = = k dr d k dx ad deote yx () = Rr () 1

Spherical Wave Fuctios (cot.) We the have 1 d R dr r + + kr = dr dr ( ) 1 0 d x = kr = k dr d dx k k ( + ) + x = 1 d x dy y dx k dx 1 0 d ( xy ) + x ( + 1) y = 0 dx 13

Spherical Wave Fuctios (cot.) or d ( xy ) + x ( + 1) y = 0 dx + + ( + 1) = 0 x y xy x y spherical Bessel equatio Solutio: b (x) Note the lower case b. 14

Spherical Wave Fuctios (cot.) ad let Hece ( ) Deote = 1/ b x x g x ( ) yx ( ) = b x ( ) 1/ 1 3/ b ( x) = x g x g 1 1 3 b ( x) = x g x g x g + x g 4 1/ 3/ 3/ 5/ 3 x g x g + x g + x g x g 4 ( ) 3/ 1/ 1/ 1/ 1/ ( ) 1/ + x + 1 x g( x) = 0 15

Spherical Wave Fuctios (cot.) Multiply by 1/ x 3 x g xg + g + ( xg g) + x ( + 1 ) g( x) = 0 4 Cobie these ters or 1 xg + xg g + xg ( + 1) g = 0 4 Use 1 1 + ( + 1) = + 4 Cobie these ters Defie γ + 1 16

Spherical Wave Fuctios (cot.) ( ) We the have γ x g + xg + x g = This is Bessel s equatio of order γ. 0 Hece g Jγ ( x) = Yγ ( x) so that b 1 J+ 1/( x) π = x Y + 1/( x) added for coveiece 17

Spherical Wave Fuctios (cot.) Defie π j( x) J+ 1/( x) x π y( x) Y+ 1/( x) x The j ( ) () ( ) kr R r = b kr = y ( ) kr 18

Suary ( ) ( ) ψ + k ψ = 0 ( φ ) ( φ ) j kr P (cos ) cos θ ψ = y kr Q (cos ) si θ π b( x) = B+ 1/( x) x I geeral, w 19

Properties of Spherical Bessel Fuctios = π b( x) = B+ 1/( x) x Bessel fuctios of half-iteger order are give by closed-for expressios. J ( x) k + k ( 1) x! ( 1) k = 0 k! ( + k)! z =Γ z+ = This becoes a closed-for expressio! = + 1/ 0

Properties of Spherical Bessel Fuctios (cot.) Exaples: J1/( x) = si x π x ( ) J 1/( x) = cos x π x ( ) ( x) si J3/ ( x) = cos x π x x ( x) ( ) cos J 3/ ( x) = + si x π x x ( ) J( x)cos( π ) J ( x) Y ( x) si( π ) 1/ 1/ ( ) Y ( x) = J x 3/ 3/ ( ) Y ( x) = J x 1

Properties of Spherical Bessel Fuctios (cot.) Start with: k ( ) ( ) Proof for = 1/ J1/( x) = si x π x ( ) k ( ) ( ) + k 1/+ k 1 x 1 x J ( x) = J1/( x) = k= 0 k! + k! k= 0 k! 1/ + k! Hece k ( 1) ( k) x x J1/( x) = k = 0 k! 1/ +! k + 1

Properties of Spherical Bessel Fuctios (cot.) Exaie the factorial expressio: Note: x = x( x ) ( ) =! 1! 1/! π / 1 + 1/! = + 1/ 1/ 3/... 3/! ( k ) ( k )( k )( k ) ( ) ( k 1/)( k 1/)( k 3/ )...( 3/ )( π /) = + 1 = + = = = ( k 1)( k 1)( k 3 )...( 3 )( π / ) + 1 4...4 ( k 1 )!( π /) ( k)( k )( k ) 1 + k 1... 1 ( k 1 )!( π /) ( k)( k )( k ) 1 + k k! ( k 1 )!( π /) k k k k 3

Properties of Spherical Bessel Fuctios (cot.) Hece ( 1) x x J1/( x) = k k = 0 ( ) ( ) 1 k + 1! π / k! k k! k k + 1 Hece, we have k ( 1) ( k + ) π x x J1/( x) = k = 0 1! k + 1 k + 1 4

Properties of Spherical Bessel Fuctios (cot.) k ( 1) ( k + ) π x x J1/( x) = k = 0 1! k + 1 k + 1 or k = 0 k ( 1) ( k + ) π x J 1/( x ) = x 1! k + 1 We the recogize that π x J x x ( ) si 1/ = ( ) 5

Properties of Legedre Fuctios Relatio to Legedre fuctios (whe w = = iteger): ( ) / d P ( x) = 1 x P( x) dx ( ) / d Q ( x) = 1 x Q( x) dx These also hold for. For w (ot a iteger) the associated Legedre fuctio is defied i ters of the hypergeoetric fuctio. 6

Properties of Legedre Fuctios (cot.) Rodriguez s forula (for = ): 1 d P x x! dx ( ) ( = 1) Legedre polyoial (a polyoial of order ) ( ) ( ) 0 P0 x = x 1 = 1 1 d P1 ( x) = ( x 1) = x dx 1 d 1 P x = x 1 = 3x 1 8 dx ( ) ( ) ( ) 7

Properties of Legedre Fuctios (cot.) P x = > Note: ( ) 0, This follows fro these two relatios: ( ( ) 1 ) / d P x = x P( x) dx 1 d ( ) ( 1) P x = x = polyoial of order! dx 8

Properties of Legedre Fuctios (cot.) w=, = ( ( ) 1 ) / d Q x = x Q( x) dx Q ( ) x = ifiite series, ot a polyoial (ay blow up) ( 1) Q ± = (see ext slide) The Q fuctios all ted to ifiity as x ± 1 θ 0, π Recall : x = cosθ 9

Properties of Legedre Fuctios (cot.) Lowest-order Q fuctios: Q 0 ( x) 1 1+ x = l 1 x Q Q 1 ( x) ( x) x 1+ x = l 1 1 x 3 x 1 1+ 3 l x = x 4 1 x 30

Properties of Legedre Fuctios (cot.) Negative idex idetities: ( ) ( ) P x P x ( + 1) = (This idetity also holds for.) ( ) π ( ) ( ) ( ) Q + x = P x + Q x ( 1) 1 31

Properties of Legedre Fuctios (cot.) P ( x) Q ( x) = (see Harrigto, Appedix E) ifiite series = ifiite series P ( x) ( ) ( ) ( ) ( ) ( ) ( ) ( ) N ( ) N 1 +! 1 x si π 1! +! 1 x = = 0!! π = + 1! N = largest iteger less tha or equal to. P P (1) = 1 ( 1) = Q ( x) ( ) cos( π ) ( ) ( π ) π P x P x = si Q ( ± 1) = Both are valid solutios, which are liearly idepedet for (see ext slide) 3

Properties of Legedre Fuctios (cot.) Proof that a valid solutio is P ( x) d dx The ( ) 1 x P ( x) + ( + 1) P ( x) 0 ( 1 x ) = Let t = x d ( 1 t )( 1 ) P ( t) + ( + 1) P ( t) 0 = dt ( 1 t ) d dx d = dt or (t x) d dx ( ) 1 x P ( x) + ( + 1) P ( x) 0 ( 1 x ) = Hece, a valid solutio is P ( x) 33

Properties of Legedre Fuctios (cot.) P ( x) ad P ( x) are two liearly idepedet solutios. Valid idepedet solutios: P Q (cos θ ) (cos θ ) or P P (cos θ ) ( cos θ ) We have a choice which set we wish to use. 34

Properties of Legedre Fuctios (cot.) = P ( x) = ( 1) P ( x) (They are liearly depedet.) I this case we ust use P Q (cos θ ) (cos θ ) 35

Properties of Legedre Fuctios (cot.) Suary of z-axis properties (x = cos (θ )) = P 1 = 1 P 1 = 1 ( ) ( ) Q ( ± 1) = Q ( ± 1) = P 1 = 1 P 1 = ( ) ( ) ( ) 36

Properties of Legedre Fuctios (cot.) z P ( ) x allowed P ( x) allowed y x x = cosθ Q ( ) ( ) x ad Q x are ot allowed o ± zaxis. 37

Properties of Legedre Fuctios (cot.) Outside or iside sphere z z x y Iside hollow coe y Oly P (x) is allowed x Both P (x) ad P (x) are allowed 38

Properties of Legedre Fuctios (cot.) z Outside coe x y Oly P (x) is allowed 39

Properties of Legedre Fuctios (cot.) z x y Oly P (x) is allowed Iside iverted coe Note: The physics is the sae as for the upright coe, but the atheatical for of the solutio is differet! 40

Properties of Legedre Fuctios (cot.) z y x Both P (x) ad P (x) are allowed Outside iverted coe 41

Properties of Legedre Fuctios (cot.) z Outside bicoe x y P (x) ad P (x) are allowed Q (x) ad Q (x) are allowed 4