Решение уравнений состояния

Методы решения
Аналитические	Численные
Операторный метод Метод матричной экспоненты	Явный (неявный) методы Эйлера Метод Рунге-Кутта Метод Адамса

Пример

Рис. 4.3.

Операторный метод:

1. Имеем систему уравнений состояния

2. Отобразим эту систему уравнений в операторную область

3. После преобразований получим

4. Переходим к оригиналам

6.

Строим график:

Рис. 4.4.

Численный метод (явный метод Эйлера).

Решение получаемое численным методом имеет вид таблицы, в которой даны дискретные значения искомой величины в фиксированные моменты времени. Ось времени при этом разбивается на интервалы длительностью в один шаг интегрирования h. Дискретные значения функции вычисляются по

Рис. 4.5.

.рекурентной формуле, которая является дискретным аналогом непрерывного дифференциального уравнения. Существуют различные способы получить рекурентное уравнение. В методе Эйлера на каждом шаге интегрирования кривую x(t) апроксимируют рядом Тейлора, содержащем первые два члена:

. (4.3)

Эта формула позволяет по одному предыдущему x_k значению рассчитать следующее x_k+1. Для вычисления значения искомой величины в момент времени t необходимо последовательно рассчитать [t/h]-1 предыдущих значений.

Численные методы расчета относятся к приближенным методам. Чем больше сложность рекурентного уравнения, тем меньше зависит точность вычисления от величины шага интегрирования. С уменьшением шага интегрирования решение получается точнее, но такой способ повышения точности расчета имеет ограничения. Во-первых, увеличивается количество операций и время вычисления, а во-вторых, накапливается другая погрешность, связанная с округлением вещественных чисел в вычислительном устройстве.

Пример

Рис. 4.6.

Метод Эйлера:

1. Имеем уравнения в форме Коши.

2. Получим рекурентные уравнения для искомых величин по методу явному Эйлера.

3. В течении 0.25 с рассчитаем переходный процесс для двух значений шага интегрирования (h=1.0 мс и h=5.0 мс) и сравним эти значения с точным решением, полученным аналитическим методом.

Таблица дискретных значений

переменных состояния (h=5.0 мс)

t=5.000 мc iL1=.3500 A Uc= 96.250 В

t=10.000 мc iL1=.7600 A Uc= 88.594 В

t=15.000 мc iL1= 1.1980 A Uc= 87.020 В

t=20.000 мc iL1= 1.6330 A Uc= 91.117 В

t=25.000 мc iL1= 2.0366 A Uc= 100.140 В

t=30.000 мc iL1= 2.3850 A Uc= 113.080 В

t=35.000 мc iL1= 2.6599 A Uc= 128.757 В

t=180.000 мc iL1= 2.0826 A Uc= 182.312 В

t=185.000 мc iL1= 2.0189 A Uc= 185.555 В

t=190.000 мc iL1= 1.9407 A Uc= 187.598 В

t=195.000 мc iL1= 1.8542 A Uc= 188.407 В

t=200.000 мc iL1= 1.7658 A Uc= 188.033 В

t=205.000 мc iL1= 1.6815 A Uc= 186.602 В

t=210.000 мc iL1= 1.6064 A Uc= 184.295 В

t=215.000 мc iL1= 1.5448 A Uc= 181.339 В

t=220.000 мc iL1= 1.4995 A Uc= 177.982 В

t=225.000 мc iL1= 1.4721 A Uc= 174.478 В

t=230.000 мc iL1= 1.4629 A Uc= 171.069 В

t=235.000 мc iL1= 1.4710 A Uc= 167.972 В

t=240.000 мc iL1= 1.4944 A Uc= 165.363 В

t=245.000 мc iL1= 1.5302 A Uc= 163.372 В

t=250.000 мc iL1= 1.5751 A Uc= 162.078 В

Рис. 4.7.

Дата добавления: 2014-01-07