用Matlab求贝塞尔函数根.docxVIP

(完满word版)用Matlab求贝塞尔函数根 (完满word版)用Matlab求贝塞尔函数根 (完满word版)用Matlab求贝塞尔函数根 用 Matlab 求贝塞尔函数的根 clear clc R=1; maxv = 5; maxs = 5; j= zeros(maxv, maxs);%用于存放第一类贝塞尔函数的根，同一行阶数相同，列数表示第几个根 incr = 4.0; for v=0:maxv-1 h = v+1.9*v^(1/3)+1; if (v==0) j(v+1,1) = fzero(@(x)besselj(v,x),2); else j(v+1,1) = fzero(@(x)besselj(v,x),h); end for s=2:maxs j(v+1,s) = fzero(@(x)besselj(v,x),j(v+1,s-1)+incr); end end j maxv1 =5; maxs1 = 5; jd= zeros(maxv1, maxs1);%用于存放导数的根 incr1 = 3.0; for v1=0:maxv1-1 h1 = v1+2; if (v1==0) jd(v1+1,1) = fzero(@(x)besseld(v1,x),3); else jd(v1+1,1) = fzero(@(x)besseld(v1,x),h1); end for