- •Определение 2. Величина называется абсолютной погрешностью представления числа X с помощью числа .
- •1.2 Распространение ошибок округления в арифметических операциях.
- •1.7.2. Простейшие свойства многочленов Чебышева.
- •11. Численное интегрирование. Использование функциональных рядов.
- •12. Квадратурные формулы на основе интерполяции. Формулы для коэффициентов и остаточного члена.
- •13. Базовые квадратурные формулы прямоугольников, трапеций и парабол. Формулы теоретической погрешности.
- •14. Обобщенные квадратурные формулы трапеций и Симпсона
- •15. Теоретические оценки погрешности обобщенных формул трапеций и Симпсона.
- •17. Общие свойства полиномов ортогональных с весом.
- •18. Алгебраическая степень точности квадратурной формулы. Квадратурные формулы Гаусса-Кристоффеля.
- •19. Теорема о необходимых условиях выбора узлов в формулах Гаусса-Кристоффеля.
- •20. Классические ортогональные многочлены. Построение ортогональных многочленов Лежандра на каноническом отрезке [-1,1].
- •2. Покажем, что - неподвижная точка. Из сходимостиследует, т.Е.
- •22. Метод простых итераций решения функциональных уравнений и систем. Условия сходимости.
- •23. Метод Ньютона. Геометрическая интерпретация. Теорема о сходимости метода в одномерном случае.
- •24. Метод Ньютона в многомерном случае. Организация итерационного алгоритма.
- •25. Численные методы решения слау. Прямые и итерационные методы – общие понятия.
- •26. Нормы вещественных квадратных матриц. Спектральные свойства матриц.
- •27. Обусловленность матриц и систем лау. Определение числа обусловленности.
- •28. Метод итераций для слау второго рода. Теорема о достаточных условиях сходимости.
- •29. Спектральный признак сходимости (теорема о необх. И достат. Условиях сходимости).
- •30. Стационарные итерационные процедуры. Приведение слау первого рода к системе второго рода.
- •31. Метод простых итераций Ричардсона. Условия сходимости.
- •32. Теорема о выборе ускоряющего множителя в методе Ричардсона.
- •33. Метод Якоби. Организация алгоритма. Теорема о достаточных условиях сходимости.
- •34. Метод Зейделя как ускорение метода Якоби. Организация алгоритма. Теорема об условиях сходимости.
- •35. Метод последовательной верхней релаксации. Организация алгоритма.
- •36. Численное дифференцирование на основе интерполяции.
- •37. Численное дифференцирование на равномерной сетке, основанное на тэйлоровском разложении. Теоремы об аппроксимации первой и второй производной.
- •38. Задача Коши. Постановки задачи в одномерном случае для первой и второй производной. Сведение к системе уравнений с первой производной.
- •39. Метод Эйлера. Алгоритм, геометрическая интерпретация, порядок точности.
- •40. Методы Рунге-Кутты повышенной точности. Метод “предиктор – корректор” и метод «средней точки».
- •41. Общая постановка краевой задачи для оду 2-го порядка. Классификация граничных условий.
- •42. Метод “стрельбы” решения краевой задачи с граничными условиями первого рода.
- •43. Метод конечных разностей решения линейной краевой задачи для оду 2-го порядка.
- •45. Устойчивость разностных схем. Спектральный признак устойчивости для уравнений с постоянными коэффициентами. Примеры для оду и уравнений в частных производных.
- •16. Правило Рунге практической оценки погрешности квадратурной формулы Симпсона.
- •10. Среднеквадратичное приближение. Метод наименьших квадратов в регрессионном анализе. Вывод уравнения прямой регрессии.
41. Общая постановка краевой задачи для оду 2-го порядка. Классификация граничных условий.
Рассмотрим вначале общую нелинейную постановку краевой задачи для ОДУ второго порядка:
(18)
В постановке задачи (18) принята следующая классификация граничных условий:
при - краевые условия –однородные; при - краевоеусловие первого рода на левом конце; при - краевое условиевторого рода на левом конце; при - краевые условиятретьего рода на левом конце. На правом конце отрезка краевые условия классифицируются аналогично.
Отметим основное отличие краевой задачи (18) от задачи Коши: в задаче Коши начальные условия задаются в одной точке (как правило, на левом конце отрезка), а в краевой задаче – на обоих концах. Естественно попытаться свести краевую задачу к задаче Коши, т.к. для нее разработаны эффективные итерационные методы решения. Таковым, например, является метод стрельбы.
42. Метод “стрельбы” решения краевой задачи с граничными условиями первого рода.
Для простоты рассмотрим задачу (18) с краевыми условиями первого рода:
(19)
Заменим задачу (19) на следующую задачу Коши:
(20)
Зададим какое-либо значение параметру и решим задачу (20) подходящим методом Рунге-Кутты. Обозначим полученное решениеи вычислим разность- погрешность решения на правом конце. Корректируем угловой коэффициентна левом конце в зависимости от знака погрешности, заменяя его на. Находим новое решениеи т.д. до тех пор, пока не выполнится условие, где- заданная погрешность. Название метода связано с его геометрической интерпретацией:стрельба из точки с координатами под углом, определяемым угловым коэффициентом.
43. Метод конечных разностей решения линейной краевой задачи для оду 2-го порядка.
Наиболее эффективно этот метод применяется для линейной краевой задачи, постановка которой для дифференциального уравнения второго порядка имеет вид:
|
(21) |
Введем равномерную сетку - шаг сетки.
Определим на этой сетке сеточные функции:
, где - точное решение на сетке.
Запишем основное уравнение системы (21) для фиксированного узла :
. (22)
В каждом внутреннем узле сетки () аппроксимируем производные по формулам второго порядка точности (см. теор. 4.1 – 4.3):
Подставим полученные соотношения в (22):
,
где обозначено -сеточный оператор, действующий на сетке. Умножая обе части уравнения на и отбрасывая остаточный член, получим приближение к точному уравнению:
.
Приводя подобные члены, получим окончательно:
, |
(23) |
Общее число неизвестных на сетке равно (), а система (23) содержитуравнение для неизвестных. Недостающие значения должны быть определены из граничных условий. Порядок аппроксимации основного уравнения системы. При постановке задачи с краевыми условиями первого рода порядок аппроксимации всей схемы не ухудшится и будет равен. В случае краевых условий второго или третьего рода необходимо аппроксимировать их также со вторым порядком, чтобы не потерять общий второй порядок аппроксимации всей схемы. Рассмотрим один из вариантов аппроксимации краевого условия третьего рода на левом конце отрезка.
Граничное условие третьего рода имеет вид (см. урав.(18)):
(24)
Разложим функцию по Тэйлору с центром в точкеи найдем значение:
найдем из уравнения (21), аппроксимируя его в точке сетки (непрерывность вторых производных на отрезкепозволяет это сделать):
.
Подставляя найденное значение второй производной в выражение для , получаем:
.
Отбрасывая остаточный член и приводя подобные в последнем уравнении, получаем систему линейных уравнений для неизвестных :
Исключая из полученной системы , получаем линейное уравнение, связывающее два сеточных значенияи:
,
где параметры легко выражаются через. Аналогично аппроксимируется третье краевое условие на правом конце отрезка, которое приводится к виду
.
В результате получаем неявную конечно-разностную систему уравнений, записанную в каноническом виде:
(25)
Система (25) - трехдиагональная, для которой разработаны специальные эффективные методы численного решения, например, метод прогонки [1,6,11].
Исследуем вопрос о сходимости решения системы (25). Одно из условий метода Якоби (и соответственно, методов Зейделя и верхней релаксации): - очевидно выполняется. Следующее условие требует «диагонального преобладания», т.е.. Подставляя сюда значения коэффициентов системы (23), получаем:
Полученное неравенство выполняется, если одновременно выполняются условия
Отсюда следует, что должны выполняться условия
.
Усиливая последнее неравенство, получаем условие на шаг сетки:
.