马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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 |