马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我自己完成的部分代码求出1阶贝塞尔函数的根,关于后面的部分不知如何进行下去为好,请指点一下
clear all; maxv=2; maxs=10;
theta=linspace(0,2*pi,100); t=0.4;
[r]=root_bes(maxv,maxs); 求bessel1阶函数的根
function [r]=root_bes(maxv,maxs)
r= zeros(maxv-1, maxs); incr = 4.0;
for v=maxv-1
h = v+1.9*v^(1/3)+1;
if (v==0)
r(v+1,1) = fzero(@(x)besselj(v,x),2);
else
r(v+1,1) = fzero(@(x)besselj(v,x),h);
end
for s=2:maxs
r(v+1,s) = fzero(@(x)besselj(v,x),r(v+1,s-1)+incr);
end
end
r
[ 本帖最后由 ChaChing 于 2010-3-30 21:25 编辑 ] |