
Математические основы теории систем
..pdf
151
В этих случаях нет необходимости использовать интеграл свёртки и записывать решение в форме уравнения (3.60). Подставив соотношение (3.61) в уравнение (3.58), получим:
Y (s) = |
B(s) |
|
N (s) |
+ |
M (s) |
= |
B(s)N (s) |
+ |
M (s) |
. |
A(s) |
P(s) |
A(s) |
D(s) |
|
||||||
|
|
|
|
|
A(s) |
Для перехода в область переменной t можно воспользоваться любым методом определения обратного преобразования Лапласа, например теоремой разложения.
······················ |
Пример 3.18 ······················ |
Решить уравнение
|
|
d2 y |
+3 |
dy |
+ 2y = |
dr |
+3r , |
|
|
dt2 |
dt |
|
|||
|
|
|
|
dt |
|||
при r(t) = t при t ≥ 0, |
и нулевых начальных условиях. |
||||||
0 при t < 0, |
|
|
|
|
|
|
|
Применим преобразование Лапласа к дифференциальному уравнению, учитывая, что L{t} =1s2 :
(s2 + 3s + 2)Y (s) = (s + 3) s12 .
Решим полученное алгебраическое уравнение относительно Y (s):
Y (s) = s2 (s2s++33s + 2) .
Перейдём к оригиналу, используя результат примера 3.48:
y(t) = |
|
− |
7 |
+ |
3 |
t + 2e |
−t |
− |
1 |
−2t |
|
|
|
|
e |
1(t). |
|||||
|
|
4 |
|
2 |
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
····································································
··························································
Контрольные вопросы по главе 3
··························································
1.В каких случаях возможна линеаризация нелинейных уравнений?
2.В чем различия и что общего между исходным нелинейным уравнением и линеаризованным?
152
3.Как записывается общее решение однородного линейного дифференциального уравнения в случае некратных и кратных корней характеристического уравнения?
4.Какова форма записи решения для комплексных корней характеристического уравнения?
5.Какие методы существуют для нахождения частного решения уравнения?
6.К каким вынуждающим функциям применим метод неопределенных коэффициентов при решении неоднородных дифференциальных уравнений?
7.В чем состоит необходимое и достаточное условие линейной независимости n решений однородного линейного дифференциального уравнения n-го порядка?
8.К каким функциям не применимо преобразование Фурье?
9.Как получить из преобразования Лапласа преобразование Фурье?
10.Что такое передаточная функция системы?
11.Назовите основные свойства преобразования Лапласа.
12.Что такое абсцисса абсолютной сходимости?
13.Как получить оригинал по изображению?
14.Чем отличается преобразование Карсона – Хевисайда от преобразования Лапласа?

153
4 Операторное описание дискретных по времени систем
Будем полагать в этом разделе, что функции r(t) и y(t), то есть входной и выходной сигналы системы, определены на счетном множестве моментов времени. Другими словами, время течет дискретно, квантами через равные промежутки, обозначаемые в этом разделе буквой Т, то есть t = kT , где k N0 . Для упрощения записи можно выбрать соответствующий масштаб по оси времени и положить Т =1, то есть считать r (k ) и y(k ) как функции, определенные только при целых значениях k.
4.1 Прямой и обратный разностные операторы
4.1.1Оператор сдвига и разностный оператор
Определим оператор сдвига E так:
|
{ |
|
( |
|
)} |
|
( |
) |
|
E |
|
y |
|
k |
|
= y |
|
k +1 . |
(4.1) |
Последовательное применение этого оператора дает в общем случае |
|
||||||||
En {y(k)} = y(k + n) , |
(4.2) |
||||||||
где n N0 . |
|
|
|
|
|
|
|
|
|
Разностный оператор ∆ можно определить как |
|
||||||||
∆y(k) = y(k +1) − y(k) . |
(4.3) |
Оператор, определяемый формулой (4.3), называют еще правым разностным оператором, и он задает так называемую первую прямую разность функции y(k ), в отличие от используемого иногда левого разностного оператора ,
определяемого выражением
y(k) = y(k) − y(k −1)
и задающего первую обратную разность функции y(k ).
Выражение (4.3) с учетом (4.1) можно записать в виде
∆y(k) = (E −1) y(k) ,
где операторы ∆ и Е связаны соотношением |
|
∆ = E −1. |
(4.4) |
Разности второго, третьего и более высокого порядков определяются по очевидным формулам:

