目标函数最优解

zzz.dan

使用fmincon求目标函数的最优解时，其中x【4】的值为12 14 16 18 20 22 25 28 32中的一个，因此我在约束函数nonlcon里设置了一行逻辑或的等式约束条件  但是最后求解出来的时候  x【4】是带小数的12.XXX  并没有满足这个约束 这时候什么情况 我该怎么解决  求大佬指导!谢谢
9 S% T. E+ h2 A7 W& c
9 [5 m# A: d: v- g$ M9 ^- `/ ^程序如下：
, M0 R1 k' @4 K9 d; e& N5 l8 L' mclc;
, L$ O. e% n! ]/ Kclear all;
close all;

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

A=[-1,0,0,0;0,-1,0,0;1,-1,0,0;]
B = [-200;0;0]
5 {* h' h3 Z" d' N2 o) BAeq = []
$ l. O: R; M+ B0 {; P  iBeq = []
# B7 ~3 K' C0 U  J( x, @LB =[200,0,4,12];
, b( \  \1 O4 w9 IUB =[1000,1000,20,32]  
[x, fval] = fmincon('fun', x0, A, B, Aeq, Beq, LB, UB,'nonlcon')
) D% K* k0 {( S  A: p( e
fprintf('the minimum value in x:\n');
fprintf('%f\n%f\n',x(1),x(2),x(3),x(4));
fprintf('the minimum value:\n');
fprintf('%f\n',fval);
: S4 p0 b6 o' C9 X! d


function [c,ceq] = nonlcon(x);
3 O1 p" C3 E! p8 V9 ?- hfai=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 ;
-x(3)*(pi*x(4)^2)/(4*x(1)*x(2))+0.008;
* V3 v. C* C" ~4.52-x(3)*(pi*(10*x(4))^2)/4];
ceq=[mod(x(3),1);


end

nkkopd

帮你顶一下

大小的小

你等式约束的方式写成(x-12)*(x-14)*(x-16)*(x-18)*(x-20)*(x-22)*(x-25)*(x-28)=0，设置下精度试试

llbnmo

二楼正解。。

shuddkk