目标函数最优解

zzz.dan

使用fmincon求目标函数的最优解时，其中x【4】的值为12 14 16 18 20 22 25 28 32中的一个，因此我在约束函数nonlcon里设置了一行逻辑或的等式约束条件  但是最后求解出来的时候  x【4】是带小数的12.XXX  并没有满足这个约束 这时候什么情况 我该怎么解决  求大佬指导!谢谢
! o2 f: [; e# g
& h# C) g) Y8 |# Z0 u: Z% w程序如下：
clc;
clear all;
close all;

x0 = [100,100,100,100];
* F- e0 b% c; g! T- s9 i' R8 s
A=[-1,0,0,0;0,-1,0,0;1,-1,0,0;]
B = [-200;0;0]
1 s: n6 [+ L- mAeq = []
Beq = []
LB =[200,0,4,12];
UB =[1000,1000,20,32]  
# f! a0 j: K$ j[x, fval] = fmincon('fun', x0, A, B, Aeq, Beq, LB, UB,'nonlcon')

fprintf('the minimum value in x:\n');
+ y0 u( H4 b  @. w: F) Ffprintf('%f\n%f\n',x(1),x(2),x(3),x(4));
fprintf('the minimum value:\n');
# d# n! p; k9 \5 y$ m9 v2 nfprintf('%f\n',fval);
, @) u& X6 Z+ W9 D9 r) R

9 E* k9 z" q$ b* ?
function [c,ceq] = nonlcon(x);
4 ]) `7 U, v1 a' D9 Rfai=1.2687-0.0291*(6500/x(1))+0.0001*(6500/x(1))^2;
$ p$ \" [) }2 g8 D4 o* Zc=[651000-(1/1.2)*(10*x(1)*x(2)+310*x(3)*(pi*x(4)^2)/4)*fai ;
1 R8 c9 N2 G% D1 _; q% `-x(3)*(pi*x(4)^2)/(4*x(1)*x(2))+0.008;
; y. d7 C9 r% w4 m( i2 Z* I4.52-x(3)*(pi*(10*x(4))^2)/4];
ceq=[mod(x(3),1);


- T- c  a% b1 ^, O" vend

. c; s- K. u  h7 _. c, [5 [3 w3 v) m) _

nkkopd

帮你顶一下

大小的小

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

llbnmo

二楼正解。。

shuddkk