Plotting the Bessel function equation

30 次查看（过去 30 天）
sha
sha 2015-1-7
编辑： Image Analyst 2015-1-7
Hello, Can anyone tell me how to plot
J_1^2(x)/J_0^2(x) as shown in the attached file?
The code I have prepared so far is as follows:
syms x
b1=(besselj(0, x).^2)
b2=(besselj(1, x).^2)
b= b2/b1;
plot(b,x)
Does it make sense? Thanks

请先登录，再回答此问题。

采纳的回答

Josh Meyer
Josh Meyer 2015-1-7
编辑：Josh Meyer 2015-1-7
I think you're pretty close, but I did the plotting numerically instead of using a symbolic variable. If you're interested in the peaks of the ratio, you'll probably want to zoom out quite a bit more (at least 1e5). Hope this helps.
x = 0:0.01:20;
b1 = besselj(0, x);
b2 = besselj(1, x);
b = b2.^2./b1.^2;
plot(x,b,x,b1,x,b2)
legend('J_1^2(x)/J_0^2(x)','J_0(x)','J_1(x)')
axis([0 20 -0.5 15])

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Bessel functions 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by