课堂授课专题3：特殊函数及可视化.pptVIP

Page * 贝塞尔函数 贝塞尔函数的表达式 第一类贝塞尔函数： 第二类贝塞尔函数： 虚宗量贝塞尔方程 将贝塞尔方程的宗量x变换为虚数ix，于是得到虚宗量贝塞方程： Page * 贝塞尔函数 特殊贝塞尔函数：虚宗量贝塞尔函数 Iv(x)为v阶第一类虚宗量贝塞尔函数（第一类修正贝塞尔函数） （1）当 (整数)时，虚宗量贝塞方程的通解为： （C和D为任意常数） （2）当 取任意值时，定义： Kv(x)为v阶第二类虚宗量贝塞尔函数（或麦克唐纳函数，或第二类修正贝塞尔函数） Page * 贝塞尔函数 当v 取任意值时，虚宗量贝塞方程的通解为： 贝塞尔函数的计算和图形绘制 j=Besselj(nu,z) nu为阶，z为贝塞尔函数的常点 （或复数变量） Besselj 第一类贝塞尔函数，简称贝塞尔函数 Bessely 第二类贝塞尔函数，又称诺依曼函数 Besselh 第三类贝塞尔函数，又称汉克尔函数 Besseli 第一类虚宗量贝塞尔函数，又称虚宗量贝塞尔函数 Besselk 第二类虚宗量贝塞尔函数，又称虚宗量汉克尔函数 Page * 贝塞尔函数的图形 绘制贝塞尔函数图形 y=besselj(0:3,(0:0.2:10)); figure(1) plot((0:0.2:10),y) legend(J_0,J_1,J_2,J_3) Page * 贝塞尔函数 寻找贝塞尔函数零点 %方法I(插值法) x=0:0.05:50; y=besselj(0,x); LD=[]; for k=1:length(y)-1 if y(k)*y(k+1)0 h = interp1(y(k:k+1),x(k:k+1),0); %插值函数 LD=[LD,h]; end end %方法II(利用Matlab内建函数) j=1; LD=[ ]; y=inline(besselj(0,x),x) %直接定义函数的表达式 for k=1:16 while ~(besselj(0,j)*besselj(0,j+1)0) j=j+1; end q=fzero(y,j); %查找一元连续函数的零点 j=j+1; LD=[LD,q]; k=k+1; end Page * 贝塞尔函数 贝塞尔函数及其零点 LD = 2.4049 5.5201 8.6537 11.7915 14.9309 18.0711 21.2116 24.3525 27.4935 30.6346 33.7758 36.9171 40.0584 43.1998 46.3412 49.4826 Page * 贝塞尔函数 绘制诺依曼函数的图形 y=bessely(0:1,(0:0.02:10)); plot((0:0.02:10),y) legend(N_0,N_1) axis ([0 10 -3.5 1]) grid on Page * 贝塞尔函数 绘制虚宗量贝塞尔函数的图形 I=besseli(0:1,(0.1:0.3:3)); plot((0.1:0.3:3),I) legend(I_0,I_1) axis ([0 3 0 5]) Page * 贝塞尔函数 绘制虚宗量汉克尔函数的图形 K=besselk(0:1,(0.1:0.1:3)); plot((0.1:0.1:3),K) legend(K_0,K_1) axis ([0 3 0 10]) 问题由来：与时间有关的三维方程的变量分离 分解为两个常微分方程： 　　　　　　　　　　　 亥姆霍兹方程 球贝塞尔方程及其解 2. 三维热传导(输运)方程 　　　　　　　　　　　　 分离变量：　　　　　　　　　　　　 分解为两个常微分方程 　　　　　　　　　　　　 亥姆霍兹方程: 球贝塞尔方程及其解 亥姆霍兹方程在球坐标系中的变量分离 分离变量形式的解：　　　　　　　　　　　　 球贝塞尔方程及其解 分解为两个常微分方程: 　 　(10.4.23) 　 (10.4.24) 球函数方程 球贝塞尔方程及其解 l阶球贝塞尔方程 阶的贝塞尔方程 。球贝塞尔方程的解称为球贝塞尔函数。 球贝塞尔方程及其解 Page * 球贝塞尔函数 第一类球贝塞尔函数： 第二类球贝塞尔函数（球诺依曼函数）： 第二类球贝塞尔函数（球汉克尔函数