Федеральное агентство железнодорожного транспорта
Омский государственный университет путей сообщения
Кафедра «Автоматика и системы управления»
РЕШЕНИЕ УРАВНЕНИЙ И ИНТЕРПОЛЯЦИЯ ФУНКЦИЙ
Пояснительная записка к курсовому проекту
по дисциплине «Численные методы»
ИНМВ.3000000.000 ПЗ
|
|
Студентка гр. ИС 20253
_________ А.С.Хананина «__»________2014 г.
Руководитель –доцент кафедры АиСУ ____________А. C. Окишев «__»________2014 г. |
Омск 2014
УДК 519.6
Пояснительная записка содержит 30 страниц, 18 рисунков, 5 таблиц, 3 источника.
Численный метод, нелинейное уравнение, корень, итерация, сходимость, аппроксимация, интерполяция, многочлен Лагранжа, задача Коши, погрешность, кубический сплайн, обыкновенное дифференциальное уравнение.
Объектом исследования являются приближенные численные методы решения некоторых математических и инженерных задач, а также программное обеспечение, реализующее эти методы.
Цель работы – ознакомиться с численными методами решения нелинейных и дифференциальных уравнений и интерполяции функций, решить предложенные типовые задачи с помощью предоставленного преподавателем программного обеспечения, сформулировать выводы по полученным решениям, отметить достоинства и недостатки методов, сравнить удобство использования и эффективность работы каждой программы.
Пояснительная записка к курсовому проекту оформлена в текстовом редакторе MicrosoftOffice 2010 с установленным интерактивным редактором формул MathType 6.6. Графики нелинейных функций построены с помощью программы AdvancedGrapher. При решении обыкновенных дифференциальных уравнений использовалась среда математического моделирования Mathcad 14.
Содержание
1 Решение нелинейных уравнений. 5
1.3 Решение нелинейного уравнения методом простых итераций. 9
1.4 Решение нелинейного уравнения методом Ньютона. 12
2.1 Интерполяционная формула Лагранжа. 16
2.2 Интерполяционные формулы Ньютона. 16
2.3 Интерполяция кубическими сплайнами. 17
2.4 Интерполяция функций сплайнами по таблице значений. 19
3 Решение обыкновенных дифференциальных уравнений. 21
3.2 Решение ОДУ методом Эйлера. 23
3.3 Решение ОДУ методом Эйлера-Коши. 25
Заключение. 28
В связи с развитием новой вычислительной техники инженерная практика наших дней все чаще и чаще встречается с математическими задачами, точное решение которых получить весьма сложно или невозможно. В этих случаях обычно прибегают к тем или иным приближенным вычислениям. Вот почему приближенные и численные методы математического анализа получили за последние годы широкое развитие и приобрели исключительно важное значение.
Новые вычислительные средства вызвали переоценку известных методов решения задач с точки зрения целесообразности их реализации на ЭВМ и стимулировали создание более эффективных, что привело к появлению новой дисциплины – вычислительной математики. Предметом изучения последней являются численные методы решения задач математического анализа: изучение алгоритмов и условий сходимости итерационных методов, определение границ применимости методов, исследования оценок погрешностей методов и вычислений. Главным разделом вычислительной математики является реализация численных методов на ЭВМ, то есть составление программы для требуемого алгоритма и решения с ее помощью конкретной задачи.
Любая прикладная задача формируется исходя из определенного физического смысла некоторого процесса (распределение тепла в стержне, описание траектории движения объектов). Прикладная математическая задача может быть сформулирована, например, из описания некоторой экономической модели (задача распределения ресурсов, задача планирования производства, транспортная задача перевозки грузов, оптимальных в заданном смысле). Следовательно, для постановки любой прикладной задачи нужна математическая модель. Поэтому, можно выделить следующие этапы решения задач на ЭВМ:
1) описание математической модели задачи на основе физической или экономической модели;
2) изучение методов решения поставленной математической модели задачи и создание новых методов;
3) выбор метода решения задачи исходя из заданной точности решения и особенностей задачи;
4) составление блок-схемы программы для решения задачи на ЭВМ;
5) отладка программы и оценка полученных результатов;
6) решение задачи на ЭВМ, построение графиков, получение оценки погрешностей, обоснование результатов.
В курсовом проекте рассматриваются не прикладные, а типовые математические задачи, которые могут возникнуть при переходе от реальных систем к их математическим моделям, поэтому основное внимание уделяется последнему этапу.
1 Решение нелинейных уравнений
Нелинейными уравнениями называются уравнения вида
,
где – нелинейная функция, которая может относиться к трем типам:
1) нелинейная алгебраическая функция вида
;
2) трансцендентные функции – тригонометрические, обратные тригонометрические, логарифмические, показательные и гиперболические функции;
3) различные комбинации этих функций, например, .
Решением нелинейного уравнения является такая точка , которая при подстановке в уравнение обращает его в тождество. На практике не всегда удается подобрать такое решение. В этом случае решение уравнения находят с применением приближенных (численных) методов. Тогда решением будет являться такая точка , при подстановке которой в уравнение последнее будет выполняться с определенной степенью точности, т.е. , где e – малая величина. Нахождение таких решений и составляет основу численных методов и вычислительной математики.
Решение нелинейных уравнений разделяется на два этапа: отделение корней уравнений и уточнение корней нелинейных уравнений.
На первом этапе необходимо исследовать уравнение и выяснить, имеются корни или нет. Если корни имеются, то узнать, сколько их, и затем определить интервалы, в каждом из которых находится единственный корень.
Первый способ отделения корней – графический. Исходя из уравнения , можно построить график функции . Тогда точка пересечения графика с осью абсцисс является приближенным значением корня. Если f(x) имеет сложный вид, то ее можно представить в виде разности двух функций . Так как , то выполняется равенство .Если построить два графика , , то абсцисса точки их пересечения будет приближенным значением корня уравнения.
Второй способ отделения корней нелинейных уравнений – аналитический. Он основывается на следующих трех теоремах.
Теорема 1. Если функция непрерывна на отрезке и меняет на концах отрезка знак (т.е. ), то на содержится хотя бы один корень.
Теорема 2. Если функция непрерывна на отрезке , выполняется условие вида и производная сохраняет знак на , то на отрезке имеется единственный корень.
Теорема 3. Если функция является многочленом n-ой степени и на концах отрезка меняет знак, то на имеется нечетное количество корней (если производная сохраняет знак на , то корень единственный). Если на концах отрезка функция не меняет знак, то уравнение либо не имеет корней на , либо имеет четное количество корней.
При аналитическом методе исследований необходимо выявить интервалы монотонности функции . Для этого необходимо вычислить критические точки , в которых первая производная равна нулю или не существует. Тогда вся числовая ось разбивается на интервалы монотонности . На каждом из них определяется знак производной , где . Затем выделяются те интервалы монотонности, на которых функция меняет знак.
На втором этапе на каждом из этих интервалов для поиска корня используются численные итерационные методы уточнения корней, например методы половинного деления, простых итераций или Ньютона.
1.1 Метод простых итераций
Пусть известно, что нелинейное уравнение имеет на отрезке единственный вещественный корень . Требуется найти этот корень с заданной точностью. Применяя тождественные преобразования, приведем уравнение к виду
.
Выберем произвольно приближенное значение корня и вычислим . Найденное значение подставим в правую часть соотношения (2) и вычислим . Продолжая процесс вычислений дальше, получим числовую последовательность . Если существует предел этой последовательности, то он и является корнем уравнения . В самом деле, пусть . Тогда, переходя к пределу в равенстве
и учитывая непрерывность функции на отрезке , получим
или .
Корень можно вычислить с заданной точностью по итерационной формуле
Достаточное условие, при котором итерационный процесс сходится, определяет следующая теорема: пусть функция определена и дифференцируема на отрезке , причем все ее значения и выполняется условие
,
тогда процесс итераций сходится независимо от начального значения и предельное значение является единственным корнем уравнения на .
Геометрическая интерпретация метода простых итераций заключается в следующем: если построить два графика функций и , то абсцисса точки их пересечения будет корнем . Построим итерационный процесс. Зададим .Вычислим – первое приближение и – второе приближение. В первом случае (рисунок 1, а) процесс сходящийся ( ), во втором (рисунок 1, б) – расходящийся ( ).
а б
Рисунок 1 – Сходящийся (а) и расходящийся (б) итерационные процессы
Часто, если итерационный процесс расходится из-за невыполнения условия , нелинейное уравнение можно привести к виду, допускающему сходящиеся итерации.
Выполнения условия сходимости можно добиться путем перехода от исходного уравнения к эквивалентному виду следующим образом: сначала умножить обе части уравнения (1) на , а затем прибавить к обеим частямx, тогда . Обозначив ,получим уравнение . Константа с выбирается так, чтобы выполнялось достаточное условие сходимости итерационного процесса , т.е.
Условие равносильно двойному неравенству
,
Поэтому константа выбирается из соотношений:
Метод простых итераций и почти все другие итерационные методы имеют два достоинства:
1) являются универсальными и самоисправляющимися, то есть любая неточность на каком-либо шаге итераций отражается не на конечном результате, а лишь на количестве итераций. Подобные ошибки устойчивы даже по отношению к грубым ошибкам (сбоям ЭВМ), если только ошибка не выбрасывает очередное приближение за пределы области сходимости;
2) позволяют достигнуть любой заданной точности при любом начальном приближении .
Недостатки методов:
– трудность приведения уравнения к виду .
– если начальное приближение находится далеко от корня, то число итераций при этом увеличивается, а объем вычислений возрастает.
Процесс итераций заканчивается при выполнении двух критериев:
1) Когда два последних приближения отличается между собой по модулю на заданную величинуe:
.
Одного критерия недостаточно, так как в случае крутизны графика, данное условие будет выполнено, но может находиться далеко от корня;
2) Когда последнее вычисленное приближение к корню удовлетворяет уравнению с заданной точностью:
.
Отдельно критерия бывает недостаточно, так как при пологой функции условие может быть выполнено, но может быть далеко от корня.
1.2 Метод Ньютона
Пусть уравнение имеет на интервале единственный корень, причем существует непрерывная на производная . Метод Ньютона служит для уточнения корней нелинейных уравнений в заданном интервале. Его можно рассматривать как частный случай метода простых итераций, если
.
Тогда итерационный процесс осуществляется по формуле:
.
Геометрически этот процесс представлен на рисунке 2. Он означает замену на каждой итерацииk графика кривой касательной к ней в точках с координатами .
Достаточное условие сходимости обеспечивается выбором начальной точки . Начальным приближением служит один из концов отрезка , в зависимости от того, в каком из них выполняется достаточное условие сходимости
.
При произвольном начальном приближении итерации сходятся, если
.
Метод Ньютона рекомендуется применять для нахождения простых действительных корней уравнения .
Достоинством метода является то, что он обладает скоростью сходимости, близкой к квадратичной.
Недостатки метода:
– не при любом начальном приближении метод Ньютона сходится, а лишь при таком, для которого ;
– если , т.е. касательная к графику почти параллельна оси абсцисс, то и метод расходится;
– если ,т.е. касательная к графику почти параллельна оси ординат, то и продвижения к корню не будет.
Последних трудностей можно избежать, применив модификацию метода Ньютона, в которой используется только касательная в точке начального приближения. Рабочая формула при этом имеет вид:
.
1.3 Решение нелинейного уравнения методом простых итераций
Решите нелинейное уравнение
x3/ 7 + 5 = 0
методом простой итерации:
- a) выполните приблизительную оценку корней предложенного уравнения с помощью системы MathCad, графической программы AdvancedGrapher или вручную;
б) преобразуйте уравнение вида f (x) = 0 к итерационному виду x = φ(x), гарантировав при этом сходимость метода простой итерации;
в) с помощью программы MeProst уточните один из корней первого уравнения вида f (x)=0 с точностью ε = 10-3, выбрав подходящий отрезок для построения графической иллюстрации работы метода простой итерации;
г) в пояснительной записке приведите график функции f (x) с указанием выбранных для уточнения корней, результаты уточнения корней программой MeProst, значения корней, требуемое число итераций, погрешность решения и графические иллюстрации процесса сходимости к корню.
Для приблизительной оценки корней нелинейного уравнения построим график функции с помощью программы AdvancedGrapher (рисунок 3).
Из рисунка 3 видно, что график функции пересекает ось абсцисс в одной точке, поэтому уравнение имеет один корень.
Преобразуем уравнение
к виду x+(x^3/7+5)/-20
т.е. для уточнения корня будем использовать итерационную последовательность:
где – номер итерации.
С помощью программы MeProst выполним уточнение корня. В разделе меню «Ввод» задаем следующие параметры (Рисунок 4):
Рисунок 4 –Ввод параметров в программе Newton
- отрезок, на котором находится уточняемый корень ;
- необходимая точность приближения к корню ;
- начальное приближение к корню
- максимальное допустимое число итераций .
Чтобы получить наглядное представление о процессе сходимости к корню, выберем пошаговый режим расчета. В открывшемся окне (рисунок 5) видим графики функций (зеленый) и (красный). Очевидно, что корень уравнения определяется как точка пересечения этих двух графиков. Точка начального приближения показана синим цветом. Чтобы вычислить следующее приближение геометрическим способом, необходимо найти точку пересечения прямой и кривой и ординату этой точки взять в качестве . Последующие приближения определяются аналогично.
Рисунок 5 – Метод простых итераций
Как следует из рисунка 5, итерационный процесс сходится (с каждый новым шагом мы приближаемся к решению).
Корнем данного уравнения является
1.4 Решение нелинейного уравнения методом Ньютона
Решите нелинейное уравнение
exp(– x) – 5,5 + x = 0
методом Ньютона:
- a) выполните приблизительную оценку корней уравнения с помощью системы MathCad, графической программы AdvancedGrapher или вручную;
б) с помощью программы Newton уточните один из корней первого уравнения вида f (x)=0 с точностью ε = 10-3, выбрав подходящий отрезок для построения графической иллюстрации работы метода Ньютона;
в) в пояснительной записке приведите график функции f(x) с указанием выбранных для уточнения корней, результаты уточнения корней программой Newton, значения корней, требуемое число итераций, погрешность решения и графические иллюстрации процесса сходимости к корню.
Для приблизительной оценки корней нелинейного уравнения построим график функции с помощью программы AdvancedGrapher (рисунок 6).
Как видно из рисунка 6, график функции пересекает ось абсцисс в двух точках, поэтому первый корень x1= -2, второй корень x2= 5,5.
С помощью программы Newton выполним уточнение выбранного корня. Для этого в разделе меню «Ввод» задаем следующие параметры (рисунок 7):
- отрезок построения касательных к графику функции ;
- необходимая точность приближения к корню ;
- начальное приближение к корню ;
- максимальное допустимое число итераций .
Рисунок 7 – Ввод параметров в программе Newton
Рисунок 8 – Четвертая итерация метода Ньютона
Как видно из рисунка 8, корень уравнения 5 найден с заданной точностью за итераций, погрешность решения составляет .
Задача интерполяции функций возникает в тех случаях, когда:
функция задана в виде таблицы, и необходимо знать значения функции для промежуточных значений аргументов, расположенных в таблице между узлами , а также для аргументов, расположенных вне таблицы;
известна лишь таблица функции и требуется определить ее аналитическое выражение;
известно аналитическое выражение функции, но оно имеет очень сложный вид, вследствие чего возникает необходимость представления этой функции в более простом виде. Например, при вычислении определенных интегралов вида можно заменить подынтегральную функцию некоторой приближенной функцией в виде многочлена. Тогда
.
Построив интерполяционный многочлен любого вида, также можно расширить таблицу как влево, так и вправо, вычисляя построенный многочлен в точках, не принадлежащих таблице (задача экстраполяции). Кроме того, построив интерполяционный многочлен, можно уплотнить таблицу, определяя значения функции для промежуточных аргументов между узловыми точками.
2.1 Интерполяционная формула Лагранжа
Пусть задана система точек , в которых известны значения функции (см. таблицу 1).
Таблица 1 – Экспериментальные значения функции
|
Номер точки |
0 |
1 |
… |
n |
|
… |
||||
|
… |
Установим зависимость одного ряда чисел от другого и построим новую функцию, которая с определенной степенью точности будет приближена к заданной.
Рассмотрим пример построения интерполяционного многочлена Лагранжа по заданной системе точек (в общем случае для неравноотстоящих аргументов). Построим некоторый многочлен таким образом, чтобы его значения совпали со значениями функции, заданными в таблице, для тех же аргументов, то есть . Лагранж предложил строить многочлен степени nв виде:
Здесь в каждом слагаемом коэффициенту ai соответствует разность .Найдем неизвестные коэффициенты , называемые коэффициентами Лагранжа, используя условие :
При :
;
.
Следовательно, коэффициент a0 вычисляется по следующей формуле:
.
При :
;
.
Следовательно, коэффициент a1 вычисляется по формуле:
.
Значения остальных коэффициентов вычисляются аналогично.
С учетом найденных коэффициентов интерполяционный многочлен Лагранжа запишется в виде
.
Погрешность формулы определяется остаточным членом:
,
где x–точка наименьшего промежутка, содержащего все узлы и точкуx.
2.2 Интерполяционные формулы Ньютона
Если таблица, для которой построена формула Лагранжа, задана для равноотстоящих узлов , то формула Лагранжа упрощается.
Обозначим .Тогда
;
.
С учетом введенных обозначений формула Лагранжа запишется в виде:
Запишем формулу Лагранжа в случае, если :
Получили формулу линейной интерполяции вида:
,
где – табличные разности первого порядка.
При получаем формулу квадратичной интерполяции вида:
,
где – табличные разности второго порядка.
Таким образом, первая интерполяционная формула Ньютона будет иметь вид:
.
Остаточный член формулы имеет вид:
,
где ,x – точка наименьшего промежутка, содержащего все узлы и точкуx.
Первая интерполяционная формула рекомендуется для интерполяции в начале таблицы.
Вторая интерполяционная формула Ньютона используется при вычислении функции для значений , близких к концу таблицы и для продолжения таблицы. Обозначим через , тогда
.
Тогда получим следующую формулу Ньютона:
.
Погрешность формулы
имеет тот же смысл, что и в первой формуле Ньютона.
2.3 Интерполяция кубическими сплайнами
Кубическим сплайном называется функция , которая:
- на каждом отрезке является многочленом третьей степени;
- имеет непрерывные первую и вторую производные на всем отрезке [a, b];
- в точках xi выполняются равенства и ;
- имеет граничные условия вида .
По количеству точек, учитываемых при аппроксимации функции, кубические сплайны можно разделить на локальные и глобальные.
Локальные сплайны используют производные , вычисляемые по трем заданным точкам, ближайшим к точке x, в которой необходимо вычислить значение функции.
Глобальные сплайны используют массив вторых производных по всем узловым точкам , который рассчитывается заранее путем решения системы линейных алгебраических уравнений (СЛАУ).
Пусть имеется таблица значений , где i изменяется от 0 до n.
Получается n+1 узел, причем сетка узлов неравномерная, тогда на интервале для нахождения значения сплайна в некоторой введенной точке необходимо проверить, принадлежит ли данная точка интервалу значений исходной таблицы.
- Если принадлежит, то формула вычисления значения сплайна в точке имеет вид:
где – расстояние между соседними узлами;
– значения функции в узловых точках;
– значения вторых производных.
Вторые производные mi находятся из трехдиагональной СЛАУ:
.
Всего в системе n–1 уравнение и n+1 неизвестная: . Система не полностью определяет mi. Сводим ее к трехдиагональной СЛАУ заданием граничных условий. Для нормального сплайна , . Теперь в СЛАУ n–1неизвестная иона является трехдиагональной. Такая система решается методом прогонки.
- Асимптотическое поведение сплайна вне интервала описывается формулами:
- линейная экстраполяция при (влево):
;
- линейная экстраполяция при (вправо):
.
2.4 Интерполяция функций сплайнами по таблице значений
Постройте график кубического сплайна и вычислите значения неизвестной функции f(x) в промежуточных точках.
Для каждого из двух предложенных наборов экспериментальных данных (см. таблицу 2):
а) с помощью программы Splacc постройте график кубического нормального сплайна и приведите его в пояснительной записке;
б) с помощью программы Splacc вычислите значения неизвестной функции f(x) во всех предложенных промежуточных точках;
в) вычислите те же значения неизвестной функции f(x) в промежуточных точках, введя их единым массивом в программе Spladd;
г) сравните программы Splacc и Spladd, укажите достоинства и недостатки каждой из них, а также отличия в полученных числовых значениях, если таковые имеются.
Таблица 2 – Экспериментальные данные для кубической интерполяции
|
10 |
Экспериментальные точки |
||||||||||
|
1 |
x |
-6,5 |
-5,5 |
-4,5 |
-2,5 |
-1 |
0 |
2 |
4,5 |
5,5 |
8,5 |
|
y |
-6 |
-1 |
1,5 |
3,5 |
2,5 |
-1,5 |
-3,75 |
-1,5 |
2 |
8,5 |
|
|
2 |
x |
-9,5 |
-6,5 |
-4 |
-2,5 |
-0,5 |
1,5 |
3 |
4,5 |
9,5 |
|
|
y |
5,5 |
-1 |
-4,5 |
-2 |
-5 |
-2,5 |
0 |
1,5 |
-1,5 |
|
|
Графики сплайнов для первой и второй функции, построенные с помощью программы Splacc, представлены на рисунке10.
Рисунок 10 – Интерполяция сплайном для первой (а) и второй (б) функции в программе Splacc
Значения сплайна в расчетных точках приведены в таблицах 3,4. Расчеты значений сплайнов программами Splacc и Spladd дали одинаковые результаты.
Таблица 3 – Значения сплайна в расчетных точках в программе Splacc
|
10 |
Точки вычисления f(x) |
||||||||
|
1 |
x |
-6 |
-5 |
-3,5 |
-0,5 |
1 |
3 |
7 |
10 |
|
y |
-3.28 |
0.524 |
2.76 |
0.56 |
-3.5 |
-3.65 |
5.99 |
10.77 |
|
|
2 |
x |
-8 |
-5 |
-3 |
-1,5 |
0,5 |
4 |
8 |
|
|
y |
2.43 |
-4.28 |
-2.59 |
-3.23 |
-4.44 |
1.14 |
0.19 |
|
|
Рисунок 11 – Значения сплайна в расчетных точках в программе Spladd для первой (А) и второй функции (Б)
3 Решение обыкновенных дифференциальных уравнений
Дифференциальные уравнения являются основным математическим инструментом моделирования и анализа разнообразных явлений и процессов в науке и технике.
Методы их решения подразделяются на два класса:
аналитические методы, в которых решение получается в виде аналитических функций;
численные (приближенные) методы, где искомые интегральные кривые получают в виде таблиц их численных значений.
Применение аналитических методов позволяет исследовать полученные решения методами математического анализа и сделать соответствующие выводы о свойствах моделируемого явления или процесса. К сожалению, с помощью таких методов можно решать достаточно ограниченный круг реальных задач. Численные методы позволяют получить с определенной точностью приближенное решение практически любой задачи.
Решить дифференциальное уравнение
численным методом означает, что для заданной последовательности аргументов и числа , не определяя аналитического вида функции , найти значения , удовлетворяющие начальным условиям:
.
Этот метод является сравнительно грубым и применяется в основном для ориентировочных расчетов. Однако идеи, положенные в основу метода Эйлера, являются исходными для ряда других численных методов.
Пусть дано обыкновенное дифференциальное уравнение (ОДУ) с начальными условиями (задача Коши):
и удовлетворяются условия существования и единственности решения.
Требуется найти решение задачи Коши на отрезке . Найдем решение в виде таблицы . Для этого разобьем отрезок наn равных частей и построим последовательность
,
где – шаг интегрирования.
Проинтегрируем исходное уравнение на отрезке :
Полученное соотношение можно переписать как
.
Если считать подынтегральную функцию постоянной на участке и равной значению в начальной точке этого интервала , то получим
.
Подставляя полученный результат в формулу , получим основную расчетную формулу метода Эйлера:
.
Вычисление значений осуществляется с использованием формулы следующим образом. По заданным начальным условиям иy0, полагая в выражении , вычисляется значение
.
Далее, определяя значение аргумента x по формуле , используя найденное значение y1 и полагая в формуле , вычисляем следующее приближенное значение интегральной кривой как
.
Поступая аналогичным образом при , определяем все остальные значенияyk, в том числе последнее значение
,
которое соответствует значению аргумента .
Таким образом, соединяя на координатной плоскости точки отрезками прямых в качестве приближенного представления искомой интегральной кривой , получаем ломанную линию с вершинами в точках .
Метод Эйлера может быть применен к решению систем дифференциальных уравнений.
Пусть задана система двух уравнений первого порядка:
с начальными условиями
.
Необходимо найти решение этой задачи Коши. Проводя аналогичные рассуждения, получаем расчетные формулы вида:
гдеh– шаг интегрирования.
При расчетах полагается, что и . В результате применения расчетной схемы получается приближенное представление интегральных кривых и в форме двух ломанных Эйлера, построенных по полученным таблицам . Точность метода Эйлера .
3.2 Решение ОДУ методом Эйлера
Решите ОДУ первого порядка (задачу Коши)
|
h=0,2 |
методом Эйлера:
а) изучите устройство и работу шаблона МЭ в среде математического моделирования MathCAD и получите методом Эйлера таблицу значений неизвестной функции y(x);
б) аппроксимируйте полученные данные кубическим сплайном и постройте приближенный график функции y(x) в программе Splacc.
Вводим исходные данные в шаблон МЭ.mcd и по формуле определяем значение y1(листинг 1).
Листинг 1 – Первый шаг метода Эйлера
Подставляя найденное значение y1 в формулу, вычисляем следующее значение функции y2 (листинг 2).
Листинг 2 – Второй шаг метода Эйлера
Остальные значения функции вычисляются аналогично до тех пор, пока не будет достигнут правый конец отрезка интегрирования (листинг 3).
Листинг 3 – Последний шаг метода Эйлера
Полученные значения неизвестной функции представлены в таблице 4.
Таблица 4 – Значения функции
|
k |
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
xk |
-1.5 |
-1,3 |
-1,1 |
-0,9 |
-0,7 |
-0,5 |
-0,3 |
-0,1 |
0,1 |
0,3 |
0,5 |
|
yk |
-2,518 |
-3,214 |
-3,715 |
-4,068 |
-4,309 |
-4,463 |
-4,548 |
-4,578 |
-4,566 |
-4,519 |
-4,443 |
Аппроксимируем таблично заданную функцию кубическим сплайном с помощью программы Splacc. График сплайна показан на рисунке 15. Как видно из рисунка, решение ОДУ практически линейно.
Рисунок 15 – Интерполяция численного решения ОДУ кубическим сплайном
3.3 Решение ОДУ методом Эйлера-Коши
Решите ОДУ первого порядка (задачу Коши)
|
h=0,3 |
методом Эйлера-Коши:
а) изучите устройство и работу шаблона МЭК в среде математического моделирования MathCAD и получите методом Эйлера-Коши таблицу значений неизвестной функции y(x);
б) аппроксимируйте полученные данные многочленом Лагранжа и постройте приближенный график функции y(x) в программе Lagr;
Вводим исходные данные в шаблон МЭК.mcd и определяем значение y1 (листинг 4).
Листинг 4 – Первый шаг метода Эйлера-Коши
Остальные значения функции вычисляются аналогично по шаблону МЭК до тех пор, пока не будет достигнут правый конец отрезка интегрирования (листинг 5).
Листинг 5 – Последний шаг метода Эйлера-Коши
Полученные значения неизвестной функции представлены в таблице 5.
Таблица 5 – Значения функции
|
k |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
xk |
0,3 |
0,6 |
0,9 |
1,2 |
1,5 |
1,8 |
2,1 |
2,4 |
2,7 |
3,0 |
|
yk |
2,079 |
1,269 |
0,83 |
0,806 |
1,212 |
2,013 |
2,465 |
1,872 |
1,141 |
0,75 |
Аппроксимируем таблично заданную функцию многочленом Лагранжа с помощью программы Lagr. График многочлена приведен на рисунке 18. Как следует из рисунка, в данном случае решение ОДУ имеет нелинейный характер.
Рисунок 18 – Интерполяция численного решения ОДУ многочленом Лагранжа
Недостаток аналитических методов – использование целого ряда допущений и предположений в процессе построения математических моделей и невозможность, в некоторых случаях, получить решение в явном виде из-за неразрешимости уравнений в аналитической форме, отсутствия первообразных для подынтегральных функций и т.п. В этих случаях широко применяются численные методы.
Численные методы основываются на построении конечной последовательности действий над числами. Применение численных методов сводится к замене математических операций и отношений соответствующими операциями над числами, например, к замене интегралов суммами, бесконечных сумм – конечными и т.п. Результатом применения численных методов являются таблицы и графики зависимостей, раскрывающих свойства объекта. Численные методы являются продолжением аналитических методов в тех случаях, когда результат не может быть получен в явном виде. Численные методы по сравнению с аналитическими методами позволяют решать значительно более широкий круг задач.
Предложенные преподавателем программы показали себя эффективно для решения конкретных задач, под которые они были написаны. Но их использование для решения различных инженерных задач нецелесообразно из-за крайне ограниченного функционала. Возможный выход – это самостоятельное написание программ под конкретные задачи, но это чревато большими трудозатрами. Альтернатива – использование коммерческих пакетов (например, системы математического моделирования Mathcad, Matlab, Maple и т.д.), позволяющих решать очень широкий круг задач. Использование различных скриптов к этим коммерческим пакетам также позволяет еще сильнее расширить их функционал.
Примеров инженерных задач, в которых могут потребоваться изученные численные методы, можно привести целое множество. Прежде всего это те направления, где требуется математическое моделирование различных процессов. Например, моделирование уличной дорожной сети (для оптимизации использования общественного и частного транспорта), моделирование использования различных инженерных сооружений для расчета эффективности их использования и расчета оптимальных, а также максимальных нагрузок на них, моделирование взаимодействия инженерных сооружений и окружающей среды для расчета возможных экологических последствий их использования.
Библиографический список
- ВержбицкийВ.М. Численные методы. Линейная алгебра и нелинейные уравнения. М.: Оникс 21 век, 2005. – 636 с.
- ВержбицкийВ.М. Численные методы. Математический анализ и обыкновенные дифференциальные уравнения. М.: Оникс 21 век, 2005. – 400 с.
- Бахвалов Н.С. Численные методы/ Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков – М.: Бином. Лаборатория знаний, 2008. – 432 с.
Скачать: