A=[5/12 -1 1 0 0 0;5/2 1 0 1 0 0;-1 0 0 0 1 0;0 1 0 0 0 1];
b=[35/12 35/2 5 5]';
c=[-1 -1 0 0 0 0]';
N=3;
M=1000;
W=M*A'*A;
X=rands(6,1);
K=3;
n=0;
for i=1:1:(N+K)^2
if abs(W(i))>n
n=abs(W(i));
else
continue;
end
end
n
m=N*n;
h1=1/m;
a=M*A'*b;
q=c-a;
f=X'*A'*A*X-2*b'*A*X+b'*b;
g=0.00000001;
e=1;
k=0;
while abs(e)>g
E1=M*X'*A'*A*X+(c'-2*M*b'*A)*X+M*b'*b;
y=X-h1*(M*A'*A*X+q);
X(1)=y(1);
X(2)=y(2);
for j=3:1:6
X(j)=max(y(j),0);
end
E2=M*X'*A'*A*X+(c'-2*M*b'*A)*X+M*b'*b;
e=E2-E1;
E1=E2;
k=k+1;
end
E1
k
X
评论0