Учебники

Системы управления — государственная космическая модель

Модель пространства состояний системы линейного инварианта времени (LTI) может быть представлена ​​как

 точкаX=AX+BU

Y=CX+DU

Первое и второе уравнения известны как уравнение состояния и выходное уравнение соответственно.

Куда,

  • X и  dotX — вектор состояния и вектор дифференциального состояния соответственно.

  • U и Y — входной вектор и выходной вектор соответственно.

  • А — системная матрица.

  • B и C — входная и выходная матрицы.

  • D — матрица прямой связи.

X и  dotX — вектор состояния и вектор дифференциального состояния соответственно.

U и Y — входной вектор и выходной вектор соответственно.

А — системная матрица.

B и C — входная и выходная матрицы.

D — матрица прямой связи.

Основные понятия государственной космической модели

Следующая основная терминология используется в этой главе.

государственный

Это группа переменных, которая суммирует историю системы, чтобы предсказать будущие значения (результаты).

Государственная переменная

Количество требуемых переменных состояния равно количеству элементов хранения, присутствующих в системе.

Примеры — ток, протекающий через индуктор, напряжение на конденсаторе

Государственный вектор

Это вектор, который содержит переменные состояния в качестве элементов.

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

Модель пространства состояний из дифференциального уравнения

Рассмотрим следующую серию схемы RLC. Он имеет входное напряжение, vi(t), а ток, протекающий по цепи, составляет i(t).

Дифференциальное уравнение RLC

В этой схеме есть два накопителя (индуктор и конденсатор). Таким образом, число переменных состояния равно двум, и этими переменными состояния являются ток, протекающий через индуктор, i(t) и напряжение на конденсаторе, vc(t).

Из схемы выходное напряжение v0(t) равно напряжению на конденсаторе, vc(t).

v0(т)=vc(т)

Примените КВЛ вокруг петли.

vi(т)=Ri(т)+L гидроразрыва текстdя(т) текстdт+vc(т)

 Rightarrow frac textdi(t) textdt= fracRi(t)L fracvc(t)L+ гидроразрываvi(т)L

Напряжение на конденсаторе составляет —

vc(t)= frac1C inti(t)dt

Дифференцируйте вышеприведенное уравнение по времени.

 гидроразрыва текстdvc(т) текстdт= гидроразрывая(т)C,

Вектор состояния, X= beginbmatrixi(t)vc(t) endbmatrix

Вектор дифференциального состояния,  dotX= beginbmatrix frac textdi(t) textdt frac textdvc(t) textdt endbmatrix

Мы можем организовать дифференциальные уравнения и выходное уравнение в стандартную форму модели пространства состояний как,

\ dot {X} = \ begin {bmatrix} \ frac {\ text {d} i (t)} {\ text {d} t} \\\ frac {\ text {d} v_c (t)} { \ text {d} t} \ end {bmatrix} = \ begin {bmatrix} — \ frac {R} {L} & — \ frac {1} {L} \\\ frac {1} {C} & 0 \ end {bmatrix} \ begin {bmatrix} i (t) \\ v_c (t) \ end {bmatrix} + \ begin {bmatrix} \ frac {1} {L} \\ 0 \ end {bmatrix} \ begin {bmatrix } v_i (t) \ end {bmatrix}

Y = \ begin {bmatrix} 0 & 1 \ end {bmatrix} \ begin {bmatrix} i (t) \\ v_c (t) \ end {bmatrix}

Куда,

A = \ begin {bmatrix} — \ frac {R} {L} & — \ frac {1} {L} \\\ frac {1} {C} & 0 \ end {bmatrix}, \: B = \ begin {bmatrix} \ frac {1} {L} \\ 0 \ end {bmatrix}, \: C = \ begin {bmatrix} 0 & 1 \ end {bmatrix} \: и \: D = \ begin {bmatrix } 0 \ end {bmatrix}

Модель пространства состояний из передаточной функции