154
∆2 y(k) = ∆(∆y(k)) = y(k + 2) − 2y(k +1) + y(k) ,
∆3 y(k ) = ∆(∆2 y(k )) = y(k + 3) − 3y(k + 2) + 3y(k +1) + y(k )
или, в общем случае,
∆n y(k) = ∑(−1)n |
n y(k + n − r), |
(4.5) |
n |
|
|
r=1 |
r |
|
где через n обозначены биномиальные коэффициенты.
r
С учетом уравнений (4.2) и (4.4) из выражения (4.5) получим:
∆n y(k ) = (E −1)n y(k) = ∑(−1)n |
n En−r y(k ). |
n |
|
r=1 |
r |
Операторы ∆ и E являются линейными операторами, то есть справедливы следующие соотношения, например, для оператора ∆:
∆cy(k) = c∆y(k),
∆n (y(k) + x(k)) = ∆n y(k) + ∆n x(k) ,
∆n∆m y(k) = ∆m∆n y(k) = ∆m+n y(k).
где с – константа, m и n – целые положительные числа.
Таким образом, оператор ∆ для функций дискретного переменного является аналогом дифференциального оператора p = ddt для непрерывных функций. Чтобы еще раз подчеркнуть эту аналогию, рассмотрим производную от непрерывной функции f (t):
df (t) = lim f (t + T )− f (t) . dt T →0 T
Если функцию f (t) рассматривать только в дискретные моменты времени t = kT (k N0), то оператор сдвига и разностный оператор дадут выражения:
Ef (t) = f (t + T ) и ∆f (t) = f (t + T )− f (t). Тогда получим:
|
|
df (t) |
= lim |
∆f (t) |
, |
|
|
|
|
|
|
|
|||
|
|
dt |
T →0 |
T |
|
||
или для случая m-й производной: |
|
|
|
|
|||
|
dm f (t) |
= lim |
∆m f (t) |
. |
|||
|
|
|
|
|
|
||
|
|
dtm |
Tm |
||||
|
|
T→0 |
|

155
Существуют и разностные формулы, аналогичные (но не идентичные) формулам дифференцирования произведения и дроби.
Например:
∆(y(k) x(k)) = ∆y(k)∆x(k) + y(k)∆x(k) + ∆y(k) x(k),
∆ |
y(k) |
= |
z(k)∆y(k)− y(k)∆z(k) |
. |
|
|
|||
z(k) |
z(k)z(k +1) |
Дифференцирование многочленов аналогично вычислению разностей факториальных многочленов [10]. Произвольный обыкновенный многочлен можно представить суммой факториальных многочленов. Факториальный многочлен m-го порядка определяется как
(k )(m) = k (k −1)(k − 2)...(k − m +1), |
(4.6) |
где m – положительное целое число. |
|
Согласно определению разностного оператора, имеем: |
|
∆(k )(m) = m(k )(m−1) = mk (k −1)(k − 2)...(k − m + 2). |
(4.7) |
4.1.2 Обратный разностный оператор
Найдем теперь обратный оператор, аналогичный интегральному оператору p−1 , при этом
t |
|
p−1 ( f (t))= ∫ f (t)dt + c = ∫ f (τ)dτ + K , |
(4.8) |
t0 |
|
где с и K – постоянные интегрирования.
Нижний предел t0 , вообще говоря, произвольный и определяется началом отсчета времени при анализе системы (обычно моментом поступления входного воздействия). Величина t0 формирует часть постоянной интегрирования, именно:
c = K − ∫ f (t)dt .
t=t0
Выражение y(t) = p−1 ( f (t)) является решением уравнения py(t) = f (t). Соответственно, выполняется соотношение pp−1 f (t) = f (t).
По аналогии обратный оператор ∆−1 должен иметь такой вид, чтобы выражение y(k) = ∆−1 f (k) являлось решением уравнения
∆y(k) = f (k) |
(4.9) |

