基于fmincon函数的最优化计算

x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon)

求解问题的标准型为
min F(X)
s.t
AX <= b（线性不等式约束）
AeqX = beq（线性等式约束）
G(x) <= 0（非线性不等式约束）
Ceq(X) = 0（非线性等式约束）
lb <= X <= ub（变量约束）
 

load case'mydata4'
P=[Pbi;Psj];
s=size(Pbi,1);
t=size(Psj,1);
A=eye(s+t);
b=zeros(s+t,1);
for i=1:s
    b(i,1)=Qli(i,1)-Qki(i,1);
end
for j=s+1:s+t
    b(j,1)=Qxj(j-s,1);
end
x0=zeros(s+t,1);%%将交易电量的初值赋值为0，其大小与买家、卖家个数一致
Eb=ones(s,1);
Es=ones(t,1)*(-1);
Aeq=[Eb',Es'];
beq=0;
vlb=x0;%%将交易电量的下限赋值为0
vub=[];
f=@(Q)P'*Q*-1;
[Q,f]=fmincon(f,x0,A,b,Aeq,beq,vlb,vub);