ftp43.m

基于matlab的约束非线性规划算法库

function [f, J] = ftp43(x, ctrl)
%Call: [f J]=ftp43(x,ctrl)
%Evaluate the residuals corresponding
%to the objective in a nonlinear lsq problem.
%If ctrl>0 the Jacobian J is computed analytically.
%Objective: 0.5*norm2(f)
%The problem is defined as
%  min 0.5*norm2(f(x))
%      s.t. c(x)>=0
%where
%c(1)=-x(1)^2-x(2)^2-x(3)^2-x(4)^2-x(1)+x(2)-x(3)+x(4)+8;
%c(2)=-x(1)^2-2*x(2)^2-x(3)^2-2*x(4)^2+x(1)+x(4)+10;
%c(3)=-2*x(1)^2-x(2)^2-x(3)^2-2*x(1)+x(2)+x(4)+5;
%
f(1)=x(1)-2.5;
f(2)=x(2)-2.5;
f(3)=sqrt(2)*x(3)-10.5/sqrt(2);
f(4)=x(4)+3.5;
f=f(:);
%The Jacobian is set to the empty matix to indicate that it should
%be computed using forward differences
J=[];