梯度投影法 matlab 源码程序
function [x,minf] = minRosen(f,A,b,x0,var,eps)
format long;
if nargin == 5
eps = 1.0e-6;
end
syms l;
x0 = transpose(x0);
n = length(var);
sz = size(A);
m = sz(1);
gf = jacobian(f,var);
bConti = 1;
while bConti
k = 0;
s = 0;
A1 = A;
A2 = A;
b1 = b;
b2 = b;
for i=1:m