Рассмотрим два типа передаточных функций, основанные на типе терминов, присутствующих в числителе.

  • Передаточная функция с постоянным членом в Числителе.
  • Передаточная функция, имеющая полиномиальную функцию ‘s’ в Числителе.

Передаточная функция с постоянным членом в Числителе

Рассмотрим следующую передаточную функцию системы

$$ \ гидроразрыва {Y (s)} {U (s)} = \ {гидроразрыва b_0} {s ^ п + a_ {п-1} s ^ {п-1} + … + a_1s + A_0} $ $

Переставить приведенное выше уравнение как

(sn+an1sn1+...+a0)Y(s)=b0U(s)

Примените обратное преобразование Лапласа с обеих сторон.

 гидроразрыва текстdпу(т) текстdтп+ап1 гидроразрыва текстdN1у(т) текстdтN1+...a1 гидроразрыва текстdу(г) текстdт+a0y(т)=B0и(т)

Позволять

у(г)=x1

 гидроразрыва текстdу(г) текстdт=x2= точках1

 гидроразрыва текстd2у(т) текстdт2=x3= точках2

.

.

.

$$ \ гидроразрыва {\ текст {d} ^ {N-1} у (г)} {\ текст {d} т ^ {N-1}} = x_n = \ точки {х} _ {п-1} $ $

 гидроразрыва текстdпу(т) текстdтп= точкахn

и u(t)=u

Затем,

 dotxn+an1xn+...+a1x2+a0x1=b0u

Из приведенного выше уравнения мы можем написать следующее уравнение состояния.

 dotxn=a0x1a1x2...an1xn+b0u

Выходное уравнение —

у(г)=у=x1

Модель пространства состояний —

 dotX= beginbmatrix dotx1 dotx2 vdots dotxn1 dotxn endbmatrix

= \ begin {bmatrix} 0 & 1 & 0 & \ dotso & 0 & 0 \\ 0 & 0 & 1 & \ dotso & 0 & 0 \\\ vdots & \ vdots & \ vdots & \ dotso & \ vdots & \ vdots \\ 0 & 0 & 0 & \ dotso & 0 & 1 \\ — a_0 & -a_1 & -a_2 & \ dotso & -a_ {n-2} & -a_ {n-1} \ end {bmatrix } \ begin {bmatrix} x_1 \\ x_2 \\\ vdots \\ x_ {n-1} \\ x_n \ end {bmatrix} + \ begin {bmatrix} 0 \\ 0 \\\ vdots \\ 0 \\ b_0 \ end {bmatrix} \ begin {bmatrix} u \ end {bmatrix}

Y = \ begin {bmatrix} 1 & 0 & \ dotso & 0 & 0 \ end {bmatrix} \ begin {bmatrix} x_1 \\ x_2 \\\ vdots \\ x_ {n-1} \\ x_n \ end {bmatrix}

Здесь D= left[0 right].

пример

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

 гидроразрываY(ы)U(ы)= гидроразрыва1S2+S+1

Переставить, приведенное выше уравнение как,

(S2+S+1)Y(s)=U(ы)

Примените обратное преобразование Лапласа с обеих сторон.

 гидроразрыва текстd2у(т) текстdт2+ гидроразрыва текстdу(г) текстdт+у(т)=U(T)

Позволять

у(г)=x1

 гидроразрыва текстdу(г) текстdт=x2= точках1

и u(t)=u

Тогда уравнение состояния

 точках2=x1x2+и

Выходное уравнение

у(г)=у=x1

Модель пространства состояний

\ dot {X} = \ begin {bmatrix} \ dot {x} _1 \\\ dot {x} _2 \ end {bmatrix} = \ begin {bmatrix} 0 & 1 \\ — 1 & -1 \ end {bmatrix} \ begin {bmatrix} x_1 \\ x_2 \ end {bmatrix} + \ begin {bmatrix} 0 \\ 1 \ end {bmatrix} \ left [u \ right]

