| 查看: 1837 | 回复: 12 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
多元函数求最值
|
||||
|
刚接触matlab，我想用matlab求下面函数的最大、最小值以及取最值的时候自变量取值，是一个有4个自变量的函数，自变量都有一定的变化范围。具体函数是： f=[d*（s+b+d）+0.858*r*r-0.858*（r+d/2）*（r+d/2）]/[s*b-0.858*r*r]；其中，2.79≤s≤2.83;7.97≤b≤8.07;0.13≤d≤0.18;0.6≤r≤1.有木有大神帮帮忙知道一下啊  |
回复此楼» 猜你喜欢
354求调剂
已经有9人回复
-
化学工程321分求调剂
已经有22人回复
-
考研调剂
已经有3人回复
-
326求调剂
已经有8人回复
-
333求调剂
已经有7人回复
-
一志愿东华大学控制学硕320求调剂
已经有3人回复
-
考研化学学硕调剂，一志愿985
已经有7人回复
-
0703化学调剂
已经有15人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有14人回复
-
326求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 连续随机过程的均值函数是连续函数吗？哪里有证明？
已经有10人回复
- 小弟新学matlab，想借用 fmincon求个函数的最小值，求助错误该法，感谢！
已经有4人回复
- 怎么得到旋轨耦合中不同J值的本征值和波函数？
已经有6人回复
- 为什么画出的对关联函数第三峰刚开始,之后的函数值就全变为零了,怎么解决这个问题
已经有5人回复
- 多元函数如何积分啊？
已经有3人回复
- 求解一个函数的最大值,谢谢!!
已经有5人回复
- 【求助】关于聚酯多元醇的羟值问题
已经有5人回复
- 【求助】相关函数归一化后，输出值怎么变成NaN了？？
已经有4人回复
- 【求助】关于fminsearch函数初值选择
已经有3人回复
- 【求助】函数最小值问题？
已经有9人回复
- 【求助】酸酐与多元醇反应产物的羟值测定
已经有4人回复
- 【求助】控制律设计Lyapunov函数的参数取值范围
已经有8人回复
- 【求助】只知道有限个函数值求积分
已经有10人回复
- 【求助】求多元函数的极值
已经有11人回复
- 【求助】多元醇羟值测定
已经有11人回复
yqx1985
木虫 (著名写手)
云中仙
- 应助: 34 (小学生)
- 贵宾: 0.032
- 金币: 1608.8
- 散金: 1182
- 红花: 33
- 帖子: 1544
- 在线: 620.3小时
- 虫号: 550008
- 注册: 2008-04-22
- 性别: GG
- 专业: 药物设计与药物信息
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
iwill1k: 金币+3, ★★★★★最佳答案 2012-04-14 11:14:40
感谢参与，应助指数 +1
iwill1k: 金币+3, ★★★★★最佳答案 2012-04-14 11:14:40
|
的却如此呢，1stOpt给的这个结果的确奇怪，我又试了次 ====== 结果 ====== 迭代数: 17 计算用时(时:分:秒:毫秒): 00:00:00:125 计算中止原因: 达到收敛判定标准 优化算法: 共扼梯度法 + 通用全局优化法 函数表达式: (d*(s+b+d)+0.858*r*r-0.858*(r+d/2)*(r+d/2))/(s*b-0.858*r*r) 目标函数值(最小): 0.0599967759074685 s: 2.82999099858591 b: 8.06997582502897 d: 0.130000407383572 r: 0.999995774967586 ====== 计算结束 ====== |
8楼2012-04-10 11:14:50
【答案】应助回帖
感谢参与，应助指数 +1
| 百度你会不？http://www.google.com.hk/#hl=zh- ... rgbld.&bav=on.2,or.r_gc.r_pw.,cf.osb&fp=b7e86341b10bb823&biw=1366&bih=667 |
2楼2012-04-09 11:43:42
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 程序强帖: 2
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
>> f=@(x)(x(3).*(x(1)+x(2)+x(3))+0.858*x(4).*x(4)-0.858*(x(4)+x(3)/2).*(x(4)+x(3)/2))./(x(1).*x(2)-0.858*x(4).*x(4)); opt=optimset('Algorithm','interior-point'); [X,FVAL,EXITFLAG]= fmincon(f,[2.8 8 0.15 0.8],[],[],[],[],[2.79 7.97 0.13 0.6],[2.83 8.07 0.18 1],[],opt) Local minimum found that satisfies the constraints. Optimization completed because the objective function is non-decreasing in feasible directions, to within the default value of the function tolerance, and constraints were satisfied to within the default value of the constraint tolerance. X = 2.8300 8.0700 0.1300 0.9998 FVAL = 0.0600 EXITFLAG = 1 >> f1=@(x)-(x(3).*(x(1)+x(2)+x(3))+0.858*x(4).*x(4)-0.858*(x(4)+x(3)/2).*(x(4)+x(3)/2))./(x(1).*x(2)-0.858*x(4).*x(4)); >> [X,FVAL,EXITFLAG]= fmincon(f1,[2.8 8 0.15 0.8],[],[],[],[],[2.79 7.97 0.13 0.6],[2.83 8.07 0.18 1],[],opt) Local minimum found that satisfies the constraints. Optimization completed because the objective function is non-decreasing in feasible directions, to within the default value of the function tolerance, and constraints were satisfied to within the default value of the constraint tolerance. X = 2.7900 7.9700 0.1800 0.6000 FVAL = -0.0853 EXITFLAG = 1 |
3楼2012-04-09 12:35:37
yqx1985
木虫 (著名写手)
云中仙
- 应助: 34 (小学生)
- 贵宾: 0.032
- 金币: 1608.8
- 散金: 1182
- 红花: 33
- 帖子: 1544
- 在线: 620.3小时
- 虫号: 550008
- 注册: 2008-04-22
- 性别: GG
- 专业: 药物设计与药物信息
楼上正解 给你贴个1stOpt的 Parameters s[2.79,2.83], b[7.97,8.07],d[0.13,0.18],r[0.6,1]; Minimum =f; Function f=(d*(s+b+d)+0.858*r*r-0.858*(r+d/2)*(r+d/2))/(s*b-0.858*r*r); ====== 结果 ====== 迭代数: 18 计算用时(时:分:秒:毫秒): 00:00:00:94 计算中止原因: 达到收敛判定标准 优化算法: 共扼梯度法 + 通用全局优化法 函数表达式: (d*(s+b+d)+0.858*r*r-0.858*(r+d/2)*(r+d/2))/(s*b-0.858*r*r) 目标函数值(最小): 0 s: 2.82999999999992 b: 8.06833196340274 d: 0 r: 0.830131999171392 ====== 计算结束 ====== |
4楼2012-04-09 13:43:37
40