24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 2037  |  回复: 8

老花

木虫 (著名写手)

老花

[求助] 帮忙给看看这个matlab优化函数 问题

主:
A=[];b=[];Aeq=[];Beq=[];Lb=[];Ub=[];
[x,f]=fmincon('myobj',[0;0],[],[],[],[],[],[],'mycon');

目标函数:function y=myobj(x)
y=(-1)*8*x(2)*(1-x(1))*x(3);

非线性约束function[c,ceq]=mycon(x)
cep=[x1*(1-x(1))*x(3)-x(2)*(1+x(2));2*acos(exp(-x(4)))-x(3);1.5*2/sin(x(5))-x(4);atan((1-x(1))/(1+x(2)))];
请帮忙看看程序的错误,谢谢!
回复此楼
平凡中不平凡。我们都希望长寿,但是只有有意义的人生和生命中那些闪光的瞬间才是最重要的。用我们的精神而不是钟表来度量时间的流逝吧。....
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

starlee

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
老花: 金币+100, ★★★★★最佳答案, 结果有点偏差,我再研究一下,谢谢啊! 2012-06-28 23:12:22
老花: 回帖置顶 2012-06-28 23:24:42
A=[];B=[];Aeq=[];Beq=[];Lb=[];Ub=[];
[x,fval]=fmincon('myobj',rand(1,3),A,B,Aeq,Beq,Lb,Ub,'mycon')

function f=myobj(x)
f=-1*8*x(2)*(1-x(1))*2/pi*acos(exp(-1.5*2/sin(x(3))));

function [c,ceq]=mycon(x)
c=0;
ceq(1)=x(1)*(1-x(1))*2/pi*acos(exp(-1.5*2/sin(x(3))))-x(2)*(1+x(2));
ceq(2)=atan((1-x(1))/(1+x(2)))-x(3);

结果:
x =

    0.3187    0.1832    0.5224


fval =

   -0.9971
结果合不合理?
7楼2012-06-28 22:59:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

starlee

金虫 (正式写手)

目标函数只到x(3)
约束条件中却有x(4),这不合理啊,楼主
2楼2012-06-28 21:27:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老花

木虫 (著名写手)

老花

引用回帖:
2楼: Originally posted by starlee at 2012-06-28 21:27:48
目标函数只到x(3)
约束条件中却有x(4),这不合理啊,楼主

问题是里面没有啊! 我自己加上啊?
平凡中不平凡。我们都希望长寿,但是只有有意义的人生和生命中那些闪光的瞬间才是最重要的。用我们的精神而不是钟表来度量时间的流逝吧。....
3楼2012-06-28 21:46:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

starlee

金虫 (正式写手)

把x(4) x(5)有关的约束条件去掉,倒是能算个结果出来
4楼2012-06-28 22:09:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老花

木虫 (著名写手)

老花

function[c,ceq]=mycon(x)
cep=[x1*(1-x(1))*2/pi*acos(exp(-1.5*2/sin(x(3))))-x(2)*(1+x(2));atan((1-x(1))/(1+x(2)))-x(3)];

function y=myobj(x)
y=(-1)*8*x(2)*(1-x(1))*2/pi*acos(exp(-1.5*2/sin(x(3)))) ;

修改后还是这样! 请大家指点
Error in myobj (line 2)
y=(-1)*8*x(2)*(1-x(1))*2/pi*acos(exp(-1.5*2/sin(x(3)))) ;

Error in fmincon (line 601)
      initVals.f = feval(funfcn{3},X,varargin{:});

Caused by:
    Failure in initial user-supplied objective function evaluation.
    FMINCON cannot continue.
平凡中不平凡。我们都希望长寿,但是只有有意义的人生和生命中那些闪光的瞬间才是最重要的。用我们的精神而不是钟表来度量时间的流逝吧。....
5楼2012-06-28 22:33:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老花

木虫 (著名写手)

