|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
使用fmincon求目标函数的最优解时，其中x【4】的值为12 14 16 18 20 22 25 28 32中的一个，因此我在约束函数nonlcon里设置了一行逻辑或的等式约束条件 但是最后求解出来的时候 x【4】是带小数的12.XXX 并没有满足这个约束 这时候什么情况 我该怎么解决 求大佬指导!谢谢
6 p# y1 @; v+ p& a9 h% j, R* Q9 J5 c/ H! h, C
程序如下：
6 N8 n+ C7 ]. ~clc;
3 v( ~* X( C: E1 G8 Iclear all;( G5 V. b4 ~5 t6 ^6 j- D. e0 \" g- ]2 b
close all;
# K" Q6 p/ I- c: Y; g% V7 C# ?6 u' P4 O, e/ A4 R5 Z
x0 = [100,100,100,100];
& {+ `9 f5 X: j5 W" |5 C4 @7 K. `* u
A=[-1,0,0,0;0,-1,0,0;1,-1,0,0;]
4 @' ?2 ~0 A; P+ L6 N1 bB = [-200;0;0]6 u! X: Z/ p- c7 [- ?
Aeq = []' p- g, H: {' z+ x+ Z) @! w
Beq = []) s+ t2 Z0 @2 D$ `! X, S; A' ?
LB =[200,0,4,12];
1 v: M/ E8 T0 M- S) M& ^UB =[1000,1000,20,32]
% b, e, l- k- s7 r' g1 C. s[x, fval] = fmincon('fun', x0, A, B, Aeq, Beq, LB, UB,'nonlcon')9 K1 L: W) c1 Y$ E) n
8 }7 F, f+ s/ Rfprintf('the minimum value in x:\n');
. ^) F0 W+ M% D* z7 j2 a cfprintf('%f\n%f\n',x(1),x(2),x(3),x(4));
7 _, z3 k& B! q) U1 ~fprintf('the minimum value:\n');+ X( g. n3 O, C L4 q) r8 _( A, }1 ]
fprintf('%f\n',fval);
, I3 |. k# p3 v+ f
+ Q2 A% K: m- W- S3 D
' p. \( C. t( ]: |7 G _% }. m$ d# J& V% ~
function [c,ceq] = nonlcon(x);$ t) t+ m6 h. f
fai=1.2687-0.0291*(6500/x(1))+0.0001*(6500/x(1))^2;
& E1 V( ]* s& N4 y2 Lc=[651000-(1/1.2)*(10*x(1)*x(2)+310*x(3)*(pi*x(4)^2)/4)*fai ;
; O; W6 \* ?' C* K+ b9 Q! v: Q-x(3)*(pi*x(4)^2)/(4*x(1)*x(2))+0.008;
( _" F `, q% L' f) N! E6 z, d4.52-x(3)*(pi*(10*x(4))^2)/4];
) q6 K5 v/ y8 d! c7 N; |, I* yceq=[mod(x(3),1);
5 {3 Y# v/ y. V T4 z
, R9 W: Y) v; @1 j! J
0 I( @. a1 q- C5 N+ [# Send2 j0 I( Y* G, Q% y+ n" q
0 @ L0 z$ F( j0 g" O
|
|
/1 
发表于 2020-7-17 10:52
收藏
淘帖
支持！
反对！