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

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

当前主题已经存档。

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

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

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

回复此楼

» 猜你喜欢

321求调剂已经有5人回复
324求调剂已经有6人回复
359求调剂已经有7人回复
一志愿华北电力大学（北京），材料科学与工程学硕265，求调剂已经有11人回复
368求调剂已经有7人回复
288求调剂，一志愿华南理工大学071005 已经有3人回复
求生物学调剂已经有13人回复
284求调剂已经有3人回复
调剂已经有8人回复
一志愿211，化学学硕，310分，本科重点双非，求调剂已经有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
专业: 能源化工

谢谢各位提示！
已经搞定了：）

赞一下

回复此楼

沉默坚守

9楼2008-04-01 11:04:02

已阅回复此楼关注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 个回答

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

[考研] 292求调剂 +11	2022080213 2026-04-04	13/650	2026-04-04 18:38 by macy2011
[考研] 调剂 +3	是可乐不是可乐 2026-04-04	3/150	2026-04-04 17:21 by wxiongid
[考研] 085600，专业课化工原理，320分求调剂 +8	大馋小子 2026-03-29	8/400	2026-04-04 16:21 by dongzh2009
[考博] 申博 +4	IQwQl 2026-04-04	4/200	2026-04-04 13:57 by 工大学长
[考研] 材料383求调剂 +4	郭阳阳阳成 2026-04-04	4/200	2026-04-04 12:07 by 无际的草原
[基金申请] 同年申报，是不是就不会派给你函评本子了？ +4	redcom 2026-03-29	4/200	2026-04-04 10:08 by kudofaye
[考研] 求调剂 +4	xzghyuj 2026-04-04	4/200	2026-04-04 09:51 by 唐沐儿
[考研] 387求调剂 +4	爱吃片豆土 2026-04-03	5/250	2026-04-04 08:10 by 岸上的一条鱼
[考研] 281求调剂 +10	aaawhy 2026-04-03	10/500	2026-04-03 21:42 by lbsjt
[考研] 285求调剂 +5	AZMK 2026-04-03	8/400	2026-04-03 18:17 by AZMK
[考研] 303求调剂 +10	DLkz1314. 2026-03-30	10/500	2026-04-03 18:03 by Jimmyandyou
[考研] 315分 085602 求调剂 +15	26考研上岸版26 2026-04-02	15/750	2026-04-03 12:45 by xingguangj
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +5	蓝笺片 2026-04-02	6/300	2026-04-02 21:37 by dongzh2009
[考研] 材料求调剂 +10	呢呢妮妮 2026-04-01	13/650	2026-04-02 09:17 by olim
[考研] 11408 321分求调剂 +3	huchun12138 2026-03-30	4/200	2026-04-01 22:48 by guanxin1001
[考研] 265求调剂 +11	yelck 2026-04-01	12/600	2026-04-01 19:12 by 549790059
[考研] 311求调剂 +10	李芷新1 2026-03-31	10/500	2026-04-01 14:38 by chenqifeng666
[考研] 0855机械初试280求调剂 +3	kazenotori 2026-03-31	3/150	2026-04-01 10:08 by JourneyLucky
[考研] 材料科学与工程求调剂 +13	深V宿舍吧 2026-03-29	13/650	2026-03-31 19:50 by Dyhoer
[考研] 269求调剂 +4	我想读研11 2026-03-31	4/200	2026-03-31 10:04 by cal0306