24小时热门版块排行榜    

查看: 2199  |  回复: 6

xiaolobster

新虫 (初入文坛)


[交流] 用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,但是都没能给出正确的结果。
求高手指点!
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

bruceleem

铜虫 (著名写手)


★ ★ ★ ★ ★ ★
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的回帖

xiaolobster

新虫 (初入文坛)


引用回帖:
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的回帖

bruceleem

铜虫 (著名写手)


祝楼主好运!
4楼2013-11-29 09:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

真心妹夫

新虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
请问楼主最后用solve怎么解决的?
5楼2014-02-18 04:21:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

九垓

新虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主 问题解决了吗/
6楼2015-05-19 14:17:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

九垓

新虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by xiaolobster at 2013-11-29 09:04:49
非常感谢。主要是我的版本有点低,升级版本之后,用solve能解决。不过仍然十分感谢您的热情回复。5个金币要给,必须的。...

对复杂贝塞尔形式的化简,,楼主可好 小弟求助
7楼2015-05-19 14:20:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaolobster 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护,明天很可能放榜! +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
请填处理意见