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

,

一．第一类贝塞尔函数第一类贝塞尔函数的代码如下：

clear all

x=(0::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 on

xlabel('X');

ylabel('J(X)');

title('第一类贝塞尔函数曲线图');

text(1,,'J0(X)');

text(2,,'J1(X)');

：

text(4,,'J2(X)');

text(12,,'J3(x)');

二．第二类贝塞尔函数第二类贝塞尔函数的代码如下：

clear all；clc；

x=(0::25);

Y0=bessely(0,x);

}

Y1=bessely(1,x);

Y2=bessely(2,x);

Y3=bessely(3,x);

plot(x,Y0,':',x,Y1,'-.',x,Y2,'-',x,Y3,'-');

axis([0,25,-1,1]);

text(2,,'Y0(X)');

text(4,,'Y1(X)');

text(11,,'Y2(X)');

；

text(13,,'Y3(X)');

grid on

xlabel('X');

ylabel('Y(X)');

title('第二类贝塞尔函数曲线图');

三．第一类修正贝塞尔函数

第一类修正贝塞尔函数的代码如下：《

clear all；clc；

x=(0::25);

I0=besseli(0,x);

I1=besseli(1,x);

I2=besseli(2,x);

I3=besseli(3,x);

plot(x,I0,':',x,I1,'-.',x,I2,'-',x,I3,'-'); axis([0,5,0,10]);

￥

text,,'I0(X)');

text,2,'I1(X)');

text,,'I2(X)');

text,5,'I3(X)');

grid on

xlabel('X');

ylabel('I(X)');

title('第一类修正贝塞尔函数曲线图');

~

四．第二类修正贝塞尔函数第二类修正贝塞尔函数的代码如下：

clear all；clc；

x=(0::20);

K0=besselk(0,x);

K1=besselk(1,x);

K2=besselk(2,x);

！

K3=besselk(3,x);

plot(x,K0,':',x,K1,'-.',x,K2,'-',x,K3,'-');

axis([0,5,0,10]);

text,,'K0(X)');

text,,'K1(X)');

text(1,2,'K2(X)');

text,1,'K3(X)');

grid on

)

xlabel('X');

ylabel('K(X)');

title('第二类修正贝塞尔函数曲线图');

五．阶跃函数阶跃函数的代码如下：

clear all

clc;

>

format long

ha=[0::12];a=[0::2];b=[0::5];c=[0::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-','linewidth',;

grid on

axis([0,12,-15,15]);

title('l=0时图解法确定LP模的ha');

text,,'V=2');text(3,5,'V=5');text(4,8,'V=8');

xlabel('ha');ylabel(' Y');

set(h,'LineSmoothing','on');