24小时热门版块排行榜    

Znn3bq.jpeg
查看: 745  |  回复: 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、这个方程组该怎么求解,真心请教大神们解答
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
2866110dai(feixiaolin代发): 金币+5 2015-01-11 10:42:46
1stOpt求解:
CODE:
Parameter x(3);
Function  x(2)-bessk(0,x(1))/bessk(1,x(1))*sin(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)));

有三角函数,不加范围约束有多组解(x3):
1:
x1: 0.13239358527979
x2: -0.282153438833906
x3: -1.82641189392894

2:
x1: 0.132393585286807
x2: -0.282153438837244
x3: 4.45677341328588
3楼2015-01-10 21:15:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

feixiaolin

荣誉版主 (文坛精英)

优秀版主

bessel函数中,变量默认是复数。
x 换成 real(x) 试一试。
2楼2015-01-10 07:42:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 恳请有学校收留 +5 柯淮然 2026-04-12 5/250 2026-04-15 21:49 by noqvsozv
[考研] 291求调剂 +6 关忆北. 2026-04-14 6/300 2026-04-15 21:32 by 笑春风☆
[考研] 通信工程求调剂!!! +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 一志愿华中农业071010,320求调剂 +11 困困困困坤坤 2026-04-14 12/600 2026-04-15 18:58 by zhuwenxu
[考研] 材料专业344求调剂 +17 hualkop 2026-04-10 22/1100 2026-04-14 16:21 by sxdj2
[考研] 305求调剂 +8 玛卡巴卡boom 2026-04-11 8/400 2026-04-14 09:04 by pengliang8036
[考研] 机械还有还有名额吗?太难了 +8 笑笑袁 2026-04-10 8/400 2026-04-14 08:44 by screening
[考研] 一志愿哈工大 085600 277 12材科基求调剂 5+5 chenny174 2026-04-10 37/1850 2026-04-14 07:39 by Abskk
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +8 慕绝cc 2026-04-09 8/400 2026-04-13 14:08 by 张zhihao
[考研] 生物学调剂 +11 小冉要努力 2026-04-10 13/650 2026-04-13 11:46 by 电化学及催化
[考研] 一志愿华中农微生物,288分,三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 326求调剂 +6 Shansyn 2026-04-10 6/300 2026-04-12 09:46 by hammer3
[考研] 药学专硕调剂 +8 ? 一路生?花? 2026-04-10 10/500 2026-04-11 21:21 by zhouxiaoyu
[考研] 求调剂 +11 翩翩一书生 2026-04-09 11/550 2026-04-11 19:57 by 逆水乘风
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 085400 328分 求调剂 +10 喂你一个大橙子 2026-04-09 14/700 2026-04-11 19:53 by lqspecial
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 调剂 化学 307 +21 73372112 2026-04-09 23/1150 2026-04-10 23:53 by wj165256
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
信息提示
请填处理意见