24小时热门版块排行榜    

查看: 2481  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

雁南飞6698

至尊木虫 (著名写手)

恶人谷谷主

[求助] 请教一个四元二次方程的极大值

Y=87.299+2.44X1+0.969X2+2.6X3+0.172X4+0.3X1X1+0.398X1X2-0.419X1X3-0.261X1X4+0.675X2X2+ 0.618X2X3+0.401X2X4-1.2X3X3-0.098X3X4-1.575X4X4

方程如上,就是要求出这个四元二次方程的极大值,X1,X2,X3,X4取值在-2~+2之间,不仅要极大值,还有方程取极大值时这四个自变量的值。听说matlab可以很方便地解出来,可惜我不会使用那个大型软件。如果能给出过程更好。结果的形式不要以图片的形式发上来,谢谢~我是化学专业的,不会呵呵
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

成功,源自执著
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

li_haizeng

银虫 (小有名气)

用Maplea求解
Y := 87.299+2.44*X1+.969*X2+2.6*X3+.172*X4+.3*X1*X1+.398*X1*X2-.419*X1*X3-.261*X1*X4+.675*X2*X2+.618*X2*X3+.401*X2*X4-1.2*X3*X3-0.98e-1*X3*X4-1.575*X4*X4
> maximize(Y, X1 = -2 .. 2, X2 = -2 .. 2, X3 = -2 .. 2, X4 = -2 .. 2, location);
print(`output redirected...`); # input placeholder
  101.4987647, {[

    {X1 = 2., X2 = 2., X3 = 1.244888877, X4 = 0.1047621873},

    101.4987647]}
输出结果是:
101.4987647, {[{X1 = 2., X2 = 2., X3 = 1.244888877, X4 = .1047621873}, 101.4987647]}
8楼2012-07-10 15:15:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

xuyx_78

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
这个用matlab计算,要选择一个好的初值不容易.不如直接用求驻点的方法,手工求解.这是二次函数,对每个变量求偏导且令偏导数等于零得到关于四个变量的四元一次方程,手工就可以求得精确的驻点.再比较驻点与边界处的函数值就可以求出要求的.
2楼2012-07-09 16:37:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雁南飞6698

至尊木虫 (著名写手)

恶人谷谷主

引用回帖:
2楼: Originally posted by xuyx_78 at 2012-07-09 16:37:11
这个用matlab计算,要选择一个好的初值不容易.不如直接用求驻点的方法,手工求解.这是二次函数,对每个变量求偏导且令偏导数等于零得到关于四个变量的四元一次方程,手工就可以求得精确的驻点.再比较驻点与边界处的函数 ...

哥哥,您说的那些我真心不懂,更不会手动计算。只希望重赏之下必有勇夫帮我代劳
成功,源自执著
3楼2012-07-09 16:42:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
雁南飞6698: 金币+10, ★★★★★最佳答案, 我喜欢直接的 2012-07-09 17:56:44
1stOpt求解,两句话:
CODE:
ParameterDomain = [-2,2];
MaxFunction 87.299+2.44*X1+0.969*X2+2.6*X3+0.172*X4+0.3*X1*X1+0.398*X1*X2-0.419*X1*X3-0.261*X1*X4+
            0.675*X2*X2+ 0.618*X2*X3+0.401*X2*X4-1.2*X3*X3-0.098*X3*X4-1.575*X4*X4;

目标函数值(最大): 101.498764681482
x1: 2
x2: 2
x3: 1.24488886605106
x4: 0.104762170420606
4楼2012-07-09 17:31:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:39 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:03 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 00:51 by 9i7PCSZyMzot
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:15 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:03 by 9i7PCSZyMzot
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 4pRyVGCVerAA 2026-09-19 4/200 2026-09-20 00:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:51 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:39 by 9i7PCSZyMzot
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[教师之家] 现代”学阀”该如何界定 +3 iaeyuan 2026-09-19 3/150 2026-09-19 20:52 by 没昵称呀
[基金申请] 国社科系统bug了,是不是要放榜了? +7 kynobel 2026-09-16 9/450 2026-09-19 14:52 by 接好运婷
[考博] 申博发邮件 +8 Lmengk 2026-09-15 9/450 2026-09-19 09:27 by 中间体超市
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 9/450 2026-09-18 18:17 by zapen
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
信息提示
请填处理意见