老花

引用回帖:
4楼: Originally posted by starlee at 2012-06-28 22:09:02
把x(4) x(5)有关的约束条件去掉,倒是能算个结果出来

你帮忙看看我下边改的啊这个!
平凡中不平凡。我们都希望长寿,但是只有有意义的人生和生命中那些闪光的瞬间才是最重要的。用我们的精神而不是钟表来度量时间的流逝吧。....
6楼2012-06-28 22:46:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老花

木虫 (著名写手)

老花

引用回帖:
7楼: Originally posted by starlee at 2012-06-28 22:59:10
A=[];B=[];Aeq=[];Beq=[];Lb=[];Ub=[];
=fmincon('myobj',rand(1,3),A,B,Aeq,Beq,Lb,Ub,'mycon')

function f=myobj(x)
f=-1*8*x(2)*(1-x(1))*2/pi*acos(exp(-1.5*2/sin(x(3))));

function =mycon(x)
c=0;
...

rand(1,3),A,B,Aeq
啥意思啊? 解释一下呗!
平凡中不平凡。我们都希望长寿,但是只有有意义的人生和生命中那些闪光的瞬间才是最重要的。用我们的精神而不是钟表来度量时间的流逝吧。....
8楼2012-06-28 23:24:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

starlee

金虫 (正式写手)

引用回帖:
8楼: Originally posted by 老花 at 2012-06-28 23:24:31
rand(1,3),A,B,Aeq
啥意思啊? 解释一下呗!...

A,B是空矩阵啊,第一行写明了
rand(1,3)指 初值为1行3列的0~1随机数
结果是否准确与初值的关系比较大的
9楼2012-06-29 12:37:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 老花 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 324求调剂 +15 想上学求调 2026-04-03 15/750 2026-04-08 14:51 by zzucheup
[考研] 求调剂 +9 月@163.com 2026-04-07 11/550 2026-04-08 14:48 by qlm5820
[考研] 生物学328分求调剂 +7 闪电kkl 2026-04-08 7/350 2026-04-08 13:54 by 蔡苏阳
[考研] 一志愿深大085601材料工程专业(专硕)300分可以调剂去哪 +12 10160315 2026-04-02 12/600 2026-04-08 11:49 by 猪会飞
[考研] 机械专硕273请求调剂 +5 庚申壬申 2026-04-07 5/250 2026-04-08 11:46 by zhyzzh
[考博] 申博 +8 IQwQl 2026-04-04 8/400 2026-04-08 09:43 by 0608104024
[考研] 266调剂 +6 daya sun 2026-04-07 7/350 2026-04-08 08:59 by 蓝云思雨
[考研] 计算机408|在校多次国家级竞赛获奖|申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 08600生物与医药-327 +9 18755400796 2026-04-05 9/450 2026-04-06 22:35 by 52305043001
[考研] 262求调剂 +7 天下第一文 2026-04-04 8/400 2026-04-05 21:31 by 激流勇渡
[考研] 085600,320分求调剂 +7 大馋小子 2026-04-01 8/400 2026-04-05 21:19 by 学员8dgXkO
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11 adaie 2026-04-04 12/600 2026-04-05 19:04 by 蓝云思雨
[考研] 考研调剂生寻找导师 +3 顾瞻考研啊 2026-04-05 3/150 2026-04-05 18:18 by 啵啵啵0119
[考研] 271分求调剂学校 +12 zph158488! 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
[考研] 复试调剂 +6 范根培 2026-04-04 6/300 2026-04-04 14:27 by 土木硕士招生
[考研] 387求调剂 +4 爱吃片豆土 2026-04-03 5/250 2026-04-04 08:10 by 岸上的一条鱼
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 【求调剂】新能源材料本科,一志愿211,初试321 +6 求调剂学校, 2026-04-02 6/300 2026-04-02 09:41 by 晴空210210
信息提示
请填处理意见