24小时热门版块排行榜    

查看: 1606  |  回复: 6
【悬赏金币】回答本帖问题,作者羽化成鹏将赠送您 370 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

羽化成鹏

新虫 (初入文坛)

[求助] matlab解三元高次方程组。400金币。 已有2人参与

求解下面的三元高次方程组的解。
使用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)
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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)
QQ:2120156492
7楼2021-01-14 22:14:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
求解结果在 [-1,-1,0] 是整体最优的。

残差较大可能是方程组本身的问题
数值计算
2楼2021-01-13 16:04:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

羽化成鹏

新虫 (初入文坛)

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

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

大佬 系数已经变了。您仔细看一下
3楼2021-01-13 16:13:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

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

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

除非修改模型
数值计算
4楼2021-01-13 16:26:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] A区一本交叉课题组,低分调剂,招收机械电子信息通信等交叉方向 +53 lisimayy 2026-03-04 66/3300 2026-03-09 18:27 by wilburt
[考研] 标题:捡漏预警|08工科/09农学调剂!英语要求低,过线即有机会! +3 马超放烟花 2026-03-07 3/150 2026-03-09 16:27 by 17862602468
[考研] 0703化学求调剂 +4 很老实人 2026-03-09 4/200 2026-03-09 13:38 by macy2011
[考研] 296求调剂 +4 Xinyu Wu311 2026-03-09 4/200 2026-03-09 11:51 by 勇敢太监王公公
[考研] 材料与化工,291分,求调剂 +13 咕噜咕噜123123 2026-03-04 15/750 2026-03-09 09:32 by 30660438
[考博] 2026博士申请 +4 起泡酒 2026-03-08 4/200 2026-03-09 08:22 by 0611517sll
[考研] 085701环境工程专硕求调剂 +5 幼稚园大班打手 2026-03-08 5/250 2026-03-09 07:07 by houyaoxu
[考研] 297求调剂 +3 胡达灵 2026-03-05 5/250 2026-03-08 23:13 by 清风月
[教师之家] 交大前校长王树国:现在最先进的科技并不在大学实验室,而是在企业研究院 +4 zju2000 2026-03-08 6/300 2026-03-08 19:15 by zju2000
[考研] 085701环境工程专业,初试305,均过国家A区线 +7 卡卡来了@ 2026-03-07 8/400 2026-03-08 06:40 by 刘兵
[考研] 材料专硕323求调剂 +5 李白26 2026-03-07 5/250 2026-03-07 16:46 by njzyff
[考研] 304求调剂 +4 52hz~~ 2026-03-05 5/250 2026-03-07 15:47 by lature00
[考研] 276求调剂 本科天津工业大学 一志愿北京化工大学0855 +4 路lyh123 2026-03-06 6/300 2026-03-07 11:37 by shdliugang
[考研] 材料277分求调剂 +13 饭饭星球 2026-03-04 14/700 2026-03-06 16:10 by @飒飒飒飒
[考研] 275求调剂 +4 大爆炸难民 2026-03-06 5/250 2026-03-06 09:21 by guoweigw
[考研] 化工282求调剂一志愿211 +5 NA0912 2026-03-05 6/300 2026-03-05 20:10 by 2735147993
[考研] 求调剂院校 +6 云朵452 2026-03-02 12/600 2026-03-04 23:17 by 云朵452
[基金申请] 没有青基直接申请面上,感觉自己疯了 +6 kevin63t 2026-03-02 12/600 2026-03-04 22:23 by kevin63t
[考研] 一志愿985材料与化工 326分求调剂 +3 Hz795795 2026-03-04 3/150 2026-03-04 20:54 by wutongshun
[考研] 复试调剂 +4 杰尼龟aaa 2026-03-03 4/200 2026-03-03 23:06 by zhukairuo
信息提示
请填处理意见