目标函数最优解

zzz.dan

使用fmincon求目标函数的最优解时，其中x【4】的值为12 14 16 18 20 22 25 28 32中的一个，因此我在约束函数nonlcon里设置了一行逻辑或的等式约束条件  但是最后求解出来的时候  x【4】是带小数的12.XXX  并没有满足这个约束 这时候什么情况 我该怎么解决  求大佬指导!谢谢

& d# H+ Z& {" f1 ~! {2 x, |程序如下：
clc;
clear all;
- X+ Y/ C+ ]% B. M! O/ Tclose all;
8 Z2 [* q, d3 _; i; E
6 M/ G3 X5 y  L# o) Nx0 = [100,100,100,100];
8 {3 J" z' z1 ]
A=[-1,0,0,0;0,-1,0,0;1,-1,0,0;]
B = [-200;0;0]
& Y, ~$ S5 |' l. A  e5 FAeq = []
Beq = []
LB =[200,0,4,12];
. g! f1 W, m) ^8 d5 L7 t  CUB =[1000,1000,20,32]  
[x, fval] = fmincon('fun', x0, A, B, Aeq, Beq, LB, UB,'nonlcon')

6 c6 E# i  R* v% ufprintf('the minimum value in x:\n');
, g+ j0 ]# K# a4 f+ }fprintf('%f\n%f\n',x(1),x(2),x(3),x(4));
fprintf('the minimum value:\n');
' I, ?+ r! Z7 Q& c# ]5 Cfprintf('%f\n',fval);
) @: a( \# @$ z. y
8 ?" S/ J% l& R

function [c,ceq] = nonlcon(x);
: h9 J3 B7 v9 Z2 V- Qfai=1.2687-0.0291*(6500/x(1))+0.0001*(6500/x(1))^2;
# W; p0 Q6 Y5 a- Yc=[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;
9 I; G) [  Y1 f9 P8 l  p4.52-x(3)*(pi*(10*x(4))^2)/4];
+ [. V/ C2 Z5 W: f, Bceq=[mod(x(3),1);

9 p/ k/ d; ^' F
' m4 P' X/ a/ jend

nkkopd

帮你顶一下

大小的小

你等式约束的方式写成(x-12)*(x-14)*(x-16)*(x-18)*(x-20)*(x-22)*(x-25)*(x-28)=0，设置下精度试试
0 K0 C2 x' m3 @3 `" o

llbnmo

二楼正解。。

shuddkk