007最速下降和牛顿法.pptVIP

　　Matlab程序示例 计算 Min ||A-X||F S.t. Rank(X)=1 X≥0 其中 A 是m*n的矩阵 function fval=opt_2(A) %直接调用约束优化函数fmincon [m n]=size(A); [U,S,V] = svd(A); S11 = S(1,1); x=FMINCON(@(x) myfun_2(x,A,S11),ones(n+m),[],[],[],[],zeros(n+m),[]); fval = myfun_2(x,A,S11); fval = sqrt(fval); function F = myfun_2(x,A,S11) [m n]=size(A); F=0; for i=1:m for j=1:n F = F+(A(i,j)-x(i)*S11*x(m+j))^2; end end m=10 n=20 残差 14.38 CPU 58.032s function fval=opt_3(A) %直接调用无约束优化函数fminunc [m n]=size(A);