| 查看: 2153 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
用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，但是都没能给出正确的结果。 求高手指点！ |
回复此楼» 猜你喜欢
系统今天又提示维护了，估计离放榜不远了
已经有5人回复
-
这个自发加氧反应的机理是什么？
已经有7人回复
-
J-J-W不为人知的一面
已经有19人回复
-
没消息就是被刷了呗
已经有14人回复
-
微信指数没变化，科研之友没阅读
已经有8人回复
-
时间戳他又来了
已经有12人回复
-
这种情况还有戏吗
已经有11人回复
-
同事接到电话了，我却没有
已经有4人回复
-
准备明年的基金了
已经有5人回复
-
今年的WR进展到哪一步了？
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
求助2013年-2011年，武汉哪所学校的实验室合成过单价3万/克（或毫升）的剧毒成品？
+5/3265
-
今年面上没有消息，想来是没了，散个金攒攒人品，为明年祈福
+1/624
-
华东理工大学 生物工程学院 生物反应器工程国家重点室 博士后招聘
+1/84
-
实验室用简易工装、夹具、模具、测试板
+1/82
-
Postdoctoral fellow in Medicinal Chemistry, McGill University
+1/72
-
西湖大学拓扑光学、非厄米光学、太赫兹方向博士后招聘（长期有效）
+5/60
-
西湖大学物理学、光学、电子信息方向博士生招生（长期有效）
+5/50
-
招收博士后（神经退行性疾病、认知障碍疾病）
+1/34
-
大连理工大学化学学院界面介导黏附与组装课题组招收2027级推免生/考研生
+1/32
-
大连理工大学化学学院界面介导黏附与组装课题组招收推免生
+1/31
-
深圳大学 招收2027级推免研究生 (金属材料/ 3D打印/ 氢能等方向)
+1/26
-
直博夏令营--北京理工大学/国家人工智能学院AI交叉学科联合招生
+1/23
-
西安交通大学补亚忠课题组招收申请考核制博士生
+2/14
-
浙江大学医学院附属妇产科医院张霈婧课题组招科研助理——生物信息方向
+2/12
-
宁波诺丁汉大学新能源与环境材料技术研究中心招收2027年春季入学博士生
+1/7
-
四川大学周加境课题组招博士研究生/研究助理（自组装材料/生物医学方向）
+1/6
-
准现货/原装全新超微整机8x NVIDIA HGX A100 80G
+1/6
-
武汉大学海优团队博士生招聘（微纳光子学与集成光子学）
+5/5
-
（Scopus期刊）Biofunctional Materials《生物功能性材料》招募青年编委及征稿
+1/5
-
中山大学招收2027级硕士研究生（增材制造/机器学习）
+1/1
5楼2014-02-18 04:21:59
★ ★ ★ ★ ★ ★
xiaolobster: 金币+5, ★有帮助 2013-11-29 09:03:39
xzhdty: 金币+1, 谢谢参与 2013-11-30 10:51:42
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
3楼2013-11-29 09:04:49
4楼2013-11-29 09:23:25