УДК 621.3.013.62
С.С. Зельманов МАТРИЧНЫЙ РЕЗОНАНС ФОРМЫ
Московский технический университет связи и информатики (Волго-Вятский филиал МТУСИ)
Показана возможность существования резонанса формы в линейной стационарной динамической системе с п входами и одним выходом при действии на её входах векторной совокупности сигналов со сложными конструкциями. Особенностью такого резонанса как свойства внешнего состояния системы является возможность его обнаружения при отсутствии признаков частотного резонанса в системе. Эта особенность весьма существенна при поиске в системе резонансных явлений, связанных с их возможными негативными последствиями.
Ключевые слова: резонанс формы, парциальные импульсные характеристики, матричный резонатор, векторная совокупность входных сигналов, внешнее состояние системы.
Постановка задачи
Динамическая система любой физической реализации и известной структуры может быть представлена как система с п входами и одним выходом. На её входы подаются внешние воздействия, а на выходе возникает экстремальная ответная реакция. При таком подходе основная задача — это поиск некоторой комбинации управляющих воздействий, обеспечивающих экстремальное реакцию системы на её выходе.
При этом резонанс рассматривается как свойство системы с внешним описанием. Это свойство альтернативно традиционному резонансу системы с внутренним описанием.
Введем обобщение, состоящее в том, что на входах такой резонансной системы могут действовать ряд сигналов с подчас сложными сигнальными конструкциями.
Результатом наших исследований будет доказательство того факта, что в зависимости от формы этих сигналов только при определенном их сочетании на выходе системы может появиться экстремальный сигнал, пиковое значение которого может существенно превышать значения выходного сигнала при любых других сигналах на входах, - это резонанс.
Практическая ценность такой задачи состоит в том, что реальные системы действительно могут содержать много точек внешнего воздействия. При действии на эти точки определенного сочетания внешних сил в системе может возникать резонанс, которого не ждут и который может привести к непредсказуемым последствиям. При этом иные сочетания автономных воздействий внешних сил на отдельные входы системы экстремального отклика на выходе системы не вызывают. При традиционном обследовании таких систем частотным методом резонансы в них не обнаруживаются. К таким системам могут относиться мостовые конструкции, отрезки железнодорожных путей с движущимися по ним составам и другие системы с ограниченной базой, находящиеся под воздействием внешних сил в разных точках.
Поэтому при вариации видов и параметров входных сигналов необходимо установить в ней явление резонанса и предупредить нежелательные его последствия. Такой вид резонанса будет резонансом формы, и мы определим его как матричный резонанс формы. При этом необходимо показать, что в случае векторной совокупности входных сигналов при определенной их параметризации на выходе системы может иметь место экстремальный отклик, величина которого в некоторый момент времени больше, чем для любой иной комбинации входных сигналов.
© Зельманов С.С., 2011.
Решение задачи
Теорема.
Для любой линейной стационарной динамической систем, имеющей п входов и один выход, существует сочетание (различных) сигналов с фиксированной энергией, подаваемых каждый на свой вход, только при одновременном возбуждении которыми п входов системы имеет место экстремальный отклик на выходе, т.е. резонанс в системе.
Доказательство.
Теорема констатирует, что резонанс в линейной системе при вариации формы единственного входного сигнала не является исключительным видом резонанса формы. Системный анализ явления резонанса предполагает дальнейшее развитие и совершенствование данного вида резонанса.
Рассмотрим явление резонанса в линейной стационарной динамической системе с постоянными параметрами, в которой имеется п входов и один выход. Пусть сигнал имеет п-мерное пространство состояний и вектор пространства состояний в момент ^ :
x(t
(1)
Тогда модель п х 1 системы в матричном виде может быть представлена уравнениями вида
^x(t) = Ax(t) + Bu(t),
y(t) = Cx(t) + Du(t); х(0) = х0, где А - n х n матрица переходов состояния, удовлетворяющих условиям устойчивости,
Ьц, bi2, ... bin
(2)
B
J21,
■>22,
■>2 n
матрица входов,
Ьп1, Ьп 2 ? ••• Ьпп С = [с, с2,...сп ] - матрица выхода.
Положим П = ё2,...ёп] = [0,0,...0], что соответствует отсутствию прямого прохождения сигнала с входа на выход для упрощения.
Известно, что (2) имеет решение в плоскости комплексной переменной:
Чр)
Х2(р)
Xn (p)_
да
где X ( p ) = J хг (t) e~ ptdt
U{p) =
Чр)
U2 (p)
Un ( p )
(3)
где и, ( р ) = | и (X) е-
о
при нулевых начальных условиях.
х(р) = (р1-А)~1-В-и(р)
где I - единичная матрица.
Для выхода будем иметь
00
¥(р) = $у(х)е-р'сИ = С(р1-Ау1-В-и(р) (4)
о
Матрица размера - 1 х п с коэффициентами, зависящими от р (комплексной переменной) С (р1 — А) 1 • В = К (р), имеет смысл матричной передаточной функции вида
К ( р ) = [ К ( р ) ,..Кп ( р )]. (5)
Легко проследить, что
и (р у
г (р )=[ к (р )...Кп (р)]
и2 (р)
ип ( р )
=£* (р )и1( р), (6)
,=1
где к (р) - передаточная функция с , входа на выход; К, (р) - представляют собой дробно-
рациональные функции многочленов, стоящих в числителе и знаменателе и зависящих от коэффициентов матриц А, В и С.
Полюса Ki (р) определяются соответствующими числами матрицы А.
Если обозначить соответствующие оригиналы
а+да
) = — / К,(р)ерф,
Л а—да
то можно считать, что этот набор функций определяет матричную импульсную характеристику системы, т.е.
О(X) = [К (X)...Ип (X)] (7)
а функции К, (X) считать парциальными импульсными характеристиками по , -входу.
Выражение (6) можно интерпретировать как декомпозицию произвольной линейной системы с постоянными параметрами, описываемую обыкновенными дифференциальными уравнениями и имеющую п - входов и один выход.
Такая система может быть представлена блок-схемой в следующем виде (рис. 1). Можно обобщить эту модель на все линейные стационарные системы, в том числе описываемые не только обыкновенными дифференциальными уравнениями, поскольку для всех них также можно представить связь п входов и одного выхода через парциальные передаточные функции, образующие матричную импульсную характеристику (7). Только в этом случае К (р) не будут дробно-рациональными функциями переменной р.
Рис. 1. Модель матричного резонатора общего вида
Покажем, что когда на входах системы действуют сигналы ^ (г (г).....^ (г), которые связаны с парциальными импульсными характеристиками соотношением
щ (г) = к(Т - г) \(т -г), I = 1,2...и, (8)
то на выходе будет иметь место экстремальный отклик, т.е. будет иметь место матричный резонанс в том смысле, который сформулирован в теореме о матричном резонансе.
Сначала рассмотрим (8) в частотной области. Для этого применим преобразование Лапласа и учтем зеркальность по времени функции в правой части. При вычислении получим следующее выражение через матричную передаточную функцию системы;
и (р) = К * (р)ерг = -С( р1 + А)-1 • ВерТ . (9)
Это решение может рассматриваться и как решение вариационной задачи на множестве всех функций входных сигналов при фиксированной энергии каждого из них. При этом вариация должна осуществляться как формой каждого сигнала, так и числом ненулевых сигналов, добиваясь уникального сочетания сигналов на п входах, обеспечивающего появление непредсказуемого резонансного значения сигнала на выходе системы. Выражение для сигнала на выходе системы при этом будет иметь вид
У(г) =[щ(г )щ(г).....ип (г)]<
'К (г У К(г)
К§ (г)
(10)
или
I =п г
У(г) = Е [ Щ (Т
1=1 п
(11)
где щ (г)щ (г).....ип (г) -входные сигналы; Ик (г) - импульсные характеристики. Приведенное
решение для многих систем может дать в качестве сигнала физически не реализуемую функцию (например, монотонно возрастающую до бесконечности). Далее мы приведем особенности решения с учетом требования финитности находимого решения.
Каждый из возможных входных сигналов щ (г) = 0 при г < 0. При образовании импульсной характеристики из такого сигнала с помощью его «зеркального» отображения и сдвига на величину Т эта характеристика к. (г) = 0 при г > Т. Эта оговорка не может заменить строго математического представления импульсной характеристики на временном ин-
тервале в бесконечных пределах. Такое представление может быть сделано с помощью «окна» вида [1(7) —— Г)], умноженного на функцию входного сигнала, т.е
«1 (X )•[ ^) — ^ — Г )].
Тогда в общем случае:
1. Функция сигнала и его изображение будут иметь вид
иМак-{р)=т -
В(р) (р — «1 )(р —«2 )""(р — «п )
(12)
где К * (р) - изображение сопряженного комплекса передаточной функции системы, которая имеет вид
К ( р ) = и ^ ( р )•
- рГ
2. Функция «окна» и его изображение:
(13)
(14)
3. Изображение произведения функций входного сигнала и «окна» с использованием вычетов в полюсах будет иметь вид
* Я+
Г (р) = Г [« (X)^(1(/) —1(/ — г))] = — | Г (г)• Г (р — г)±
3 Л—Л
— г) аг =
Л+
т!
А (г)
1 — е
-(р—г)Г
_ рМ
2к]¿Л1 + «1 )(г + «2+ «п) (р — г) ~ 2к] Лч(г)
1 Л+]т
= !
9тт i !
аг , (15)
где р( г ) = А( г )
1 — е
-( р—г )Г '
; ч (г )=(г+«1)(г+«2 )••••(г+«п)(р—г).
Корни знаменателя: г = —«, г2 = —«2, • • • • = —ап, гп+1 = р. Тогда изображение Лапласа для финитной функции входного сигнала:
и (р) = = ъх) + Яев^ = ъ2) +-----Яев^ = ъп) + Яев^ = ъп+х)
Р( г )
Ч( г )
+
Р( г )
г = г,
1 я (г)
+
р( г)
г =
Ч ( г )
+
р( г)
г = г
Ч (г)
(16)
г = г
п+1
В частном случае, входной сигнал «(г) при ограничении его длительности в интервале [о,Г] будет иметь вид Щ ~ ~ > примем= = е~ш и —-
р + а
/2 (0 = [1(0 -1(/ - Т)] □ 1 -±е-т =
1_ /_1 р р р
-1 Л+jco
f(p)-^ \;
1 1 - е
-(p-z )T
. dz -—f dz.
2nj j (z + a) (p - z) 2nj ÄJJ6} q (z)
p(z)-1 -е (p z)T;q(z)-(z + a)(p-z)- zp-z2 +ap-az, где корни z1 - -a, z2 - p .
Fi ( p ) - F [ fi (t) • Л (t)] - ^е^ ( z - zi) + Res ( z - z2) -
(18) (19)
_ P(z - zi ) P(z - z2 ) _ P(z)
q'(z - zi) q'(z - z2) q'(z)
+
p( z)
z - z,
i q'(z)
(20)
z - z„
q'(z)-p-2z-a; q'(z-z --a)-p + 2a-a-(p + a)
q'(z - z2 -p)-p-2p -a--(p + a)
=o
i - e-(p+a)T F ( p )-F [ fi (t )• f2 (t)]-Li-+
p + a
-(p + a)
p + a
(21)
где
F [ fi (t )• f2 (t )]-
i - e
-(p+a)T
р + а
Таким образом, для системы с заданной передаточной функцией резонансный финитный сигнал будет иметь вид а) в общем случае:
U (p)=P& U (p) q'(z)
+
P( z )
z - zi q'(z)
+
p( z)
z - zo
?'( z )
+
P( z)
z - z_
q'( z)
z - z
n+i
При этом будет выполняться К (р) = и \ (р) • е
pT ■
(22)
(23)
б) в частном случае, в качестве примера покажем, что системы с нерациональными передаточными функциями могут иметь простые резонансные сигналы.
Пусть
K ( p )-
е-aT - е-pT
p + a
(24)
Тогда резонансный сигнал щ (г) = е а. Покажем это:
1 - е~(р+а)Т и1 ( р )=.
р + а
Заменив р на — р, в силу сопряженности комплекса получим
(25)
u;( p )
i _ е-(p+a)T i _ е (p-a)T
p + a
K (p)-U * (p )
е-pt -
i - е
(p-a)T
-е-pt -
-p + a
pT r-aT r-aT pT
е ^ - е
p + a
p + a
е - е
p-a
(26) (27)
Покажем теперь, что найденное решение для сигнала, названного резонансным, имеет экстремальное свойство. Для любой парциальной системы можно записать
У (X) = к!« (т)н8 (X—т)ат,5 0 =0 при г> X=т^ у) (28)
о
где Ви7г (X) - взаимнокорреляционная функция сигнала и (X) и импульсной характеристики
любого канала системы, характеризующая их взаимную энергию; к > 1 - это коэффициент пропорциональности.
Коэффициент пропорциональности является фиксированной величиной для линейной системы и связан с определением импульсной характеристики, размерность которой 1/с. Когда мы говорим, что подаем на вход сигнал, равный зеркальной импульсной характеристике, то имеем в виду, что совпадают формы, а равенство существует с некоторым коэффициентом размерности (с). Физический смысл этого коэффициента обычно связывают с временем интегрирования в системе, т.е. с длительностью импульсной характеристики. Этот коэффициент будет присутствовать и в равенствах, содержащих ВКФ и АКФ
Максимальное значение сигнала на выходе линейной системы с одним входом и одним выходом при резонансе формы будет
7 T
y(T) = — \ü*(rn)drn = kEs >y(t) для t Ф T (29)
2л
т.е. в том случае, когда сигнал описывается функцией, являющейся зеркальной во времени функции импульсной характеристики парциального канала. Тогда парциальный выходной сигнал воспроизводит во времени сдвинутую АКФ входного сигнала, причем так, что максимум на выходе достигается в момент T. По определению АКФ сигнала конечной длительности имеет размерность энергии в точке T и равна энергии сигнала. Это определяет максимально возможное пиковое значение в канале.
Если входные сигналы окажутся такими, что для каждого парциального канала будут выполнены условия (29) для одного и того же момента времени Т, то на выходе образуется сигнал с пиковым значением, равным сумме пиковых значений каждого из парциальных каналов, пропорциональный сумме энергий входных сигналов:
I n T
y(T) = £Bmax = jj«I(t)hsdt > max[ukmax(T)], (30)
I=1 0
где индекс kmax относится к максимальному входному сигналу одного из каналов. В матричной форме выражение для выходного сигнала будет иметь вид
Y(p) = K* (p)epT • U(p) = -C(pI - A)~1 • BepTU(p) / (31)
Очевидно, что большего пикового значения, чем указанное d выражении (30) на выходе данной линейной системы при фиксированной энергии входных сигналов, получить невозможно.
Выводы
В линейной стационарной системе с n входами (при n>1) и одним выходом возможен резонанс формы сигнала, при котором значение выходного сигнала в момент отсчета будет больше или равно значению максимального из входных сигналов, т.е. матричный резонанс в линейной системе с многополюсным входом и одиночным выходом возможен.
0
1. Зельманов, С.С. Исследование явления резонанса формы сигнала в согласованном фильтре / // Электросвязь. 2011. №1.
Дата поступления в редакцию 11.10.2011
S.S. Zelmanov
FROM MATRIXER RESONANCE
The possibility of the resonance form existence in the liner stationary dynamic system with n inputs and one output by the action of the vector signal totality with the composite designs is presented in this article. The peculiarity of such resonance as the characteristic of the external system state is its reality by the absence of the frequency resonance signs in the system. This peculiarity is extremely sensational by the search in the system of the resonance phenomenon's which are connected with its possible negative consequences. Key words: resonance forms, partial impulse characteristics matrix resonator, vector total combination of input signals, external state of the system.
Key words: resonance forms, partial impulse characteristics matrix resonator, vector total combination of input signals, external state of the system.