GOS for Iphone / mobile / ТАУ / Книги / ТАУ учебник
.pdfПолученный результат показывает, что lim B |
= B , что и требова- |
||
|
T → |
q |
q |
лось доказать. |
0 |
|
|
|
|
|
|
3.3.Синтез импульсных систем управления
соднозначными кусочно-линейными характеристиками
произвольного вида
Вработах, посвященных решению задачи параметрического синтеза нелинейных непрерывных и импульсных САУ обобщенным методом Галеркина [89, 90, 206 и др. ], приводятся рекуррентные аналитические выражения, полученные для типовых нелинейных элементов, характеристики которых допускают кусочно-линейную аппроксимацию.
Следует отметить целесообразность такой аппроксимации, поскольку точный учет характеристик нелинейных звеньев существенно усложняет расчеты и вместе с тем не приводит к возрастанию точности получаемых результатов. Это обусловлено тем, что применяемый для решения задачи синтеза подход является приближенным в силу обстоятельств, изложенных ранее. При использовании кусочно-линейной аппроксимации (характеристики нелинейных звеньев заменяются прямолинейными отрезками) кусочно-линейные САУ рассматриваются как линейные модели с предысторией [209] и для получения надежных результатов применяется теория обобщенных функций [205].
Втех случаях, когда рассматривается нелинейная САУ с характеристикой нелинейного элемента отличной от типовой возможны следующие пути решения данной проблемы:
– получение рекуррентных аналитических выражений Bq, Bq*, определяющих интегралы Галеркина для рассматриваемой характеристики нелинейного элемента;
– получение рекуррентных аналитических выражений Bq, Bq*, определяющих интегралы Галеркина для характеристики произвольного вида, представленной набором произвольного (в общем случае n) числа отрезков прямых на плоскости; данный подход предложен и реализуется
Л.А. Осиповым и Т. Г. Поляковой;
– применение эквивалентного представления характеристик нели-
нейных звеньев набором типовых, для которых соответствующие соотношения получены ранее [187, 208].
В первом случае требуется выполнять достаточно сложные операции интегрирования и дифференцирования обобщенных функций, что само по себе является весьма трудоемкой задачей.
81
Во втором случае рекуррентное аналитическое выражение необходимо получать лишь один раз, поскольку характеристика нелинейного звена представлена в общем виде произвольным набором отрезков прямых на плоскости. Существенный недостаток такого подхода заключается в том, что задача из математической области переходит в алгоритмическую. Это обусловлено тем, что получить рекуррентное соотношение в общем виде для произвольной нелинейной характеристики и произвольного программного движения на ее входе (представленного суммой комплексных экспонент) не сложно. Вместе с тем составить алгоритм и написать соответствующую программу (так как метод ориентирован на использование вычислительной техники) весьма непросто из-за того, что в них необходимо предусматривать процедуры определения моментов переключения нелинейности, которые зависят как от вида характеристики, так и от вида процесса на ее входе (заранее неизвестны).
Поэтому представляется целесообразным применять на практике третий подход, сводящий решение данной задачи к эквивалентным преобразованиям нелинейности произвольного вида (эквивалентная замена нелинейности произвольного вида набором типовых характеристик).
Рассмотрим это положение на примере ряда нелинейных элементов, которые соответствуют характеристикам устройств, широко применяемых в непрерывных и импульсных системах управления [197]. На рис. 3.1, 3.2 представлены эквивалентные преобразования некоторых нелинейных характеристик: характеристики золотников гидро- и пневмосервомоторов (рис. 3.1, а); зависимости коэффициента сцепления от величины относительного проскальзывания (рис. 3.1, б); статической многоступенчатой характеристики квантователя-экстраполятора (рис. 3.1, в); нелинейности типа холостой ход (рис. 3.2, а); нелинейности типа начальный скачок и ограничение (рис. 3.2, б); нелинейности типа отрицательный дефект (рис. 3.2, в). Как видно из рис. 3.1, 3.2, исходные характеристики структурно преобразованы с использованием параллельного включения комбинаций типовых нелинейных и линейных элементов. При таком структурном преобразовании процесс на входе исходного нелинейного звена совпадает с процессом на входах нескольких параллельно включенных элементов, т. е. не возникает необходимости пересчета процесса на входы типовых нелинейных звеньев, что представляет собой достаточно сложную задачу. Так, например пересчет процесса, неизбежен при последовательном включении нелинейных звеньев или
82
F(x)
c
arctgk
–b2 –b1 |
x |
x(t) |
0 b1 b2
|
|
|
|
|
|
|
|
|
|
F(x) |
|
|
|||||||
|
|
|
|
|
c |
|
|
|
|
|
|
|
|
arctgk2 |
|
|
|||
–b |
|
|
|
|
|
|
|
|
x |
x(t) |
|
||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
b |
|
|
|||||||
|
|
|
|||||||||||||||||
|
|
|
|
||||||||||||||||
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
arctgk1 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
–c |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
F(x)
–b/2 |
c |
|
x |
|
b/2 |
||
|
|
|
|
|
|
–c |
|
x(t)
F1
–b1 x
0 b1 arctgk
bF2
2x
–b2 0 arctgk
F1 c
–b
x
0b arctgk1
F2
–b
x
0 b
arctgk2
|
F1 |
|
c |
x |
|
–b/2 |
||
b/2 |
||
|
||
–c |
|
|
–3b/2 c |
F2 |
|
x |
||
|
3b/2 |
|
|
–c |
c Fg
x
–c
Рис. 3.1
F1(x)
F[x(t)]
F2(x)
F1(x)
F[x(t)]
F2(x)
F1(x)
F2(x) F[x(t)]
83
охвате нелинейного звена обратной связью другим нелинейным звеном.
Обобщая показанные на рис. 3.1, 3.2 эквивалентные преобразования, можно записать следующее аналитическое выражение [187]:
l
F x0 (t) = ∑ Fg x0 (t) , (3.45)
g =1
либо при наличии импульсного процесса на входе нелинейного элемента
|
|
0 |
|
= ∑ |
l |
|
|
0 |
|
|
x |
|
x |
(3.46) |
|||||||
F |
|
(t) |
|
Fg |
|
(t) , |
||||
|
|
|
|
g =1 |
|
|
|
|
|
|
где Fg[x0(t)], Fg[x0*(t)] – типовые нелинейные характеристики: переменный коэффициент усиления; зона нечувствительности; насыщение (ограничение); релейная характеристика, в том числе с зоной нечувствительности.
С учетом (3.45), (3.46) интегралы Bqi, Bqi* принимают вид
|
|
|
∞ |
|
|
|
l |
|
|
|
|
|
|
−ρqt |
|
|
|
l |
|
|
|
|
B |
= |
∫ |
D |
i |
∑ |
|
F x |
0 |
(t) |
|
|
e |
dt = |
|
B , |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
qi |
|
|
|
|
g |
|
|
|
|
|
|
|
|
∑ |
|
qig |
|
|
|||
|
|
|
0 |
|
g=1 |
|
|
|
|
|
|
|
|
|
|
g=1 |
|
|
||||
|
i = 0,1,…,u, |
q = |
|
1,2,… m , g, = |
|
|
1,2,… l |
, |
|
|||||||||||||
|
|
|
∞ |
|
|
|
l |
|
0 |
|
|
|
|
|
−ρqt |
|
|
l |
|
|
|
|
B |
= |
∫ |
D |
i |
∑ |
|
F x |
(t) |
|
e |
dt |
= |
|
|
|
|||||||
|
|
|
|
|
|
|
|
B , |
|
|
||||||||||||
qi |
|
|
|
|
g |
|
|
|
|
|
|
|
|
|
∑ |
|
qig |
|
|
|||
|
|
|
0 |
|
g=1 |
|
|
|
|
|
|
|
|
|
|
g=1 |
|
|
||||
i = 0,1,…,u , |
q = |
|
1,2,… m , g, = |
|
1,2,… l |
, |
(3.47) |
|||||||||||||||
где Bqig, Bqig* – рекуррентные аналитические выражения, определяющие интегралы Галеркина типовых нелинейных элементов, представленные в таблицах приложения; g – число типовых нелинейных звеньев, включенных параллельно друг другу.
Таким образом, соотношения вида (3.45)–(3.47), представляющие собой аналитическое выражение принципа эквивалентных преобразований нелинейных характеристик, позволяют использовать обобщенный метод Галеркина к системам, содержащим нелинейные зве-
84
|
|
F(x) |
||
–b c |
|
|
|
arctgk |
|
|
|
x |
|
|
|
|
||
0–cb
|
|
|
F(x) |
|||
–b |
c |
|
|
B |
||
|
|
|||||
|
|
|||||
|
|
|||||
|
|
|||||
|
|
|||||
|
||||||
|
||||||
|
|
|
x |
|||
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
b |
|||
|
|
|||||
|
|
|||||
|
|
|||||
|
|
|
–c |
|||
|
|
|
||||
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
arctgk |
|
|
|
|
||
|
|
|
|
|
|
|
F(x)
c arctgk
x
0
–c
F1
–b
x
0 b
c
x(t)
F2
–b
x
0 b
arctgk
|
|
c |
F1 |
|
|
|
|
|
|
|
|
|
|
x |
|
0 |
|||
|
- c |
|||
x(t) |
|
|
|
|
|
|
|
F2 |
|
|
|
|
||
|
|
|
|
B1 |
|
|
|
|
x |
|
|
|
||
|
|
|
b |
|
|
|
arctgk |
B1 = B–c |
|
|
|
c |
F1 |
|
|
|
|
|
|
|
|
|
|
x |
|
0 |
–c |
||
x(t) |
|
|
|
|
|
|
arctgk |
x1 |
|
|
|
|||
|
|
|
|
|
|
|
|
0 |
x |
|
|
|
||
|
|
|
|
|
Рис. 3.2 |
|
|
||
F1(x)
F[x(t)]
F2(x)
F1(x)
F[x(t)]
F2(x)
F1(x)
F[x(t)]
x1
нья с характеристикой произвольного вида, допускающей кусочнолинейную аппроксимацию.
85
3.4.Синтез импульсных САУ
сучетом конечной длительности замыкания импульсного элемента
Проектирование реальных систем с амплитудно-импульсным модулятором часто приводит к необходимости учета влияния длительности замыкания ключа на динамику синтезируемой системы. Поэтому ниже рассматриваются некоторые особенности метода ортогональных проекций применительно к параметрическому синтезу реальных САУ с АИМ. Необходимо отметить, что данные особенности связаны, прежде всего, с получением рекуррентных соотношений, определяющих соответствующие интегралы вида (3.6), поскольку, изложенная общая схема решения поставленной задачи остается неизменной для САУ любого класса из числа детерминированных.
В соответствии с терминологией, принятой во второй главе, исследуемая координата, относительно которой записано уравнение движения САУ (в случае нелинейных систем – координата входа нелинейного элемента) и входное воздействие на выходе импульсного элемента типа I описываются следующим образом:
|
∞ |
|
|
|
|
|
|
||||
x (t ) = ∑ |
; |
||||
x (nT ) 1(t − nT ) −1(t − (nT − τ)) |
|||||
|
n=0 |
|
|
|
|
|
∞ |
|
|
|
|
|
|
||||
f (t ) = ∑ |
, |
||||
f (nT ) |
1(t − nT ) −1(t − (nT − τ)) |
||||
|
n=0 |
|
|
|
|
либо в случае импульсного элемента типа II
x (t ) = x (t )∑∞ |
1(t − nT ) −1(t − (nT − τ)) ; |
n=0 |
|
f (t ) = f (t )∑∞ |
1(t − nT ) −1(t − (nT − τ)) , |
n=0
(3.48)
(3.49)
(3.50)
(3.51)
где τ = γT – ширина импульса на выходе модулятора, здесь γ < 1. Рассмотрим получение аналитических выражений, определяющих
интегралы Aqi*, Bqi*, Cqi* для импульсных элементов типа I и II. Следует отметить, что данные интегралы вычисляются в соответствии с пра-
86
вилами действия над обобщенными функциями [205] и функциональными рядами [210]. Воспользуемся следующими соотношениями:
b |
∞ |
|
∑ |
∞ b |
|
|
|
|
|
||
∫ |
|
∑ |
k |
|
∫ |
|
k |
(x )dx, |
[a ≤ x ≤b]; |
|
|
|
|
f |
|
|
f |
|
(3.52) |
||||
|
|
|
(x ) dx = |
|
|
||||||
a |
k =0 |
|
|
k =0 a |
|
|
|
|
|
||
D |
|
∞ |
f |
|
|
= |
∞ |
|
(x )dx}, |
|
∑ |
k |
(x )dx |
D{f |
k |
||||
|
|
|
|
∑ |
|
|
|||
|
|
k =0 |
|
|
|
|
k =0 |
|
|
|
|
|
|
|
|
|
|
∞
∫ f (t )δ(k ) (t − τ)dt = (−1)k f (k ) (τ),
0
k = 0,1,… ; |
(3.53) |
k = 0,1,… , |
(3.54) |
последнее из которых описывает фильтрующее свойство δ -функции, существующей в момент времени t = τ .
Если на вход синтезируемой системы действует скачкообразное воздействие амплитудой H, то на выходе импульсного элемента как типа I, так и типа II, в соответствии с (3.49), (3.51), имеем процесс вида
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
f |
(t ) = ∑ |
|
|
|
|
|
|
|
|
(3.55) |
||
|
|
|
|
|
|
|
|
||||||
|
|
H 1(t − nT ) − 1(t − (n + γ)T ) . |
|
||||||||||
|
|
|
|
n=0 |
|
|
|
|
|
|
|
|
|
Тогда с учетом (3.55) выражение C |
* принимает вид |
|
|||||||||||
|
|
|
|
|
|
|
qi |
|
|
|
|
|
|
|
∞ |
|
|
∞ |
|
|
|
|
|
|
−ρ t |
|
|
Cqi = |
∫ |
|
|
∑ |
H |
1(t |
− nT ) − 1(t |
− (n + γ)T ) |
|
e |
dt. |
||
Di |
|
|
q |
||||||||||
|
|
|
n=0 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
Далее, используя соотношения (3.52) – (3.54), получаем следующее при i = 0:
Cq0
= H
|
∞ |
|
|
∞ |
|
|
|
|
|
−ρqt |
|
= |
∫ |
D |
i |
∑ |
|
H |
1(t − nT ) − 1(t − (n + γ)T ) |
|
e |
dt = |
|
|
|
|
|
||||||||
|
|
|
n= |
0 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
∞ |
|
|
|
|
∞ |
|
|
|
|
|
|
|
|||
|
|
|
|
∫ |
|
−ρqt |
|
|
|
∫ |
|
−ρqt |
|
|
|
|
|
|||
|
|
|
= H ∑ |
|
e |
|
dt − |
|
e |
|
dt |
= |
|
|
||||||
|
|
|
n=0 nT |
|
|
|
|
(n+γ)T |
|
|
|
|
|
|
|
|||||
∞ |
|
|
−ρqnT |
|
|
−ρq (n+γ)T |
|
H (1 − e−ρqγT )∞ |
|
|
||||||||||
∑ |
|
e |
|
|
− |
e |
|
|
|
|
|
= |
|
|
|
|
∑ |
e−ρqnT . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.56) |
|||
n=0 |
|
ρq |
|
|
|
ρq |
|
|
|
|
|
ρq |
|
|
n=0 |
|
|
|||
87
Сомножитель, стоящий под знаком суммы, представляет собой бесконечно убывающую геометрическую прогрессию со знаменателем q = e–ρqT, сумма членов которой может быть определена по формуле
∑∞ |
aqn = |
|
a |
|
, |
|
при |
|
q |
|
< 1, |
(3.57) |
|||||||
|
|
||||||||||||||||||
|
|
|
|
|
|||||||||||||||
n=0 |
|
1 |
− q |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
т. е. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
−ρqnT |
|
|
|
1 |
|
|
|
|
|
|
|
|
||||
|
∑ |
e |
= |
|
|
|
|
. |
|
(3.58) |
|||||||||
|
|
|
|
|
|
− e |
−ρqT |
|
|||||||||||
|
n=0 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
С учетом (3.58) окончательно получаем |
|
|
|||||||||||||||||
|
|
= |
H (1 − e−ρqγT ) |
|
− 1 |
|
|
||||||||||||
|
Cq0 |
(1 − e−ρqT ) |
|
ρq |
; |
|
|||||||||||||
при i = 1, рассуждая аналогично, получаем следующее:
∞ |
δ (t − nT ) − δ (t − (n + γ)T ) e−ρqt dt = |
||
Cq1 = ∫ H ∑∞ |
|||
0 n=0 |
|
H (1 − e−ρqγT ) |
|
|
|
||
|
= |
(1 − e−ρqT ) |
ρ0q ; |
………………………………………………………………………
при i = v* имеем
∞ |
∞ |
|
v −1 |
v |
−1 |
|
|
Cqv = ∫ |
H ∑ |
δ( |
) (t − nT ) − δ( |
|
) (t − (n + γ)T ) e−ρqt dt = |
||
0 |
n=0 |
|
|
H (1 − e−ρqγT ) |
|
|
|
|
|
|
|
|
−1 |
||
|
|
|
|
|
|
v |
|
|
|
|
= |
|
|
ρ(q |
). |
|
|
|
(1 − e−ρqT ) |
|
|||
88
Обобщая полученные результаты, окончательно получаем следующее:
|
|
∞ |
|
|
∞ |
|
|
|
|
|
−ρ t |
|
|
|
Cqi |
= |
∫ |
|
H ∑ |
|
1(t − nT ) − (t − (n + γ)T ) |
|
e |
dt = Cqρ−iq |
1, |
|
|||
Di |
|
q |
(3.59) |
|||||||||||
|
|
|
|
n=0 |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
где C |
= |
H (1 − e−ρqγT ) |
. |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
||||||||
q |
|
|
1 − e |
−ρqT |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Теперь перейдем к вычислению интегралов A*qi. В случае импульсного элемента типа I процесс на выходе модулятора (при сигнале вида (3.2) на входе) будет иметь вид
|
∞ |
|
|
|
|
− |
αnT |
|
|
|
|
|
1(t− |
|
|
|
||
|
|
|
|
|
|
|
|
|
||||||||||
x (t ) = ∑ |
|
+ H e |
|
nT )− |
(n+ |
. (3.60) |
||||||||||||
xу |
|
cos (βnT − ϕ 0 ) |
1(t− |
γ)T ) |
||||||||||||||
|
n=0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Тогда, с учетом (3.60) выражение A*qi принимает вид |
|
|
||||||||||||||||
|
|
|
|
|
∞ |
|
|
∞ |
− αnT |
|
|
|
|
|
|
|
|
|
|
|
= |
∫ |
D |
i |
∑ |
cos (βnT − ϕ 0 ) |
|
|
|
||||||||
|
|
Aqi |
|
xу + H e |
|
|
||||||||||||
|
|
|
|
|
|
|
n=0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−ρqt |
|
|
|
|
||
|
|
|
|
|
|
1(t− |
nT )− 1 (t− (n+ |
γ)T ) |
|
e |
dt. |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Далее в результате выкладок, аналогичных вычислению данного интеграла в случае идеального импульсного элемента [90, 215, 231], получаем рекуррентное аналитическое соотношение следующего вида:
|
∞ |
|
|
∞ |
− αnT |
|
|
|
|
|
∫ |
D |
i |
∑ |
cos(βnT − ϕ 0 ) |
|
1(t − nT ) − 1(t − (n + γ)T ) |
|
|
||
Aqi = |
|
xу + H e |
|
|||||||
|
|
|
n=0 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
−ρqt |
|
1− e−ρqγT |
|
i |
|
|
|
e |
|
dt= |
|
A ρ |
q |
, |
i= |
0,1,…,n ; q= 1,2,…,m, |
|
|
|||||||
|
|
|
ρq |
q |
|
|
|
|
|
|
|
|
|
|
|
|
где Aq*– было определено ранее (3.8).
В случае импульсного элемента типа II имеем
|
|
|
|
− αt |
∞ |
1(t− nT )− 1(t− |
|
|
(t ) = |
cos (βt − ϕ 0 ) ∑ |
(n+ |
||||
x |
|
xу + H |
e |
||||
|
|
|
|
|
n=0 |
|
|
(3.61)
γ)T ) . (3.62)
89
С учетом (3.62) выражение A |
* принимает вид |
|
|||||
|
|
|
|
|
qi |
|
|
|
∞ |
|
|
|
− αt |
|
|
∫ |
D |
i |
cos (βt − ϕ 0 ) |
|
|||
Aqi = |
|
xу + H |
e |
||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
−ρqt |
|
∑ |
|
|
1(t − nT ) − 1(t − (n + γ)T ) |
|
|
e |
dt, |
|
|
|
|
||||||
n= |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
либо
|
∞ |
|
|
∞ |
|
|
|
|
−ρqt |
|
∫ |
D |
i |
xу ∑ |
|
1(t − nT ) − 1(t − (n + γ)T ) |
|
e |
dt + |
||
Aqi = |
|
|
|
|||||||
|
|
|
n=0 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
∞ |
|
∞ |
|
|
|
|
|
|
−ρqt |
|
|
+ |
∫ |
|
H e− αt cos(βt − ϕ 0 )∑ |
|
1(t− nT )− 1(t− |
(n+ |
γ)T ) |
|
e |
dt. |
(3.63) |
|
Di |
|
|
||||||||||
|
|
n=0 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
Следует отметить, что первый интеграл формулы (3.63) аналогичен соотношению (3.59), поэтому рекуррентное выражение для его определения будет следующим:
∞ |
|
|
∞ |
|
|
|
|
|
|
|
|
−ρqt |
|
∫ |
D |
i |
xу ∑ |
|
1(t − nT ) − (t − (n + |
γ)T ) |
|
e |
dt = |
||||
|
|
|
|||||||||||
|
|
n=0 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
xу (1 − e−ρqγT ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
ρiq−1. |
|
|
|
(3.64) |
||
|
|
|
|
|
1 − e−ρqT |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||
Для вычисления второго интеграла соотношения (3.63) необходимо определить обобщенную производную i-го порядка от непрерывной на интервале от (t–nT) до (t–(n+γ )T) функции x(t). Обобщенная производная i-го порядка от данной функции определяется следующим выражением:
Di {x (t )1(t − nT )}= x(i) (t )1(t − nT )+ x(i−1) (nT )δ(t − nT )+ … |
|
+x (nT )δ(i−1) (t − nT ). |
(3.65) |
90
