版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jlzeng

木虫 (正式写手)

应助: 10 (幼儿园)
金币: 2413.7
散金: 5
红花: 4
帖子: 471
在线: 214小时
虫号: 427388
注册: 2007-07-29
性别: GG
专业: 能源化工

[交流] 请帮我看一段matlab程序

问题：
Min  10X2+10X3+40
Sub. to  1≤X1≤4
1≤X2≤4
1≤X3≤4
1≤X4≤4
-(6.9402131+0.252469*x(1)-0.276517*x(2)+0.207394*x(3)+0.314560*x(4)+0.039335*x(1)*x(1)+0.039335*x(2)*x(2)-0.003563*x(3)*x(3)-0.006647*x(4)*x(4)-0.006296*x(1)*x(2)-0.034093*x(1)*x(3)-0.059668*x(1)*x(4)+0.010462*x(2)*x(3)+0.005134*x(2)*x(4)-0.008807*x(3)*x(4))+8.62≤0
6.9402131+0.252469*x(1)-0.276517*x(2)+0.207394*x(3)+0.314560*x(4)+0.039335*x(1)*x(1)+0.039335*x(2)*x(2)-0.003563*x(3)*x(3)-0.006647*x(4)*x(4)-0.006296*x(1)*x(2)-0.034093*x(1)*x(3)-0.059668*x(1)*x(4)+0.010462*x(2)*x(3)+0.005134*x(2)*x(4)-0.008807*x(3)*x(4)-9.17≤0
0.6730551+0.301146*x(1)-0.078675*x(2)+0.011658*x(3)+0.178452*x(4)+0.009664*x(1)*x(1)+0.009664*x(2)*x(2)+0.001284*x(3)*x(3)-0.002713*x(4)*x(4)-0.012422*x(1)*x(2)-0.023157*x(1)*x(3)-0.036011*x(1)*x(4)+0.024744*x(2)*x(3)-0.020623*x(2)*x(4)-0.002149*x(3)*x(4)-1.5≤0
-(0.6730551+0.301146*x(1)-0.078675*x(2)+0.011658*x(3)+0.178452*x(4)+0.009664*x(1)*x(1)+0.009664*x(2)*x(2)+0.001284*x(3)*x(3)-0.002713*x(4)*x(4)-0.012422*x(1)*x(2)-0.023157*x(1)*x(3)-0.036011*x(1)*x(4)+0.024744*x(2)*x(3)-0.020623*x(2)*x(4)-0.002149*x(3)*x(4)) ≤0
-(51.8277931+6.142814*x(1)+1.235285*x(2)-1.799889*x(3)+4.118358*x(4)-0.434893*x(1)*x(1)-0.434893*x(2)*x(2)+0.033447*x(3)*x(3)+0.076468*x(4)*x(4)-0.201259*x(1)*x(2)-0.652022*x(1)*x(3)+0.820377*x(1)*x(4)+1.839250*x(2)*x(3)-1.451273*x(2)*x(4)-0.023900*x(3)*x(4))+95≤0
51.8277931+6.142814*x(1)+1.235285*x(2)-1.799889*x(3)+4.118358*x(4)-0.434893*x(1)*x(1)-0.434893*x(2)*x(2)+0.033447*x(3)*x(3)+0.076468*x(4)*x(4)-0.201259*x(1)*x(2)-0.652022*x(1)*x(3)+0.820377*x(1)*x(4)+1.839250*x(2)*x(3)-1.451273*x(2)*x(4)-0.023900*x(3)*x(4)-100≤0

我的程序如下：
m文件
function [c,ceq]=mycon(x)
c=[-(6.9402131+0.252469*x(1)-0.276517*x(2)+0.207394*x(3)+0.314560*x(4)+0.039335*x(1)*x(1)+0.039335*x(2)*x(2)-0.003563*x(3)*x(3)-0.006647*x(4)*x(4)-0.006296*x(1)*x(2)-0.034093*x(1)*x(3)-0.059668*x(1)*x(4)+0.010462*x(2)*x(3)+0.005134*x(2)*x(4)-0.008807*x(3)*x(4))+8.62;
6.9402131+0.252469*x(1)-0.276517*x(2)+0.207394*x(3)+0.314560*x(4)+0.039335*x(1)*x(1)+0.039335*x(2)*x(2)-0.003563*x(3)*x(3)-0.006647*x(4)*x(4)-0.006296*x(1)*x(2)-0.034093*x(1)*x(3)-0.059668*x(1)*x(4)+0.010462*x(2)*x(3)+0.005134*x(2)*x(4)-0.008807*x(3)*x(4)-9.17;
0.6730551+0.301146*x(1)-0.078675*x(2)+0.011658*x(3)+0.178452*x(4)+0.009664*x(1)*x(1)+0.009664*x(2)*x(2)+0.001284*x(3)*x(3)-0.002713*x(4)*x(4)-0.012422*x(1)*x(2)-0.023157*x(1)*x(3)-0.036011*x(1)*x(4)+0.024744*x(2)*x(3)-0.020623*x(2)*x(4)-0.002149*x(3)*x(4)-1.5;
-(0.6730551+0.301146*x(1)-0.078675*x(2)+0.011658*x(3)+0.178452*x(4)+0.009664*x(1)*x(1)+0.009664*x(2)*x(2)+0.001284*x(3)*x(3)-0.002713*x(4)*x(4)-0.012422*x(1)*x(2)-0.023157*x(1)*x(3)-0.036011*x(1)*x(4)+0.024744*x(2)*x(3)-0.020623*x(2)*x(4)-0.002149*x(3)*x(4));
-(51.8277931+6.142814*x(1)+1.235285*x(2)-1.799889*x(3)+4.118358*x(4)-0.434893*x(1)*x(1)-0.434893*x(2)*x(2)+0.033447*x(3)*x(3)+0.076468*x(4)*x(4)-0.201259*x(1)*x(2)-0.652022*x(1)*x(3)+0.820377*x(1)*x(4)+1.839250*x(2)*x(3)-1.451273*x(2)*x(4)-0.023900*x(3)*x(4))+95;
51.8277931+6.142814*x(1)+1.235285*x(2)-1.799889*x(3)+4.118358*x(4)-0.434893*x(1)*x(1)-0.434893*x(2)*x(2)+0.033447*x(3)*x(3)+0.076468*x(4)*x(4)-0.201259*x(1)*x(2)-0.652022*x(1)*x(3)+0.820377*x(1)*x(4)+1.839250*x(2)*x(3)-1.451273*x(2)*x(4)-0.023900*x(3)*x(4)-100];
ceq=[];

主程序：
>> fun='x(1)*0+x(2)*10+10+x(3)*10+30+x(4)*0';
>> x0=[1 1 1 1];
>> A=[];
>> b=[];
>> Aeq=[];
>> beq=[];
>> lb=[1 1 1 1];
>> ub=[4 4 4 4];
>> [x,fval,exitflag,output,lambda,grad,hessian]=fmincon(fun,x0,A,b,Aeq,beq,lb,ub,@mycon)

运行结果：
Warning: Large-scale (trust region) method does not currently solve this type of problem,
switching to medium-scale (line search).
> In fmincon at 260
Maximum number of function evaluations exceeded;
increase OPTIONS.MaxFunEvals.
x =

4.2343 3.1404 4.3459 4.3459

fval =

  114.8638

exitflag =

   0

output =

   iterations: 66
      funcCount: 401
      stepsize: 1
      algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
firstorderopt: 10.0000
   cgiterations: []
      message: [1x79 char]

lambda =

      lower: [4x1 double]
      upper: [4x1 double]
      eqlin: [0x1 double]
   eqnonlin: [0x1 double]
   ineqlin: [0x1 double]
ineqnonlin: [6x1 double]

grad =

      0
10.0000
10.0000
      0

hessian =

0.5552 0.1233 0.0784 -0.4841
0.1233 0.0670 0.0061 -0.1023
0.0784 0.0061 0.0198 -0.0586
-0.4841 -0.1023 -0.0586 0.4462

请问是什么原因出错了？谢谢

回复此楼

» 猜你喜欢

面上意见出来了已经有20人回复
能否申诉？已经有8人回复
梦想已经有5人回复
中青基了要发朋友圈吗？已经有7人回复
有多少人是今天查系统知道结果的？已经有19人回复
29号明天会评吗已经有4人回复
国社科又开始会评了，不知道这次命运如何已经有10人回复
基金不中，共勉已经有12人回复
为什么到现在没收到通知？已经有5人回复
有没有仍没收到信息的已经有8人回复

沉默坚守

1楼 2008-03-31 10:01:58

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jlzeng

木虫 (正式写手)

应助: 10 (幼儿园)
金币: 2413.7
散金: 5
红花: 4
帖子: 471
在线: 214小时
虫号: 427388
注册: 2007-07-29
性别: GG
专业: 能源化工

csfn(金币+0,VIP+0):可以到此帖交流：http://emuch.net/bbs/viewthread.php?tid=759829&fpage=1

高手能否将我的程序进行适当修改？
我没学过matlab，很多东西都不明白：）

赞一下(10人)

回复此楼

沉默坚守

7楼2008-03-31 11:59:08

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 9 个回答

cuihang

铁杆木虫 (正式写手)

VASPer

应助: 0 (幼儿园)
金币: 7507.1
散金: 17
红花: 2
帖子: 697
在线: 78.5小时
虫号: 447551
注册: 2007-11-01
性别: GG
专业: 凝聚态物性I:结构、力学和

★ ★
spur(金币+2,VIP+0):感谢参与！

提示说你选用的算法规模太大不适合这个问题，建议换成较小的线性搜索
只是个警告，其实不必理会它

[ Last edited by cuihang on 2008-3-31 at 10:05 ]

赞一下(10人)

回复此楼

聪明人知道自己是傻子，傻子以为自己是聪明人

2楼2008-03-31 10:04:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jlzeng

木虫 (正式写手)

应助: 10 (幼儿园)
金币: 2413.7
散金: 5
红花: 4
帖子: 471
在线: 214小时
虫号: 427388
注册: 2007-07-29
性别: GG
专业: 能源化工

那是不是要换其他函数啊？
用哪个比较好呢？

赞一下

回复此楼

沉默坚守

3楼2008-03-31 10:06:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)

应助: 61 (初中生)
贵宾: 0.33
金币: 5959.3
散金: 1
红花: 21
帖子: 1884
在线: 491.1小时
虫号: 142265
注册: 2005-12-21
专业: 凝聚态物性 II ：电子结构

★ ★ ★
spur(金币+3,VIP+0):感谢应助！欢迎常来！！

Maximum number of function evaluations exceeded;
increase OPTIONS.MaxFunEvals.
___________________________
看英文，加大OPTIONS.MaxFunEvals

Warning: Large-scale (trust region) method does not currently solve this type of problem,
switching to medium-scale (line search).
提示说你选用的算法规模太大不适合这个问题，已经自动换成较小的线性搜索

exitflag =

0
这个是不成功的意思

赞一下(10人)

回复此楼

4楼2008-03-31 11:04:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 9 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 面上意见出来了 +11	黄鸟于飞Chao 2026-08-29	20/1000	2026-08-31 17:36 by dxcharlary
[基金申请] 能否申诉？ +7	echo8914667 2026-08-30	8/400	2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗？ +7	349506619 2026-08-28	7/350	2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 投票: 有多少人是今天查系统知道结果的？ +17	爱看书的可乐 2026-08-26	19/950	2026-08-31 11:51 by xiangy672
[基金申请] 29号明天会评吗 +4	笨笨唐 2026-08-28	4/200	2026-08-31 09:30 by huixian257
[基金申请] 国社科又开始会评了，不知道这次命运如何 +6	雨打竹帘 2026-08-30	10/500	2026-08-31 09:29 by huixian257
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +6	gdd2018 2026-08-28	11/550	2026-08-30 08:58 by 超级无敌华子
[基金申请] 怎么查啊 +6	huang1991js 2026-08-26	6/300	2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5	yuleib84 2026-08-26	6/300	2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5	ssyjh 2026-08-27	5/250	2026-08-27 20:50 by gdfollow
[基金申请] 基金未中，这种答复是模板吗？ +5	zhaosm1982 2026-08-27	6/300	2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5	叶九微 2026-08-26	5/250	2026-08-27 10:35 by l_zh2008
[基金申请] 我不理解！ +15	Edward_pc 2026-08-26	23/1150	2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25	木水思豆 2026-08-25	28/1400	2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4	bjdxyxy 2026-08-26	4/200	2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7	一怀馨秋 2026-08-26	7/350	2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18	Ldrop2023 2026-08-26	18/900	2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4	yanglien 2026-08-26	5/250	2026-08-26 11:10 by wenfengw83
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6	ziyangfang 2026-08-25	9/450	2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9	yuleib84 2026-08-25	10/500	2026-08-25 17:14 by alexon