|
|
156 |
|
или чтобы удовлетворялось равенство |
|
||
|
|
∆∆−1 f (k) = f (k). |
(4.10) |
|
Так как |
|
|
k−1 |
|
= ( f (k) + f (k −1) + ...+ f (0) + K)− ( f (k −1) |
+ ...+ f (0) + K ) = f (k), |
∆ ∑ f (n)+ K |
|||
n=0 |
|
|
|
то обратный оператор, удовлетворяющий уравнениям (4.9) и (4.10), имеет вид
|
k−1 |
|
∆−1 f (k) = ∑ f (n)+ K . |
(4.11) |
|
|
n=0 |
|
Уравнение (4.11), определяющее обратный оператор, можно переписать |
||
так: |
|
|
n=k−1 |
n=k |
|
∆−1 f (k) = ∑ f (n)+ c = ∑ f (n −1)+ c, |
(4.12) |
где суммирование производится по фиктивной переменной n.
Нижний предел в уравнении (4.12) не указан, так как можно объединить произвольное число членов f (0), f (1), f (2)... в уравнении (4.11) с постоянной суммирования K и образовать новую постоянную с. Таким образом, произвольный предел в уравнении (4.12) является аналогом постоянной t0 в уравнении (4.8), и его выбор определяется наиболее выгодным образом для каждого конкретного случая.
В отличие от интегралов, точное вычисление которых требует определенного искусства, а иногда и невозможно, для операторов ∆−1 таких сложностей не существует, однако вычисление ∆−1 f (k) непосредственно по формуле (4.12) довольно утомительно, особенно при больших k. Поэтому желательно выражать ∆−1 f (k) в замкнутом свернутом виде. Суммирование конечных рядов (как и вычисление интегралов) подчиняется определенным правилам. Существуют таблицы формул суммирования, правило суммирования по частям (аналог интегрированию по частям), используются многочлены Бернулли, разложение функций на простые дроби и т. д.
······················· Пример 4.1 ······················
k |
1 |
|
|
Найти сумму так называемого телескопического ряда ∑ |
. |
||
n(n −1) |
|||
n=2 |
|
157
Сумму такого ряда нетрудно найти, если представить выражение под знаком суммы в виде простых дробей:
|
|
|
|
|
|
|
|
1 |
|
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
− n . |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
n(n −1) |
n −1 |
|
|
|
|
|
|
|
|
|
|||||||||
Тогда результат суммирования будет: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
k |
1 |
k 1 |
|
1 |
|
|
|
1 |
1 |
|
1 |
|
1 |
|
1 |
|
|
1 |
|
|||||||
∑ |
|
= ∑ |
|
|
− |
|
|
= 1 |
− |
|
+ |
|
|
− |
|
|
+ ...+ |
|
|
− |
|
|
=1− |
|
. |
|
n(n −1) |
|
|
2 |
3 |
|
k |
k |
|||||||||||||||||||
n=2 |
n=2 n −1 |
|
n |
|
|
2 |
|
|
|
k −1 |
|
|
|
|
·······································································
Но не всегда конечные суммы можно свернуть и выразить в замкнутой форме. В некоторых случаях можно пользоваться верхней и нижней оценкой таких сумм.
При использовании факториальных многочленов из уравнения (4.7) можно получить с учетом формулы (4.12):
∆−1 (k)(m) = |
1 |
|
(k)(m+1) + K . |
|
m +1 |
||||
|
|
|||
Аналогами дифференциальных уравнений для дискретной переменной |
являются разностные уравнения или, как их еще называют, уравнения в конечных разностях.
4.2 Разностные линейные уравнения динамики
4.2.1 Общие свойства разностных уравнений
Общий вид разностного уравнения, связывающего выход y(k) с входом
r(k) системы с дискретным временем, следующий: |
|
a0 y(k + n)+ a1y(k + n −1)+ ...+ an y(k) = b0r(k + m)+ ...+ bmr(k). |
(4.13) |
Это же уравнение (4.13) можно представить в другом виде: |
|
(c0∆n + c1∆n−1 + ...+ cn )y(k ) = (d0∆m + d1∆m−1 + ...+ dm )r(k). |
(4.14) |
Переход от одной формы уравнения к другой очевиден, если иметь в виду соотношения (4.1)–(4.5). Уравнение (4.14) – более близкий аналог уравнению (3.5), а уравнение (4.13) легче решать, и поэтому оно более распространено.
Для линейных систем коэффициенты уравнений (4.13) и (4.14) не зависят от y или r, а для стационарных систем они независимы и от k, то есть являются постоянными величинами.
Применяя оператор сдвига E, уравнение (4.13) перепишем в виде
158 |
|
(a0En + a1En−1 + ...+ an ) y(k ) = (b0Em + b1Em−1 + ...+ bm )r(k ). |
(4.15) |
Поскольку вход r(k) считается известным, правую часть уравнения (4.13) (или (4.15)) можно обозначить как известную вынуждающую функцию F (k ) и записать уравнение (4.13) более компактно:
a0 y(k + n) + a1 y(k + n −1) + ...+ an y(k) = F (k). |
(4.16) |
Разностными уравнениями описываются системы, в которых процессы являются функциями дискретного переменного. Чаще всего эта дискретная переменная – время, но это может быть и положение, пространственные координаты, например в периодических структурах.
Уравнение (4.13), если a0 ≠ 0 и an ≠ 0 , является разностным уравнением n- го порядка. Если an = 0, а a0 ≠ 0 и an−1 ≠ 0 , то получим уравнение n −1-го порядка. То есть в отличие от дифференциального уравнения, порядок разностного уравнения определяется разностью высшей и низшей степеней Е. При использовании оператора ∆, например в уравнении (4.14), установить порядок уравнения непосредственно по его виду невозможно. Уравнение (4.15) называется неоднородным разностным уравнением, в отличие от однородного уравнения
(a0En + a1En−1 + ...+ an ) y(k) = 0 . |
(4.17) |
|||
Разностные уравнения, по сути, являются рекуррентными формулами. |
||||
Уравнение (4.16) можно решить относительно y(k + n) : |
|
|||
y(k + n) = |
1 |
(F(k)− a y(k + n −1) |
−...− a y(k)). |
(4.18) |
|
||||
1 |
n |
|
||
|
a0 |
|
|
При известных значениях y(0)÷ y(n −1) (начальные условия) y(k ) можно непосредственно найти для всех k ≥ n путем последовательного применения соотношения (4.18).
Таким образом, в отличие от дифференциального уравнения, y(k ) можно найти непосредственно по разностному уравнению для любых значений k. Но обычно не прибегают к итерационной процедуре, описываемой уравнением (4.18), а находят решение в замкнутой форме.
4.2.2 Решение однородных разностных уравнений
Однородное разностное уравнение n-го порядка содержит n линейно независимых решений. Обозначим n решений уравнения (4.17) при a0 ≠ 0 и an ≠ 0
|
|
|
159 |
|
|
|
|
|
через |
y1 (k), y2 (k),..., yn (k). Тогда условием (необходимым и достаточным) ли- |
|||||||
нейной независимости этих решений будет: |
|
|
|
|
|
|||
|
|
y1 |
y2 ... |
yn |
|
|
|
|
|
|
|
|
|
||||
|
C (k) = |
Ey1 |
Ey2 ... |
Eyn |
|
≠ 0. |
(4.19) |
|
|
|
. |
. |
. |
. |
|
|
|
|
|
En−1 y |
. |
.. |
En−1 y |
n |
|
|
|
|
1 |
|
|
|
|
|
|
|
Определитель C (k ) называется определителем Касорати. |
|
||||||
|
Поскольку уравнение (4.17) линейное, то его решением будет и линейная |
|||||||
комбинация независимых решений |
yi (k), то есть общее решение уравнения |
|||||||
(4.17) можно записать как |
|
|
|
|
|
|
|
|
|
yо (k) = c1y1 (k)+ c2 y2 (k)+ ...+ cn yn (k), |
(4.20) |
||||||
где ci |
– постоянные, не зависящие от k. |
|
|
|
|
|
||
|
Решение уравнения (4.17) можно по аналогии с дифференциальными |
|||||||
уравнениями искать в форме |
|
|
|
|
|
|
|
y(k) = euk ,
где u – неизвестная постоянная величина, подлежащая определению.
Но удобнее ввести обозначение z = eu и предполагаемое решение записать в виде
|
y(k) = zk . |
|
|
|
(4.21) |
|
Подставляя решение (4.21) в (4.17) и учитывая соотношение Enzk = znzk , |
||||||
получим характеристическое уравнение: |
|
|
|
|
||
a zn |
+ a zn−1 |
+ ...+ a |
n |
= 0. |
|
(4.22) |
0 |
1 |
|
|
|
|
|
При различных корнях характеристического уравнения z1, z2 ,..., zn |
общее |
|||||
решение уравнения (4.17) получит вид |
|
|
|
|
|
|
y (k) |
= c zk + c zk +...+ c zk . |
|
(4.23) |
|||
о |
1 1 |
2 2 |
|
n n |
|
|
Можно показать, что при различных zi |
отдельные решения y = zk |
удо- |
||||
|
|
|
|
i |
i |
|
влетворяют условию (4.19) и, следовательно, независимы. Если же, например, корень z1 имеет кратность m, то составляющая общего решения, соответствующая этому корню, равна
yо (k) = c1z1k + c2kz1k + ...+ cmkm−1z1k .

