| 查看: 1871 | 回复: 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，但是都没能给出正确的结果。 求高手指点！ |
回复此楼» 猜你喜欢
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有9人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有11人回复
-
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了
已经有7人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有8人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有9人回复
-
有院领导为了换新车，用横向课题经费买了俩车
已经有10人回复
-
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有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/987
-
湖南大学材料院陶益杰老师招收2026年联合培养硕士学生一名
+1/182
-
西南科技大学曹克课题组招收2026级申请考核制有机化学博士研究生
+1/176
-
2025难忘的时刻
+1/82
-
欢迎报考南京农业大学植物环境适应课题组课题组2026级博士生。
+1/76
-
中国科学院深圳先进技术研究院——招聘博士后
+3/56
-
有没有人做过这种结构的顺式体向反式体的转化？
+1/36
-
同济大学脑机智能团队脑机接口方向招生招聘
+1/31
-
同济大学脑机智能团队脑机接口方向招生招聘
+1/30
-
澳大利亚麦考瑞大学(Macquarie University)国际博士硕士全额奖学金-计算机-26年中开学
+1/19
-
新加坡 南洋理工大学- 智能光子/ 传感 PHD 全奖一名 2026 - 8 月入学
+1/17
-
澳门大学生物医学影像实验室诚招博士生（2026秋季入学）
+1/10
-
26储能博士申请自荐
+1/6
-
美国密苏里大学“柔性电子”课题组诚招博士研究生
+1/3
-
美国密苏里大学“柔性电子”课题组诚招博士研究生
+1/2
-
【经验分享】CRISPR基因敲除细胞系构建全流程踩坑指南——从递送方式选择到克隆筛选
+1/2
-
澳科大招收2026年秋季入学生物材料方向全奖博士研究生（3月5日截止）
+1/1
-
自荐：大模型ai辅助论文阅读软件：EasyReader论文易读
+1/1
-
岭南大学（香港）诚招材料方向研究助理/RA/博士后
+1/1
-
2026申博
+1/1
7楼2015-05-19 14:20:12
★ ★ ★ ★ ★ ★
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