| 查看: 1784 | 回复: 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，但是都没能给出正确的结果。 求高手指点！ |
回复此楼» 猜你喜欢
自荐读博
已经有4人回复
-
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有6人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有4人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
需要合成515-64-0，50g，能接单的留言
已经有3人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
-
临港实验室与上科大联培博士招生1名
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
» 抢金币啦！回帖就可以得到:
-
虚位以待，共探前沿：热烈欢迎报考北京工业大学纳米多孔材料课题组2026级博士研究生
+1/257
-
【宁德时代招聘】AI 物理学家
+1/172
-
中国科大化学与材料科学学院/苏州高研院刘东/熊宇杰教授团队诚聘博士后
+5/125
-
天津大学化学系吴立朋课题组申请考核制博士招生/博后招聘
+1/79
-
加拿大卡尔加里大学 量子通信和信息方向 硕士/博士招生
+1/60
-
深圳理工大学梁国进课题组招聘研究助理教授、博后多名（电化学储能方向）
+1/49
-
Analytical Science Advances（Wiley出版社）长期征稿中...
+1/37
-
盐湖所镁基储氢材料课题组招聘
+1/32
-
新年快乐
+2/32
-
江西师范大学化学与材料学院2026年博士研究生招生
+1/26
-
SCI，计算机相关可以写
+1/26
-
重庆大学前沿院，黄小洋教授课题组，招收2026年非均相催化方向学术博士2名
+1/15
-
天津大学化学系吴立朋课题组申请考核制博士招生/博后招聘-有机化学，金属有机
+1/14
-
北京理工大学计算材料化学课题组招聘2026年博士生/博后
+3/12
-
天津大学化学系吴立朋课题组申请考核制博士招生/博后招聘-有机化学，金属有机
+1/11
-
SCI计算机
+1/6
-
南京林业大学”申请-考核”制学术学位博士研究生招生
+1/5
-
中科院和北京工商大学招收2026博士/化学或生物背景
+1/2
-
河南师范大学植物生殖生物学科研团队博士招聘
+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