这个是第一类贝塞尔函数,想知道ym咋求出来

img

这个是第一类贝塞尔函数,想知道ym咋求出来!
这个是第一类贝塞尔函数,想知道ym咋求出来!

f = @(x)0.61*besselj(0,x) - 0.39 * besselj(2,x);
x = 0:0.01:1000;
plot(x,f(x))

img
看图,因为与x轴有无数交点,所以求解数值解很多个,我这里给出0到1000里面的数值解:

f = @(x)0.61*besselj(0,x) - 0.39 * besselj(2,x);
x = [];
for i = 1:1000
    x = [x; fsolve(f, i)];
end
x = uniquetol(x, 1e-3)
syms x
eq = 0.61*besselj(0,x) - 0.39 * besselj(2,x);
solve(eq, x)

结果为

ans =
 
-193.98996985366718682905757542382