Построение математической модели, описывающей процесс решения дифференциального уравнения
Основные правила расчета значений дифференциального уравнения. Изучение выполнения оценки погрешности вычислений, осуществления аппроксимации решений. Разработка алгоритма и написание соответствующей программы. Построение интерполяционного многочлена.
| Рубрика | Математика |
| Вид | курсовая работа |
| Язык | русский |
| Дата добавления | 11.12.2013 |
| Размер файла | 212,6 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Лист задания
Начальные условия: , , , , ,
Постановка задачи:
1. Рассчитать погрешности уточненных значений
1.1 Рассчитать уравнения с шагом h
1.2 Рассчитать уравнения с шагом h/2
1.3 Оценить погрешности вычислений при решении задачи
1.4 Рассчитать уточненные решения yут.
1.5 Составить таблицу данных 1
1.6 Построить график 1 - значений yh, yh/2, yут
2. Аппроксимирование квадратного трехчлена методом наименьших квадратов
2.1 Составить таблицу 2 - рассчитанных значений для расчета коэффициентов
2.2 Составить систему уравнений
2.3 Решить систему уравнений методом Гаусса
2.4 Составить таблицу 3 - данных для расчета погрешности аппроксимации
2.5 Построить график 2 - значений yh и F(x).
3. Интерполирование
3.1 Построить интерполяционный многочлен Лагранжа и сгустить значения
3.2 Рассчитать погрешность интерполяции
3.3 Составить таблицу данных 4
3.4 Построить график 3 - значений y(x), F(x), P(x)
4. Проанализировать полученные результаты
5. Составить программу для проверки правильности расчетов
Содержание
Введение
1. Расчет погрешностей и уточненных значений
1.1 Расчет уравнений с шагом h
1.2 Расчет уравнений с шагом h/2
1.3 Оценка погрешности вычислений при решении задачи
1.4 Расчет уточненных решений yут
1.5 Таблица данных 1
1.6 График 1 - значений yh, yh/2, yут
2. Аппроксимирование квадратного трехчлена методом наименьших квадратов
2.1 Таблица 2 - рассчитанных значений для расчета коэффициентов
2.2 Составление системы уравнений
2.3 Решение системы уравнений методом Гаусса
2.4 Таблица 3 - данных для расчета погрешности аппроксимации "оапп"
2.5 График 2 - значений yh и F(x)
3. Интерполяция
3.1 Построение интерполяционного многочлена Лагранжа
3.2 Расчет погрешности интерполяции
3.3 Таблица данных 4
3.4 График 3 - значений y(х), F(X) и P(x)
4. Анализ полученных результатов
Заключение
Введение
Настоящее время характеризуется резким расширением математики, что связано с созданием и развитием средств вычислительной техники. В результате появления ЭВМ с программным управлением, менее чем за 50 лет скорость выполнения операций возросла от 0,1 операции в секунду при ручном счете до 1012 операций на современных серийных ЭВМ.
Мнение о всемогуществе современных ЭВМ порождает впечатление, что разработка численных методов не столь важна. В действительности же, расширение возможностей приложения математики обусловило математизацию химии, экономики, биологии, географии, геологии и д.р. суть математизации состоит в построении новых математических моделей явлений и процессов, а также разработке методов их исследования.
1. Расчет погрешностей и уточненных значений
1.1 Расчет уравнений с шагом h
Начальные условия:
, , , ,
метод Эйлера-Коши, аппроксимирование методом квадратного трехчлена, составление многочлена Лагранжа.
Необходимо вычислить табличные значения решений квадратного уравнения, используя формулы:
; ; ;
1. х1 = 0,1
2. х2 = 0,2
3. х3 = 0,3
4. х4 = 0,4
5. х5 = 0,5
6. х6= 0,6
7. х7 = 0,7
8. х8 = 0,8
9. х9 = 0,9
1.2 Расчет уравнений с шагом h/2
1. х1 = 0,1
2. х2 = 0,2
5. х5 = 0,5
9. х9 = 0,9
1.3 Оценка погрешности вычислений при решении задачи
Допустимая погрешность на шаге определяется его максимальной величиной. Шаг интегрирования должен быть таким, чтобы быть существенно меньше интеграла, на котором решение дифференциального уравнения значительно изменяется. Удобный способ оценки погрешности и правильного выбора шага дает метод Рунге.
Погрешность для метода второго порядка точности можно вычислить по следующей формуле:
В результате расчетов пришлось уменьшить шаг до , чтобы выполнить условие вычисления погрешности.
1.4 Расчет уточненных решений yут
Если оh/2 по модулю не больше допустимой погрешности, то шаг выбран правильно и находятся уточненные решения:
В противном случае шаг уменьшается вдвое и все повторяется.
1. х1 = 0,1
2. х2 = 0,2
9. х9 = 0,9
.
1.5 Таблица данных 1
|
оhi |
|||||
|
0,100000 |
1,470387 |
1,471685 |
0,001731 |
1,478388 |
|
|
0,200000 |
2,173681 |
2,176455 |
0,003698 |
2,193262 |
|
|
0,300000 |
3,205241 |
3,210517 |
0,007035 |
3,244853 |
|
|
0,400000 |
4,709109 |
4,718558 |
0,012599 |
4,782654 |
|
|
0,500000 |
6,894874 |
6,911161 |
0,021716 |
7,024763 |
|
|
0,600000 |
10,066320 |
10,093667 |
0,036462 |
10,288284 |
|
|
0,700000 |
14,663307 |
14,708347 |
0,060053 |
15,033829 |
|
|
0,800000 |
21,322495 |
21,395582 |
0,097449 |
21,930238 |
|
|
0,900000 |
30,965262 |
31,082471 |
0,156278 |
31,948581 |
1.6 График 1 - значений yh, yh/2, yут
2. Аппроксимирование квадратного трехчлена методом наименьших квадратов
2.1 Таблица 2 - рассчитанных значений для расчета коэффициентов
|
0 |
1 |
1 |
1 |
1 |
1 |
1 |
||
|
0,100000 |
1,470387 |
1,155 |
1,21 |
1,331 |
1,4641 |
1,2705 |
||
|
0,200000 |
2,173681 |
1,318604 |
1,44 |
1,728 |
2,0736 |
1,582325 |
||
|
0,300000 |
3,205241 |
1,490628 |
1,69 |
2,197 |
2,8561 |
1,937817 |
||
|
0,400000 |
4,709109 |
1,670904 |
1,96 |
2,744 |
3,8416 |
2,339266 |
||
|
0,500000 |
6,894874 |
1,859283 |
2,25 |
3,375 |
5,0625 |
2,788925 |
||
|
0,600000 |
10,066320 |
2,055629 |
2,56 |
4,096 |
6,5536 |
3,289006 |
||
|
0,700000 |
14,663307 |
2,259817 |
2,89 |
4,913 |
8,3521 |
3,841689 |
||
|
0,800000 |
21,322495 |
2,471733 |
3,24 |
5,832 |
10,4976 |
4,449119 |
||
|
0,900000 |
30,965262 |
2,691274 |
3,61 |
6,859 |
13,0321 |
5,113421 |
||
|
У |
18,6 |
15,092739 |
24,044054 |
30,26 |
51,336 |
90,1814 |
40,069714 |
2.2 Составление системы уравнений
Сумма квадрата разности соответствует значению функции ѓ(x) и F(x):
Данная функция - сумма трех переменных (a, b, c), задача сводится к отысканию ее минимума. Функция трех переменных имеет минимум, когда все ее частные производные равны нулю.
Получим:
1) - первое уравнение системы
2) - второе уравнение системы
, где
Таким образом, получена система линейных уравнений следующего вида (после замены сумм коэффициентов):
,где , ,
, , , ,
Решение полученной системы дает значения параметров a, b и с для приближенной функции вида .
2.3 Решение системы уравнений методом Гаусса
; ;
2.4 Таблица 3 - данных для расчета погрешности аппроксимации "оапп"
|
0 |
1 |
2,282894 |
1,282894 |
1,645818 |
|
|
0,1 |
1,470387 |
1,206049 |
-0,264338 |
0,069874 |
|
|
0,2 |
2,173681 |
1,146702 |
-1,026979 |
1,054685 |
|
|
0,3 |
3,205241 |
2,104853 |
-1,100388 |
1,210853 |
|
|
0,4 |
4,709109 |
4,080502 |
-0,628607 |
0,395146 |
|
|
0,5 |
6,894874 |
7,073649 |
0,178775 |
0,031961 |
|
|
0,6 |
10,066320 |
11,084294 |
1,017974 |
1,036271 |
|
|
0,7 |
14,663307 |
16,112437 |
1,449130 |
2,099978 |
|
|
0,8 |
21,322495 |
22,158078 |
0,835583 |
0,698199 |
|
|
0,9 |
30,965262 |
29,221217 |
-1,744045 |
3,041694 |
|
|
У = 11,284479 |
|||||
|
д = 1,062284 |
2.5 График 2 - значений yh и F(x)
3. Интерполяция
3.1 Построение интерполяционного многочлена Лагранжа и сгущение значений
Очевидный прием решения данной задачи - вычисление значений ѓ(x), воспользовавшись аналитическими значениями функции ѓ. Для этого - по исходной информации, заданной табличными значениями, отыскивается приближенное значение функции F(x), которая в некотором смысле близка к ѓ(x). И аналитическим выражением которой можно воспользоваться для вычисления, считая, что . Классический подход к решению данной задачи, а именно - построение приближенной функции основано на требовании строгого совпадения значений ѓ(x) и F(x) в xi, при данная функция - интерполируема, а - узлы интерполяции. Как правило, интерполяция функции F(x) отыскивается в виде полинома в степени не больше n.
Значения совпадают со значениями заданной функции в узлах интерполяции.
Интерполирование геометрически обозначает необходимость нахождения алгебраической кривой вида yn, проходящую через систему точек. Значения и значения исходной функции должны быть одинаковы в узлах интерполяции.
Пусть: , , тогда многочлен Лагранжа будет иметь следующий вид:
В соответствии с заданием получим:
Подставив значения , получим:
Теперь подставим значения в поучившийся многочлен. Значения «x» возьмем из таблицы, взяв шаг равный 0,02.
3.2 Расчет погрешности интерполяции
Если известно выражение интерполируемой функции ѓ(х), то можно применить формулу для оценки погрешности интерполирования, как погрешности метода. Остаточный член интерполяционного многочлена имеет следующий вид:
.
Предполагая, что ѓ(х) имеет производную до n+1 порядка включительно, можно предположить, оценивая формулу:
В связи с тем, что формулы Ньютона и Лагранжа - различны способы записи одного и того же интерполяционного многочлена, то для подсчета погрешности интерполяции может быть использована одна и та же формула:
Исходя из задания, формула погрешности будет выглядеть так:
Подставив данные, получим следующие значения погрешностей:
3.3 Таблица данных 4
|
0 |
1 |
0 |
|
|
0,02 |
1,064286 |
0,000668 |
|
|
0,04 |
1,141542 |
0,001142 |
|
|
0,06 |
1,231044 |
0,001437 |
|
|
0,08 |
1,332232 |
0,001566 |
|
|
0,1 |
1,444705 |
0,001558 |
|
|
0,12 |
1,568225 |
0,001408 |
|
|
0,14 |
1,702716 |
0,001150 |
|
|
0,16 |
1,848264 |
0,000811 |
|
|
0,18 |
2,005116 |
0,000418 |
|
|
0,2 |
2,173681 |
0 |
|
|
0,22 |
2,354529 |
0,000414 |
|
|
0,24 |
2,548393 |
0,000793 |
|
|
0,26 |
2,756166 |
0,001110 |
|
|
0,28 |
2,978905 |
0,001339 |
|
|
0,3 |
3,217826 |
0,001455 |
|
|
0,32 |
3,474310 |
0,001442 |
|
|
0,34 |
3,749897 |
0,001289 |
|
|
0,36 |
4,046289 |
0,000992 |
|
|
0,38 |
4,365351 |
0,000556 |
|
|
0,4 |
4,709109 |
0 |
|
|
0,42 |
5,079751 |
0,000648 |
|
|
0,44 |
5,479626 |
0,001347 |
|
|
0,46 |
5,911245 |
0,002041 |
|
|
0,48 |
6,377281 |
0,002662 |
|
|
0,5 |
6,880570 |
0,003131 |
|
|
0,52 |
7,424106 |
0,003358 |
|
|
0,54 |
8,011049 |
0,003247 |
|
|
0,56 |
8,644717 |
0,002705 |
|
|
0,58 |
9,328593 |
0,001644 |
|
|
0,6 |
10,066320 |
0 |
|
|
0,62 |
10,861702 |
0,002259 |
|
|
0,64 |
11,718706 |
0,005109 |
|
|
0,66 |
12,641461 |
0,008445 |
|
|
0,68 |
13,634256 |
0,012046 |
|
|
0,7 |
14,701543 |
0,015541 |
|
|
0,72 |
15,847937 |
0,018350 |
|
|
0,74 |
17,078210 |
0,019625 |
|
|
0,76 |
18,397302 |
0,018174 |
|
|
0,78 |
19,810310 |
0,012362 |
|
|
0,8 |
21,322495 |
0 |
|
|
0,82 |
22,939279 |
0,021793 |
|
|
0,84 |
24,666245 |
0,056758 |
|
|
0,86 |
26,509139 |
0,109695 |
|
|
0,88 |
28,473868 |
0,186701 |
|
|
0,9 |
30,566501 |
0,295448 |
3.4 График 3 - значений y(х), F(X) и P(x)
4. Анализ полученных результатов
Интегрирование методом Эйлера-Коши подтвердило то, что данный метод обладает достаточно высокой точностью. Погрешность вычислений метода Эйлера-Коши зависит от шага h и пропорциональна шагу h во второй степени. Поэтому метод Эйлера-Коши является методом второго порядка. Погрешность не превысила указанной величины.
В результате аппроксимирования были вычислены значения и построена приближающая функция F(x). Значения полученной функции, как и следует из задачи аппроксимирования, получились достаточно близки к табличным. Найденная погрешность аппроксимации была в пределах нормы.
В ходе интерполирования с использованием многочлена Лагранжа была выполнена задача по построению приближенной функции к физической величине f(x) по известным значениям в некоторых точках.
Заключение
дифференциальный уравнение алгоритм многочлен
При выполнении данной курсовой работы осуществлено построение математической модели, описывающей процесс решения дифференциального уравнения. На её основе разработан алгоритм и написана программа.
При построении математической модели произведен расчет значений дифференциального уравнения с шагом h и с шагом h/2, выполнена оценка погрешности вычислений, осуществлена аппроксимация решений дифференциального уравнения. По полученным данным составлена система уравнений.
Используя метод Гаусса, найдены решения данной системы. Затем построен интерполяционный многочлен.
Размещено на Allbest.ru
Подобные документы
Порядок решения дифференциального уравнения 1-го порядка. Поиск частного решения дифференциального уравнения, удовлетворяющего указанным начальным условиям. Особенности применения метода Эйлера. Составление характеристического уравнения матрицы системы.
контрольная работа [332,6 K], добавлен 14.12.2012Проверка непрерывности заданных функций. Интегрирование заданного уравнения и выполние преобразования с ним. Интегрирование однородного дифференциального уравнения. Решение линейного дифференциального уравнения. Общее решение неоднородного уравнения.
контрольная работа [65,3 K], добавлен 15.12.2010Дифференциальные уравнения Риккати. Общее решение линейного уравнения. Нахождение всех возможных решений дифференциального уравнения Бернулли. Решение уравнений с разделяющимися переменными. Общее и особое решения дифференциального уравнения Клеро.
курсовая работа [347,1 K], добавлен 26.01.2015Последовательность решения линейной краевой задачи. Особенности метода прогонки. Алгоритм метода конечных разностей: построение сетки в заданной области, замена дифференциального оператора. Решение СЛАУ методом Гаусса, конечно-разностные уравнения.
контрольная работа [366,5 K], добавлен 28.07.2013Построение таблицы и графика решения линейного дифференциального уравнения. Зависимость погрешности решения от выбора шага интегрирования. Метод Адамса-Башфорта и его применение. Основные функции и переменные, использованные в реализованной программе.
контрольная работа [2,0 M], добавлен 13.06.2012Понятие дифференциального уравнения. Нахождение первообразной для заданной функции. Нахождение решения дифференциального уравнения. Выделение определенной интегральной кривой. Понятие произвольных независимых постоянных. Уравнение в частных производных.
презентация [42,8 K], добавлен 17.09.2013Вычисление и исследование предела и производной функции, построение графиков. Вычисление неопределенных интегралов, площади фигуры, ограниченной графиками функций. Нахождение решения дифференциального уравнения и построение графиков частных решений.
контрольная работа [153,6 K], добавлен 19.01.2010Порядок и принципы составления дифференциального уравнения, методика нахождения неизвестных значений. Замена исходного дифференциального уравнения на систему n-линейных уравнений относительно n-неизвестных. Формирование и решение системы уравнений.
задача [118,8 K], добавлен 20.09.2013Решение дифференциального уравнения методом численного интегрирования Адамса. Методы, основанные на применении производных высших порядков. Формулы, обеспечивающие более высокую степень точности, требующие вычисления третьей производной искомого решения.
курсовая работа [81,9 K], добавлен 29.08.2010Общий вид линейного однородного уравнения. Нахождение производных, вещественные и равные корни характеристического уравнения. Пример решения дифференциального уравнения с постоянными коэффициентами. Общее и частное решение неоднородного уравнения.
презентация [206,3 K], добавлен 17.09.2013
