/
aytsan_ns
/
Num_meth
Обзор
Документация
Войти
/
aytsan_ns
/
Num_meth
Код
Запросы
0
Задачи
Пакеты
0
Релизы
0
Аналитика
Безопасность
master
Lab14/Lab_14.m
75 строк
1 KB
aytsan-ns
Исправила ошибки в 14 работе. Изменила метод Гаусса для решения СЛАУ на метод прогонки. Обновила отчет.
10 сен 2024, 15:01
10 сен 2024, 15:01
468dd84
Код
Авторство
О чём код?
f = @(x) cos(x); result1 = load('result1.txt'); result2 = load('result2.txt'); error1 = load('error1.txt'); error2 = load('error2.txt'); x_result1 = result1(:,1) y_result1 = result1(:,2) x_result2 = result2(:,1) y_result2 = result2(:,2) x_error1 = error1(:,1) y_error1 = error1(:,2) x_error2 = error2(:,1) y_error2 = error2(:,2) figure; hold on; fplot(f, [0 pi/2], 'b'); plot(x_result1, y_result1, 'y'); plot(x_result2, y_result2, 'r'); grid on legend('������ �������', 'h = pi/6', 'h = pi/12', 'Location', 'Best'); hold off; figure; hold on; plot(x_error1, y_error1, 'b'); plot(x_error2, y_error2, 'r'); grid on legend('h = pi/6', 'h = pi/12', 'Location', 'Best'); ylabel('������', 'Interpreter', 'none'); hold off; result = load('result.txt'); error = load('error.txt'); x_result = result(:,1) y_result = result(:,2) x_error = error(:,1) y_error = error(:,2) figure; hold on; fplot(f, [0 pi/2], 'b'); plot(x_result, y_result, 'r'); plot(x_result, y_result, 'or'); grid on legend('������ �������', 'N = 20', 'Location', 'Best'); hold off; figure; hold on; plot(x_error, y_error, 'b'); grid on ylabel('������', 'Interpreter', 'none'); hold off; result_error_eps = load('result_error_eps.txt'); x_result_error_eps = result_error_eps(:,1) y_result_error_eps = result_error_eps(:,2) figure; hold on; plot(x_result_error_eps, y_result_error_eps, 'r'); grid on set(gca, 'XScale', 'log'); set(gca, 'YScale', 'log'); ylabel('����������� �����������', 'Interpreter', 'none'); xlabel('��������', 'Interpreter', 'none'); hold off;