ЧИСЛЕННЫЕ МЕТОДЫ РЕШЕНИЯ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ. - презентация

1 ЧИСЛЕННЫЕ МЕТОДЫ РЕШЕНИЯ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

2 В настоящее время разработано большое число методов численного интегрирования систем дифференциальных уравнений. К их числу можно отнести метод Рунге- Кутта, явный и неявный методе Эйлера, метод Милна и т. д. Однако, несмотря на большое разнообразие этих методов, алгоритм программ для всех их примерно одинаков и состоит из следующих блоков, а именно:

3 Алгоритм программ блока исходных и расчета дополнительных данных; блока формирования начальных условий и итерационных циклов; блока формирования итерационных уравнений в зависимости от принятого метода численного интегрирования дифференциальных уравнений; блока формирования решения дифференциальных уравнений и обработки полученных результатов.

4 ЭЛЕМЕНТЫ ЧИСЛЕННЫХ МЕТОДОВ Основным элементом численных методов является производная функции. Производная функции - есть предел отношения приращения функции к приращению независимой переменной при стремлении к нулю приращения независимой переменной

5 При численном нахождении производной заменяют отношение бесконечно малых приращений функций и аргумента отношением конечных разностей. Очевидно, что чем меньше будет приращение аргумента, тем точнее численное значение производной.

6 Методы графического представления производной В основе методов графического представления производной лежит геометрический смысл производной. Для вычисления первой производной разработаны двухточечные методы численного дифференцирования.

7 Двухточечные методы Для двухточечных методов при вычислении производных используется значение функции в двух точках. Приращение аргумента задается тремя способами, откладывая Δ x = h вправо, влево и в обе стороны от исследуемой точки. Соответственно получается три двухточечных метода численного дифференцирования

8

9 Метод 1

10 Метод 2

11 Метод 3

12 Численное решение дифференциальных уравнений Дифференциальным уравнением первого порядка называется уравнение вида или

13 Функция y(x), при подстановке которой уравнение обращается в тождество, называется решением дифференциального уравнения.

14 Метод Эйлера В основе метода Эйлера лежит идея графического построения решения дифференциального уравнения c начальными условиями y 0 =y(x 0 )

15

16 Варианты вывода формул Существуют аналитический и графический варианты вывода расчетных формул метода Эйлера Представим это уравнение в виде

17 Тогда можно записать:

18 Расчетные формулы 1-го шага Тогда расчетные формулы для первого шага можно представить в виде:

19 Расчетные формулы i -го шага Расчетные формулы i-го шага по аналогии с первым шагом можно записать в виде:

20 Если обозначить то расчетные формулы можно записать в виде

21 Численное решение системы дифференциальных уравнений Системой дифференциальных уравнений называется система вида

22 или где x – независимый аргумент, y i – зависимая функция, y i|x=x0 =y i0 – начальные условия. Функции y i (x), при подстановке которой система уравнений обращается в тождество, называется решением системы дифференциальных уравнений.

23 Метод Эйлера Итеррацонные уравнения для численного решения системы дифференциальных уравнений можно записать в виде:

24 Математическая модель двигателя постоянного тока Математическое описания процессов электромеханического преобразования энергии в двигателе постоянного тока (ДПТ) содержит: систему уравнений равновесия напряжений, в которых в качестве переменных приняты значения токов в обмотках ОВ и ОЯ; уравнение механического равновесия;

25 выражение для электромагнитного момента. В результате определенного вида преобразований получают систему уравнений, описывающую электромеханические процессы в двигателе постоянного тока, в удобном для математического моделирования виде, в форме уравнений Коши:

26 Система уравнений ДПТ

27 Для реализации такой модели в среде Mathcad с использованием метода Эйлера необходимо: сформировать исходные данные, которые включают в себя параметры двигателя, напряжения ОВ и ОЯ, суммарный момент инерции двигателя и момент статической нагрузки на валу двигателя; определить и записать начальные условия для исследуемых переменных ДПТ; определяют число итераций.

28 Система итерационных уравнений