160
······················· Пример 4.2 ······················
Найти общее решение разностного уравнения
y(k + 2) + 0,3y(k +1) + 0,02y(k) = 0. Составим характеристическое уравнение
z2 + 0,3z + 0,02 = 0
и найдем его корни z1 = −0,2; z2 = −0,1. Осталось записать решение в форме (4.23):
yо (k) = c1 (−0,2)k + c2 (−0,1)k .
·······································································
Для любого комплексного корня уравнения (4.22) с действительными коэффициентами должен существовать и комплексно сопряженный корень. Решение разностного уравнения, соответствующее паре комплексно сопряжённых корней:
|
|
|
|
|
|
z |
= ρe± jθ , |
(4.24) |
||||||
|
|
|
1,2 |
|
|
|
|
|
|
|
|
|||
записывается в форме |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y (k) = ρk |
(Acosθk + Bsinθk) = Cρk cos(θk + ϕ) , |
(4.25) |
||||||||||||
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где A, B, C и φ – действительные постоянные, связанные друг с другом извест- |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
ными формулами приведения |
|
тригонометрических функций C = |
A2 + B2 , |
|||||||||||
ϕ = − arctan B A . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
······················· |
|
|
|
|
|
|
|
Пример 4.3 ······················ |
||||||
Решить уравнение y(k + 2)+ y(k +1)+ y(k) = 0 . |
|
|
||||||||||||
Корни характеристического уравнения z2 + z +1= 0 равны |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
π |
|
|
z = − 1 ± |
|
|
3 |
|
1 |
+ 3e± jarctan |
|
|
||||||
|
|
|
|
= e± j 3 . |
|
|
||||||||
|
|
j = |
3 |
|
|
|||||||||
|
|
|
|
|
|
|
||||||||
1,2 |
2 |
|
2 |
|
|
4 |
4 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
Таким образом, в выражении (4.24) ρ = 1, ϕ = π3 и решение согласно (4.25) равно y(k) = Acos π3k + Bsin π3k .
·······································································
Особое внимание нужно уделять нулевым корням характеристического уравнения (4.22). Если an = 0, а a0 ≠ 0 и an−1 ≠ 0 в уравнении (4.22), то характе-