查看: 411|回复: 3|关注: 0

[已答复] 多元多次方程组求数值解

[复制链接]

新手

5 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
发表于 2019-3-14 21:33:32 | 显示全部楼层 |阅读模式
function M=answer(x)

        A=1;B=1/3;C=5.77;E1=210;E2=13.8;E3=210;v1=0.3;v2=0.25;v3=0.3;ri=54.3;r1=63.5;r2=77.39;r0=88.9;pj=69;p0=0;f=pi/6;
        M=[(1+v1)*2*(1-v1)*ri^2*r1*pj/E1/(r1^2-ri^2)-(1+v1)*(r1^3*(1-2*v1)+ri^2*r1)*x(1)/E1/(r1^2-ri^2)-((1+v2)*(1-2*v2)*x(1)*r1/E2+x(5)/r1);
            (1+v2)*(1-2*v2)*C*cot(f)/E2*(1-(1+x(1)/C/cot(f))*(x(2)/r1)^(B-1))*x(2)+x(5)/x(2)-((1+v2)*(x(2)^3*(1-2*v2)+r2^2*x(2))*x(3)/E2/(r2^2-x(2)^2)-(1+v2)*2*(1-v2)*r2^2*x(2)*x(4)/E2/(r2^2-x(2)^2));
            (1+v2)*2*(1-v2)*x(2)^2*r2*x(3)/E2/(r2^2-x(2)^2)-(1+v2)*(r2^3*(1-2*v2)+x(2)^2*r2)*x(4)/E2/(r2^2-x(2)^2)-((1+v3)*(r2^3*(1-2*v3)+r0^2*r2)*x(4)/E3/(r0^2-r2^2)-(1+v3)*2*(1-v3)*r0^2*r2*p0/E3/(r0^2-r2^2));
            C*cot(f)*((1+x(1)/C/cot(f))*(x(2)/r1)^(B-1)-1)-x(3);
            ((A*r2^2+x(2)^2*sin(f))*x(3)-(r2^2-x(2)^2)*C*cos(f))/r2^2/(A+sin(f))-x(4)];

  end
x0=[25;70;18;25;20];
    options=optimset('Display','iter');
    [x,fval,exitflag]=fsolve(@answer,x0,options)


请问解和初值有关系吗,我的出的解和预期结果不符,有其他的计算方法吗

论坛优秀回答者

权威

9151 麦片

财富积分



11

主题

1万

帖子

1994

最佳答案
  • 关注者: 518
发表于 2019-3-15 07:55:19 | 显示全部楼层
与初值有很大的关系,fsolve只给出初值附近的解。
你的预期是啥?
[x,fval,flag]=fsolve(@answer,[10,65,10,8,3700])试试,
x =

   1.0e+03 *

    0.0100    0.0638    0.0099    0.0078    3.6800


fval =

   1.0e-10 *

         0
    0.1046
   -0.0653
    0.0045
   -0.0001


flag =

     1

新手

5 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
 楼主| 发表于 2019-3-15 09:43:32 | 显示全部楼层
x2应该是70多

论坛优秀回答者

中级

932 麦片

财富积分


5001500


0

主题

2356

帖子

202

最佳答案
  • 关注者: 134
发表于 2019-3-15 11:59:23 | 显示全部楼层
参考下1stOpt计算结果(不需初值):

x1: 10.0117823187245
x2: 63.8456756868464
x3: 9.93950652996427
x4: 7.81727357258045
x5: 3679.98588141515
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

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