Question

matlab 中贝塞尔函数计算出现NAN

matlab 中贝塞尔函数计算出现NAN
例：
rho = logspace(0,4,100);
nl = length(rho);
k = sqrt( -1i * miu ./ rho' * w(N) + repmat(w(N) * w(N) * epu * miu, nl, 1) ); % 不同频率时、不同电阻率条件下的波数

tmp=-1i * k * r;
tmp_VMD = 1i * k * r/2;
tmp2_VMD = besseli(1,tmp_VMD) .* besselk(1,tmp_VMD); %I1*K1
tmp4_VMD = besseli(2,tmp_VMD) .* besselk(2,tmp_VMD);
w = 8000Hz时就没值了，出现NAN。求助

Answer 1

MATLAB提供了计算贝塞尔函数的函数，具体包括：

• besselj - 第一类贝塞尔函数，或简称贝塞尔函数；

• bessely - 第二类贝塞尔函数，又称诺伊曼函数（Neumann function）；

• besseli - 第一类修正贝塞尔函数；

• besselk - 第二类修正贝塞尔函数；

• besselh - 第三类贝塞尔函数，又称汉克尔函数（Hankel function）。

这几个函数的调用语法基本相同，例如

J = besselj(nu,Z)

J = besselj(nu,Z,1)

[J,ierr] = besselj(nu,Z)

其中，nu为贝塞尔函数的阶数，Z为函数自变量。阶数必须为实数，但Z可以是复数。

值得一提的是，上述函数是MATLAB基本模块（也就是说不需要任何附加的工具箱）提供的特殊函数，采用数值方法计算；而符号数学工具箱则提供了第一和第二类的4个贝塞尔函数，名称和调用方式都与MATLAB基本系统的4个函数完全一致，但支持微分、积分等符号运算。