Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

A Bessel decomposition for arbitrary function

  1. Apr 19, 2016 #1
    Orthogonality condition for the 1st-kind Bessel function J_m
    $$\int_0^R J_m(\alpha_{mp})J_m(\alpha_{mq})rdr=\delta_{pq}\frac{R^2}{2}J_{m \pm 1}^2(\alpha_{mn}),$$
    where α_{mn} is the n^{th} positive root of J_m(r), suggests that an original function f(r) could be decomposed into a series of 1-st kind Bessel functions of order m:
    $$f(r)=\sum_n a_n\cdot J_m(\alpha_{m,n}*r/R),$$ with $$a_n=\frac{2}{R^2J_{m \pm 1}^2(\alpha_{mn})}\int_0^R J_m(\alpha_{mn}*r/R) f(r) r dr$$. However what if the original function is something like:
    $$f(r)=\sum_n a_n\cdot J_m(\alpha_{m,n}*r/R)+b_n\cdot J_{m+1}(\alpha_{m+1,n}*r/R),$$
    which combines Bessel functions of different orders. The final question is can any function be decomposed into Bessel functions someshow, like FFT/IFFT(fast fourier transform / inverse fast Fourier transform)?
     
  2. jcsd
  3. Apr 19, 2016 #2

    Orodruin

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper
    Gold Member

    A priori, you could expand any Bessel function ##J_n(\alpha_{nk'}r/R)## in terms of the Bessel functions ##J_m(\alpha_{mk}r/R)## for any fixed ##n## and ##m##. It is just a matter of using a different basis for your function space.
     
  4. Apr 19, 2016 #3
    Thanks for reply. I did NOT know Bessel function of order n can be expanded into Bessel functions of order m. Could you please give me an identity?
     
    Last edited: Apr 19, 2016
  5. Apr 19, 2016 #4

    Orodruin

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper
    Gold Member

    While it can be done, it is not clear why you would want to do this. In general, Bessel function expansions will arise naturally as Bessel functions of different orders are eigenfunctions of different Sturm-Liouville operators. It is therefore usually quite clear what Bessel functions you should use for your expansion.
     
  6. Apr 19, 2016 #5
    In my case, I do not know which order of Bessel functions should be used as the basis functions.

    In many other cases, the original profile ##f(r)## contains Bessel functions of different orders. I am wondering how can I decompose it with ##J_m(r)## with fixed ##m =1##?
     
  7. Apr 19, 2016 #6

    Orodruin

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper
    Gold Member

    Perhaps it would be a good idea to state the type of problem you are looking to solve. I think it would help us understand what you need better.
     
  8. Apr 19, 2016 #7
    What I am doing now is simple. I have a random series ##f(r)##, with boundary condition ##f(0)=0## and ##f(1)=0##. I want to approximate it using ##J_1(\alpha_{1n}r)##, which satisfies the B.Cs.

    I had a plan to solve a system of ##nr## linear equations numerically, where ##nr## is the number of grid points along ##r## direction, as follows:
    $$\sum_{i=1}^{nr} a_i J_1(\alpha_{1,i}\cdot r_i)=f(r_i),$$ for ##i=1,2,3...,nr##
    Not sure if it is what you are going to suggest?

    Not sure how do it analytically using orthogonality and other identities of Bessel functions.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted



Similar Discussions: Bessel decomposition for arbitrary function
  1. Bessel function ... (Replies: 10)

Loading...