ТЕМА 1. СКАЛЯРНЫЕ НЕЛИНЕЙНЫЕ УРАВНЕНИЯ
Задание 1
Исследовать скорость сходимости (требуемое количество итераций ) метода а) дихотомии и хорд; б) Ньютона; в) секущих и г) простых итераций для численного решения уравнения Кеплера
с точностью до
в зависимости от параметров уравнения
и
на сетке
и
. Представить результаты графически как поверхность (или карту линий уровней) зависимости
.
Рекомендации: ;
(начальное приближение); для метода секущих
.
Задание 2
С точностью до методом I) дихотомии; II) секущих найти все корни многочлена а) Лежандра; б) Чебышева
степени
на отрезке
. Многочлены задаются рекуррентно: для варианта а)
для варианта б)
Рассмотреть случаи . Сравнить (для варианта б) вычисленные корни полинома Чебышева с их точными значениями
Рекомендации: I) для метода дихотомии выполнять поиск корней на начальных отрезках c граничными значениями
и
; II) для метода секущих принимать начальные приближения
и
; для каждой пары приближений найти соответствующий корень многочлена. Из всех найденных корней выбрать
различных.
Задание 3
Определить область сходимости метода Ньютона с точностью до к нулевому корню функции а)
; б)
(константы произвольные).
Рекомендации: рассмотреть начальные приближения на сетке
; для каждого найти решение соответствующего уравнения; первое наименьшее из начальных приближений, для которого численное решение ненулевое, будет задавать правую границу искомой области сходимости метода к нулевому решению. Для определения левой границы выполнить то же самое на сетке
.
Задание 4
Исследовать зависимость получаемого методом Ньютона численного решения уравнения а)
; б)
, а также скорости сходимости (количество итераций
) с точностью до
от начального приближения
Представить зависимости
и
графически.
ТЕМА 2. СИСТЕМЫЛИНЕЙНЫХ УРАВНЕНИЙ И
СОБСТВЕННЫЕ ЗНАЧЕНИЯ
Задание 1
Программно реализовать метод Гаусса для численного решения систем линейных уравнений произвольного порядка. Опробовать метод на примере системы уравнений четвертого порядка
, где
Сравнить численное решение с точным. Оценить вычислительные ошибки численного решения.
Задание 2
Программно реализовать метод Гаусса для вычисления определителя матрицы произвольного порядка. Опробовать метод на примере матрицы четвертого порядка
Сравнить численное решение с точным значением определителя. Оценить ошибку численного решения.
Задание 3
Программно реализовать метод Гаусса для вычисления обратной матрицы произвольного порядка. Опробовать метод на примере матрицы четвертого порядка
Сравнить численно полученную обратную матрицу с точной
Оценить вычислительные ошибки численного обращения матрицы.
Задание 4
Программно реализовать метод простых итераций и метод Зейделя для численного решения систем линейных уравнений произвольного порядка. Опробовать методы на примере системы уравнений четвертого порядка, где
Сравнить скорость сходимости методов при достижении точности решения . В качестве начального приближения выбирать нулевое
.
Задание 5
Программно реализовать метод Зейделя для численного решения систем линейных уравнений произвольного порядка, где предполагается симметризация Гаусса:
. Опробовать метод на примере нормальной системы уравнений четвертого порядка, где
Сравнить скорость сходимости метода при достижении точности решения с симметризацией и без. В качестве начального приближения выбирать нулевое
.
Задание 6
Степенным методом определить максимальное собственное число и соответствующий собственный вектор
для нормальных матриц
Решение находить с точностью до . По количеству выполненных итераций оценить скорость сходимости метода. Вычислить число обусловленности матрицы как
Задание 7
Программно реализовать метод вращений Якоби для вычисления собственных чисел . Решение находить с точностью до
. Опробовать метод на примере матрицы четвертого порядка
По количеству выполненных итераций оценить скорость сходимости метода. Вычислить число обусловленности матрицы.
Задание 8
Программно реализовать степенной метод для вычисления максимального собственного числа (нормальной) матрицы произвольного порядка
с точностью до
. Найти норму матрицы четвертого порядка
как , где
― максимальное собственное число нормальной матрицы
. Вычислить норму иным способом, как
на множестве случайных векторов , равномерно распределенных внутри гиперкуба
. Выполнить
испытаний. Сравнить численные значения норм.
Задание 9
Программно реализовать метод вращений Якоби для вычисления обратной нормальной матрицы произвольного порядка. Решение находить с точностью вычисления собственных чисел до
. Опробовать метод на примере матрицы четвертого порядка
Сравнить численно полученную обратную матрицу с точной
Оценить вычислительные ошибки численного обращения матрицы.
Задание 10
Исследовать скорость сходимости (количество итераций ) метода Зейделя при решении системы линейных уравнений второго порядка
с точностью до . В качестве начального приближения выбирать нулевое
. Рассмотреть варианты
. Графически представить зависимость скорости сходимости
от числа
.
ТЕМА 3. АППРОКСИМАЦИЯ ФУНКЦИЙ
Задание 1
Составить программу для полиномиальной интерполяции функции
по ее
узловым значениям на отрезке
. Исследовать поведение ошибки
на отрезке интерполяции при
для равномерной сетки:
и неравномерной сетки Чебышева:
Варианты полинома : а) канонический; б) Лагранжа; в) Ньютона; г) Эйткена–Невилла.
Задание 2
Составить программу для вычисления первой производной от функции I) и II)
, используя формулу дифференцирования интерполяционных многочленов: а) Лагранжа и б) канонического на равномерной сетке узлов: I)
и II)
, где
― количество узлов. Представить графически зависимость ошибки вычисления производной
от числа
.
Задание 3
Построить кубический сплайн для функции Рунге
на равномерной сетке с
узловыми значениями:
. Оценить поведение ошибки
внутри отрезка интерполяции
. Представить результаты графически.
Задание 4
Используя узловые значения функции на равномерной сетке
методом наименьших квадратов построить аппроксимирующий канонический полином
пятой степени. Вычислить среднеквадратическую ошибку полинома
при
. Сравнить коэффициенты полинома с коэффициентами ряда Тейлора для
.
Задание 5
Методом наименьших квадратов представить приближенно функцию квадратичной функцией
, используя значения
на сетке узлов
при
. Представить аппроксимирующую функцию
для случаев
и вычислить соответствующие среднеквадратические ошибки.
Задание 6
Методом наименьших квадратов представить приближенно функцию квадратичной функцией
, используя значения
на сетке узлов
. Здесь
— случайная величина, распределенная равномерно на отрезке
. Показать, как зависит точность вычисления коэффициентов аппроксимации от количества измерений
.