| 查看: 807 | 回复: 2 | |||
2866110dai金虫 (小有名气)
|
[求助]
包含贝塞尔函数的非线性方程组 已有1人参与
|
|
方程组及求解程序如下： f11=@(x)[x(2)-besselk(0,x(1))./besselk(1,x(1)).*sin(x(3)); %满足的方程组， x=x(1),y=x(2),phi0=x(3) x(1)-0.5*sin(pi+x(3)-pi/3); x(2)+(0.3-0.5*(1-cos(pi+x(3)-pi/3)))]; options = optimset('Display','final','TolFun',1e-10,'TolX',1e-10); X=fsolve(f11,[2,2,pi/3],options) ********************************************************************************************** 求出来的结果是虚数，Columns 1 through 2 0.146326396146141 + 0.572064042319110i 0.953776042482729 - 0.111051645262746i Column 3 0.855457842453923 - 0.993689273363200i 然而这是个实际问题，并且x(3)是角度值，解是存在的，而且一定是实数。 然后我改变初始点，如[1,1,pi/3],提示如下： No solution found. fsolve stopped because the last step was ineffective. However, the vector of function values is not near zero, as measured by the selected value of the function tolerance. <stopping criteria details> X = 0.313681166680663 0.435001517684952 0.535653117999125。 ************************************************************************************************* 问题：1、fsolve求解的是多元方程的 实数解，而我在初值[2,2,pi/3]求出来的却是虚数，这是怎么回事 2、在初值[1,1,pi/3],既然提示No solution found，怎么最后又有X值， 0.313681166680663 0.435001517684952 0.535653117999125。 3、这个方程组该怎么求解，真心请教大神们解答 |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有43人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有5人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有13人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有11人回复
-
读博的好处
已经有4人回复
-
filecode
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关于mathematics的问题，用mathematics对一些点拟合出一个椭圆方程
已经有24人回复
- [求助]求解极座标上的偏微分方程(Laplace’s Equation)!!!!
已经有6人回复
- 求教关于origin8.0 自定义函数拟合方面的问题
已经有5人回复
- 急！在用origin进行非线性拟合怎样创建新的拟合函数？
已经有14人回复
- 求助用origin自带的ExpAssoc函数进行非线性拟合
已经有8人回复
- 用matlab求解含有bessel函数的方程
已经有5人回复
- fortran中计算线性积分
已经有7人回复
- 求帮忙，关于含有bessel函数的方程的编程求解问题
已经有11人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2015-01-10 07:42:02
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
3楼2015-01-10 21:15:06
