| 查看: 1087 | 回复: 7 | ||
liuhongbo木虫 (著名写手)
|
[求助]
两个可变参数 T, E; 两个自变量x,y;的高次多项式极值怎么求？ 已有1人参与
|
|
f=1/2*a*(x^2+y^2)+1/4*b*(x^4+y^4)+1/6*c*(x^6+y^6)+g*x*y-(x+y)*E a=g+a1*(T-T0) a1=1 b=c=1/3 g=0.5 T0=500 T,E均为可变参数，T的变化范围是0到700；E的变化范围为0到0.6E0，E0 = (6b^2/(25c))*sqrt(3*abs(b)/(10*c)). 需要求函数的极值，请问怎么用matlab实现呢？ 如果太复杂，自变量可限制x>0; -x<y<x 感谢感谢！ |
回复此楼» 本帖已获得的红花（最新10朵）
liuhongbo» 猜你喜欢
最失望的一年
已经有17人回复
-
拟解决的关键科学问题还要不要写
已经有8人回复
-
为什么nbs上溴 没有产物点出现呢
已经有6人回复
-
求推荐博导
已经有4人回复
-
存款400万可以在学校里躺平吗
已经有34人回复
-
求助一下有机合成大神
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
疑惑？
已经有5人回复
2楼2016-05-19 13:48:47
liuhongbo
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 19985.5
- 散金: 8
- 红花: 4
- 帖子: 2122
- 在线: 751.9小时
- 虫号: 638484
- 注册: 2008-10-27
- 专业: 功能陶瓷
3楼2016-05-19 18:50:30
4楼2016-05-19 20:18:58
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liuhongbo: 金币+15, ★★★★★最佳答案, merci! 2016-05-21 07:27:58
liuhongbo: 金币+5, ★★★★★最佳答案 2016-05-21 07:33:05
感谢参与，应助指数 +1
liuhongbo: 金币+15, ★★★★★最佳答案, merci! 2016-05-21 07:27:58
liuhongbo: 金币+5, ★★★★★最佳答案 2016-05-21 07:33:05
|
function [] = exercise() options=optimset('LargeScale','off','Display','iter','TolX',1e-4,'TolFun',1e-4,'MaxFunEvals',1e6,'GradObj','off','GradConstr','off','Algorithm','interior-point'); global T E A=[];b=[]; Aeq=[]; beq=[]; lb=[0 -Inf]; ub=[Inf Inf]; T=100;E=0.02; xy0=rand(2,1); [XYopt,fval,exitflag]=fmincon(@obj,xy0,A,b,Aeq,beq,lb,ub,[],options); disp(XYopt); function f=obj(xy0) global T E x=xy0(1);y=xy0(2); g=0.5;b=1/3;c=1/3;a1=1; T0=500; a=g+a1*(T-T0); f=1/2*a*(x^2+y^2)+1/4*b*(x^4+y^4)+1/6*c*(x^6+y^6)+g*x*y-(x+y)*E; 貌似你的问题里x和y是对称的 |
5楼2016-05-19 20:32:27
liuhongbo
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 19985.5
- 散金: 8
- 红花: 4
- 帖子: 2122
- 在线: 751.9小时
- 虫号: 638484
- 注册: 2008-10-27
- 专业: 功能陶瓷
6楼2016-05-20 06:45:32
7楼2016-05-20 11:24:57
8楼2016-05-20 11:25:47