|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
8 w8 n3 [4 d2 x( s6 v
7 y$ m+ x3 N/ F/ `2 o& kfunction f=fun1(x)9 w! ^1 }6 W, n7 N! d5 Q) u3 D
f=-5*sin(x(1))*sin(x(2))*sin(x(3))*sin(x(4))*sin(x(5))-sin(5*x(1)*5*x(2)*5*x(3)*5*x(4)*5*x(5))+8;2 `5 U, j4 O" ~. y' w
上述代码已经保存为fun1.m文件。( ]. x: ^3 r4 d' k6 [* P4 ~
* H8 C# y, d0 ^+ `7 B' V0 E, I
再执行下述代码:
, D X( S+ k/ F+ o' l# GA=[];b=[];Aeq=[];beq=[];lb=[0,0,0,0,0];ub=[0.9*pi,0.9*pi,0.9*pi,0.9*pi,0.9*pi];: d/ e5 ~+ s$ p4 C5 f
x0=[0,0,0,0,0];
' X4 t) O% Y8 }[x,fval]=fmincon('fun1',x0,A,b,Aeq,beq,lb,ub)
- q6 e" X/ [6 m1 B- i/ h, d- Z, }
2 Q0 r1 \ K A. U; ?4 @. ~; z% N9 ~ R6 G! W3 j0 Z2 o( c0 v
运行结果为:
: q3 u) ^& ?6 ux =
; g% K4 d" K% P. K7 \0 f9 M/ h2 N/ q& d
1.4298 1.4298 1.4298 1.4298 1.4298
5 }9 L5 G7 o$ L; c' z/ k2 {/ u7 w7 _9 N$ { F
/ V* a2 h7 }1 p/ j" u0 g) b8 O; g
fval = [* K% ^. _2 X ]: |
4 ~! f7 d7 r0 l) e3 l' w7 g" h
2.24313 B# E* A3 B( V4 A+ n- l4 J! {
显然结果精度不够。
1 o ~+ h) G; L# D2 B中间出现了) M7 A. K0 e3 W
[backcolor=rgba(255, 255, 255, 0.8)]fmincon停止是因为当前步骤的大小小于. ^) o2 @7 B1 r4 y! r5 L, @; J
[backcolor=rgba(255, 255, 255, 0.8)]默认值的步骤尺寸公差和约束
/ n' R6 Z3 c7 n[backcolor=rgba(255, 255, 255, 0.8)]满意的默认值约束宽容。
: C& _6 F7 Z o0 H/ E' L
% r) t3 K3 X% B' W7 o9 {请教大神,这个错误需要如何解决?非常感谢! |
|