Файл: Пугачев, В. С. Основы статистической теории автоматических систем.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 190

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

K„(t,t)=DyUi

 

 

 

 

2

 

 

0 , ю

 

 

 

 

 

 

6„ = 0,05

 

 

Л

А

4

 

 

1

\ г

 

4

8

12

16

t,c

Рис. 5.5. Дисперсия выходной координаты

 

представлены графики

зависимости

дисперсии

выходной коорди­

наты от времени при

различных

значениях

интенсивности по­

мехи.

Анализ результатов вычисления показывает, что в данной системе после включения имеет место переходный процесс. Время переход­ ного процесса дисперсии и корреляционной функции существенно зависит от интенсивности помехи. При отсутствии возмущения (GN = = 0) переходный процесс по дисперсии затухает, а весовая функция

g n

0) описывает гармонические колебания частоты со = 1/7\

что

соответствует периоду колебаний Тп = 3,14 с; это непосред­

ственно следует из рис. 5.4. Таким образом, корреляционная функ­ ция К 11 {t, 0) при отсутствии случайного возмущения характеризует процесс развития автоколебаний в данной системе.

При действии случайного возмущения весовая функция g ^ (t, 0) затухает до нуля, причем время затухания тем меньше, чем больше интенсивность возмущения. На первый взгляд кажется странным что весовая функция статистически линеаризованной автоколебательной системы затухает, поскольку она является частным интегралом си­ стемы (5.51), не содержащей входного возмущения. Объяснение данному факту можно дать, если учесть, что коэффициент статисти­ ческой линеаризации зависит от математического ожидания и дис­ персии выходной координаты. Изменение этого коэффициента во времени в соответствии с изменением дисперсии выходной перемен­ ной (при расчете принято mN = 0) приводит к демпфированию коле­ баний в системе, и весовая функция g ^ (t, 0) стремится к нулю при t —>оо. Демпфирование тем сильнее, чем больше интенсивность воз­ мущения, поскольку это увеличивает дисперсию 0ц.

Рассмотрим далее установившийся режим работы системы для случая, когда математическое ожидание входного возмущения равно нулю. Определим дисперсию и спектральную плотность выходной координаты.

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

к

(5.52)

Ф(«) = Т0т*& + 2Д7У- + T0s + kkx

151


где коэффициент статистической линеаризации k 1 в

соответствии

со второй формулой (5.47) при /?гх = 0 имеет вид

 

кх = 21/]^2поу.

(5.53)

Вычислим дисперсию выходной переменной, воспользовавшись формулой

СО

 

Оу = J | Ф (гео) |2 SN (со) dco.

(5.54)

Подставим в эту формулу значение частотной

характеристики

из выражения (5.52) и постоянную спектральную плотность белого шума SN — GN/2n. Воспользовавшись значениями интегралов от дробно-рациональных функций, приведенных в приложении 2, получаем

kG,\r

(5.55)

2T0kl Tkk\

 

Решая уравнения (5.53), (5.55) относительно среднего квадрати­ ческого отклонения и коэффициента статистической линеаризации, находим

_

kTl (1 - г ц) .

, __

2Т0

(5.56)

а

T0 V Ш

1

