| 查看: 2719 | 回复: 16 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
yijianmeixyt铁杆木虫 (著名写手)
一剪梅
|
[求助]
Matlab中的fmincon功能求助 已有2人参与
|
|
|
本人在非线性约束优化问题中遇到难题，请高人帮忙，谢谢！ 问题描述： 求解函数 f=w1+w2+40*cos(a1)*cos(b1)+30* sin(a1)cos(b1)-40 在非线性约束条件 0.05≤w1+20*sin(b1)+20* sin(a1)*cos(b1)≤0.05 0.05≤w1-20*sin(b1)+20* sin(a1)*cos(b1)≤0.05 0.05≤w1-20*sin(b1)-20* sin(a1)*cos(b1)≤0.05 0.05≤w1+20*sin(b1)-20* sin(a1)*cos(b1)≤0.05 0.05≤40* sin(a1)*cos(b1)≤0.05 0.05≤w2-20*sin(b2)+50* sin(a2)*cos(b2)≤0.05 0.05≤w2-20*sin(b2)-50* sin(a2)*cos(b2)≤0.05 0.05≤w2+20*sin(b2)+50* sin(a2)*cos(b2)≤0.05 0.05≤w2+20*sin(b2)-50* sin(a2)*cos(b2)≤0.05 下的优化解。约束条件中的角度a和b都非常小，一般只有0.006*180/π左右，而w的范围在±0.05。 我使用fmincon()命令来计算，程序代码如下： %/////主程序 clc; clear; % w1=x(1); % w2=x(2); % cos(a1)=x(3) % cos(b1)=x(4) % sin(b2)=x(5) % sin(a2)=x(6) x0=[0;0;0;0;0;0]; ub=[0.05;0.05;1;1;0.007;0.007]; lb=-ub; options = optimset('display','iter','Algorithm','active-set'); [x,fval,exitflag]=fmincon(@myfun,x0,[],[],[],[],lb,ub,@nonlcon1,options); disp(fval); %/////////目标函数 function f= myfun(x) f=x(1)+40*x(3)*x(4)+30*x(4)*sqrt(1-x(3)^2)-40+x(2)+0*x(5)+0*x(6); end %////////约束 function [c ceq]=nonlcon1(x) c=[x(1)-20*sqrt(1-x(4)^2)-20*x(4)*sqrt(1-x(3)^2)-0.05;-(x(1)-20*sqrt(1-x(4)^2)-20*x(4)*sqrt(1-x(3)^2)-0.05)-0.1;... x(1)+20*sqrt(1-x(4)^2)-20*x(4)*sqrt(1-x(3)^2)-0.05;-(x(1)+20*sqrt(1-x(4)^2)-20*x(4)*sqrt(1-x(3)^2)-0.05)-0.1;... x(1)+20*sqrt(1-x(4)^2)+20*x(4)*sqrt(1-x(3)^2)-0.05;-(x(1)+20*sqrt(1-x(4)^2)+20*x(4)*sqrt(1-x(3)^2)-0.05)-0.1;... x(1)-20*sqrt(1-x(4)^2)+20*x(4)*sqrt(1-x(3)^2)-0.05;-(x(1)-20*sqrt(1-x(4)^2)+20*x(4)*sqrt(1-x(3)^2)-0.05)-0.1;... x(2)-20*x(5)-50*sqrt(1-x(5)^2)*x(6)-0.05;-(x(2)-20*x(5)-50*sqrt(1-x(5)^2)*x(6)-0.05)-0.1;... x(2)-20*x(5)+50*sqrt(1-x(5)^2)*x(6)-0.05;-(x(2)-20*x(5)+50*sqrt(1-x(5)^2)*x(6)-0.05)-0.1;... x(2)+20*x(5)+50*sqrt(1-x(5)^2)*x(6)-0.05;-(x(2)+20*x(5)+50*sqrt(1-x(5)^2)*x(6)-0.05)-0.1;... x(2)+20*x(5)-50*sqrt(1-x(5)^2)*x(6)-0.05;-(x(2)+20*x(5)-50*sqrt(1-x(5)^2)*x(6)-0.05)-0.1]; ceq=[]; end 结果明显不对，请求高人指点，谢谢！ |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab中fmincon的fun是矩阵运算怎么办？
已经有5人回复
- matlab线性规划
已经有5人回复
- 求教matlab线性优化求最小值的问题
已经有7人回复
- matlab 绘图功能
已经有7人回复
- 求助matlab解非线性规划。求代码。
已经有13人回复
- matlab里有没有类似于origin合并图形的功能
已经有8人回复
- matlab fmincon优化函数 怎样知道每次得到的迭代点
已经有3人回复
- 小弟新学matlab，想借用 fmincon求个函数的最小值，求助错误该法，感谢！
已经有4人回复
- fmincon优化问题
已经有5人回复
- 【资源】《matlab辅助激光光学分析与应用》 pdf 里面含有很多matlab程序
已经有81人回复
yijianmeixyt
铁杆木虫 (著名写手)
一剪梅
- 应助: 158 (高中生)
- 金币: 16506.5
- 散金: 59
- 红花: 5
- 帖子: 1880
- 在线: 304.6小时
- 虫号: 1161506
- 注册: 2010-12-03
- 性别: GG
- 专业: 机构学与机器人
3楼2014-03-26 09:12:37
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2014-03-26 08:58:09
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
yijianmeixyt: 金币+5, ★有帮助, 提供了一个解决思路 2014-03-26 09:32:30
感谢参与，应助指数 +1
yijianmeixyt: 金币+5, ★有帮助, 提供了一个解决思路 2014-03-26 09:32:30
|
fmincon用的方法本质上是非线性目标函数的迭代方法优化; 大部分情况下得不到好的结果. 尤其是初值不理想的情况下; 所以,最好用专门的全局优化算法. 1.模拟退火, 2. 遗传算法; 3.粒子群; 4. 差分演化; 5.其它蚁群蜂群鱼群等算法... 得到的结果再作为fmincon的初值. 此外,含有正余弦函数的优化问题,最好变成不含正余弦的两个变量形式 然后增加一个两变量平方和为1的非线性约束. |
4楼2014-03-26 09:24:51
yijianmeixyt
铁杆木虫 (著名写手)
一剪梅
- 应助: 158 (高中生)
- 金币: 16506.5
- 散金: 59
- 红花: 5
- 帖子: 1880
- 在线: 304.6小时
- 虫号: 1161506
- 注册: 2010-12-03
- 性别: GG
- 专业: 机构学与机器人
5楼2014-03-26 09:32:04
40