MATLAB Bessel函数的零点计算

hinihaoaaa

已于 2023-11-28 23:12:17 修改

阅读量1.7k

点赞数 6

CC 4.0 BY-SA版权

文章标签： matlab

于 2023-11-28 23:10:25 首次发布

本文链接：https://blog.csdn.net/w178191520/article/details/134679216

本文介绍了如何使用编程方法计算9阶贝塞尔函数的零点，特别关注了理论上的精度限制，并探讨了Airy光束、高斯光束和Bessel光束在经过大气湍流相位屏后的传输特性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

0～9阶的贝塞尔函数零点的计算公式，其中理论上计算零点个数N在50以内时较为精确：

Airy光束，基模高斯光束 Bessel光束经过大气湍流相位屏的传输特性

function rootBessel = besal_pi0(n,N)    

% 求解n阶贝塞尔函数的零点(0-9)
% n为贝塞尔函数阶数 
% N为要求的零点数量

j = zeros(n+1, N);    % 贝塞尔函数的根
incr = 4.0;
for v = 0 : n
   h = v + 1.9*v^(1/3)+1;
   if (v == 0)             % 0阶贝塞尔函数的第一个零点
       j(v+1,1) = fzero(@(x)besselj(v,x),2);
   else                    % 1阶及以上阶贝塞尔函数的第一个零点
       j(v+1,1) = fzero(@(x)besselj(v,x),h);
   end
   for s = 2 : N           % 贝塞尔函数的第2个及后面的零点
       j(v+1,s) = fzero(@(x)besselj(v,x),j(v+1,s-1)+incr);
   end    
end

rootBessel = j;