| 查看: 1755 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
matlab对含贝塞尔函数的超越方程求根 已有1人参与
|
||
|
有一个含有贝塞尔函数的超越方程，用solve函数求解，输入matlab直接给出了一个复数结果，但是实际上超越方程是有很多个根的，不知道给出的这个解是第几阶的？怎么样才能得到它的第n（任意的）阶根呢？采用牛顿迭代法给定不同的初值可以计算得到不同的根，但是也不能确定是第几阶。不知道论坛里有没有朋友有更好的计算方法？函数就不贴了，下面是牛顿迭代法的matlab程序，函数见里面的f。请了解的朋友给点意见，感谢！ clc;clear;close all; r0=0.0125; r=r0; fre=1000; z=0; l=0.122; r0=0.0125; r=r0; rou0=1000; a0=1500; pi=3.1415926; k0=2*pi*f/a0; Zw=1800000*(1+i); %牛顿迭代法 syms f kr; f=subs('besselj(0,kr*r0)-kr*r0*i*Zw/(rou0*a0*k0*r0)*besselj(1,kr*r0)'); df=diff(f); kr0=0.1-0.1*i; e1=0.1e-10; e2=0.000005; N=100; k=1; while (k<N) kr=kr0; if abs(eval(f))<e1 fprintf('奇异!\nkr=%.6f\n迭代次数为:%d\n',kr0,k) break else kr1=kr0-eval(f)/eval(df); if abs(kr1-kr0)<e2 fprintf('kr=%.6f\n迭代次数为:%d\n',kr1,k) break else kr0=kr1; k=k+1; end end end if k>=N fprintf('失败\n') end |
回复此楼» 猜你喜欢
一志愿哈尔滨工业大学材料与化工方向336分
已经有5人回复
-
0856材料化工调剂 总分330
已经有9人回复
-
085602 化工专硕 338分 求调剂
已经有9人回复
-
352分 化工与材料
已经有4人回复
-
求调剂
已经有8人回复
-
材料求调剂 一志愿哈工大总分298分，前三科223分
已经有3人回复
-
材料求调剂一志愿哈工大324
已经有4人回复
-
085600 286分 材料求调剂
已经有4人回复
-
求调剂
已经有7人回复
-
275求调剂
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助matlab写变限积分
已经有4人回复
- 求助：如何求含贝塞尔函数的超越方程的复数解
已经有3人回复
- 有关于mathematics的问题，用mathematics对一些点拟合出一个椭圆方程
已经有24人回复
- [求助]求解极座标上的偏微分方程(Laplace’s Equation)!!!!
已经有6人回复
- 求教关于origin8.0 自定义函数拟合方面的问题
已经有5人回复
- 急！在用origin进行非线性拟合怎样创建新的拟合函数？
已经有14人回复
- 师兄师姐们，帮小弟看看这个 Laplace变换后的二阶常微分方程怎么解啊
已经有23人回复
- VS2012如何快速输入函数？
已经有6人回复
- origin 里 “三次样条函数拟合” 与“三次多项式拟合”是一回事吗？如图
已经有3人回复
- 用matlab求解含有bessel函数的方程
已经有5人回复
- 求解高次复杂函数
已经有11人回复
- 微分方程求解（通解中含hypergeom函数）
已经有17人回复
- 用matlab求bessel函数的根
已经有6人回复
- 一个高一函数题求解
已经有7人回复
- 【求助】复超越方程的解法
已经有19人回复
3楼2015-03-23 14:15:33
mathstudy
金虫 (正式写手)
- 数学EPI: 2
- 应助: 161 (高中生)
- 金币: 4320.1
- 散金: 38
- 红花: 16
- 帖子: 446
- 在线: 141.1小时
- 虫号: 2515489
- 注册: 2013-06-20
- 性别: GG
- 专业: 计算数学与科学工程计算
2楼2015-03-23 09:33:06
mathstudy
金虫 (正式写手)
- 数学EPI: 2
- 应助: 161 (高中生)
- 金币: 4320.1
- 散金: 38
- 红花: 16
- 帖子: 446
- 在线: 141.1小时
- 虫号: 2515489
- 注册: 2013-06-20
- 性别: GG
- 专业: 计算数学与科学工程计算
4楼2015-03-24 09:44:38
5楼2015-03-24 14:20:53
