目标函数最优解

zzz.dan

使用fmincon求目标函数的最优解时，其中x【4】的值为12 14 16 18 20 22 25 28 32中的一个，因此我在约束函数nonlcon里设置了一行逻辑或的等式约束条件  但是最后求解出来的时候  x【4】是带小数的12.XXX  并没有满足这个约束 这时候什么情况 我该怎么解决  求大佬指导!谢谢
2 w% }3 x$ w) ~- r' b, p! F# Q( X
/ ]% c4 A" `; Y/ k: T+ h程序如下：
clc;
clear all;
: P6 ^: ~6 z- Xclose all;

x0 = [100,100,100,100];
* I; B9 `8 K- o5 @) C
A=[-1,0,0,0;0,-1,0,0;1,-1,0,0;]
B = [-200;0;0]
1 O5 S5 r  s% dAeq = []
" B! {/ |/ D% F( MBeq = []
" Y+ t" u6 o, e% pLB =[200,0,4,12];
6 k# ]3 Q2 n: q. dUB =[1000,1000,20,32]  
[x, fval] = fmincon('fun', x0, A, B, Aeq, Beq, LB, UB,'nonlcon')

) p% y& N8 p* o9 L; C& Lfprintf('the minimum value in x:\n');
; g* G1 C0 ?3 J+ S+ Kfprintf('%f\n%f\n',x(1),x(2),x(3),x(4));
4 [$ r' ^( I% E. M2 }/ W8 cfprintf('the minimum value:\n');
$ F1 A! P5 |! ?- l. j- Pfprintf('%f\n',fval);
' d1 F6 t! D' }+ @) g- H


5 d& e4 l6 X# Q1 k, o4 vfunction [c,ceq] = nonlcon(x);
fai=1.2687-0.0291*(6500/x(1))+0.0001*(6500/x(1))^2;
c=[651000-(1/1.2)*(10*x(1)*x(2)+310*x(3)*(pi*x(4)^2)/4)*fai ;
: M- V! n4 `9 R. Q- Q8 T" r5 T# b-x(3)*(pi*x(4)^2)/(4*x(1)*x(2))+0.008;
0 ~+ X: k3 B- H# {# w9 @( N4.52-x(3)*(pi*(10*x(4))^2)/4];
& |! c) t$ Z& U# `+ e, Gceq=[mod(x(3),1);
9 C7 r2 H; s) F6 G5 F9 x" b4 s

end

$ R4 H: T, ?8 Q" |5 I6 A3 T

nkkopd

帮你顶一下

大小的小

你等式约束的方式写成(x-12)*(x-14)*(x-16)*(x-18)*(x-20)*(x-22)*(x-25)*(x-28)=0，设置下精度试试
, M3 F; Q% |8 L& T% d5 q$ i5 m

llbnmo

二楼正解。。

shuddkk