QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 9408|回复: 7
打印 上一主题 下一主题

求助: fmincon的参数问题

[复制链接]
字体大小: 正常 放大

3

主题

3

听众

20

积分

升级  15.79%

该用户从未签到

自我介绍
对新事物永远热爱
跳转到指定楼层
1#
发表于 2010-5-20 10:10 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
t=2.3263;
2 V% I4 q4 O  q3 Q. k: ~& vk=0.5244;
* V7 [- N+ {+ Z8 Y3 m, g* BR=0.15;+ \2 n5 N3 u" E1 A
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
7 ]$ U3 l  n+ [# J# S* C4 oE=[1234,864,-617,-246,987;864,642,586,-198,-815;-617,586,1433,469,-963;-246,-198,469,346,25;987,-815,-963,25,1827];
3 q1 c; [7 z, z6 ?7 f. CE=E/10000;
: B- z- ^( }, P9 }  n' l6 T8 O" ]c=0.0075;$ P( G5 z) B3 w. F5 X' g9 a: E6 L
options = optimset('Algorithm','active-set'); % run active-set algorithm  \0 S8 g$ Y2 Z; y$ I( Y
aeq=[1+c,1+c,1+c,1+c,1+c];5 w. U6 _: u7 w; M1 M
beq=1;. [8 B$ ^0 l% J# H, k) N
lb=[0,0,0,0,0]';
* |7 r. ]+ u/ o5 `3 i: J" H: q/ xub=[1,1,1,1,1]';
( w7 I! v. Y* ux0=[0.02,0.03,0.05,0.11,0.12]';
% x' s  y2 i; V! w[x,fval,exitflag,output] = fmincon(@(x) myfun( x,t,E,u,c),x0,[],[],aeq,beq,lb,ub,@(x) mycon( x,t,k,E,u,c,R),options)# ^' C7 s8 O: f" E* p5 {' }
7 l" c. i1 B0 o7 x) {% z2 n; J
function f = myfun( x,t,E,u,c )    %目标函数2 V% R2 x& K% Y9 f, U
%UNTITLED3 Summary of this function goes here
; }' Z/ h$ ]* V. t, D1 m4 n% S( k%   Detailed explanation goes here3 W) \# w( \9 j4 M/ ]4 p
f=t*sqrt(x'*E*x)-(u'-c)*x;& n/ e+ I! j3 A- X1 N5 b
end4 N4 X- e& W- A" q' l1 o6 V
+ W2 k2 F3 a6 W5 C' E" c
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
) c, m! M3 Q5 q* I- J2 r  x%UNTITLED4 Summary of this function goes here9 n  D/ F" x2 f7 }1 }, ]
%   Detailed explanation goes here+ e# M. b$ \( K! n
C=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
& L4 D  }% b8 [8 ]Ceq=[];% a: Z7 G3 a( Z6 p
end5 [6 R6 \5 Z* d; B) {, p
- i; ]7 c8 Q# L0 }" P2 G! ?
% S) {# B  U: F: a
运行结果:
- c7 r& ]8 _& {" j2 [* JSolver stopped prematurely.
$ O2 \! M9 Q8 i- A% I" O, Jfmincon stopped because it exceeded the function evaluation limit,
9 w" e" a% [; K  X0 Hoptions.MaxFunEvals = 500 (the default value).
+ ~' e4 u: l7 ^4 ^% a) D
! ^8 `7 g2 k7 ~( }x =
. O0 j7 g6 q0 x# }& W  k5 U   -0.01754 {+ B% I) A" e& A
    0.4875' B& w, B  j. _
   -0.0175
! W# b9 }6 u+ y. v" V5 T4 o    0.3074
5 l# H% P8 u/ P' o  r/ m    0.2328& d" I) m# c; Z

* f& T6 g8 S3 }: Qfval =/ v! T6 z  I2 F4 x# T: G
   -0.03903 f- d  b; a. V2 U

" r* V: ^$ h+ gexitflag =& h; t3 E  w9 V1 c. t$ t+ g+ a
     0
* [. \2 r1 e4 k! p) Q
- e. u5 _' _$ W0 B! t. z  a+ Koutput = 7 D5 ]5 Q0 [" j
         iterations: 51
8 s& k; L# E" Q4 ]8 M: T- i          funcCount: 5058 w" |" b# z( F9 i! ?+ U$ S( L
       lssteplength: 0.1250
9 V2 f3 \, a! p; S           stepsize: 0.0812" M% }/ D( }% V) T
          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
! W: _3 c( J) X8 Z. v! U" ~      firstorderopt: 2.47322 j( w+ H: ?: y! A- y4 G. ^
    constrviolation: 0.0509
/ k( m: c9 p, |- r3 y6 H9 l, _- l            message: [1x145 char]
- @2 ~8 J, g* ~* x! I2 I9 p各位大虾,程序如上所述,为什么函数提前就结束了,是不是要改默认参数设置,小弟也试过,但是捣鼓不出来,希望各位高人不吝赐教,小弟先拜谢了  Orz
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
zfq12 实名认证       

2

主题

6

听众

981

积分

升级  95.25%

该用户从未签到

自我介绍
我很喜欢数学,希望在数学中国能结识更多的朋友。。。。。

新人进步奖

群组西南大学建模组

群组数学建摸协会

群组Matlab讨论组

回复

使用道具 举报

zfq12 实名认证       

2

主题

6

听众

981

积分

升级  95.25%

该用户从未签到

自我介绍
我很喜欢数学,希望在数学中国能结识更多的朋友。。。。。

新人进步奖

群组西南大学建模组

群组数学建摸协会

群组Matlab讨论组

回复

使用道具 举报

zhangal 实名认证       

5

主题

3

听众

45

积分

升级  42.11%

该用户从未签到

自我介绍
冰雪聪明 热情大方 诚心诚意
不会啊,帮不上》》》》》》》》》》》》》》》》》》》》》》》》》》》》
回复

使用道具 举报

3

主题

3

听众

20

积分

升级  15.79%

该用户从未签到

自我介绍
对新事物永远热爱
回复 3# zfq12
  l4 _* p8 ]% V# z# U
. @: X: ~1 l5 W# N0 k" F
0 V* i, o3 i1 L) c" \9 B$ G    我也试过把lb,ub上下界去掉,结果程序跑了一个上午都还在busy。。。不知何故啊
回复

使用道具 举报

0

主题

0

听众

2

积分

升级  40%

该用户从未签到

回复

使用道具 举报

alair002        
头像被屏蔽

1

主题

4

听众

328

积分

升级  9.33%

  • TA的每日心情
    擦汗
    2012-2-6 07:40
  • 签到天数: 6 天

    [LV.2]偶尔看看I

    提示: 作者被禁止或删除 内容自动屏蔽
    回复

    使用道具 举报

    0

    主题

    3

    听众

    7

    积分

    升级  2.11%

    该用户从未签到

    自我介绍
    888888
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-7-22 18:59 , Processed in 0.475633 second(s), 96 queries .

    回顶部