Lab11/Lab_11.m - master - aytsan_ns/Num_meth

Lab11/Lab_11.m
56 строк1 KB
aytsan-ns
Нашла и исправила ошибки в 11 работе.
09 сен 2024, 19:39
09 сен 2024, 19:3958a57a6
error_rado = dlmread('error_rado.txt');
error_rich = dlmread('error_rich.txt');

X = [0.1 0.01 0.001 0.0001 0.00001 0.000001 0.0000001]

figure;
hold on;
plot(X, error_rado, 'b');
plot(X, error_rich, 'y');
grid on
set(gca, 'XScale', 'log');
set(gca, 'YScale', 'log');
x = linspace(0.0000001, 0.1, 100);
plot(x, x, '--r');
legend('����', '���������', '����������a', 'Location', 'Best');
xlabel('��������', 'Interpreter', 'none');
ylabel('������', 'Interpreter', 'none');
hold off;

iter_rado = dlmread('iter_rado.txt');

figure;
hold on;
plot(X, iter_rado, 'b');
grid on
set(gca, 'XScale', 'log');
set(gca, 'YScale', 'log');
xlabel('��������', 'Interpreter', 'none');
ylabel('����� ��������', 'Interpreter', 'none');
hold off;

error_h_rado = load('error_h_rado.txt');

h_rado = error_h_rado(:,1)
error_rado = error_h_rado(:,2)
log_h_rado = log2(h_rado);
log_error_rado = log2(error_rado);

figure;
hold on;
plot(log_h_rado, log_error_rado, 'b');
grid on;

coeffs = polyfit(log_h_rado, log_error_rado, 1);
p = coeffs(1);  % ������� �������� (������ �����)
log_C = coeffs(2);  % �������� ��������� C
C = 2^log_C;  % ���������

fit_error = polyval(coeffs, log_h_rado);
plot(log_h_rado, fit_error, 'r--', 'LineWidth', 1);

legend('����', '�������� �������������', 'Location', 'best');
title(sprintf('������� ��������: %.2f, ���������: %.2e', p, C));
xlabel('log_2(h)', 'Interpreter', 'none');
ylabel('log_2(����������� ������)', 'Interpreter', 'none');
hold off;