^ ( 1 + Ю

 

где параметр ц характеризует интенсивность входного случайного возмущения,

р = nGN/2l2T.

.

(5.57)

Спектральная плотность выходной переменной

Sy (со) = |Ф (ш) |2Sn ( с о ) .

(5.58)

Подставляя в формулу для частотной характеристики значение коэффициента статистической линеаризации (5.56) и вычисляя квад­ рат модуля, получаем следующее выражение для спектральной плот­ ности выходного сигнала:

Sy (со) =

кЧ°-Т3

 

(5.59)

2Г 2

H -7

 

Я Го (соГ)6 + 2 (соГ)4 +

(соТ)г +

 

 

р+1

(1 +Юа

Исследуем зависимость спектральной плотности от частоты и па­ раметров системы. При со = 0 спектральная плотность

Sy (0) = /e W p (1 + p)2/4jx27l.

При со = оо Sy {оо) = 0. Спектральная плотность имеет мак­ симум на частоте сот , определяемой из уравнения (d,Sy/d(d)am = 0.

Выполняя дифференцирование и решая уравнение, получаем сле­ дующую формулу для частоты точки экстремума:

соm

+1

Ч

^

 

2

1

- i / 2 5 + ц

при р,<;7,

 

 

 

У

 

L

з

3

1 + ц .

 

(5.60)

 

0

 

 

 

при р, >

7.

152


Максимальная величина спектральной плотности

S y (®ш) ---

 

 

27fe27'3/2p

 

(5.61)

71 — р

_

р \ 7

54

2л27’п

1+ Р

/ 25 + ц У

+ (1 + Р )2

 

 

Vн - р )

 

При

увеличении интенсивности помехи (р —>оо) максимум ча­

стотной

характеристики смещается в область более

низких частот,

а величина максимума увеличивается. При р > 7

максимум спек­

тральной плотности находится в нуле.

При уменьшении интенсивности помехи (р —>0) максимум спек­ тральной плотности смещается вправо и уменьшается до некоторого значения, а затем вновь увеличивается. При отсутствии случайного возмущения (р = 0) спектральная плотность превращается в 6-функ- цню в точках со = ±1/7\ В этом можно убедиться, анализируя выражения (5.59), (5.60), (5.61) при р —»0. Вычисляя неопределен­ ность в выражении (5.59) при р —>0, получаем

limS„(w)

6 W 2 б ( с о - ^ ) + б ( с о + ^ )

(5.62)

Ц->0

4лТ1

 

График спектральной плотности в относительных единицах при различных значениях параметра р в области положительных частот представлен на рис. 5.6. При интенсивности помехи = 0,1 и ука­ занных выше значениях параметров имеем р = 0,323. По формуле (5.60) резонансная частота а тТ = ±0,891. Таким образом, сдвиг частоты по отношению к частоте собственных автоколебаний, имею­ щим место при отсутствии шума, составляет А = +0,109.

153


Вырождение спектральной плотности в

б-функцию на

частоте

со = 1/Г означает,

что в системе при отсутствии случайных

возму­

щений существуют

автоколебания

с частотой

со =

1/Г. Чтобы

убе­

диться в этом, определим частоту

и амплитуду автоколебаний

си­

стемы с помощью гармонической линеаризации. Линеаризуя

нели­

нейность, получим следующую частотную характеристику:

 

 

сю (('со) = Г(1-р ((м):1+ 2Т0Т (ко)'- 4- Г '

(

(5.63)

где /ег — коэффициент

гармонической линеаризации,

 

 

 

 

/гг = 4Ипа.

 

 

(5.64)

В этом выражении

а — амплитуда гармоники.

В соответствии

с методом гармонической линеаризации приравниваем к нулю дей­ ствительную и мнимую части характеристического уравнения и полу­ чаем два уравнения относительно частоты и амплитуды автоколеба­ ний:

— Т 0Г2со3 + Г 0со = 0;

 

— 2Г0Гсо2 + Akllna = 0.

(5.65)

Решая эти уравнения, получаем

 

о = 1/7; а = 2klTlnT0.

(5.66)

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

гармонической линеаризации, получаем

 

kT = 2T jk T .

(5.67)

Сравнивая это выражение с выражением для коэффициента ста­ тистической линеаризации при р = 0, убеждаемся, что эти коэф­ фициенты совпадают. Это подтверждает известный факт, что гармо­ ническая линеаризация является частным случаем статистической линеаризации,.

Равенство коэффициентов статистической и гармонической лине­ аризации при отсутствии случайных возмущений позволяет выразить амплитуду автоколебаний через среднее квадратическое отклонение. Приравнивая коэффициенты статистической и гармонической лине­ аризации, определяем амплитуду автоколебаний как функцию среднего квадратического отклонения:

 

а = Аау/]/г2л.

(5.68)

Как известно,

дисперсия

гармонического

сигнала,

тривать его как процесс со случайной фазой, распределенной в ин­ тервале 0 — 2я, связана с амплитудой сигнала соотношением

а2

=

а2/2.

(5.69)

Это выражение отличается

от

соотношения

(5.68),чтообъясн

погрешностью метода статистической линеаризации, постулирую­ щего нормальный закон распределения входного сигнала нелиней­ ности. В данном случае закон распределения гармонического сигнала

154