example6_8.m

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

%编写目标函数
function y=target67(x)
y=-x(3)/(x(3)+x(4));

%编写非线性约束函数
function [c,ceq]=fun67(x)
%无非线性不等式约束
c=[];
   
r0=3;
r1=10*x(1)/(0.3333+x(1));
r2=4.1667*x(2)/(0.6667*(1+x(3)/(1+0.5*x(1)))+x(2));
r3=6.5*x(2)/(0.6429*(1+x(4)/(7*(1+0.6667*x(3)))))+...
    x(2)*(1+x(4)/(3.5+2.3333*x(3)));
r4=6.25*x(3)/(0.1875+x(3));
r5=3.516*x(4)/(0.3333+x(4));
%非线性等式约束
ceq(1)=r0-r1*x(5);
ceq(2)=r1*x(5)-r2*x(6)-r3*x(7);
ceq(3)=r2*x(6)-r4*x(9);
ceq(4)=r3*x(7)-r5*x(9);

%优化变量边界约束
lb=[0.1285 0.1603 0.2699 0.4668 0.2 0.2 0.2 0.2 0.2]';
ub=[0.1571 0.1959 0.3199 0.5706 5.0 5.0 5.0 5.0 5.0]';
%无限型约束
A=[];
b=[];
Aeq=[];
beq=[];
%设定初始值
x0=[0.1428 0.1781 0.2999 0.5187 0.5 0.5 0.5 0.5 0.5]';
%调用优化函数
[x,fval]=fmincon('target67',x0,A,b,Aeq,beq,lb,ub,'fun67')