| 查看: 2202 | 回复: 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，但是都没能给出正确的结果。 求高手指点！ |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有78人回复
-
filecode=后面第一个是大写字母
已经有8人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有13人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有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/396
-
北京信息科技大学机电工程学院硕士研究生招生线上宣讲会-8月16日
+5/250
-
[求助]请问水下焊接技术有啥进步？
+1/190
-
西安交大马伟/王亚洲课题组诚聘联培博士和硕士学生多名-有机生物电子学方向
+1/180
-
山东征女友，坐标济南
+1/161
-
北京信息科技大学机电工程学院硕士研究生招生线上宣讲会——8月16日
+1/81
-
中国科学院空天院精密干涉测量团队招聘光学/算法/机械方向博士后/特别研究助理
+1/80
-
中国科学院空天院精密干涉测量团队招聘光学/算法/机械方向博士后/特别研究助理/工程师
+1/78
-
广西大学 化学化工学院 新入职教师招生 化学、化学工程与技术、材料化工专业
+1/47
-
广东医科大学博士后招聘-基础医学，生信等方向
+1/38
-
天津大学环境学院童银栋教授招收博士后（生源要素循环/温室效应/新污染物方向）
+1/30
-
华南理工大学招2027年生物材料、药物递送申请考核制博士生
+1/24
-
西北工业大学 生科院 秦绪军课题组诚招博士后、博士/硕士研究生【长期有效】
+1/13
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/4
-
UTS 招收2027年入学博士生
+1/3
-
中国科学院生物物理研究所田勇课题组招聘联合培养生
+1/3
-
南京医科大学-冯璋课题组-招收2027级生化与分子生物学专业博士和硕士研究生
+1/3
-
河海大学海洋环境污染与生态毒理学研究组招收2027级硕士生 及博士研究生
+1/3
-
海南大学化学院—功能分子器件团队2027博士招生+博士后招聘
+1/2
-
国家级人才团队诚聘生物信息相关专业的科研助理
+1/1
6楼2015-05-19 14:17:14
★ ★ ★ ★ ★ ★
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
