%% Gradiente Coniugato MODIFICATO, pcg, direzioni, angoli, discesa
% Gradiente Coniugato MODIFICATO
x = x0;
r = b - A*x;
p = r;
for iter = 1:3
alfa = (p' * r) / (p' * A*p);
x = x + alfa * p;
r = r - alfa * A*p;
beta = (p' * A * r) / (p' * A*p);
p_Old = p;
p = r - beta * p;
% Bonus: Trovo angoli fra direzioni successive
theta = acos((p_Old' * p) / (norm(p) * norm(p_Old)));
theta = rad2deg(theta); % Converto in deg°
fprintf('Iterazione: %d; norm(p) = %.4f; theta[°] = %.4f; \n', ...
iter, norm(p), theta);
end
%% Metodo del gradiente a mano
phi = @(y) 1/2*y'*A*y-y'*b;
x0 = b;
r0 = b-A*x0;
alpha0 = (r0'*r0)/(r0'*A*r0);
x1 = x0 + alpha0*r0;
disp(phi(x0)); disp(phi(x1));
