目标函数最优解

zzz.dan

使用fmincon求目标函数的最优解时，其中x【4】的值为12 14 16 18 20 22 25 28 32中的一个，因此我在约束函数nonlcon里设置了一行逻辑或的等式约束条件  但是最后求解出来的时候  x【4】是带小数的12.XXX  并没有满足这个约束 这时候什么情况 我该怎么解决  求大佬指导!谢谢
+ L7 y1 q. P5 H% G; T2 {
) ]7 o0 j3 i& V程序如下：
clc;
clear all;
8 u6 I; _8 x% z& G1 `close all;

x0 = [100,100,100,100];

A=[-1,0,0,0;0,-1,0,0;1,-1,0,0;]
B = [-200;0;0]
Aeq = []
Beq = []
LB =[200,0,4,12];
9 M2 A1 {9 [- b9 R0 K: I; B. q5 hUB =[1000,1000,20,32]  
[x, fval] = fmincon('fun', x0, A, B, Aeq, Beq, LB, UB,'nonlcon')

9 M$ U: [& J" t7 G0 X; Wfprintf('the minimum value in x:\n');
fprintf('%f\n%f\n',x(1),x(2),x(3),x(4));
- h" R7 E/ C5 I) ~5 H  \8 Lfprintf('the minimum value:\n');
1 V. X4 f1 ^$ }0 T5 F4 P% Wfprintf('%f\n',fval);
3 I6 G! Q9 _$ W) T' V
- K' z- u; J' O5 {
$ M% X% B6 f$ y. z% ?$ i
function [c,ceq] = nonlcon(x);
  S3 y5 a) m' B9 R& E$ U: j4 Afai=1.2687-0.0291*(6500/x(1))+0.0001*(6500/x(1))^2;
: M3 u) Z+ V$ K5 k% Gc=[651000-(1/1.2)*(10*x(1)*x(2)+310*x(3)*(pi*x(4)^2)/4)*fai ;
-x(3)*(pi*x(4)^2)/(4*x(1)*x(2))+0.008;
  K2 O. [' G; ~/ G5 k4.52-x(3)*(pi*(10*x(4))^2)/4];
ceq=[mod(x(3),1);


" j: V, x: g) O, n3 [4 h6 u7 wend
; e( S2 S) Y8 Q9 L1 `7 w

nkkopd

帮你顶一下

大小的小

你等式约束的方式写成(x-12)*(x-14)*(x-16)*(x-18)*(x-20)*(x-22)*(x-25)*(x-28)=0，设置下精度试试
" e( L/ G9 m0 A% c* E

llbnmo

二楼正解。。

shuddkk