声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1740|回复: 4

[综合讨论] 解非线性方程组收敛不了怎么解决。

[复制链接]
发表于 2011-7-13 16:37 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
function x=maafun1
format long
x0 = [0;0;0;20000;14000;0;0;0];
options=optimset('MaxFunEvals',200000,'MaxIter',600000);
[x,fval] = fsolve(@f,x0,options) ;
function F=f(x)
k=5000
ga=19475
ei=415019052
p=20000
h=172-20000/k
ae=0
ea=200000*796.497
a=(p*(p+ga)/ei/ga)^0.5
b=ga/(ga+p)
c=1000
hc=2
F =[k*x(2)-ea*(((h-ae-x(2))^2+(x(1))^2)^0.5-h)/h*(1/(1+(x(3))^2))^0.5;
x(1)-(h-ae-x(2))*x(3);
(x(4)*(1-cos(a*h))-c*h)/p+(p*(c-k*x(2)*x(3))+ga*c)*sin(a*h)/(ga*a*p)-x(1);
(c*cos(a*h)/p+a*b*x(5)*sin(a*h)/p-c/p)-x(7);
c*cos(a*h)/p+(k*x(2)*x(3)*(1-cos(a*h))*p-(ga+p)*c+a*ga*x(4)*sin(a*h))/(ga+p)/p-x(7);
-x(5)*cos(a*h)/p+c*sin(a*h)/a/b/p+x(5)/p-c*h/p-x(6);
x(4)+x(5)-c*(h+hc+h)-p*x(8);
x(1)+x(6)+x(7)*hc-x(8)]


解一个多元非线性方程组   'MaxFunEvals',200000,'MaxIter',600000  迭代次数已经取得非常大了,但是还是收敛不了。我这方程检查过很多遍不会错的。请问要怎么解决。希望高手给予帮助,谢谢。
结果如下:F =

  -0.00196998854290
  -0.11934246867301
   2.67324300248162
  -0.04637127539383
  -0.04497784503437
   2.38019764407663
  -0.01351242000237
  14.04077335942650

Maximum number of function evaluations reached:
increase options.MaxFunEvals.

ans =

  1.0e+005 *

  -0.00013232140533
   0.00000519156024
  -0.00000078294316
  -2.46425742461584
  -2.75264386385792
  -0.00015561707175
  -0.00000074942350
  -0.00042984506407
回复
分享到:

使用道具 举报

发表于 2011-7-13 22:29 | 显示全部楼层
1.以前没有相关经验, 又没太多时间去试, 抱歉
2.大约瞄了下LZ的主题, 好像都是类似的东西
3.感觉LZ没把相关的式子/背景列出, 可能有经验的人较不愿意出手
 楼主| 发表于 2011-7-14 08:55 | 显示全部楼层
回复 2 # ChaChing 的帖子

式子背景较为复杂,一时半会也说不清。是一个力学问题
就是根据本构关系,几何关系,平衡关系。得出来的。
其中的3,4,5,6式是由微分方程解出来的。已验证过不会有错。

帮我想想看还有什么办法

LS很热心,表示感谢  
发表于 2011-7-14 09:08 | 显示全部楼层
用1stOpt得到的一组结果:
x1        -0.541506717377644
x2        337.780291226606
x3        0.00318946358946136
x4        468459.791578956
x5        1080718.99003213
x6        60.6845552283401
x7        0.207945340129349
x8        60.5589390805543

评分

1

查看全部评分

 楼主| 发表于 2011-7-14 09:50 | 显示全部楼层
回复 4 # dingd 的帖子

估计是有两组解吧,条件没限制好   h-ae-x(2)  这项必须是大于0的。

1stOpt 软件我以前下过,盗版的好像用不了。。。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 16:24 , Processed in 0.053334 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表