| 查看: 1958 | 回复: 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，但是都没能给出正确的结果。 求高手指点！ |
回复此楼» 猜你喜欢
327求调剂
已经有11人回复
-
08工科求调剂290分
已经有11人回复
-
一志愿985初试354分生物调剂
已经有3人回复
-
一志愿2110，化学学硕310分，本科重点双非求调剂
已经有13人回复
-
353求调剂
已经有8人回复
-
308求调剂
已经有16人回复
-
334求调剂
已经有14人回复
-
一志愿0807 数一英一 313 有没有二轮调剂
已经有7人回复
-
复试调剂，一志愿郑州大学材料与化工289分
已经有14人回复
-
086000调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
» 抢金币啦！回帖就可以得到:
-
安徽农业大学李亚辉教授诚招 化学、农药学、植保类研究生，接收调剂
+5/589
-
攀枝花学院大量招收调剂：材料科学与工程(学硕：0805)和材料工程（专硕：0856）。
+1/187
-
重庆长江师范学院材料工程方向接收材料类和机械类考生调剂硕士研究生
+1/145
-
湖北大学生科院邢琼老师课题组招硕士/博士,目前还有调剂名额！
+1/41
-
安徽工业大学工程研究院接受调剂
+1/39
-
武汉纺织大学学硕调剂
+2/36
-
西安工程大学材料学院杨杰课题组招收考研调剂生
+1/35
-
药物化学、材料与化工专业还有大量名额，速报名！！！
+1/13
-
长沙理工大学土木水利招收专硕调剂生
+1/11
-
速报名！！！材料、化学等相关专业均可调剂
+1/10
-
【上海调剂】985联合培养！别错过能送你去丹麦/交大/同济/中科院的神仙导师！
+1/9
-
济南大学物理科学与技术学院粒子物理方向招收调剂研究生
+1/7
-
河南理工大学化学化工学院招收材料、化工类专业调剂研究生
+1/7
-
山东交通学院航运学院关于2026年硕士研究生招生考试调剂复试有关事项的公告
+1/5
-
欢迎调剂攀枝花学院化学工程专业，资源与环境专业，线上复试
+1/5
-
济南大学前沿交叉科学研究院王金刚教授课题组接受化学、化学工程调剂
+1/5
-
湖南大学白玉罡课题组招收2026年入学博士生（第二批补录）
+1/4
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/4
-
河北大学招收调剂生
+1/4
-
武汉大学国家特聘青年教授侯绍聪课题组招收多名物理/化学/材料/电子/电气等专业博生生
+1/2
5楼2014-02-18 04:21:59
★ ★ ★ ★ ★ ★
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
