Plotting Bessel function in MATLAB

math matlab bessel-functions
4alala picture 4alala · May 5, 2010 · Viewed 12.9k times · Source

In MATLAB how do you plot

f(r) = { 2*J1(a*r) / r }^2

where a = 2*pi and J1 is Bessel function of the 1st kind and r = sqrt(x^2 + y^2)

This should plot in 3D, i.e. kind of be like a bubble (not sure how to do this)

Answer

Jacob picture Jacob · May 5, 2010

Use besselj --- the Bessel function of first kind --- to generate J1. I suppose you have to vary a and r to generate the "bubble".

I generated the following by varying x and y from -1:0.01:1 and plotting meshing points (x,y,f), I don't know if this is what you want.

Code

a = 2*pi;
[X Y] = meshgrid(-1:0.01:1,-1:0.01:1);
R = sqrt(X.^2+Y.^2);
f = (2*besselj(1,a*R(:))./R(:)).^2;
mesh(X,Y,reshape(f,size(X)));
axis vis3d;

Log plot

Doresdoom suggestion, I replaced axis vis3d; with set(gca,'Zscale','Log').

alt text

Mesh

alt text

Related questions

Optional args in MATLAB functions

How can I declare function in MATLAB with optional arguments? For example: function [a] = train(x, y, opt), where opt must be an optional argument.

Read More
Using fourier analysis for time series prediction

For data that is known to have seasonal, or daily patterns I'd like to use fourier analysis be used to make predictions. After running fft on time series data, I obtain coefficients. How can I use these coefficients for prediction? …

Read More
How do I use the least squares approximation in MATLAB?

For a homework assignment in linear algebra, I have solved the following equation using MATLAB's \ operator (which is the recommended way of doing it): A = [0.2 0.25; 0.4 0.5; 0.4 0.25]; y = [0.9 1.7 1.2]'; x = A \ y which produces the following answer: x = 1.7000 2.0800 For the next part …

Read More