x0 = [3.0000 5.0000 0.0707 7.0000 0 0.9293 0 0 3.9293 0 6.0000 ...
10.0707 6.3875 4.3943 5.7511 7.1867]';
A=[1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0;0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0];
b=[20;20];
Aeq=[1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0;...
0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0;...
0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0;...
0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0;...
0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0;...
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0];
beq=[3;5;4;7;6;11];
VLB=[0;0;0;0;0;0;0;0;0;0;0;0];
VUB=[];
[x,fval]=fmincon('f1',x0,A,b,Aeq,beq,VLB,VUB)
function [x,val,k]=bfgs(fun,gfun,x0,varargin)
%功能: 用BFGS算法求解无约束问题: min f(x)
%输入: x0是初始点, fun, gfun分别是目标函数及其梯度;
% varargin是输入的可变参数变量, 简单调用bfgs时可以忽略它,
% 但若其它程序循环调用该程序时将发挥重要的作用
%输出: x, val分别是近似最优点和最优值, k是迭代次数.
maxk=500; %给出最大迭代次数
rho=0.55; sigma1=0.4; epsilon1=1e-5;
k=0; n=length(x0);
Bk=eye(n); %Bk=feval('Hess',x0);
while(k<maxk)
gk=feval(gfun,x0,varargin{:}); %计算梯度
if(norm(gk)<epsilon1), break; end %检验终止准则
dk=-Bk\gk; %解方程组, 计算搜索方向
m=0; mk=0;
while(m<20) % 用Armijo搜索求步长
newf=feval(fun,x0+rho^m*dk,varargin{:});
oldf=feval(fun,x0,varargin{:});
if(newf<oldf+sigma1*rho^m*gk'*dk)
mk=m; break;
end
m=m+1;
end
%BFGS校正
x=x0+rho^mk*dk;
sk=x-x0; yk=feval(gfun,x,varargin{:})-gk;
if(yk'*sk>0)
Bk=Bk-(Bk*sk*sk'*Bk)/(sk'*Bk*sk)+(yk*yk')/(yk'*sk);
end
k=k+1; x0=x;
end
function g=df1(x)
g = [sqrt((x(13)-1.25)^2+(x(14)-1.25)^2), sqrt((x(13)-8.75)^2+(x(14)-0.75)^2),...
sqrt((x(13)-0.5)^2+(x(14)-4.75)^2), sqrt((x(13)-5.75)^2+(x(14)-5.0)^2),...
sqrt((x(13)-3.0)^2+(x(14)-6.5)^2), sqrt((x(13)-7.25)^2+(x(14)-7.25)^2),...
sqrt((x(15)-1.25)^2+(x(16)-1.25)^2), sqrt((x(15)-8.75)^2+(x(16)-0.75)^2),...
sqrt((x(15)-0.5)^2+(x(16)-4.75)^2), sqrt((x(15)-5.75)^2+(x(16)-5.0)^2),...
sqrt((x(15)-3.0)^2+(x(16)-6.5)^2), sqrt((x(15)-7.25)^2+(x(16)-7.25)^2)]';
g(13)=x(1)*(x(13)-1.25)/g(1)+x(2)*(x(13)-8.75)/g(2)+x(3)*(x(13)-0.5)/g(3)+...
x(4)*(x(13)-5.75)/g(4)+x(5)*(x(13)-3.0)/g(5)+x(6)*(x(13)-7.25)/g(6);
g(14)=x(1)*(x(14)-1.25)/g(1)+x(2)*(x(14)-0.75)/g(2)+x(3)*(x(14)-4.75)/g(3)+...
x(4)*(x(14)-5.0)/g(4)+x(5)*(x(14)-6.5)/g(5)+x(6)*(x(14)-7.25)/g(6);
g(15)=x(7)*(x(15)-1.25)/g(7)+x(8)*(x(15)-8.75)/g(8)+x(9)*(x(15)-0.5)/g(9)+...
x(10)*(x(15)-5.75)/g(10)+x(11)*(x(15)-3.0)/g(11)+x(12)*(x(15)-7.25)/g(12);
g(16)=x(7)*(x(16)-1.25)/g(7)+x(8)*(x(16)-0.75)/g(8)+x(9)*(x(16)-4.75)/g(9)+...
x(10)*(x(16)-5.0)/g(10)+x(11)*(x(16)-6.5)/g(11)+x(12)*(x(16)-7.25)/g(12);
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- 10.
- 11.
- 12.
- 13.
- 14.
- 15.
- 16.
- 17.
- 18.
- 19.
- 20.
- 21.
- 22.
- 23.
- 24.
- 25.
- 26.
- 27.
- 28.
- 29.
- 30.
- 31.
- 32.
- 33.
- 34.
- 35.
- 36.
- 37.
- 38.
- 39.
- 40.
- 41.
- 42.
- 43.
- 44.
- 45.
- 46.
- 47.
- 48.
- 49.
- 50.
- 51.
- 52.
- 53.
- 54.
- 55.
- 56.
- 57.
- 58.
- 59.
- 60.