Lab12/Lab_12.m - master - aytsan_ns/Num_meth

Lab12/Lab_12.m

60 строк1 KB

aytsan-ns

Внесла изменения в код 12 работы. Сделала шаг адаптивным и уменьшающимся в 2 раза, если погрешность больше точности и увеличивающимся в 2 раза, если погрешность в 2 раза меньше точности. Также изменился график изменения шага. Прикрепила актуальный отчет.

23 сен 2024, 21:47

23 сен 2024, 21:47028bc1a

f = @(x) x.*sin(x);

result_h1 = load('result_h1.txt');
result_h2 = load('result_h2.txt');
error_result_h1 = load('error_result_h1.txt');
error_result_h2 = load('error_result_h2.txt');

x_result_h1 = result_h1(:,1)
y_result_h1 = result_h1(:,2)
x_result_h2 = result_h2(:,1)
y_result_h2 = result_h2(:,2)

x_error_result_h1 = error_result_h1(:,1)
y_error_result_h1 = error_result_h1(:,2)
x_error_result_h2 = error_result_h2(:,1)
y_error_result_h2 = error_result_h2(:,2)

figure;
hold on;
fplot(f, [pi/2 2*pi], 'b');
plot(x_result_h1, y_result_h1, 'y');
plot(x_result_h2, y_result_h2, 'r');
grid on
legend('������ �������', 'h = pi/2', 'h = pi/6', 'Location', 'Best');
hold off;

figure;
hold on;
plot(x_result_h1, y_error_result_h1, 'b');
plot(x_result_h2, y_error_result_h2, 'r');
grid on
legend('h = pi/2', 'h = pi/6', 'Location', 'Best');
ylabel('������', 'Interpreter', 'none');
hold off;

result_h = load('result_h.txt');
x_result_h = result_h(:,1)
y_result_h = result_h(:,2)

figure;
hold on;
plot(x_result_h, y_result_h, 'b');
grid on
%set(gca, 'YScale', 'log');
ylabel('���', 'Interpreter', 'none');
hold off;

error_eps = load('error_eps.txt');
x_error_eps = error_eps(:,1)
y_error_eps = error_eps(:,2)

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