Bessel problem "inf+ infi"

Question

Salwa Ben Mbarek 2021-4-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/802076-bessel-problem-inf-infi

评论： Salwa Ben Mbarek 2021-4-15

采纳的回答： Alan Stevens

Hello,

I'm getting a "inf+ infi" error or nan when I compute the following code with bessel functions:

freq = linspace(4*1e6, 31*1e6, 0.5*1e6) ;

w = 2*pi.*freq; % angular frequency

%r = d/2; % wire radius

p= 6.5 *1e-3;

d= 0.6*1e-3;

mu0 = 4*pi*1e-7;

mur= 200;

eta0 = 377;

sigw= 1.32*1e6;

Rw = (pi*sigw*(r^2))^(-1);

%tau = mu0*sigw*(r^2);

tau = mu0*mur*sigw*(r^2);

num = sqrt(1i.*w.*tau).*besselj(0, sqrt(1j.*w.*tau));

The program is simple tough, I have numerical value for sqrt(1j.*w.*tau) but when I apply Bessel to it I don't get a numerical value.

Could you please help me ?

Thank you.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Alan Stevens 2021-4-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/802076-bessel-problem-inf-infi#answer_675126

在 MATLAB Online 中打开

The following works with no inf or nan problems

freq = linspace(4*1e6, 31*1e6, 100) ;
w = 2*pi.*freq; % angular frequency
p= 6.5 *1e-3;
d= 0.6*1e-3;
r = d/2; % wire radius
mu0 = 4*pi*1e-7;
mur= 200;
eta0 = 377; 
sigw= 1.32*1e6;
Rw = (pi*sigw*(r^2))^(-1);
tau = mu0*mur*sigw*(r^2);
num = sqrt(1i.*w.*tau).*besselj(0, sqrt(1j.*w.*tau));