24小时热门版块排行榜    

查看: 1620  |  回复: 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的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +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的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
独孤神宇: 金币+5, 鼓励交流 2021-01-13 20:41:29
羽化成鹏: 金币+30, 有帮助, 答案正确的 2021-01-14 15:14:08
使用fsolve函数,可以得到
X =       5.7508
Y =       4.6978
Z =       35.247
残值
Xerr=   0.00048828
Yerr=-0.00097656
Zerr=-0.00048828
QQ:2120156492
5楼2021-01-13 20:04:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

羽化成鹏

新虫 (初入文坛)

引用回帖:
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

答案是对的。可否附上代码?是否设置呢初始值?
谢谢。
6楼2021-01-14 15:13:26
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 羽化成鹏 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿安徽大学材料工程专硕313分,求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 5/250 2026-03-13 22:06 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 3/150 2026-03-13 15:14 by userper
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 纺织、生物、化学、材料相关专业招生了 +4 耶耶业 2026-03-09 7/350 2026-03-12 19:05 by Equinoxhua
[考研] 081200计算机科学与技术数一英一306求调剂 +3 intankt 2026-03-08 3/150 2026-03-11 10:50 by dandan413
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子!
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考研] 308求调剂 +4 是Lupa啊 2026-03-08 7/350 2026-03-10 18:25 by 清风月
[考研] 085602化工求调剂 +7 董boxing 2026-03-10 7/350 2026-03-10 17:07 by BruceLiu320
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
[考研] 320求调剂 +4 魏zy 2026-03-08 4/200 2026-03-09 16:14 by ruiyingmiao
[考研] 070300化学求调剂 +5 扑风铃的猫 2026-03-08 10/500 2026-03-09 13:58 by macy2011
[考研] 337求调剂 +3 睡醒,。 2026-03-09 3/150 2026-03-09 10:02 by 求调剂zz
[考研] 083000环境科学与工程调剂 +5 加油呀fxy 2026-03-07 6/300 2026-03-08 21:07 by sdgdsdgd
[考研] 347求调剂 +4 浮云满足 2026-03-07 4/200 2026-03-08 16:46 by 星空星月
信息提示
请填处理意见