| 查看: 1913 | 回复: 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,但是都没能给出正确的结果。 求高手指点! |
» 猜你喜欢
281求调剂(0805)
已经有15人回复
085600材料与化工调剂 324分
已经有8人回复
能源材料化学课题组招收硕士研究生8-10名
已经有10人回复
267一志愿南京工业大学0817化工求调剂
已经有8人回复
274求调剂
已经有4人回复
328求调剂,英语六级551,有科研经历
已经有10人回复
328求调剂,英语六级551,有科研经历
已经有7人回复
材料专业求调剂
已经有5人回复
一志愿天津大学化学工艺专业(081702)315分求调剂
已经有10人回复
344求调剂
已经有7人回复
» 本主题相关价值贴推荐,对您同样有帮助:
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人回复
» 抢金币啦!回帖就可以得到:
持续创业者,A8离异男征婚,事业遇到瓶颈,寻找创业生活伴侣
+1/496
山东征女友,坐标济南
+1/165
化学化工学院 招收化工、化学、材料等相关方向研究生(学硕、专硕都有调剂名额)
+1/91
重庆医科大学神经科学中心/病生与生理系诚邀海内外优秀博士加盟博士后队伍(长期)
+1/81
浙江农林大学招收调剂
+1/80
博士后出站经验请教
+1/78
山东师范大学有机化学专业胡忠燕老师课题组招收2026届硕士研究生以及调剂生
+1/40
河南工学院诚聘优秀凝聚态物理和电子科学与技术相关专业硕博研究生
+2/40
0854电子信息调剂,闽南师大光电芯片研发实验室,柯少颖教授课题组,集成电路
+1/37
找工作经验求助
+1/32
Call for paper: [Fermentation – IF 3.3]
+2/22
山东省植物基蛋白技术创新中心课题组联合培养硕士研究生招募启事
+1/12
2026年西南科技大学功能涂层课题组简介
+1/11
【2026博士招生/博后招聘】北京航空航天大学潘彪课题组——AI芯片设计方向
+1/11
【博士招生】武汉科技大学招收材料、化工、环境、冶金类2026年“申请-考核”制博士生
+1/8
新加坡国立大学药学系化学生物学课题组招PhD
+1/8
招收化工与材料学科点研究生
+1/4
南林化工院招申请审核制博士
+1/3
福建师范大学海峡柔性电子学院招收2026调剂硕士
+1/2
红外发射率测试(液体,2mm的小块体)/红外热成像发射率设定
+1/1
★ ★ ★ ★ ★ ★
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
5楼2014-02-18 04:21:59
6楼2015-05-19 14:17:14
7楼2015-05-19 14:20:12













回复此楼
