Методы численного решения дифференциальных уравнений. Усовершенствованный метод Эйлера. Метод Адамса - Бэшфорта

Страницы работы

Содержание работы

       МОДЕЛИРОВАНИЕ СИСТЕМ         АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ         СКОРОСТЬЮ ЭЛЕКТРОПРИВОДОВ         ПОСТОЯННОГО ТОКА

 МЕТОДИЧЕСКИЕ УКАЗАНИЯ к контрольным работам

 по дисциплине «Моделирование в электроприводах»

для студентов специальности 1-53 01 05

          «Автоматизированные электроприводы»            заочной формы обучения

Введение

По дисциплине “Моделирование в электроприводах”, согласно учебного плана, выполняется контрольная работа.

Цель работы – закрепить полученные теоретические знания и научиться применять их для практических расчетов.

Контрольная работа выполняется в тетради. Все расчеты, рисунки должны быть оформлены в соответствии с требованиями ЕСКД. При ответе на теоретические вопросы обязательно давать ссылки на использованную литературу. Все выполненные в ходе работы распечатки вклеиваются в тетрадь. Контрольная работа должна содержать:

–  титульный лист;

–  ответы на 2 теоретических вопроса;

–  задание и исходные данные;

–  все выполненные в ходе работы расчеты;

–  структурную       схему          системы           управления,        которая

рассчитывается;

–  диаграммы рассчитанных переходных процессов;

–  анализ динамических показателей (перерегулирование и длительность процесса);

1.  Теоретические сведения

1.1.  Методы численного решения дифференциальных уравнений

Расчетные соотношения и параметры основных численных методов приведены в табл. 1.1. В расчетных соотношениях приняты обозначения:

Y – вектор переменных модели, i – номер текущего шага расчета, для которого рассчитываются значения переменных модели (соответствует моменту времени t),

Δt – шаг расчета (интервал времени между соседними шагами расчета),

f – вектор значений правых частей дифференциальных уравнений, записанных в форме Коши. ν – номер текущей итерации.


Таблица 1.1. Расчетные соотношения и параметры численных методов.

Название метода, порядок точности,

условие численной устойчивости

Расчетное соотношение

Эйлера. p =1, Δt < T .

Yi = Yi−1 + Δt fi−1

Неявный метод Эйлера. p =1.

Yi(ν) = Yi−1 + Δt fi(ν−1)

Усовершенствованный метод Эйлера. p = 2, Δt < 2⋅T .

Yi(1) = Yi−1 + Δt fi−1

Yi( )

,

 .

Метод РунгеКутта. p = 4,

Δt < 2,7853⋅T .

Yi−0,5( )1 =Yi−1 + Δt fi−1,

2

Yi−0,5( )2 =Yi−1 + Δt fi−0,5( )1 ,

2

Yi( )1 =Yi−1 + Δt fi−0,5( )2 ,

Yi( )

 .

Метод АдамсаБэшфорта. p = 2.

Yi  fi

−1 i−2) − f

Метод АдамсаМултона. p = 4.

Yi( )1

Yi( )

.

,

Метод Хэмминга.

Yi( )1

Yi( )

Yi( )

( )

,

,

 .

Метод

Башарина. p = 2, Δt < 2⋅T .

Yi( )ν = Yi−1

                     ⎛                                       (ν−1)

+ Δt f ⎜⎜ti Δ2t , X⎜⎝ti Δ2t ⎟,Yi1 +2Yi     ⎟⎠⎞ ⎝

Таким образом Yi – вектор значений переменных на текущем шаге расчета, т. е в момент времени t = i⋅Δt , Yi1 – на предыдущем шаге расчета ( t − Δt ) и т.д. Значения Yi0,5 соответствуют времени посредине между текущим и предыдущим шагами расчета.

Некоторые методы имеют несколько итераций, номер которой записывается в верхнем индексе. Т.е. значения переменных модели на текущем шаге рассчитываются несколько раз, постепенно уточняясь.

Если в расчетном соотношении метода отсутствует выражение для первой итерации (Yi( )1 ), то его следует рассчитывать по методу

Эйлера.

Индексы значений правых частей уравнений показывают какие значения переменных Y и времени подставляются в выражения правых частей f (t, X,Y) . Например: fi1 = f (t − Δt, X(t − Δt),Yi1) , fi(ν−1) = f (t, X( )t ,Yi(ν−1)), fi−0,5(1) = f (t − 0,5⋅Δt, X(t − 0,5⋅Δt),Yi−0,5( )1 ).

1.2.  Модели элементов и систем управления электроприводами

С моделями элементов систем управления электроприводами можно ознакомиться в [1], [2].

Модели систем, в целом, и методики расчета параметров приведены в [3].

Похожие материалы

Информация о работе