| 查看: 1902 | 回复: 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人回复
-
267一志愿南京工业大学0817化工求调剂
已经有6人回复
-
293求调剂
已经有4人回复
-
285化工学硕求调剂（081700）
已经有7人回复
-
302求调剂
已经有8人回复
-
08工科 320总分 求调剂
已经有3人回复
-
290求调剂
已经有9人回复
-
271求调剂
已经有13人回复
-
本人考085602 化学工程 专硕
已经有16人回复
-
304求调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
» 抢金币啦！回帖就可以得到:
-
测试█TEM/ EPR/ XPS/PY-GCMS/TG-IR/XRF/BET/MIP/核磁/EA/ICP，VX: 761711562。
+1/85
-
山东师范大学有机化学专业胡忠燕老师课题组招收2026届硕士研究生以及调剂生
+1/43
-
华南师范大学（211）博士招生- 电子、自动化、机械、生物学、物理相关专业
+2/42
-
中科院生态环境研究中心国重实验室招聘客座研究生1-2名
+1/31
-
北京航空航天大学招生微流控与智能打印技术方向博士研究生
+1/29
-
纺织科学与工程、材料化工方向招收研究生
+1/28
-
厦门大学航空航天学院智能制造课题组招2026年申请审核制博士生1名
+1/22
-
南京林业大学-国家级青年人才团队 招2026级申请考核制博士（合成化学、植物化学方向）
+1/15
-
【2026博士招生/博后招聘】北京航空航天大学潘彪课题组——AI芯片设计方向
+1/13
-
【全奖招生】北师港浸大ESLAS实验室招收密码工程/网络安全/计算机视觉博士/博士后
+1/11
-
【博士招生】武汉科技大学招收材料、化工、环境、冶金类2026年“申请-考核”制博士生
+1/10
-
香港大学David Srolovitz教授课题组招聘「香江学者」博士后
+1/9
-
【第三轮招生，31号截止】澳科大诚招2026秋全奖博士研究生（纳米医学/生物材料方向）
+1/6
-
上海第二工业大学-朱大海课题组招生（过线就能调剂）
+1/5
-
PD-L1抗体与PD-1抗体有何差异？
+1/4
-
连发两篇 Nature！大样本单细胞测序解析认知、衰老与灵长类进化
+1/3
-
海南大学徐月山老师招生第二批博士名额2~3个，2026年9月份入学（高端设备开发方向）
+1/2
-
欢迎报考浙江工业大学埃克塞特联合学院硕士研究生
+1/1
-
澳洲维多利亚大学计算机类全额奖学金博士招生 （邮箱+微信可联系）
+1/1
-
中科院上海高等研究院光源科学中心孙天啸团队诚招2位助理研究员
+1/1
4楼2013-11-29 09:23:25
★ ★ ★ ★ ★ ★
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
5楼2014-02-18 04:21:59
