24小时热门版块排行榜    

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

羽化成鹏

新虫 (初入文坛)

引用回帖:
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的回帖
查看全部 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程307,求调剂 +3 我要燃烧你的梦 2026-03-08 3/150 2026-03-08 16:03 by 斩魂滴兔子!
[教师之家] 交大前校长王树国:现在最先进的科技并不在大学实验室,而是在企业研究院 +3 zju2000 2026-03-08 4/200 2026-03-08 14:55 by Equinoxhua
[考研] 【求调剂】293分环境工程求调剂材料/化工,服从调剂,抗压能力强! +11 xiiiia 2026-03-04 11/550 2026-03-08 14:28 by 30660438
[考研] 诚求调剂,323分 有专利/科研/实习经历 +6 橙子cyx 2026-03-06 8/400 2026-03-07 19:06 by 新篇章DFSS
[考研] 一志愿211,0860总分286食品工程专业求调剂 +6 林林Winnie 2026-03-05 6/300 2026-03-06 20:06 by 我想上研究生.
[考研] 266材料化工求调剂 +5 哇塞王帅 2026-03-03 5/250 2026-03-06 16:44 by 学员FnSWZj
[考研] 0703 东华大学 理学化学均过 a 线,总分 281 求调剂 +6 香香 Lu 2026-03-03 11/550 2026-03-06 15:37 by @飒飒飒飒
[考研] 不限学校专业的调剂同学看过来 +5 啊摆啊摆 2026-03-05 9/450 2026-03-06 12:06 by 啊摆啊摆
[考博] 2026申博自荐 六级440电催化方向 +4 樱落成影花成双 2026-03-05 4/200 2026-03-06 07:57 by L135790
[考研] 材料085600 303求调剂 +7 1bygone 2026-03-04 7/350 2026-03-05 09:39 by oxidpl
[论文投稿] 主编终审卡者不动 100+4 Stray2021 2026-03-03 4/200 2026-03-05 08:04 by bdhaoxiang
[考研] 322分 085600求调剂,有互联网+国金及主持省级大创经历 +6 熊境喆 2026-03-04 6/300 2026-03-04 20:32 by kakakapanpan
[考研] 0857求调剂 +4 CDDD1013 2026-03-04 4/200 2026-03-04 16:08 by zhukairuo
[考研] 264求调剂 +3 thext 2026-03-03 3/150 2026-03-04 09:29 by 热情沙漠
[考研] 289求调剂 +7 BrightLL 2026-03-02 9/450 2026-03-03 16:02 by tgxtgxtgx9
[考研] 一志愿东北大学化学314分求调剂 +5 lr1212.. 2026-03-02 5/250 2026-03-03 15:26 by xingtian2025
[考研] 299求调剂 +5 kkcoco25 2026-03-02 9/450 2026-03-03 12:55 by 公瑾逍遥
[考研] 材料工程专硕283求调剂 5+8 ,!? 2026-03-02 10/500 2026-03-03 12:01 by EBSD
[考研] 338求调剂 +5 18162027187 2026-03-02 6/300 2026-03-03 07:52 by njzyff
[考研] 化工京区271求调剂 +7 11ing 2026-03-02 7/350 2026-03-03 07:30 by 利好利好.
信息提示
请填处理意见