
книги из ГПНТБ / Солодовников В.В. Статистическая динамика линейных систем автоматического управления
.pdf1 8 0 |
ОПРЕДЕЛЕНИЙ СТАТИСТИЧЕСКИХ ХАРАКТЕРИСТИК |
[ГЛ. IV |
В том случае, если аппроксимируемая кривая на значительном протяжении имеет наклон, не представляющий собой целое кратное 6 дб и лежащий в пределах от 6 ( к — 1) до 6/е на октаву, то этот
|
|
|
|
Рис. 4.22. |
|
|
ее участок |
можно |
заменить |
ломаной |
кривой (рис. 4.23), |
состоящей |
|
из сопрягающихся |
прямолинейных |
отрезков, из которых первый |
||||
имеет |
наклон в 6k |
децибел, |
второй |
в 6 ( k — 1) децибел, третий — |
||
опять |
в 6& |
децибел и т. д., |
причем |
точки сопряжения |
необходимо |
Рис. 4.23.
выбирать попеременно расположенными на 3 дб выше и ниже аппрок симируемой кривой.
Применение только что изложенного метода практически значи тельно проще, чем применение остальных методов, рассмотренных в этой главе (по крайней мере, если задана кривая спектральной плотности). Следует, однако, иметь в виду, что в случае необходи мости приближения логарифмических кривых спектральной плотности, имеющих резкие спады и подъемы, метод может привести к дробно
I ll ПРИМЕНЕНИЕ ДРОБНО-РАЦИОНАЛЬНЫХ ФУНКЦИЙ 181
рациональной функции, имеющей числитель и знаменатель очень высокого порядка.
В такого рода случаях для определения 5 (ш) целесообразно
использовать |
метод, |
применяемый для |
определения передаточных |
функций 1). |
|
|
|
П р и м е р. |
Предположим, что логарифмическая кривая спектральной |
||
плотности S (ш) |
имеет |
вид, Изображенный |
на рис. 4.24. Эта кривая может |
быть аппроксимирована пятью полубесконечными характеристиками со сле дующими сопрягающими частотами и крутизной:
“1 = |
0,03, |
0)2 = |
3,0, |
“3 = |
7,0, |
“ 1 = |
12,0, |
“?>= |
50,0, |
*1 = k o =
*3 =
л- II
*5 =
1.
1,
2, со
1.
причем начальная ордината равна 62 дб. |
|
|
выражение |
|||
Пользуясь изложенным выше |
методом, получим следующее |
|||||
для спектральной |
плотности: |
|
|
|
|
|
S (a ) = |
1,6- 10в |
|
МтЛИйЛСО \2 |
(4.84) |
||
|
и |
|
чы- |
з |
+ Ш У |
|
|
|
1 0,03) |
|
|
Истинные значения спектральной плотности для различных значений ш приведены в прилагаемой ниже таблице.
Таблица значений S (<о)
(1) |
S H |
О) |
S (« ) |
0 |
1,6 -100 |
8 |
5,89 |
0,1 |
13,25 • 104 |
10 |
4,175 |
0,5 |
5,57 • 103 |
20 |
0,136 |
1,0 |
211,75 |
30 |
0,045 |
4,0 |
32,67 |
40 |
0,0095 |
6,0 |
10,39 |
100 |
1,69- IQ-4 |
Подставив значения ш, приведенные в таблице, в (4.84), легко убедиться, что точность аппроксимации заданной кривой спектраль ной плотности при помощи выражения (4.84) имеет порядок 1%.
!) Основы автоматического регулирования, под ред. В. В. Солодовни-
кова, Машгиз, 1959, т. II, ч. 2, гл. XVI, § 3, стр. 111—121.
70
60
50
10
30
20
to
О
■10
-20
-30
-40
-50
Рис. 4.24.
12] ПРИМЕНЕНИЕ ФУНКЦИЙ ЛЯГЕРРА 183
12. Метод приближения кривых спектральной плотности1) и спектральный анализатор, основанные на применении функций Лягерра
Перейдем к изложению метода приближения, основанного на при менении функций Лягерра.
Этот метод приближения основан на допущении, что в качестве исходных данных используется не кривая спектральной плотности 5 (ш), полученная в результате обработки записи стационарного случайного
процесса |
x(t), |
а сама |
запись, |
например, в виде осциллограммы. |
||||||
Для пояснения идеи метода необходимо привести некоторые све |
||||||||||
дения о полиномах Лягерра. |
|
|
|
|
|
|
||||
Функцией Лягерра |
v-ro порядка |
называется2*) |
|
|
|
|
||||
|
|
1 |
|
|
|
|
|
|
|
|
|
V+ — |
|
|
|
|
|
|
|
|
|
|
1 |
2 / |
|
Г - ’ + |
. . . + |
|
|
|
|
|
U 0 = e - ‘ { |
|
|
|
|
|
|
||||
|
> - 1 ) ! |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
||
|
+ |
2'|~Х+Т |
у| (— 1)х |
. . . + 2 |
а (— l)v| . |
(4.85) |
||||
|
|
|
|
у.! [ ( v — к . ) ! ] 5 |
|
|
|
|
|
|
Преобразование Фурье Z.v(y'u>) для функции Лягерра имеет весьма |
||||||||||
простой вид |
|
|
|
|
|
|
|
|
|
|
|
L M |
= V 2 |
|
|
(v = 0 , 1 , 2 , . . . ) . |
|
(4.86) |
|||
Это легко показать следующим образом. |
|
|
|
|
||||||
Рассмотрим |
функцию от /, |
равную нулю при |
t < 0 |
и |
равную |
|||||
Ре~1 при |
^ > 0. Ее преобразование |
Фурье |
|
|
|
|
||||
|
|
|
00 |
|
|
|
|
|
|
|
|
|
|
J p |
e-‘e - l « d t = |
|
|
|
|
(4.87) |
|
|
|
|
О |
|
|
|
|
|
|
|
Поэтому, |
принимая |
во |
внимание (4.87), из (4.85) |
получим: |
|
|
||||
СО |
|
|
|
2Ч |
|
2v-iv |
|
|
|
|
J Д(Ое~/ш' dt |
|
|
|
|
• + |
|
|
|||
|
|
( 1 + » ,+1 |
( l + » v |
|
|
|
||||
о |
|
|
|
________ v!________ |
|
|
i z i i n |
|||
|
|
+ |
2“ х(— 1)х к.1 |
( v — |
к . ) ! ( 1 - ) - jut)'*-* |
+ |
• •• + |
0 + У « ) Г |
Но выражение в фигурных скобках представляет собой не что иное,
как |
|
|
|
[ 2 - ( ! + » ] » |
. ... |
( 1 - > Д |
|
|
|
|
|
|
|
|
|||
|
|
|
|
(1 _|_ушр+1 |
~ |
(1 + » ^ + 1 |
|
|
1) См. W i e n e r |
N., |
Extrapolation, |
interpolation and |
smoothing |
of statio |
|||
nary time series, |
J. Wiley, N. Y., 1949. |
|
|
|
stationary |
|||
2) W i e n e r |
N.. |
Extrapolation, interpolation and smoothing of |
||||||
time series, |
John |
Wiley, |
N. Y., 1949; Д ж е й м с X. и др„ |
Теория |
следящих |
|||
Систем, ИЛ, |
Москва, |
195J. |
|
|
|
|
184 |
ОПРЕДЕЛЕНИЕ СТАТИСТИЧЕСКИХ ХАРАКТЕРИСТИК |
[ГЛ. IV |
и, следовательно,
СО
О
что и требовалось доказать.
Применяя теорему о вычетах, легко показать, что
|
|
|
|
|
|
СО |
|
|
|
оо |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 4 |
/ ^ |
= |
1 |
|
|
(4.88) |
И |
при V |
ф |
|
СО |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
Г |
• |
* |
|
1 |
Г ( |
|
— t |
rfu) = 0. |
(4.89) |
||
|
|
2^ |
J |
A* C/w) L, (У’ш) (1ш — — |
|
|
|
|
|||||||
|
|
|
— СО |
|
|
|
|
£ ( l + > f +1“ 4 |
|
|
|||||
|
Имея |
в виду данное выше определение функций Лягерра, |
заметим, |
||||||||||||
что если |
некоторая |
функция л'(-с) удовлетворяет |
условиям |
|
|||||||||||
|
|
|
х (оо) = |
0, |
J |
[х (т)]2 d-z < |
оо, |
|
|
|
|
|
|||
|
|
|
|
|
|
|
о |
|
|
|
|
|
|
|
|
то |
ее можно |
разложить |
в |
ряд |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
ОО |
|
|
|
|
|
|
|
|
|
|
|
|
|
* ( х) = |
2 |
|
|
|
|
|
(4.90) |
|
где |
|
|
|
|
|
|
|
(1=0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
c* — f |
(х) V (х) rfx- |
|
|
|
(4.91) |
|||
|
|
|
|
|
|
|
|
о |
|
|
|
|
|
|
|
|
Следовательно, учитывая (4.90) и (4.86), |
можем |
написать: |
||||||||||||
°° |
|
|
|
со |
|
“ |
|
|
|
|
|
|
|
|
|
[ |
х (т) е ~ ^ |
dT = |
^ |
f |
l„ (т) |
dT = / 2 |
У] с, |
~r J®) |
• (4-92) |
||||||
о |
|
|
|
|
(i=o |
о |
|
|
|
|
ii =o |
|
|
||
|
Предположим теперь, что нам задан |
график |
некоторой |
стацио |
|||||||||||
нарной |
случайной функции х(() |
на интервале [0, 2Т]. Введем в рас |
|||||||||||||
смотрение |
функцию хг(/), |
определив |
ее следующим |
образом: |
|||||||||||
|
|
|
|
|
* г (0 = |
|
х (t) |
при |
0 < |
t < 2Г, |
|
|
|||
|
|
|
|
|
|
0 |
при |
остальных |
t. |
|
|
||||
|
|
|
|
|
|
|
|
|
|
||||||
|
Тогда |
ее |
преобразование Фурье будет |
иметь |
вид |
|
|
||||||||
|
|
|
|
|
|
|
|
|
ат |
|
|
|
|
|
|
|
|
|
|
|
|
Хт (Уш) = j |
x T{t) e~jwt dt, |
|
|
(4.93) |
12] |
ПРИМЕНЕНИЕ |
ФУНКЦИЙ ЛЯГЕРРА |
|
185 |
||
и в качестве |
спектральной плотности этой функции |
при достаточно |
||||
больших Т приближенно (см. § 12 гл. III) можно рассматривать вы |
||||||
ражение |
|
|
|
|
|
|
|
Sxt W = |
?t \X t U*) I2- |
|
(4.94) |
||
Разложим |
функцию х •/•(/) |
в |
ряд по |
функциям |
Лягерра, т. е. |
|
предположим, |
что |
|
N |
|
|
|
|
|
|
|
|
|
|
|
x T(t) = |
y i cJv(t). |
|
|
(4.95) |
|
|
|
|
»=о |
|
|
|
Так как функции Лягерра |
образуют |
полную |
ортогональную си |
|||
стему функций, обращающихся в нуль при У < 0 , |
то при достаточно |
большом А/ равенство (4.95) будет удовлетворяться с любой сте
пенью точности, причем коэффициенты |
ряда (4.95) |
должны быть |
вычислены по формуле1) |
|
|
27 |
|
|
с , = f |
■ |
(4.96) |
О |
|
|
Интеграл вида (4.96) может быть вычислен графически путем перемножения ординат кривых /v (У) и хт{1) и планиметрирования
площади, ограниченной |
получившейся результирующей кривой. Для |
|||||||
облегчения |
такого |
рода |
вычислений в |
приложении |
II дана таблица |
|||
функций Лягерра для значений v от 0 до 5. |
|
|
||||||
Если же коэффициенты cv вычислены, то представление спектраль |
||||||||
ной плотности в виде дробно-рациональной функции |
уже не пред |
|||||||
ставляет большого |
труда. |
|
|
|
|
|||
Действительно, |
взяв |
преобразование |
Фурье от |
обеих частей ра |
||||
венства |
(4.95), получим: |
|
N |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
^ |
г ( » = 2 с А ( » |
|
(4.97) |
||
следовательно, |
|
|
v = О |
|
|
|
||
|
|
|
|
|
|
|||
|
|
N |
|
|
N |
|
|
|
5. » = |
? f |
5 / vL .,( » |
2Т |
V |
|
|
(4.98) |
|
jLi сА 0 ’ш) |
|
|||||||
|
|
v =0 |
|
L v = 0 |
|
|
|
|
Так |
как |
каждая |
из функций. |
L4(jш) |
представляет |
собой дробно |
рациональную функцию вида (4.86), то и SXT(u>) будет представлять
собой дробно-рациональную |
функцию. Заметим, что первый множи |
||||
тель |
в правой части |
(4.98), |
заключенный в |
квадратные скобки, |
со |
держит все полюсы |
на мнимой оси в верхней полуплоскости, а вто |
||||
рой |
множитель — в нижней полуплоскости.' |
|
|
||
Изложенный метод приближения, как это было указано, обладает |
|||||
тем |
преимуществом, |
что |
он позволяет |
найти приближение |
для |
!) См., например, цитированную на стр. 1S3 книгу X. Джеймса и др.
1 8 6 ОПРЕДЕЛЕНИЕ СТАТИСТИЧЕСКИХ ХАРАКТЕРИСТИК [ГЛ. IV
функции спектральной плотности S (ш) непосредственно по записи соответствующей стационарной случайной функции без необходимости предварительной обработки этой записи. Но процесс вычисления коэффициентов cv является довольно кропотливым.
Коэффициенты cv могут быть определены не только при помощи графического вычисления интегралов (4.96), но и при помощи спе циального вычислительного устройства. Такого рода вычислитель
ное устройство, представляющее |
собой, |
по существу, спектральный |
||||||
анализатор, основанный на формулах (4.98), (4.91), |
должно состоять |
|||||||
из генератора |
функций Лягерра, умножителя и интегратора. |
|||||||
Генератор |
функций Лягерра |
можно |
выполнить в виде програм |
|||||
много механизма или в виде специального фильтра'). |
|
|
||||||
Если к фильтру приложено единичное напряжение, то токи, про |
||||||||
текающие по |
последовательным |
ветвям |
фильтра, |
будут |
изменяться |
|||
по функции Лягерра. |
|
|
|
|
|
|
||
Преимущество |
этого |
спектрального |
анализатора |
по |
сравнению |
|||
с рассмотренными- |
выше |
заключается в |
том, что, |
определив при |
помощи него коэффициенты 5 (ш), можно сразу написать аналити
ческое |
выражение |
для спектральной плотности, в |
то |
время |
как |
|
другие |
анализаторы |
позволяют получить лишь кривую спектральной |
||||
плотности, |
которую |
обычно нужно потом аппроксимировать каким- |
||||
либо известным аналитическим выражением. |
|
|
|
|||
J) См., |
например, |
в а н д е р П о л ь Б. и Б р е м м е р |
X., |
Операцион |
||
ное исчисление на основе двустороннего преобразования Лапласа, ИЛ, |
1952. |
ГЛАВА V
ПРОХОЖДЕНИЕ СЛУЧАЙНОГО СИГНАЛА ЧЕРЕЗ ДИНАМИЧЕСКУЮ СИСТЕМУ И ОПРЕДЕЛЕНИЕ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК СИСТЕМ РЕГУЛИРОВАНИЯ ПО ЭКСПЕРИМЕНТАЛЬНЫМ ДАННЫМ
1. Введение
Перейдем теперь к изложению способа применения теории слу
чайных процессов, |
вкратце |
изложенной |
в предыдущих параграфах, |
|
к анализу |
точности |
работы |
динамических систем, находящихся под |
|
влиянием |
стационарных случайных воздействий. |
|||
Нами |
уже указывалось, |
что точность |
работы системы в указан |
ных условиях не может быть полно охарактеризована при помощи обычного понятия качества, применяемого для характеристики дина мических систем, находящихся под влиянием воздействий, представ ляющих собой заданные функции времени, ввиду того, что в этом случае такие понятия, как, например, установившееся состояние, время переходного процесса и т. д., теряют смысл.
Если воздействия являются случайными функциями, то ошибка
воспроизведения е(£), т. е. разность |
|
e(t) = m(t) — x(t) |
(5.1) |
между полезным сигналом на входе m(t) и выходом x(t), также представляет собой случайную функцию. Поэтому в этих условиях речь может идти об определении не мгновенных, а лишь некоторых средних значений этой величины.
Таким средним значением обычно служит среднеквадратическое значение ошибки т . е. квадратный корень из величины
т
1
<5.2)
которой пользуются как критерием, определяющим точность или качество работы системы при наличии стационарных случайных воз действий.
188 ПРОХОЖДЕНИЕ СЛУЧАЙНОГО СИГНАЛА ЧЕРЕЗ ДНИАМИЧ. СИСТЕМУ [ГЛ. V
Оптимальными условиями работы системы согласно этому кри терию являются такие условия, при которых величина среднеквад ратической или стандартной ошибки имеет минимум.
В чем смысл критерия средиеквадратической ошибки |
(или сокра |
|||
щенно критерия СКО)? |
|
|
|
|
Согласно этому критерию нежелательность ошибки |
пропорцио |
|||
нальна квадрату |
ее величины вне |
зависимости от |
момента времени, |
|
в который имеет |
место ошибка. |
Вообще говоря, |
такого рода кри |
терий вполне логичен в тех случаях, когда нежелательность, ошибки
возрастает вместе с ее величиной. |
Однако возможны случаи, когда |
все достаточно большие ошибки |
в одинаковой мере нежелательны. |
В этих случаях следовало бы применять уже другой критерий,
учитывающий это |
обстоятельство. |
|
• Кроме того, критерий СКО обеспечивает |
малость лишь среднего, |
|
а не мгновенного |
значения ошибки, что также не всегда является |
|
достаточным. |
|
|
Наконец, применение критерия СКО может оказаться нерацио |
||
нальным в тех случаях, когда в некоторые |
моменты времени вели |
чина ошибки может быть сравнительно большой, а в некоторые моменты времени необходимо, чтобы она была, наоборот, мини мальной.
Таким образом, критерий СКО, как и всякий другой критерий точности, конечно, не является универсальным.
В предыдущей главе было показано, что среднее значение квад рата стационарной случайной переменной x ( t ) может быть весьма просто найдено, если известна корреляционная функция Rx (t) пли спектральная плотность Sx (w), соответствующая х (/).
Действительно,
ОО
(5.3)
Если ошибка воспроизведения e(t) является стационарной слу чайной функцией, то точно так же
ОО
(5.4)
где через /?£(-) и *$Е(ш) обозначены соответственно корреляционная функция и спектральная плотность ошибки s(t).
Таким образом, для того чтобы вычислить среднеквадратическую ошибку с, необходимо предварительно найти либо /?6(т), либо 5,(ш). Поэтому задача этой главы состоит в том, чтобы:
1) показать, каким образом может быть вычислена корреляционн функция или спектральная плотность величины на выходе системы,
2] |
СВЯЗЬ |
СТАТИСТИЧЕСКИХ ХАРАКТЕРИСТИК |
ВХОДА И ВЫХОДА |
189 |
|
если |
известна |
корреляционная функция |
или |
спектральная плотность |
|
сигнала на ее входе; |
|
|
|
||
2) |
изложить способ вычисления интеграла (5.4) в том случа |
||||
когда |
функция 5 6(ш) представляет собой |
дробно-рациональную функ |
цию от ш.
2. Связь между корреляционными функциями и спектральными плотностями величин на входе и на выходе линейной динамической системы *)
Рассмотрим линейную динамическую систему (рис. 5.1), имею щую импульсную переходную функцию k(t) и передаточную функ цию Ф (у'ш).
Предположим, что на вход этой системы подан стационарный
случайный |
сигнал m(t), |
имеющий |
корреляционную функцию |
/?т (т) |
||
и спектральную плотность S m(ш). |
|
|
|
|||
|
m l t l |
|
|
xiv |
|
|
|
|
|
Рис. |
5.1. |
|
|
Найдем |
корреляционную |
функцию |
R r(~) и спектральную |
плот- |
||
ость 5^.(ш) величины x ( t ) на выходе. |
|
|
||||
Имеем2): |
|
|
|
|
|
|
|
|
|
|
г |
|
|
|
RjA т) = |
lim |
o f |
f x ( t - \ - - ) x ( t ) d t . |
(5.5) |
|
|
|
Г > С О |
“ |
|
|
|
Но |
|
|
—T |
|
|
|
|
|
оо |
|
|
|
|
|
|
|
|
|
|
|
|
x ( l ) = |
J m |
— |
(л) d\. |
(5.6) |
|
|
|
|
—со |
|
|
|
- Подставляя (5.6) в (5.5) и, меняя порядок интегрирования, будем иметь:
оо оо ( Т 1
ЯжС О = lim |
J |
ей. |
J dr\ k(-r\)k(k){ |
j ' |
— т])m(t—X) dt |
|||
|
^ > 00 — СО |
|
— о о |
[ |
~ Т |
|
J |
|
|
|
|
|
|
|
|
|
(5.7) |
]) |
Д ж е й м.с |
X., |
Ни к о л ь е |
Н., |
Ф и л л и п с Р., |
Теория |
следящих |
|
систем, ИЛ, 1951. |
|
|
|
|
|
|
||
а) |
Здесь и |
везде |
далее знак |
предела при 7"-»-оо |
означает |
лишь, что |
при достаточно большом Т вероятность того, что погрешность определения R (т) интегралом (5.5) будет больше данной сколь угодно малой величины, —
бесконечно мала.