example6_9.m

MATLAB7辅助控制系统设计与仿真 第6章

%参数a
a=1;
%目标函数系数
f=[0 0 0 0 0 0 0 0 0 0 0 -6 0 0]';
%等式约束矩阵系数
Aeq=[0 0 0 0 0 0 0 0 0 0 0 6 2 1;
      1 1 0 0 0 0 0 0 0 0 0 0 0 0;
      1 0 -1 -1 -1 0 0 0 0 0 0 0 0 0 ;
      0 1 1 -1 0 0 -1 0 0 0 0 0 0 0;
      0 0 0 0 1 -1 1 0 0 0 0 0 0 0 ;
      0 0 0 0 0 1 0 0 1 1 0 0 0 0 ;
      0 0 1 0 0 0 0 -1 0 0 0 0 -1 0;
      0 0 0 0 0 0 -1 1 1 0 0 0 0 0 ;
      0 0 0 0 0 0 0 0 -1 -1 -1 0 0 0;
      0 0 0 1 0 0 0 0 0 0 -1 -1 0 0 ;
      -1 1 0 0 0 0 0 0 0 2 0 0 0 0 ];
beq=[3*a a 0 0 0 0 0 0 0 0 0 ]';
%无不等式约束
A=[];
b=[];
%边界约束
lb=[0 0 0 0 0 0 0 0 0 0 0 0 0 0]';
ub=[];
%初始点
x0=[0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ]';
%调用函数linprog
[x,fval]=linprog(f,A,b,Aeq,beq,lb,ub,x0)