贝塞尔函数MATLAB仿真(含程序)

1 .第一类贝塞尔函数第一个贝塞尔函数的代码如下：clear allx=(0:0.01:25 )J0=besselj(0，x )J1=besselj(1，x )J2=besselj(2，x )J3=besselj(3，x )plot(x，j0，x，J1，- .x，J2，-，x，J3，- )axis ( 0，25，- 1，1 )grid onxlabel(X )ylabel(J(X ) )；title (第一类贝塞尔函数的图表)文本(1，0.8，J0(X ) )文本(2，0.6，J1(X ) )文本(4，0.4，J2(X ) )文本(12，0.2，J3(x ) )2 .第二类贝塞尔函数第二个贝塞尔函数的代码如下：clear all； clc；x=(0:0.01:25 )Y0=bessely(0，x )Y1=bessely(1，x )Y2=bessely(2，x )Y3=bessely(3，x )plot(x，y 0，x，Y1，- .x，Y2，-，x，Y3，- )axis ( 0，25，- 1，1 )文本(2，0.6，Y0(X ) )文本(4，0.4，Y1(X ) )文本(11，0.24，Y2(X ) )文本(13，0.24，Y3(X ) )grid onxlabel(X )ylabel(Y(X ) )；title (第二类贝塞尔函数的图表)3 .第一类修正贝塞尔函数第一个修改贝塞尔函数的代码如下所示clear all； clc；x=(0:0.01:25 )I0=besseli(0，x )I1=besseli(1，x )I2=besseli(2，x )I3=besseli(3，x )plot(x，I 0，x，I1，- .x，I2，-，x，I3，- )axis ( 0，5，0，10 )文本(0. 5，1.3，I0(X ) )文本(2. 5，2，I1(X ) )文本(3. 5，3.1，I2(X ) )文本(4. 4，5，I3(X ) )grid onxlabel(X )ylabel(I(X ) )；title (第一类修正贝塞尔函数图)4 .第二类修正贝塞尔函数第二个修改贝塞尔函数的代码如下：clear all； clc；x=(0:0.01:20 )K0=besselk(0，x )K1=besselk(1，x )K2=besselk(2，x )K3=besselk(3，x )plot(x，k 0，x，K1，- .x，K2，-，x，K3，- )axis ( 0，5，0，10 )文本(0. 5，0.5，K0(X ) )文本(0. 3，3.5，K1(X ) )text (1，2，K2(X ) )；文本(1. 8，1，K3(X ) )grid onxlabel(X )ylabel(K(X ) )；title (第二类修正贝塞尔函数图)5 .阶跃函数步进函数的代码如下所示clear allclc；format long (格式长)ha=0:0.01:12； a=0:0.001:2； b=0:0.001:5； c=0:0.001:8；y=ha.*besselj(1，ha)./besselj(0，ha )y1=sqrt(4-a.2).*besselk(1，sqrt(4-a.2)./besselk(0，sqrt(4-a.2)；y2=sqrt(25-b.2).*besselk(1，sqrt(25-b.2)./besselk(0，sqrt(25-b.2)；y3=sqrt(64-c.2).*besselk(1，sqrt(64-c.2)./besselk(0，sqrt(64-c.2)；h=plot(ha，y，k-，a，y1，r-，b，y2，b-，c，y3，g-，线宽，1.2 )；grid onaxis ( 0，12，- 15，15 )title (决定l=0时的图解法LP模型的ha