《最优化方法及其Matlab程序设计》马昌凤 部分习题答案

第12章 序列二次规划法

6.用SQP方法的Matlab程序求解下列优化问题：

（1）

%%% 目标函数 f(x) %%%%%%%%%%%
function f=f1(x)
% f=-pi*x(1)^2*x(2);
f = (x(1)-2)^4 + (x(1)-2*x(2))^2;
%%%% 目标函数 f(x) 的梯度%%%%%
function df=df1(x)
% df=[-2*pi*x(1)*x(2), -pi*x(1)^2]';
df = [4*(x(1)-2)^3+2*(x(1)-2*x(2)); -4*(x(1)-2*x(2))];

%%%% 约束函数 %%%%%%
function [h,g]=cons(x)
% h=[pi*x(1)*x(2)+pi*x(1)^2-150];
h=[];
% g=[x(1);x(2)];
g = [x(1)^2 + x(2)];

%% 约束函数 Jacobi矩阵%%%%
function [dh,dg]=dcons(x)
% dh=[pi*x(2)+2*pi*x(1), pi*x(1)];
dh = [];
% dg=[1 0; 0 1];
dg = [2*x(1), 1];

调用：

x0 = [0 0]';
mu0 = [ ];
lam0 = [0];
[x,mu,lam,val,k] = sqpm(x0,mu0,lam0);

结果

（2）

%%% 目标函数 f(x) %%%%%%%%%%%
function f=f1(x)
% f=-pi*x(1)^2*x(2);
% f = (x(1)-2)^4 + (x(1)-2*x(2))^2;
s = -x(1)-x(2);
f = exp(s) + x(1)^2 + 2*x(1)*x(2) + x(2)^2 + 2*x(1) + 6*x(2);

%%%% 目标函数 f(x) 的梯度%%%%%
function df=df1(x)
s = -x(1)-x(2);
df = [-exp(s)+2*x(1)+2*x(2)+2; -exp(s)+2*x(1)+2*x(2)+6];

%%%% 约束函数 %%%%%%
function [h,g]=cons(x)
% h=[pi*x(1)*x(2)+pi*x(1)^2-150];
h=[];
g = [2-x(1)-x(2),x(1),x(2)];

%% 约束函数 Jacobi矩阵%%%%
function [dh,dg]=dcons(x)
% dh=[pi*x(2)+2*pi*x(1), pi*x(1)];
dh = [];
dg = [-1 -1;1 0;0 1];

调用

x0 = [2 0]';
mu0 = [ ];
lam0 = [0 0 0];
[x,mu,lam,val,k] = sqpm(x0,mu0,lam0);

结果