- •Рабочая программа учебной дисциплины
- •Объем дисциплины, виды учебной работы и формы контроля в соответствии с учебным планом специальностей по формам обучения
- •Тема 7. Численное решение краевых задач для обыкновенных дифференциальных уравнений
- •Тема 8. Обзор методов решения уравнений в частных производных
- •Раздел 2. Численные методы решения дифференциальных уравнений
- •Тестовые задания
- •Раздел 1. Методы решения задач алгебры и математического анализа
- •1. Чем вызвана неустранимая погрешность?
- •3. Пусть а - точное, а - приближенное значение некоторого числа. Дайте опре деление относительной погрешности.
- •7. Чем обусловлено появление погрешности округления при численном реше нии поставленной задачи?
- •8. Дайте определение сплайн-функции.
- •9. Сформулируйте постановку задачи интерполирования функции.
- •15. Назовите достоинства и недостатки интерполяционных формул Лагранжа.
- •16. В чем состоит сущность метода наименьших квадратов?
- •17. Когда удобно пользоваться интерполяционной схемой Эйткена?
- •19. Назовите области применения интерполирования функций.
- •С какой точностью можно вычислить по интерполяционной формуле Ла гранжа 1п 100,5 по известным значениям 1п 100, 1п 101, 1п 102 и 1п 103. А)4,5-10"5; б)6,7-10"7; в)2,3-10"9.
- •23. Опишите методику вычисления определенного интеграла по формулам прямоугольников.
- •25. Определить величину шага к по оценке остаточного члена для вычисления
- •26. Назовите области применения формул численного интегрирования.
- •29. Выбор шага интегрирования для обеспечения заданной точности вычисле ния интеграла с помощью метода двойного пересчета.
- •30. Проведите сравнение формул численного интегрирования по точности на основании остаточных членов формул.
- •32. Отличие метода Гаусса с выбором главного (ведущего) элемента от метода Гаусса решения системы линейных алгебраических уравнений.
- •33. В чем преимущество метода Зейделя для решения системы линейных ал гебраических уравнений перед методом простой итерации?
- •34. Для решения систем линейных алгебраических уравнений какого вида раз работан метод прогонки?
- •35. Опишите метод Гаусса решения системы линейных алгебраических уравнений.
- •36. Почему метод простой итерации решения систем линейных алгебраических уравнений называется самоисправляющимся?
- •37. Каковы недостатки решения системы уравнений по правилу Крамера?
- •38. Опишите метод Якоби (простой итерации) решения системы линейных ал гебраических уравнений.
- •39. Опишите метод деления отрезка пополам.
- •41. В чем достоинство и недостаток метода Ньютона нахождения корней нели нейного уравнения?
- •43. Решение нелинейного уравнения методом простой итерации.
- •44. Проведите сравнение методов деления отрезка пополам (доп) и Ньютона по различным критериям (универсальность, скорость сходимости).
- •45. Назовите основные этапы процесса нахождения корня нелинейного уравнения.
- •Раздел 2. Численные методы решения дифференциальных уравнений 1. В чем достоинство неявных методов решения дифференциальных уравнений?
- •3. Оценить погрешность аппроксимации правой разностной производной
- •4. Численное решение методом Эйлера задачи Коши для обыкновенных диф ференциальных уравнений.
- •5. Почему метод Рунге-Кутта называется самостартующим?
- •6. Опишите построение разностной схемы для численного решения обыкно венного дифференциального уравнения.
- •7. Разностная аппроксимация дифференциальных операторов.
- •8. Оценить погрешность аппроксимации центральной разностной производной
- •9. Какой метод численного решения дифференциального уравнения называет ся многошаговым?
- •10. Опишите сущность разностной аппроксимации задачи Коши для обыкно венного дифференциального уравнения первого порядка.
- •11. Оценить погрешность аппроксимации левой разностной производной
- •14. Какая конечно-разностная схема, аппроксимирующая дифференциальное уравнение в частных производных, называтся согласованной?
- •15. Какая задача для уравнений в частных производных называется корректно поставленной?
- •16. Физический смысл условия Куранта-Фридлихса-Леви.
- •17. Какая конечно-разностная схема называется слабо неустойчивой (устойчивой)?
- •18. Какие физические процессы описывают уравнения в частных производных эллиптического типа?
- •19. Укажите методы построения конечно-разностных схем, аппроксимирующих дифференциальное уравнение в частных производных.
- •20. Дайте определение маршевой задачи для уравнений в частных производных.
8. Оценить погрешность аппроксимации центральной разностной производной
У(ХМ ) - У(Х;_Л ) УМ -У 1-[ гг ~ , ,
— =^±н—^_!_9 разложив в ряд Тейлора решение дифференциальной
задачи в окрестности узла х,-.
а) О(п3), где п - шаг разностной сетки. б) О(п2), где п - шаг разностной сетки, в) О(п/3), где п - шаг разностной сетки.
9. Какой метод численного решения дифференциального уравнения называет ся многошаговым?
а) Метод, дающий формулу для вычисления значения функции ум = у(хм), являю щегося решением дифференциального уравнения, по к предыдущим значениям
уг, уг_х,..., уг_к+1, называется многошаговым, если к=\.
б) Метод, дающий формулу для вычисления значения функции у1+х =у(хм), являю щегося решением дифференциального уравнения, при известном значении у( и из вестном шаге к.
в) Метод, дающий формулу для вычисления значения функции ум = у(х1+х), являю щегося решением дифференциального уравнения, по к предыдущим значениям уг, уг_х,..., уг_к+х, называется многошаговым, если к>\.
10. Опишите сущность разностной аппроксимации задачи Коши для обыкно венного дифференциального уравнения первого порядка.
а) Исходным пунктом при разностной аппроксимации является вычисление погреш-
30
ности 8. Полагается, что погрешность 8 можно представить в виде интерполяционного полинома и-ой степени. Если разностная схема устойчива, то рост любого возмущения ограничен, следовательно, можно взять интеграл Фурье от данного полинома. В таком случае погрешность аппроксимации стремится к нулю при измельчении разностной сетки.
б) Исходным пунктом при разностной аппроксимации является замена области не прерывного изменения аргумента некоторым конечным множеством точек, лежа щих в этой области (разностная сетка) Эта дискретная модель среды описывается сеточными функциями, которые определены в узлах сетки. Дифференциальное уравнение заменяются соответствующими конечно-разностными соотношениями. В итоге исследуемая задача Коши заменяется, или, как говорят, аппроксимируется системой разностных уравнений - разностной схемой.
в) Заданный отрезок [а, Ь] заменяется системой частичных отрезков [хг, х1+\\ равной длины, называемой разностной сеткой. Расстояние между концами интервала хг-+1 - х{, = И есть единичная длина сетки. На каждом отрезке [хг, хг-+1] осуществляется чис ленное решение дифференциального уравнения. Общее решение на [а, Ь] вычисля ется как сумма частных решений по всем частичным отрезкам.
11. Оценить погрешность аппроксимации левой разностной производной
Xх;) ~ Ххг-1 ) У г ~ Уг-\ гг о , ,
— =±г—■/-±^-^ разложив в ряд Тейлора решение дифференциальной за-
дачи в окрестности узла Х{ (к - шаг разностной сетки).
а) О(п). б) О(п2). в) О(п3).
12. Применяя метод Эйлера, найти решение обыкновенного дифференциально го уравнения у'=у-2х/у на интервале [0; 1] с начальным условием у(0) = 1, выбрав шаг к = 0,2.
а)Х0,2) = 1,2000; у(0,4) = 1,4205; у(0,6) = 1,9562; ХОД = 2,3646;Х1,0) = 3,0644. б)Х0,2) = 0,9200; у(0,4) = 0,9040; у(0,6) = 0,8612; ХОД) = 0,7942;XI,0) = 0,7321. в)Х0,2) = 1,2000; ^0,4) = 1,3733; у(0,6) = 1,5294; ХОД) = 1,6786;Х1,0) = 1,8237.
13. Сформулируйте краевую задачу для обыкновенного дифференциального уравнения второго порядка у" = /(х,у,у').
а) Краевая задача для обыкновенного дифференциального уравнения второго поряд ка у" = /(х,у,у') заключается в нахождении функции /(х,у,у'), определеной во мно жестве В = {о < х < X, | у -у01 < I/} и удовлетворяющей т в области В условию Лип шица: \/(х,уъу')-/(х,у2,у')\<К\у1-у2\, где^=СОП81>0.
б) Найти функцию у = у(х), которая внутри отрезка [а, Ь] удовлетворяет уравнению у" =/(х,у,у'), а на концах отрезка - краевым условиям: ф\\у{а), у'(а)] = 0, ср2[уф),у\Ъ)\=0.
в) Краевая задача для обыкновенного дифференциального уравнения второго по рядка у" = /(х,у,у') заключается в отыскании функции у = у(х), удовлетворяющей этому уравнению и начальным условиям: Х-^о) = .Уо, У'(хо) = Уо ■
