24小时热门版块排行榜    

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

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的回帖

whyhow

铁杆木虫 (著名写手)

带你飞翔

引用回帖:
3楼: Originally posted by 雁南飞6698 at 2012-07-09 16:42:27
哥哥,您说的那些我真心不懂,更不会手动计算。只希望重赏之下必有勇夫帮我代劳...

你是有两千多个币币的富人, 拿十个出来好意思说重赏嘛!
青春有千万种,却没有一种可以重来
5楼2012-07-09 17:42:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雁南飞6698

至尊木虫 (著名写手)

恶人谷谷主

引用回帖:
5楼: Originally posted by whyhow at 2012-07-09 17:42:40
你是有两千多个币币的富人, 拿十个出来好意思说重赏嘛!...

2000个在论坛只能算中下等呵呵
成功,源自执著
6楼2012-07-09 17:51:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meihongjiang

金虫 (正式写手)

请以矩阵表示来求解
7楼2012-07-09 19:46:04
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 雁南飞6698 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 今天基金会出结果吗?20260819 +13 kkkl_v 2026-08-19 13/650 2026-08-19 22:37 by 宝贝虫子
[基金申请] 时间戳变了,能看出什么问题? +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见