目标函数最优解

zzz.dan

使用fmincon求目标函数的最优解时，其中x【4】的值为12 14 16 18 20 22 25 28 32中的一个，因此我在约束函数nonlcon里设置了一行逻辑或的等式约束条件  但是最后求解出来的时候  x【4】是带小数的12.XXX  并没有满足这个约束 这时候什么情况 我该怎么解决  求大佬指导!谢谢
; {1 J2 M3 P& A; V2 t7 H. D, @3 H
* }. a) s( |; L' {- \- R4 Q程序如下：
clc;
clear all;
1 n) |. Z% o+ @. y3 E/ Wclose all;

+ L8 W& b. I, a! h+ mx0 = [100,100,100,100];
/ k% l9 |% T' b! r2 r5 f! O
A=[-1,0,0,0;0,-1,0,0;1,-1,0,0;]
- H* |) Q  [" WB = [-200;0;0]
Aeq = []
$ ^( M/ k& b+ i: o) T- ^Beq = []
  H! ?! ^0 `1 v5 ]LB =[200,0,4,12];
$ j$ [% e5 k3 @& R2 f- R! wUB =[1000,1000,20,32]  
! {8 T. x$ n* E" v' Z[x, fval] = fmincon('fun', x0, A, B, Aeq, Beq, LB, UB,'nonlcon')
$ T8 R0 V% s2 w, L/ Y5 s4 y4 e5 z
2 v/ w7 L) @) ~& Z: _9 ]7 Gfprintf('the minimum value in x:\n');
fprintf('%f\n%f\n',x(1),x(2),x(3),x(4));
1 ]& Q1 Q8 E/ a8 Y7 q8 u0 U' O& Cfprintf('the minimum value:\n');
/ Z( D5 ~6 B  ?5 d, w7 J  N; L" Cfprintf('%f\n',fval);
6 P, o" g9 @3 n9 U


function [c,ceq] = nonlcon(x);
fai=1.2687-0.0291*(6500/x(1))+0.0001*(6500/x(1))^2;
3 F' h' C+ E& T/ f8 I9 @c=[651000-(1/1.2)*(10*x(1)*x(2)+310*x(3)*(pi*x(4)^2)/4)*fai ;
, G* ?# e' b* @$ |* H/ p9 J: C-x(3)*(pi*x(4)^2)/(4*x(1)*x(2))+0.008;
4.52-x(3)*(pi*(10*x(4))^2)/4];
ceq=[mod(x(3),1);
4 Q3 `. H/ `6 j0 P% T

end
, K- x, O  f; W5 K5 x+ O- U

nkkopd

帮你顶一下

大小的小

2 I( \! I3 l5 z/ k! p0 l* O9 I你等式约束的方式写成(x-12)*(x-14)*(x-16)*(x-18)*(x-20)*(x-22)*(x-25)*(x-28)=0，设置下精度试试

llbnmo

二楼正解。。

shuddkk