- •Численное дифференцирование
- •Научиться выполнять численное дифференцирование функций, основанное на интерполяционных многочленах Ньютона и Лагранжа.
- •Экстремумы
- •Научиться находить экстремум функции одной переменной методом золотого сечения;
- •Научится находить экстремум функции одной переменной методом Ньютона.
- •Численное решение дифференциальных уравнений
- •Научиться решать дифференциальные уравнения усовершенствованным методом Эйлера;
- •Научиться решать дифференциальные уравнения методом Рунге-Кутта.
Численное решение дифференциальных уравнений
Цель работы
Научиться решать дифференциальные уравнения усовершенствованным методом Эйлера;
Научиться решать дифференциальные уравнения методом Рунге-Кутта.
Оборудование
IBM PC, Mathcad 7 Pro, Microsoft Excel.
Ход работы
3. 1 Вариант
Решите дифференциальное уравнение
усовершенствованным методом Эйлера на отрезке [ ]
при заданном начальном условии у( )= и шаге интегрирования h=0,1.
Решите дифференциальное уравнение
методом Рунге-Кутта на отрезке [ ]
при заданном начальном условии у( )= в десяти равноотстоящих точках промежутка.
3. 2 Допуск к работе
Чем численное решение дифференциальных уравнений отличается от аналитического?
__________________________________________________________________________________________________________________________________________________________________
Как вычисляется шаг интегрирования h?
__________________________________________________________________________________________________________________________________________________________________
В каком виде должно быть записано дифференциальное уравнение для численного решения его методами Эйлера и Рунге-Кутта
__________________________________________________________________________________________________________________________________________________________________
Запишите формулу для вычисления xi+1/2 в усовершенствованном методе Эйлера.
__________________________________________________________________________________________________________________________________________________________________
Запишите формулу для вычисления уi+1/2 в усовершенствованном методе Эйлера?.
__________________________________________________________________________________________________________________________________________________________________
Запишите формулу для вычисления уi+1 в усовершенствованном методе Эйлера.
__________________________________________________________________________________________________________________________________________________________________
Запишите формулы для вычисления k1, k2, k3, k4 в методе Рунге-Кутта.
__________________________________________________________________________________________________________________________________________________________________
__________________________________________________________________________________________________________________________________________________________________
Запишите формулу для вычисления уi+1 в методе Рунге-Кутта.
__________________________________________________________________________________________________________________________________________________________________
Результаты работы
4.1
усовершенствованный метод Эйлера |
||||||
xi |
yi |
f(xi;yi) |
xi+1/2 |
yi+1/2=yi+h/2*f(xi;yi) |
f(x+1/2;yi+1/2) |
yi+1=yi+h*f(xi+1/2;yi+1/2) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Первые два столбца таблицы являются решением данной задачи
хi |
|
|
|
|
|
|
|
|
|
|
|
|
yi |
|
|
|
|
|
|
|
|
|
|
|
|
4.2
метод Рунге- Кутта |
|||||||||||
xi |
yi |
k1 |
xi+h/2 |
yi+k1/2 |
k2 |
yi+k2/2 |
k3 |
xi+h |
yi+k3 |
k3 |
yi+1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Первые два столбца таблицы являются решением данной задачи
хi |
|
|
|
|
|
|
|
|
|
|
|
|
yi |
|
|
|
|
|
|
|
|
|
|
|
|
5. Вывод
В ходе выполнения данной работы ________________________________________________________________________________________________________________________________________________________________________________________________________________________