QQ登录

只需要一步,快速开始

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

求助: fmincon的参数问题

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

3

主题

3

听众

20

积分

升级  15.79%

该用户从未签到

自我介绍
对新事物永远热爱
跳转到指定楼层
1#
发表于 2010-5-20 10:10 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
t=2.3263;+ o0 T& B# ^+ e% b
k=0.5244;
( {& P4 e$ u' UR=0.15;* |& f5 X: M% e  ?5 O& \0 |7 ^/ L
u=[0.1221,0.1786,0.1452,0.1054,0.1132]';
' s& V7 J% R' i" n% {: pE=[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];
, \+ C3 [: R' k. RE=E/10000;
$ x. n( L) ], T: R9 ]c=0.0075;
. N4 p$ `. D' G% m% ~, Hoptions = optimset('Algorithm','active-set'); % run active-set algorithm! @8 _: H; L6 G
aeq=[1+c,1+c,1+c,1+c,1+c];
! e) m0 D8 E* l: l1 Kbeq=1;& h! b7 ^7 ?/ I( Q
lb=[0,0,0,0,0]';* Z( S6 {: d: [
ub=[1,1,1,1,1]';
+ R5 V9 F: Z" I' y4 ~x0=[0.02,0.03,0.05,0.11,0.12]';1 C; _8 I8 G' v' D
[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)8 N) h, b! W( U! d/ p. R

- b0 ~8 q: R( Bfunction f = myfun( x,t,E,u,c )    %目标函数
  E% ^* p% n7 H+ t6 ^2 H%UNTITLED3 Summary of this function goes here
1 w5 D8 d" }. k%   Detailed explanation goes here+ x7 F( {/ o' I4 Y: z7 x
f=t*sqrt(x'*E*x)-(u'-c)*x;
: G. w  _, A& ~! q# Wend
3 }4 t2 C) T0 Y! O$ G8 l, R) d- @( Z% z. [, U5 G% L2 d
function [C,Ceq] = mycon( x,t,k,E,u,c,R )  %非线性约束
; K2 _; s; Q/ L9 V%UNTITLED4 Summary of this function goes here9 J3 q( Q, j. `% z
%   Detailed explanation goes here
' g6 t# p; H  f2 CC=t*k*sqrt(x'*E*x)/(t-k)-t*(u'-c)*x/(t-k)+t*R/(t-k);
8 o. @) L/ ^) x& {Ceq=[];) a) f. e$ G: {0 M; H
end+ v/ a0 h: {. K" }! A. M4 s
5 U) T/ l& u* V6 d

$ J! y" l  H2 H- V  R  S  B运行结果:/ }# E& Y3 K4 T
Solver stopped prematurely.* |) K: l- h0 q* l( G0 I' y& |
fmincon stopped because it exceeded the function evaluation limit,5 S0 b* C# A9 l
options.MaxFunEvals = 500 (the default value).
( D. i! N# K5 l8 p& i+ ^* `5 Q3 @1 \' y
x =
0 n) y8 e* G& l1 {% R- D# T   -0.0175
& G- c/ a- ^3 H9 k5 ~    0.4875$ U% h1 D4 g( Q& o
   -0.01752 B8 s; e$ D# F. C9 a# o; ~
    0.30744 j8 d% R& I2 V) O' M
    0.2328
2 o* o7 q! b* d5 a+ j5 r0 W9 ~+ I9 u  c
fval =
7 s: S/ L: I+ i, ]! ^9 E- m   -0.0390
3 j) O& k4 P  v( h3 X$ S8 t
7 o2 s3 {% c7 v: sexitflag =
- ]( ^- X/ E1 x1 k& F3 q! w     0* r% A0 l. I7 b5 U0 \# s
* \  K4 k! o: I" I- x% k- y# v/ ~
output =
, O6 v! Q; [( i' K" m! Z4 Y         iterations: 51) V- b- H" q# f4 G4 Y1 G  _
          funcCount: 505
# [( y  R: G. F4 K8 |3 o6 z       lssteplength: 0.1250* \9 {* y' Q2 B: M
           stepsize: 0.0812
# }$ H. }) a2 F: o9 @9 \7 B2 [- Z          algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
! ]0 v1 F  }. E# V      firstorderopt: 2.47320 s5 ^9 U0 D: u9 w$ V; R2 P
    constrviolation: 0.0509
3 E% V: o0 L1 T            message: [1x145 char]
6 H9 G- j6 [; J6 k, F各位大虾,程序如上所述,为什么函数提前就结束了,是不是要改默认参数设置,小弟也试过,但是捣鼓不出来,希望各位高人不吝赐教,小弟先拜谢了  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 . F8 U2 e: Z- `7 f. A& g4 u) ~
' o# D4 `( W4 l

$ @; C* Y$ Q$ h& q* _5 a    我也试过把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-6-4 05:10 , Processed in 0.986897 second(s), 96 queries .

    回顶部