How to rotate a function about the y-axis

17 views (last 30 days)
I wonder if anyone can help.
I have a Bessel function defined from x=0:20. I am trying to rotate it about the y-axis, so I can get a plot which resembles an Airy function.
I am struggling to write a for loop which will do this. Does anyone have a suggestion or pointers?
Thank you
  6 Comments
Jan
Jan on 5 Oct 2011
You can plot a matrix. But rotating the Bessel function leads to a circular shape, not to a rectangular matrix.

Sign in to comment.

Accepted Answer

Teja Muppirala
Teja Muppirala on 5 Oct 2011
Something like this maybe?
[X,Y]=meshgrid(-20:0.1:20);
R = sqrt(X.^2 + Y.^2);
J=besselj(1,R);
Z=(J./R);
Z(isnan(Z)) = 0.5; %Fix the divide by zero problem
surf(X,Y,Z);
shading flat;
camlight;
  1 Comment
Katy Pointer
Katy Pointer on 5 Oct 2011
That's amazing! Thank you. I had tried something similar, but didn't quite get the radius bit right. Thanks so much :)

Sign in to comment.

More Answers (2)

Patrick Kalita
Patrick Kalita on 5 Oct 2011

Jan
Jan on 5 Oct 2011
What about using AIRY?
v = -20:0.1:20;
m = bsxfun(@plus, v, i*v');
Z = airy(m);

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by