目标函数最优解

zzz.dan

使用fmincon求目标函数的最优解时，其中x【4】的值为12 14 16 18 20 22 25 28 32中的一个，因此我在约束函数nonlcon里设置了一行逻辑或的等式约束条件  但是最后求解出来的时候  x【4】是带小数的12.XXX  并没有满足这个约束 这时候什么情况 我该怎么解决  求大佬指导!谢谢
1 Z2 }! C2 J# e4 A- \7 B
* M1 Q1 T* U6 a/ Q程序如下：
clc;
clear all;
close all;

x0 = [100,100,100,100];
" O2 Q; B* g( j6 E! _: b  d" M$ b
: m3 N3 y9 z& W9 |7 RA=[-1,0,0,0;0,-1,0,0;1,-1,0,0;]
$ O  ^  h, p1 D1 W6 L2 aB = [-200;0;0]
9 a2 x( E8 [4 b( GAeq = []
) W, s2 D2 r' P  [) kBeq = []
LB =[200,0,4,12];
: _1 a3 t. |" j% BUB =[1000,1000,20,32]  
- z; y2 V: \; S7 E+ D" u[x, fval] = fmincon('fun', x0, A, B, Aeq, Beq, LB, UB,'nonlcon')
7 _$ }' f' u( ^; f- N
4 c% y4 W8 d3 m) u' Gfprintf('the minimum value in x:\n');
- ~  |5 k/ j" Xfprintf('%f\n%f\n',x(1),x(2),x(3),x(4));
fprintf('the minimum value:\n');
$ v  {- H. d/ s; Z7 h) @! {! Gfprintf('%f\n',fval);
5 g) t. E: i* y! F) {  o* r

- z* ^! F) Q& F, f* }" \% E7 c8 r
7 H# ~  W7 p' b, U" tfunction [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 ;
9 t! c: C  b& ]) _7 ~7 ?5 K4 \-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);
- e7 T2 ?/ i' e. c2 c3 x# m

% q: f5 `! q- n0 H. J! J  l% N/ Hend

nkkopd

帮你顶一下

大小的小

/ K: N- p4 a- h2 l* c你等式约束的方式写成(x-12)*(x-14)*(x-16)*(x-18)*(x-20)*(x-22)*(x-25)*(x-28)=0，设置下精度试试
- t2 _: t8 A% G& d% p3 i% ]

llbnmo

二楼正解。。

shuddkk