Commit 4e9cef1b authored by Gidon Lucian Bauer's avatar Gidon Lucian Bauer 🏳
Browse files

MSM Bonus 03 Final

parent 2635ab7f
......@@ -77,8 +77,12 @@ for i = 2:max_iter
% b)
if discrepancy(i-1) <= eps_alpha && k_star == -1
k_star = i;
disp(k_star);
disp(norm(A*iter(:,i) - u_f_meas, 2));
% Ausgabe von k*
disp(['k* = ', num2str(k_star)]);
% Ausgabe der Discrepancy
disp(['||A*iter(:,', num2str(k_star), ') - u_f_meas||_2 = ',...
num2str(discrepancy(i-1)), ' < ', num2str(eps_alpha), ...
' = eps_alpha']);
end
end
......@@ -104,8 +108,10 @@ ylabel('temperature');
legend('Interpreter', 'latex');
hold off;
startingpoint = 1;
figure(4);
loglog(discrepancy(1:end), xi_norm(1:end), '.', 'DisplayName', 'Discrepancy');
loglog(discrepancy(startingpoint:end), xi_norm(startingpoint:end), '.',...
'DisplayName', 'Discrepancy');
hold on;
xlabel('$||A \tilde{\xi_k} - \tilde{y}||_2$', 'Interpreter', 'latex');
ylabel('$||\tilde{\xi_k}||_2$', 'Interpreter', 'latex');
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment