- •Методическое пособие
- •Введение
- •Теория вопроса
- •Погрешности исходных данных
- •Погрешности вычислений
- •Погрешности ограничения
- •Посмотрите, сколько значащих цифр осталось?
- •Индивидуальные задания и последовательность выполнения работы
- •Исходные данные к выполнению работы Таблица 1
- •Контрольные вопросы
- •Табулирование функций
- •Теория вопроса
- •Итерационный процесс продолжается до тех пор, пока значения х1(k), х2(k), х3(k) не станут близкими с заданной погрешностью к значениям х1(k-1), х2(k-1), х3(k-1).
- •Изменим знаки у элементов с нечетной суммой индексов, тогда
- •Зададим начальное приближение
- •Исключим переменную х2 из третьего уравнения системы. Для этого умножим на 0,5 и вычтем из третьего уравнения. Получим новую систему эквивалентную двум предыдущим:
- •Вычислим первые приближения
- •Индивидуальные задания и последовательность выполнения работы
- •Варианты для лабораторных работ
- •Контрольные вопросы
- •Теория вопроса
- •Индивидуальные задания и последовательность выполнения работы
- •Теория вопроса
- •Индивидуальные задания и последовательность выполнения работы
- •Контрольные вопросы
- •Теория вопроса
- •Индивидуальные задания и последовательность выполнения работы
- •Исходные данные к выполнению работы Таблица 5
- •Контрольные вопросы:
- •Теория вопроса
- •Решением дифференциального уравнения является неизвестная функция
- •Индивидуальные задания и последовательность выполнения работы
- •Контрольные вопросы
- •Рекомендуемая литература
Индивидуальные задания и последовательность выполнения работы
Используя возможности ЭТ Excel решить дифференциальное уравнение методом Рунге-Кутта по предложенному преподавателем варианту. Результаты работы сохранить в своей рабочей папке и показать преподавателю.
Исходные данные для выполнения задания Таблица 6
№ вар-та |
Дифференциальное уравнение |
Шаг, Н |
Кол-во шагов, n |
Начальное х0 |
Начальное y0 |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 |
y'=ye-x+0,7 y'=ecos(yx)-3sin(x) y'=cos(x)ye-x+1,2 y'=ecos(10xy)-0,5 y'=sin(xy)e-x+0,5 y'=ecos(xy)-1 y'=sin(xy+1) y'=ecos(0,25xy)-1 y'=cos(2,1x)y+0,8 y'=ecos(0,25y)-0,2 y'=e-x+0,2y y'=esin(0,25y)-0,2 y'=0,5x2+0,05ey y'=e-x+sin(0,25y) y'=0,5y+x/y y'=esin(10y)+cos(5x) |
0,1 0,1 0,1 0,1 0,1 0,1 0,1 0,1 0,1 0,1 0,1 0,1 0,1 0,1 0,1 0,1 |
20 10 10 10 20 10 20 10 10 10 10 10 20 10 20 10 |
0 0 0 0 1 0 1 0 0 0 0 0 1 0 2 0 |
1 1 1 1 1 1 2 1 1 1 1 1 2 1 4 1 |
Контрольные вопросы
-
Что является решением дифференциального уравнения?
-
Почему для решения дифференциального уравнения нужно задаваться начальным условием?
-
Зачем дифференциальное уравнение преобразуют к виду y'=f(x,y)?
-
Почему метод Рунге-Кутта точнее метода Эйлера?
-
За счет чего возникает погрешность в методе Эйлера? Как ее уменьшить?
-
Какова последовательность вычисления очередной точки в методе Рунге-Кутта?
-
Поясните алгоритм вычисления решения дифференциального уравнения методом Рунге-Кутта.
Рекомендуемая литература
-
Фигурнов В.Э. IBM PC для пользователя. –М.: Финансы и статистика, 1998 –288 с.
-
Турчак Л.И. Основы численных методов. –М.: Наука, 1997-320 с.
-
Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся ВТУЗов. –М.: Наука, 1986 –544 с.
-
Самарский А.А. Введение в численные методы. –М.: Наука, 1982–272 с.