函数：[x, fval] = fmincon(FUN, X0, A, B, Aeq, Beq, LB, UB, NONLCON)
返回的x：是一个向量——在取得目标函数最小时各个xi的取值；
返回的fval：目标函数的最小值；
参数FUN：目标函数；
参数X0：向量x的初始值；
参数A：线性不等式约束的系数矩阵，若没有线性不等式约束，则A = []；
参数B：线性不等式约束右端的常数列，若没有线性不等式约束，则B = []；
参数Aeq：线性等式约束的系数矩阵，若没有等式约束，则Aeq = []；
参数Beq：线性等式约束右端的常数列，若没有等式约束，则Beq = []；
参数LB：x的下界，常遇到的x1, x2, x3 >= 0，0就是下界，可用zeros(3, 1)生成一个3行1列的向量来表示，其中向量的每个元素的值为0；
参数UB：x的上界；
参数NONLCON：用.m文件定义的非线性向量函数C(x), Ceq（x)。
-------------------------------------------------------------------------------------------------------------
求f(x)的最小值及f(x)取得最小值时x的取值：
min f(x) = x1^2 + x2^2 + 8
约束条件：
x1^2 - x2 >= 0
-x1 - x2^2 + 2 = 0
x1, x2 >= 0
-------------------------------------------------------------------------------------------------------------
编写目标函数.m文件（不要写在主.m文件里）FUN.m（文件名要与函数名相同）
%写下面两行就可以了 function f = FUN(x); f = x(1)^2 + x(2)^2 + 8;
保存。
编写非线性约束向量函数.m文件（不要写在主.m文件里）NONLCON.m（文件名要与函数名相同）
%写下面3行就可以了 function [g, ceq] = NONLCON(x); g = -x(1)^2 + x(2); %非线性不等式约束 ceq = -x(1) - x(2)^2 + 2; %非线性等式约束
保存。
最后编写主体.m文件
X0 = rand(2, 1);
A = [];
B = [];
Aeq = [];
Beq = [];
LB = zeros(2, 1);
UB = [];
[x, fval] = fmincon('FUN', X0, A, B, Aeq, Beq, LB, UB, 'NONLCON')
保存运行，得
x =
1.0000
1.0000
fval =
10