Matlab编程贝塞尔函数图片

姓名：张俊红学号：S150431033 专业：集成电路工程

1．第一类贝塞尔函数曲线图

clear

format long

x=(0:0.1:20);

J_0=besselj(0,x);

J_1=besselj(1,x);

J_2=besselj(2,x);

y=0*x;

figure(1)

h1=plot(x,J_0,'k-',x,J_1,'k-.',x,J_2,'k--','linewidth',1.8);

hold on;h2=plot(x,y,'k-');

box off;

grid on;

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

xlabel('Variable X');ylabel('Variable Y');

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

text(1.2,0.8,'J_0(x)');text(2.3,0.6,'J_1(x)');text(4.5,0.3,'J_2(x)');

legend('零阶','一阶','二阶',0);

set(h1,'LineSmoothing','on');set(h2,'LineSmoothing','on');

export_fig bess01.png;

2.第二类贝塞尔函数曲线图

clear

format long

x=(0:0.1:20);

Y_0=bessely(0,x);

Y_1=bessely(1,x);

Y_2=bessely(2,x);

y=0*x;

figure

h1=plot(x,Y_0,'k-',x,Y_1,'k-.',x,Y_2,'k--','linewidth',1.8);

hold on;h2=plot(x,y,'k-');

box off;

grid on;

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

xlabel('Variable X');ylabel('Variable Y');

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

text(1.8,0.6,'Y_0(x)');text(3.4,0.5,'Y_1(x)');text(6.3,0.2,'Y_2(x)'); legend('零阶','一阶','二阶',0);

set(h1,'LineSmoothing','on');set(h2,'LineSmoothing','on'); export_fig bess02.png;

clear

format long

x=(0:0.1:20);

I_0=besseli(0,x);

I_1=besseli(1,x);

I_2=besseli(2,x);

figure

h=plot(x,I_0,'k-',x,I_1,'k-.',x,I_2,'k--','linewidth',1.8);

box off;

grid on;

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

xlabel('Variable X');ylabel('Variable Y');

title('零阶、一阶、二阶修正第一类贝塞尔函数曲线图'); text(0.7,1.6,'I_0(x)');text(2.7,2.4,'I_1(x)');text(4.1,6,'I_2(x)'); legend('零阶','一阶','二阶',0);

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

export_fig bess03.png

clear

format long

x=(0:0.1:20);

K_0=besselk(0,x);

K_1=besselk(1,x);

K_2=besselk(2,x);

figure

h=plot(x,K_0,'k-',x,K_1,'k-.',x,K_2,'k--','linewidth',1.8);

box off;

grid on;

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

xlabel('Variable X');ylabel('Variable Y');

title('零阶、一阶、二阶修正第二类贝塞尔函数曲线图'); text(0.1,0.5,'K_0(x)');text(0.3,4.5,'K_1(x)');text(0.55,8,'K_2(x)'); legend('零阶','一阶','二阶',0);

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

export_fig bess04.png;

5.L=0时图解法确定LP模的ha函数图

clear all

format long

ha=[0:0.0001: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)); y4=0*ha;

figure(1)

h1=plot(ha,y,'k-',a,y1,'r-',b,y2,'b-',c,y3,'g-','linewidth',1.8);

hold on;h2=plot(ha,y4,'k--');

box off;

grid on;

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

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

text(0.5,3.2,'V=2');text(3,5.2,'V=5');text(4,8.1,'V=8');

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

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

export_fig lp.png;