Article 679SJ Sinc approximation to Bessel function

Sinc approximation to Bessel function

by
John
from John D. Cook on (#679SJ)

The Bessel functions Jn for even n look something like the sinc function. How well can you approximate the former by sums of the latter? To make things concrete, we'll approximate J2. Here's a plot of J2.

bessel2_plot.png

And here's a plot of sinc(x) = sin(x)/x.

sinc_plot2.png

The sinc approximation for a function f(x) is given by

sinc_approx_eqn.svg

Sinc approximation can be remarkably accurate, nearly optimal in some sense.

The accuracy of the approximation increases as n gets larger and h gets smaller. We will fix n = 10. How should we pick h? The paper cited in this post suggests using

bessel_sinc_step.svg

Let's try that and see what happens.

bessel_sinc1.png

The approximation isn't very good overall, though it's excellent near 0.

Before making plots, I had a plausible argument for why the value of h suggested above might be optimal. I also had an argument for why a much larger value of h, something on the order of 8 might be optimal. Turns out both are wrong. You can get a good approximation over a larger range by choosing h around 2.6.

bessel_sinc2.png

Related postsThe post Sinc approximation to Bessel function first appeared on John D. Cook.
External Content
Source RSS or Atom Feed
Feed Location http://feeds.feedburner.com/TheEndeavour?format=xml
Feed Title John D. Cook
Feed Link https://www.johndcook.com/blog
Reply 0 comments