应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 用matlab求bessel函数的根
71/1返回列表
查看: 2201  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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，但是都没能给出正确的结果。
求高手指点！
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2013-11-28 10:33:28
已阅   回复此楼   关注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 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 11:46 by 4GBAYCdVQoK3
[基金申请] 今天系统多次维护，明天很可能放榜！ +7 zju2000 2026-08-16 7/350 2026-08-17 11:31 by Iktsuarpok9
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +10 Ldrop2023 2026-08-13 10/500 2026-08-17 09:59 by 小豌豆_发芽
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6GojgJvkDudM 2026-08-16 3/150 2026-08-17 09:49 by xLVPIRuSvCUe
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 74/3700 2026-08-16 21:05 by 医学老男孩
[论文投稿] 岩土工程学报什么时候才能终审完呐 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
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
关闭
请填处理意见
关闭