| 查看: 1951 | 回复: 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，但是都没能给出正确的结果。 求高手指点！ |
回复此楼» 猜你喜欢
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有9人回复
-
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有17人回复
-
070300化学学硕311分求调剂
已经有15人回复
-
材料调剂
已经有13人回复
-
材料与化工专硕306分找合适调剂
已经有19人回复
-
化学工程调剂289
已经有19人回复
-
化工学硕 285求调剂
已经有22人回复
-
298求调剂
已经有6人回复
-
368化学求调剂
已经有3人回复
-
326求调剂
已经有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/1080
-
安徽宿州学院 085600材料与化工(专业学位) 2026年研究生招生调剂信息
+1/85
-
延安大学长江学者团队招收材料、化工相关领域的调剂研究生数名
+1/43
-
欢迎生物医学工程、药学、生物、冶金、计算机等专业有班长团支书经历的同学-上海
+1/28
-
杭州医学院张逢质教授课题组招收药学硕士研究生若干
+1/26
-
陕西科技大学轻工学院张老师接收化工、材料、轻工类或相近专业的研究生调剂
+5/25
-
北京理工大学（珠海）招收2026级考核制博士
+1/21
-
吉大生科博一博士联培自荐-已发表+在投一作SCI3篇-踏实肯干能加班
+1/21
-
福建农林大学计算机科学与技术学硕招生
+1/16
-
★上海★上海工程技术大学·环境与资源创新中心（ICER）接收化工材料环境等学科调剂生
+1/12
-
欢迎调剂佛山大学环境工程（专硕）专业研究生
+1/6
-
国家大人才课题组招收化学、材料化工类2026年调剂生
+1/5
-
北京高校副校长团队招收机械类，环境类学硕和专硕
+1/5
-
杭州师范大学-浙江省湿地智慧监测与生态修复重点实验室团队硕士调剂公告
+1/4
-
大连工业大学“轻工技术与工程”石海强教授课题组招收硕士研究生
+1/4
-
福建理工大学材料学院招收专业代码08开头学术型硕士研究生
+1/3
-
齐齐哈尔大学栗洪彬老师招收2026年纺织工程调剂研究生(学术型硕士2名，专业硕士1名）
+1/2
-
河北大学招收等离子体物理调剂生
+1/2
-
南开大学环境学院 环境微生物方向 2026年博士招生（已招满）
+1/1
-
广东省五邑大学MOFs催化功能材料课题组招收化学/材料背景调剂生
+1/1
3楼2013-11-29 09:04:49
★ ★ ★ ★ ★ ★
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
4楼2013-11-29 09:23:25
5楼2014-02-18 04:21:59
