| 查看: 1904 | 回复: 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本,11408一志愿中科院277分,曾在中科院自动化所实习
已经有7人回复
26申博
已经有3人回复
求调剂,总分315,考的生物医药,一志愿湖南师范大学。调剂到任何专业都可以
已经有3人回复
材料专硕306英一数二
已经有7人回复
材料工程专硕274一志愿211求调剂
已经有6人回复
一志愿南京大学,080500材料科学与工程,调剂
已经有4人回复
332求调剂
已经有5人回复
271求调剂
已经有14人回复
275求调剂
已经有4人回复
302求调剂
已经有8人回复
» 本主题相关价值贴推荐,对您同样有帮助:
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位(化学化工相关)
+5/445
物理学 调剂
+1/85
物理学 调剂
+1/84
浙江农林大学招收调剂
+1/81
上海交大化院功能大分子团队招2026年硕士、联培、博士后
+1/39
2026年重庆工商大学管理科学与工程专业学术硕士招生调剂,风景总在柳暗花明处
+1/36
关于碳量子点抗菌的研究
+1/35
中科院生态环境研究中心国重实验室招聘客座研究生1-2名
+1/31
招硕士调剂生
+2/26
国家级领军人才团队高分子材料方向2026年博士研究生招生
+1/16
211/双一流石河子大学化学化工学院--电催化,CO2转化
+1/16
博士后招聘(高薪40万+)
+1/7
教育部长江学者和创新团队发展计划”入选团队招收 材料,化学与化工博士研究生
+1/7
论文投出,祈福,求好运
+1/7
课题组招收环境及相关专业调剂硕士研究生(欢迎优秀学生加入)
+1/6
南方科技大学生命科学学院基础免疫与微生物学系招收博士后
+1/4
层流压差式流量测控产品在大气运维采样流量校准中的应用与优势
+1/1
085404 270求调剂,b区求导师收留接受跨专业
+1/1
中国科学院地球化学研究所招科研助理
+1/1
层流压差式流量测控产品在环保与分析仪器领域的应用与优势
+1/1
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













回复此楼