| 查看: 1897 | 回复: 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，但是都没能给出正确的结果。 求高手指点！ |
回复此楼» 猜你喜欢
材料与化工304求B区调剂
已经有7人回复
-
0703化学336分求调剂
已经有4人回复
-
[复试调剂]西南科技大学国防/材料导师推荐
已经有6人回复
-
化学工程321分求调剂
已经有12人回复
-
211本，11408一志愿中科院277分，曾在中科院自动化所实习
已经有4人回复
-
材料专硕326求调剂
已经有5人回复
-
东南大学364求调剂
已经有5人回复
-
国自科面上基金字体
已经有7人回复
-
药学383 求调剂
已经有4人回复
-
286求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
» 抢金币啦！回帖就可以得到:
-
西湖大学工学院智能微型机器人实验室博士后招聘启事
+4/236
-
哈工大本部-电化学新能源方向博士捡漏名额
+1/77
-
何时使用 CODA™ 科里奥利质量流量仪表- 艾里卡特（Alicat）
+2/52
-
教研论文SCI期刊投稿选刊
+1/34
-
0854电子信息调剂，集成电路，芯片方向，闽南师范大学光电芯片研发实验室
+1/31
-
澳大利亚西澳大学（UWA）张金强课题组招收全奖博士生（化工新能源方向）
+1/27
-
南开大学罗景山教授课题组招收2026年入学博士生（光/电催化方向）
+2/20
-
南京林业大学-国家级青年人才团队 招2026级申请考核制博士（合成化学、植物化学方向）
+1/16
-
【实战型】【生物医药】2026青岛大学招博士生 含少数民族骨干计划2名！！！
+1/11
-
2026年上海大学王亮课题组招收电催化或光催化研究背景的博士生
+2/10
-
西京学院土木水利 2026 级研究生招生相关说明
+1/8
-
聊城大学化院青椒招收化学相关专业调剂学生
+1/7
-
华中科技大学管理学院招聘社会用工（科研助理）1名
+1/7
-
211/双一流石河子大学化学化工学院-电池/电催化研究方向
+1/6
-
发paper
+1/5
-
北京大学环境学院张子帅研究员招收科研助理
+1/3
-
药学专硕一志愿苏大本科苏大寻求B区高校调剂名额
+1/2
-
085404 270求调剂，b区求导师收留接受跨专业
+1/2
-
西安交通大学生命学院智能化诊疗装备研究中心2026年博士生招生
+1/2
-
中国科学院地球化学研究所招科研助理
+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
