- •Сравнение естественного и стандартного чисел обусловленности матрицы а также - точного значения стандартного числа обусловленности с его оценкой, вычисленной процедурой decomp:
- •Исследовать возможность улучшения обусловленности задачи посредством внесения малого случайного возмущения в матрицу системы.
- •7 Повторить эксперимент п.6 для 2-3 задач с плохо обусловленной матрицей.
- •8 Выполняя п.П. 6 и 7 , исследовать работоспособность различных методов оценки ошибок решения ( выражения (7), (12), (13) ) при наличии возмущения левой части системы.
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные данные:
- •Исходные данные:
- •Провести исследование влияния вида доминирования матрицы задачи на сходимость процедур Якоби и Гаусса-Зейделя.
- •Исходные данные:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Исследовать устойчивость решения задачи среднеквадратичного приближения к погрешности исходных данных.
- •Исходные параметры:
- •Убедиться в справедливости условий чебышевского интерполирования.
- •Исходные параметры:
- •При интегрировании нескольких систем нелинейных уравнений обратить внимание на эффективность различных методов реализации схемы прогноз-коррекция.
- •Исходные параметры:
- •Исходные параметры:
- •При интегрировании жёстких задач:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Исходные параметры:
- •Дополнительное исследование метода релаксации:
Исходные параметры:
Интервал:
-
№ задачи
Шаг
Ошибка численного решения
Явн. Эйлера
Неявн. Эйлера
РК2
РК3
РК4
2
1e-1
3е-1
1е-1
8е-2
2е-2
3е-3
1e-2
1е-2
1е-2
4е-4
9е-6
2е-7
13
1e-1
2е-1
2е-1
7е-3
2е-4
3е-6
1e-2
3е-3
3е-3
1е-4
2е-7
5е-10
Погрешности методов Эйлера ведут себя как , погрешности методов Рунге-Кутты - как ( - порядок метода). Поэтому методы Рунге-Кутты показали большую, чем методы Эйлера, точность. Также, можно заключить, что при неизменном порядке метода с уменьшением шага (при ) увеличивается точность решения; а при неизменном шаге (снова при ) большую точность показывает метод большего порядка.
Качественное поведение погрешности решений для методов не отличается – различие лишь в гладкости (в зависимости от шага) и величине ошибки (в зависимости от порядка метода) –
об этом свидетельствуют приложения 2а и 2б.
Полученные результаты можно объяснить схожестью вычислительных схем: методы явные и неявные Рунге-Кутты фактически представляют собой «модернизацию» явных и неявных методов Эйлера (в эксперименте исследовались лишь явные методы Рунге-Кутты). Действительно, явную схему Эйлера можно представить как , где ; аналогичное можно утверждать и для неявного метода Эйлера, заменив на . Таким образом, методы Эйлера можно условно назвать методами Рунге-Кутты первого порядка, чем и объясняются полученные результаты.
При интегрировании жёстких задач:
а) получить экспериментальные характеристики эффективности явных методов (РК4, явный Эйлера);
б) установить возможность и условия интегрирования задачи неявными методами с (метод трапеции);
в) сравнить эффективность применения метода Гира второго порядка и метода трапеции
а) Исследование эффективности явных методов.
Исходные параметры:
Задача №4 (линейная)
Интервал:
;
-
Eps
Шаг
Ошибка численного решения
Явн. Эйлера
РК4
1е-4
3e-4
Останов*
Останов
2e-4
1е+0
1е-1
1e-5
8e-3
1e-7
1е-6
8е-4
4е-11
1е-5
3e-5
Останов
Останов
2e-5
1е+0
3е-8
1е-5
3е-9
8е-12
1е-6
3е-10
3е-11
* - Неустойчивость вычислительной схемы.