非线性规划.doc

文档介绍：实验报告10
一、实验项目:非线性规划
二、实验目的和要求
掌握使用MATLAB优化工具箱求解非线性规划的方法;
练****建立实际问题的非线性规划模型。
三、实验内容

Min f=exp(x1 x2 x3 x4 x5)
. ∑(i=1~5)xi 2=10;
x2x3-5x4x5=0;
x1 3+x2 3+1=0;
<=xi<=,i=1,2;
<=xi<=,i=3,4,5;
fun函数:
function f=fun(x)
f=exp(x(1)*x(2)*x(3)*x(4)*x(5));
cont函数:
function [G,Geq]=cont(x);
G=[];
Geq(1)=x(1)^2+x(2)^2+x(3)^2+x(4)^2+x(5)^2-10;
Geq(2)=x(2)*x(3)-5*x(4)*x(5);
Geq(3)=x(1)^3+x(2)^3+1;
代码:
x0=[1;1;1;1;1];
A=[];b=[];
Aeq=[];beq=[];
lb=[-;-;-;-;-];
ub=[;;;;];
[x,fval,exitfllag,output]=fmincon(***@fun,x0,A,b,Aeq,beq,lb,ub,***@cont)
结果:
x =-




fval =
exitfllag =5
output =
iterations: 18
ount: 136
stepsize: 1
algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
firstorderopt: -005
cgiterations: []
message: [1x172 char]

初始点x0=[20,20],设(ai,bi)(i=1~20)为第i栋住宅楼的坐标;
a=[ ];
b=[ ];
fun函数:
function f=fun(x)
a=[ ];
b=[ ];
for i=1:15
f=abs(x(1)-a(i))+abs(x(2)-b(i));
end
count函数:
function [G,