Y = \ begin {bmatrix} 1 & 0 \ end {bmatrix} \ begin {bmatrix} x_1 \\ x_2 \ end {bmatrix}

Передаточная функция, имеющая полиномиальную функцию ‘s’ в Числителе

Рассмотрим следующую передаточную функцию системы

 fracY(s)U(s)= fracbnsn+bn1sn1+...+b1s+b0sn+an1sn1+...+a1s+a0

 Rightarrow fracY(s)U(s)= left( frac1sn+an1sn1+...+a1s+a0 right)(bnsn+bn1sn1+...+b1s+b0)

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

 fracY(s)U(s)= left( fracV(s)U(s) right) left( fracY(s)V(s) right)

Вот,

 fracV(s)U(s)= frac1sn+an1sn1+...+a1s+a0

Переставить приведенное выше уравнение как

(ып+ап1сN1+...A0)V(S)=U(ы)

Примените обратное преобразование Лапласа с обеих сторон.

 гидроразрыва текстdпу(т) текстdтп+ап1 гидроразрыва текстdN1v(т) textdtn1+...+a1 frac textdv(t) textdt+a0v(t)=u(т)

Позволять

V(T)=x1

 гидроразрыва текстdv((т) текстdт=x2= точках1

 гидроразрыва текстd2v(т) текстdт2=x3= точких2

.

.

.

$$ \ гидроразрыва {\ текст {d} ^ {N-1} V (T)} {\ текст {d} т ^ {N-1}} = x_n = \ точка {х} _ {п-1} $ $

 гидроразрыва текстdпу(т) текстdтп= точкахn

и u(t)=u

Тогда уравнение состояния

 точкахn=a0x1a1x2...ап1xn+U

Рассматривать,

 гидроразрываY(s)V(S)=bnsп+Bп1сN1+...+b1sB0

Переставить приведенное выше уравнение как

Y(s)=(bnsп+Bп1сN1+...+b1sB0)V(S)

Примените обратное преобразование Лапласа с обеих сторон.

у(г)=bn гидроразрыва текстdпу(т) текстdтп+Bп1 гидроразрыва текстdп1V(T) текстdтN1+...b1 гидроразрыва текстdv(т) текстdт+b0v(т)

Подставляя переменные состояния и y(t)=y в вышеприведенном уравнении, получим выходное уравнение как,

у=bn точкахn+bп1xn+...+b1x2+b0x1

Замените значение  dotxn в приведенном выше уравнении.

у=bn(a0x1a1x2...ап1xn+и)+Bп1xn+...+b1x2b0x1

y=(b0bna0)x1+(b1bna1)x2+...+(bn1bnan1)xn+bnu

Модель пространства состояний

 dotX= beginbmatrix dotx1 dotx2 vdots dotxn1 dotxn endbmatrix

= \ begin {bmatrix} 0 & 1 & 0 & \ dotso & 0 & 0 \\ 0 & 0 & 1 & \ dotso & 0 & 0 \\\ vdots & \ vdots & \ vdots & \ dotso & \ vdots & \ vdots \\ 0 & 0 & 0 & \ dotso & 0 & 1 \\ — a_0 & -a_1 & -a_2 & \ dotso & -a_ {n-2} & -a_ {n-1} \ end {bmatrix } \ begin {bmatrix} x_1 \\ x_2 \\\ vdots \\ x_ {n-1} \\ x_n \ end {bmatrix} + \ begin {bmatrix} 0 \\ 0 \\\ vdots \\ 0 \\ b_0 \ end {bmatrix} \ begin {bmatrix} u \ end {bmatrix}

Y=[b0bna0 quadb1bna1 quad... quadbn2bnan2 quadbn1bnan1] beginbmatrixx1x2 vdotsxn1xn endbmatrix

Если bn=0, то

Y=[b0 quadb1 quad... quadbn2 quadbn1] beginbmatrixx1x2 vdotsxn1xn endbmatrix