MATLAB仿真贝塞尔函数

最新推荐文章于 2024-03-24 15:21:42 发布

matlab练习生

最新推荐文章于 2024-03-24 15:21:42 发布

阅读量2.1k

点赞数

分类专栏： # MATLAB基础操作文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/fendoudecheng/article/details/130404716

版权

MATLAB基础操作专栏收录该内容

27 篇文章 11 订阅

订阅专栏

贝塞尔函数来源于贝塞尔方程：

$z^2\frac{d_y^2}{dz^2}+z\frac{dy}{dz}+\left(z^2-n^2\right)y=0$

第一二类贝塞尔函数分别是该方程的两个线性无关解 $J_n\left(z\right)$ 、 $Y_n\left(z\right)$
贝塞尔方程如下：

n阶第一类贝塞尔函数：（积分表达式和级数表达式）

$J_n\left(z\right)=\frac{1}{2\pi}\int_{-\pi}^{\pi}e^{-i\left(zsin\theta-n\theta\right)}d\theta$

$J_n\left(z\right)=\sum_{k=0}^{\infty}\frac{{(-1)}^k{(z/2)}^{n+2k}}{k!(n+k)!}$

第一类Hankel函数：

$H_n^{\left(1\right)}\left(z\right)=J_n\left(z\right)+iY_n\left(z\right)$

第二类Hankel函数:

$H_n^{\left(2\right)}\left(z\right)=J_n\left(z\right)-iY_n\left(z\right)$

可以看出，Hankel函数是贝塞尔函数特定的线性组合。

贝塞尔方程的一般解是上述任意两个的线性组合。

对于整数阶的贝塞尔函数，可以用FFT（快速傅里叶变换）来计算。

MATLAB包含有贝塞尔函数的表达，可以直接调用函数：

%J = besselj(nu,z);

第一个参数是贝塞尔函数的阶数，第二个参数是自变量
绘制前五阶的贝塞尔函数如下：

clear
M = 5;
N = 1001;%每一条曲线上的点数
Xmin = 0;
Xmax = 15;

x = linspace(Xmin,Xmax,N);
y = zeros(N,M);%y的每一列表示一阶贝塞尔函数
for m = 0:M-1
    y(:,m+1) = besselj(m,x);%贝塞尔函数有0阶
end

plot(x,y);
[ymax pos] = max(y);
for i = 1:M
    text(x(pos(i)+1),ymax(i)-0.05,['J_' num2str(i-1) '(x)']);
end
grid on
xlabel('x')
ylabel('J_n(x)')

结果：