目标函数最优解

zzz.dan

使用fmincon求目标函数的最优解时，其中x【4】的值为12 14 16 18 20 22 25 28 32中的一个，因此我在约束函数nonlcon里设置了一行逻辑或的等式约束条件  但是最后求解出来的时候  x【4】是带小数的12.XXX  并没有满足这个约束 这时候什么情况 我该怎么解决  求大佬指导!谢谢
! u6 V# i. X5 I4 D1 w! U0 R/ T6 _) y
程序如下：
clc;
( N+ p. [. l& V4 pclear all;
close all;

x0 = [100,100,100,100];
) l2 y2 j$ R$ ^
, \9 O% U3 R0 x( V9 a, zA=[-1,0,0,0;0,-1,0,0;1,-1,0,0;]
B = [-200;0;0]
Aeq = []
Beq = []
LB =[200,0,4,12];
UB =[1000,1000,20,32]  
[x, fval] = fmincon('fun', x0, A, B, Aeq, Beq, LB, UB,'nonlcon')

- j$ o- _& W: B1 yfprintf('the minimum value in x:\n');
fprintf('%f\n%f\n',x(1),x(2),x(3),x(4));
fprintf('the minimum value:\n');
8 Y( W* D7 @; L8 S" }+ Y( M" Kfprintf('%f\n',fval);
2 a! ^& B& f1 F9 |$ }4 N! g
1 @: y( W0 O1 F1 s$ U! m
2 d$ @* n2 \0 I( R' N+ S
& a3 Q* J! I' f' Pfunction [c,ceq] = nonlcon(x);
0 u. ~7 r) c  `5 }, @8 e5 Sfai=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 ;
% Q8 l9 B/ e, s" e  @9 y' l-x(3)*(pi*x(4)^2)/(4*x(1)*x(2))+0.008;
5 ?' N9 h/ l: t/ N4.52-x(3)*(pi*(10*x(4))^2)/4];
& R0 `- y  b9 `7 e' R7 Xceq=[mod(x(3),1);

; U7 t, ~4 z( Q1 [5 P9 C
2 E: U* g$ B1 n9 R( ]end

nkkopd

帮你顶一下

大小的小

- M( g5 o) N( E; i6 l你等式约束的方式写成(x-12)*(x-14)*(x-16)*(x-18)*(x-20)*(x-22)*(x-25)*(x-28)=0，设置下精度试试
* Z' g2 m: {# l; d

llbnmo

二楼正解。。

shuddkk