Technical Notes

Technical Notes

Technical Notes Converting Erdélyi [1] to Bracewell [2, 3] Erdélyi defined his Hankel transform of order v as [1]: ∞ 1 gy(;)υ = fxJ () xyxy2 dx ∫0 υ ( )( ) In this project, we just consider all the Hankel transforms of order 0. We unify the Erdélyi version into Bracewell’s version [2, 3], which is more easily related to standard Fourier transform notation. Ø For f(x), we divided the functions by x12, and change x to 2π r . Ø For the Hankel transform, we divided the functions by 2π y12, and then change y to q directly. Classification For each function we made, we also marked its classification for the following criteria as true or false: 1) Oscillating 2) Finite domain 3) Singularity 4) Non negative 5) Monotonic 6) Eigen function You can use our search engine to find the functions which fit your own criteria. Notes 1) The range we chose for graphing each function is rqor= 2× 5n , nZ ∈, based on the extent of each function. For example, when n =−1 , r =25×−1 = 0.4; when n = 0 , r =25×0 = 2 So in our project, you can only find the range equal to 0.4, 2, 10, and 50 (maximum). 2) All the functions were made using Mathematica (Wolfram Research, Champaign, IL, USA), and for plotting purposes the Dirac Delta function is given a width of 0.2 and height of 1000. 3) All the functions are verified using numerical integration except for a few which are impossible to verify by numerical integration (for example the Dirac Delta function (bracewell_pg249_3)). 4) Π function’s definition ⎧ 1 1 x < ⎪ 2 Π=⎨ 1 ⎪0 x ≥ ⎩⎪ 2 5) Orthogonal polynomials Legendre polynomial n 1 d 2 n Pxn ( ) =nn x−1 2!ndx( ) Jacobi polynomial n n −1 −−αβd nn + α + β (αβ, ) ( ) ⎡⎤ Pxn ( ) =nn(11− x) ( + x) ( 1− x) ( 1+ x) 2!ndx⎣⎦ Laguerre polynomial zn−α ααez d −zn+ Lzn ( ) = n ez ndz! ( ) 0 Lznn( ) = Lz( ) 6) Legendre functions 1 2µ µ 11⎛⎞z + 11 Pzv ( ) =⎜⎟21 F− vv,1;1;+−−µ z Γ−(11µ) ⎝⎠z − ( 22) 1 µπ i 2 evπΓµ++1 1 µ µµ++vv1231 ++ Qzµµ( ) z−−−v 12 z1,;;2 F⎛⎞ v v ( ) =v+12( −×) 21⎜⎟+ 232Γ(vz+ ) ⎝⎠ 222 z in the complex plane cut along the real axis from -1 to 1. 1 2µ µ 11⎛⎞+ x 11 Pxv ( ) =⎜⎟21 F− vv,1;1;+−−−µ x ,1<< x 1 Γ−(11µ) ⎝⎠ −x ( 22) 11 − iiµπµπ µµ1 −i π⎡⎤22µ µ Qxvvv( ) =++ e⎢⎥ e Qxi( 00,11) e Qxi( −−) << x 2 ⎣⎦ 00 Pzvvv( ) == P( zQz),.( ) Q v( z) 7) Bessel functions and related functions vm+2 m ⎛⎞1 (−1) z ∞ ⎜⎟2 Jz ⎝⎠ v ( ) = ∑ m=0 mvm!1Γ( ++) Yzcosec v Jz cos v J z vvv( ) =ππ⎣⎦⎡⎤( ) −− ( ) (1) Hvvv( z) =+ J( z) iY( z) (2) Hvvv( z) = J( z) − iY( z) Modified Bessel functions 1 vm+2 ⎛⎞z ∞ ⎜⎟2 Iz ⎝⎠ v ( ) = ∑ m=0 mvm!1Γ( ++) π IzIz−vv( ) − ( ) Kzv ( ) = 2 sin vπ Kelvin’s and related functions 1 πi berzi+= bei z Jze4 vvv( ) ( ) ( ) − 1 πi berzi− bei z= Jze4 vvv( ) ( ) ( ) 1 πi kerzi+= kei z Kze4 vvv( ) ( ) ( ) − 1 πi kerzi− kei z= Kze4 vvv( ) ( ) ( ) ber(zzzz) == ber00( ) , bei( ) bei ( ) ker(zzzz) == ker00( ) , kei( ) kei ( ) Struve’s functions m z vm++21 ∞ (−1) H z ( 2) v ( ) = ∑ m=0 mvm33 Γ( +22)Γ( ++) z v+1 ( 2) 33z2 =+12Fv1; , ; − 33v ( 224) ΓΓ( 22) ( +) −1 =2112−−v π ⎡⎤Γvsz+1 ⎣⎦( 2) vv, ( ) 8) Hypergeometric functions ∞ aa... k ( 1 )kk( m ) z mnFa( 11,..., a m ;γγ ,..., n ; z) = ∑ ...k ! k=0 (γγ1 )kk( n ) 21Fabcz( ,;;) is Gauss’denoted by hypergeometric series and is often denoted by Fabcz( ,;;). 11Fabcz( ,;;) is Kummer’s confluent hypergeometric series and is sometimes denoted by Φ(acz;;). 9) Confluent hypergeometric functions Whittaker’s function 11+µ −z Mzze=+22 F1 µµ− k;2+ 1; z k,11µ ( ) ( 2 ) Γ−( 22µµ)Mzkk,,µµ( ) Γ( ) M− ( z) Wzk,µ ( ) =+ 11kk Γ−−( 22µµ) Γ( +−) Parabolic cylinder functions 11v+− 1 Dz2 24 z 2 W1 z2 v ( ) = 111v , 244+ ( 2 ) n n 11zz22 42d − Dzn ( ) =(−1) en e dz ( ) ∞ sint 1 si(x) =− dt=⎡⎤ Ei( ix) − Ei( − ix) ∫x ti2 ⎣⎦ ∞ sint Si(xdtx) ==+1 π si( ) ∫x t 2 ∞ cost 1 Ci(x) =− dt=− ci( x) =⎡⎤ Ei( ix) + Ei(− ix) ∫x ti2 ⎣⎦ 10) Elliptic functions and integrals 1 π − 1 Kk=2 1− k22 sinφφπ2 d=111 F , ;1; k 2 ( ) ∫0 ( ) 22221( ) 1 π 1 Ek=2 1− k22 sinφφ2 d=111 π F− , ;1; k 2 ( ) ∫0 ( ) 22221( ) Errors/Typos in books 1) erdelyi_8.2_17: Added to the 2µ power to the end, this has been verified by numerical integration. 2) erdelyi_8.2_38: Changed a2 to a . There is no reason for the squared there. Added a>0 constraint to compensate. 3) erdelyi_8.2_39: Changed a2 to a . There is no reason for the squared there. Added a>0 constraint to compensate. Added a negative sign to Fq(), numerical integration showed that it was missing the negative sign 4) erdelyi_8.2_46: b>0 is required, numerical integration fail for b=-1, (a<0 is fine since all a term is square of symmetric) 5) erdelyi_8.2_47: Changed a2 to a . There is no reason for the squared there. Added a>0 constraint to compensate. 6) erdelyi_8.3_4: Replace x2 in Hankel transform function to y 2 . 7) erdelyi_8.3_7: Replace x2 in Hankel transform function to y 2 . 8) erdelyi_8.3_43: For Fq(), multiply factor 2, verified by numerical integration. ππ 9) Bracewell Page335 Fig. 13.3 replace ∏(v) by δπ( −vv)+ δπ( +). 22 References 1. Erdélyi, A. and H. Bateman, Tables of integral transforms. vol. 2. 1954, McGraw-Hill: New York. 451. 2. Bracewell, R.N., The Fourier transform and its applications. 3rd ed. 2000, Boston: WCB/McGraw-Hill. 3. Andrews, L.C., B.K. Shivamoggi, and Society of Photo-optical Instrumentation Engineers., Integral transforms for engineers. 1999, Bellingham, Wash.: SPIE Optical Engineering Press. x, 353 p. 4. Andrews, L.C. and R.L. Phillips, Mathematical techniques for engineers and scientists. 2003, Bellingham: SPIE-The International Society for Optical Engineering. XV, 797 p. 5. Poularikas, A.D., The handbook of formulas and tables for signal processing. The electrical engineering handbook series. 1999, Boca Raton, Fla. u.a.: CRC Press u.a. Getr. Zählung [ca. 800 S.]. 6. Parker, K.J., Lecture Notes, Medical Imaging: Theory and Implementation. 2013, University of Rochester. 7. Bracewell, R.N., Two-dimensional imaging. Prentice-Hall signal processing series. 1995, Englewood Cliffs, N.J.: Prentice Hall. xiv, 689 p. 8. Wikipedia. Fourier transform. 9 August 2013 21:20 UTC Available from: http://en.wikipedia.org/w/index.php?title=Fourier_transform&oldid=567869073. 9. Wikipedia. Hankel transform 9 August 2013 12:02 UTC Available from: http://en.wikipedia.org/w/index.php?title=Hankel_transform&oldid=567808780. 10. Waag, R.C., J.A. Campbell, J. Ridder, and P.R. Mesdag, Cross-sectional Measurements and Exploitations of ultrasonic fields. 1985, IEEE Transactions on Sonics and Ultrasonics 32(1)26-35. Acknowledgements We are grateful to the Hajim School of Engineering and Applied Sciences and the Wadsworth C. Sykes Engineering Faculty Award for supporting this work. The excellent work of University of Rochester students Shuopeng Deng and Yuchuan Zhuang is gratefully noted. We are also indebted to Professor Al Clark for his patient advice. .

View Full Text

Details

  • File Type
    pdf
  • Upload Time
    -
  • Content Languages
    English
  • Upload User
    Anonymous/Not logged-in
  • File Pages
    5 Page
  • File Size
    -

Download

Channel Download Status
Express Download Enable

Copyright

We respect the copyrights and intellectual property rights of all users. All uploaded documents are either original works of the uploader or authorized works of the rightful owners.

  • Not to be reproduced or distributed without explicit permission.
  • Not used for commercial purposes outside of approved use cases.
  • Not used to infringe on the rights of the original creators.
  • If you believe any content infringes your copyright, please contact us immediately.

Support

For help with questions, suggestions, or problems, please contact us