24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xiaolobster

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 140
帖子: 10
在线: 37.3小时
虫号: 964557

[交流] 用matlab求bessel函数的根

问题看上去很简单，方程如下所示：
v*sqrt（1-b）*J1（v*sqrt(1-b)）/J0(v*sqrt(1-b))=-v*sqrt(b)*K1(v*sqrt(b))/K0(v*sqrt(b))
其中b为未知数。v为常数，J1我第一类第1阶贝塞尔函数，J0为第一类第0阶贝塞尔函数；K1为第二类第1阶贝塞尔函数，K2为第二类第0阶贝塞尔函数。

问题：通过解上述方程，求解b的值。
我在matlab中尝试了fzero，和solve，但是都没能给出正确的结果。
求高手指点！

回复此楼

» 猜你喜欢

寻求一种能扛住强氧化性腐蚀性的容器密封件已经有7人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab如何求解含有bessel方程中的参数已经有5人回复
求助：Matlab中如何绘制三维隐函数图像已经有20人回复
谁给解释一下扩散系数的单位问题？？已经有5人回复
数值求解N-S方程一般用什么编程啊？？？已经有9人回复
Matlab画隐函数，ezplot 已经有3人回复
求助matlab的发fgoalattain函数已经有4人回复
哪些书里面有讲解关于含有积分方程的未知数求解问题啊？非常感谢！已经有9人回复
函数有第一类间断点，则不存在原函数已经有3人回复
求matlab画一个函数的图已经有6人回复
如何用Matlab画z关于x和y的函数图？谢谢大家！已经有8人回复
求助：能不能帮忙写一个用matlab求一个函数的傅里叶变换的程序已经有16人回复
如何根据已知函数画函数曲线？谢谢已经有4人回复
【求助】请问怎么求函数方程的函数表达式？已经有11人回复
【求助】关于matlab中的函数代换已经有10人回复

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2013-11-28 10:33:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xiaolobster

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 140
帖子: 10
在线: 37.3小时
虫号: 964557

引用回帖:

2楼: Originally posted by bruceleem at 2013-11-28 19:12:56
function bessel_problem
%-----------------------------------------
function f=besselfunc(b,v)
f=v.*sqrt(1-b).*besselj(1,v.*sqrt(1-b))./besselj(0,v.*sqrt(1-b))-(-v.*sqrt(b).*bessely(1,(v.*sqrt(b))) ...

非常感谢。主要是我的版本有点低，升级版本之后，用solve能解决。不过仍然十分感谢您的热情回复。5个金币要给，必须的。

赞一下

回复此楼

3楼2013-11-29 09:04:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 7 个回答

bruceleem

铜虫 (著名写手)

应助: 3 (幼儿园)
金币: 3710
帖子: 1033
在线: 284.7小时
虫号: 946615

★ ★ ★ ★ ★ ★
xiaolobster: 金币+5, ★有帮助 2013-11-29 09:03:39
xzhdty: 金币+1, 谢谢参与 2013-11-30 10:51:42

function bessel_problem
%-----------------------------------------
function f=besselfunc(b,v)
f=v.*sqrt(1-b).*besselj(1,v.*sqrt(1-b))./besselj(0,v.*sqrt(1-b))-(-v.*sqrt(b).*bessely(1,(v.*sqrt(b)))./bessely(0,(v.*sqrt(b))));
end
%-----------------------------------------
v=10;
b0=0.1;
[b,yb,exitflag]=fsolve(@(b)besselfunc(b,v),b0);
b
yb
exitflag
end

将代码保存为m文件，在命令窗口中输入bessel_problem即可求解。给v赋不同的值得到的结果不大一样，这里给v赋值为10。

赞一下

回复此楼

2楼2013-11-28 19:12:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