当前位置: 首页 > 计算模拟 >matlab解三元高次方程组。400金币。

matlab解三元高次方程组。400金币。

作者 羽化成鹏
来源: 小木虫 300 6 举报帖子
+关注

求解下面的三元高次方程组的解。
使用matlab进行汇编。

要求:没有初始值(需要输入初始值的情况不符合要求!)!!! 程序运行,即要直接得到最优结果!!!具有鲁棒性(方程中系数564080065.628041000,609249680.389719000,649348646.964450000发生调整,依然可以正常运行  )!!!


%%%
需要附上:
1,代码
2,结果 X Y Z
3,残差(验证结果)

400金币。
谢谢。


%%需求解方程:

  (X^2+Y^2+Z^2)^4 - 564080065.628041000 * (X^2+ Y^2+ 4* Z^2)  = 0
  ( (X+2)^2+Y^2+Z^2)^4 - 609249680.389719000  * ( (X+2)^2+ Y^2+ 4* Z^2) =0
  ( (X+2)^2+(Y+2)^2+Z^2)^4 - 649348646.964450000  * ( (X+2)^2+ (Y+2)^2+ 4* Z^2) =0


%%残差验证

%  残差;
a = (X^2+Y^2+Z^2)^4 - 564080065.628041000 * (X^2+ Y^2+ 4* Z^2)
b = ( (X+2)^2+Y^2+Z^2)^4 - 609249680.389719000  * ( (X+2)^2+ Y^2+ 4* Z^2)
c = ( (X+2)^2+(Y+2)^2+Z^2)^4 - 649348646.964450000  * ( (X+2)^2+ (Y+2)^2+ 4* Z^2)

% 残差比值
ratioa = 1/(X^2+Y^2+Z^2)^4 * 564080065.628041000 * (X^2+ Y^2+ 4* Z^2)
ratiob = 1/( (X+2)^2+Y^2+Z^2)^4 * 609249680.389719000  * ( (X+2)^2+ Y^2+ 4* Z^2)
ratioc = 1/( (X+2)^2+(Y+2)^2+Z^2)^4 * 649348646.964450000  * ( (X+2)^2+ (Y+2)^2+ 4* Z^2) 返回小木虫查看更多

今日热帖
  • 精华评论
  • 独孤神宇

    求解结果在 [-1,-1,0] 是整体最优的。

    残差较大可能是方程组本身的问题

  • 羽化成鹏

    引用回帖:
    2楼: Originally posted by 独孤神宇 at 2021-01-13 16:04:29
    求解结果在  是整体最优的。

    残差较大可能是方程组本身的问题

    大佬 系数已经变了。您仔细看一下

  • 独孤神宇

    引用回帖:
    3楼: Originally posted by 羽化成鹏 at 2021-01-13 16:13:39
    大佬 系数已经变了。您仔细看一下...

    这个不管怎么变,残差还是很大,你可以试一下。

    除非修改模型

  • hzlhm

    使用fsolve函数,可以得到
    X =       5.7508
    Y =       4.6978
    Z =       35.247
    残值
    Xerr=   0.00048828
    Yerr=-0.00097656
    Zerr=-0.00048828

  • 羽化成鹏

    引用回帖:
    5楼: Originally posted by hzlhm at 2021-01-13 20:04:56
    使用fsolve函数,可以得到
    X =       5.7508
    Y =       4.6978
    Z =       35.247
    残值
    Xerr=   0.00048828
    Yerr=-0.00097656
    Zerr=-0.00048828

    答案是对的。可否附上代码?是否设置呢初始值?
    谢谢

  • hzlhm

    引用回帖:
    6楼: Originally posted by 羽化成鹏 at 2021-01-14 15:13:26
    答案是对的。可否附上代码?是否设置呢初始值?
    谢谢。...

    x0=rand(3,1)*1000;
    [y,fval,h]=fsolve(@func,x0)
    X=y(1)
    Y=y(2)
    Z=y(3)

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