应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » Matlab中的fmincon功能求助
51/1返回列表
查看: 2345  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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

结果明显不对，请求高人指点，谢谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
Rejoicinginhope,patientintribulation.
1楼 2014-03-26 08:43:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yijianmeixyt

铁杆木虫 (著名写手)

一剪梅
引用回帖:
6楼: Originally posted by dingd at 2014-03-26 09:45:58
试下1stOpt，全局优化求解：

MinFunction 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)*c ...

不太看得懂，我将w1,w2,a1,b1,a2,b2的函数值作为初始值，带入fmincon函数中，迭代了4步，电脑就算不下去了，得到的优化值为0.1，我觉得离目标很近了，我估计最终的优化值是0.07左右，谢谢！
不知道我的做法对不对，如果不对，请你指出，另外，如果最终算出来且结果正确的话，我将兑现金币，谢谢！
赞一下 回复此楼
Rejoicinginhope,patientintribulation.
8楼2014-03-26 10:08:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
yijianmeixyt: 金币+5, ★有帮助, 有助于问题的解决！ 2014-03-26 09:13:02
约束写的有问题吧，比如“0.05≤40* sin(a1)*cos(b1)≤0.05”，上下界都是0.05？
赞一下 回复此楼
2楼2014-03-26 08:58:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yijianmeixyt

铁杆木虫 (著名写手)

一剪梅
引用回帖:
2楼: Originally posted by dingd at 2014-03-26 08:58:09
约束写的有问题吧，比如“0.05≤40* sin(a1)*cos(b1)≤0.05”，上下界都是0.05？

不好意思，左边是-0.05
赞一下 回复此楼
Rejoicinginhope,patientintribulation.
3楼2014-03-26 09:12:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
yijianmeixyt: 金币+5, ★有帮助, 提供了一个解决思路 2014-03-26 09:32:30
fmincon用的方法本质上是非线性目标函数的迭代方法优化;
大部分情况下得不到好的结果. 尤其是初值不理想的情况下;
所以,最好用专门的全局优化算法.
1.模拟退火,
2. 遗传算法;
3.粒子群;
4. 差分演化;
5.其它蚁群蜂群鱼群等算法...
得到的结果再作为fmincon的初值.

此外,含有正余弦函数的优化问题,最好变成不含正余弦的两个变量形式
然后增加一个两变量平方和为1的非线性约束.
赞一下(1人) 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
4楼2014-03-26 09:24:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思：在爱情中有过遗憾吗？ +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错，点赞！ +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题，请指教！ 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
关闭
请填处理意见
关闭