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

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

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

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

Добавлен: 16.10.2024

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

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

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

В частном случае, когда Тк < Т, можно пренебречь произведе­ нием постоянных времени. Тогда дисперсию вычисляют по более простой формуле

D N a [Т ы (“ + nk\) + Т -f тм]

(5.119)

~ (Т + Тм) {n/h + а [1 + а + Тм)]}

Поскольку k y есть функция сто [см. (5.113)], то данное соотноше­ ние представляет собой нелинейное уравнение относительно вели­ чины D<>. Преобразуем соотношение (5.119) к виду

d2

 

 

9

9

 

 

 

з +,

2/1

2аЬк2

2

Да/ (Т + Тм+ аТр

CT#

k2 У 2 я а [I

а (Г +

Тм)] ■о#

(Т + Тм)[1 + а ( Т + Т ы)]

 

 

 

 

 

d 2

 

 

 

 

 

 

= 0.

(5.120)

* 2 / 2 я ( Г + .П З Г Ц - а С Г + Г м )]

Если d С 2kla\, то уравнение (5.120) превращается в кубиче­ ское уравнение относительно среднего квадратического отклонения температуры:

3 I

2п1

2

 

k2^2яа[1 +о(Г + Гм)]

 

 

DN (Т + Ты + аТм2)

<%—

 

(Т + Гм)[1

+ а (Г + Гм)]

 

D NT м'

 

 

=0 .

(5.121)

k2 У2л (Т+Тм)[1 - а ( Т + Т м)]

При известных параметрах данное уравнение легко решается. При условии Тм С Т можно пренебречь постоянной времени привода. В этом случае уравнение (5.121) превращается в следующее

соотношение:

 

з

2nl

 

 

Ау

 

 

(5.122)

 

По

/?2 У2л. (х (1 -|- аТ)

 

1 + а Т Oq

= 0.

 

Корни этого уравнения равны

о® = 0, и

 

 

 

 

Од. —

til

 

 

1

(1 -)- аТ)-2пЩа?

(5.123)

2 У 2л а (1

+ аТ)

 

+ Д 'N '

IO

 

 

 

 

 

 

 

 

 

Для

значений

параметров

п = 0,01 В ^ - с -1,

/ =

25

В,

k 2 —

= 0,01

рад-град-1, Т = 10 с,

DN =

100 град2,

а =

0,1

с-1

[62]

среднее квадратическое значение колебаний температуры объекта Сто = 0,5 град. При а = 0,01 с-1 и тех же значениях остальных пара­

метров

среднее квадратическое отклонение температуры

==

= 0,05

град.

 

165


Рассмотренный выше частный случай соответствует нелинейному элементу без зоны нечувствительности. Если рассмотреть другой крайний случай, когда зона нечувствительности равна бесконеч­ ности, то вместо уравнения (5.120) будет справедливо следующее уравнение:

з

DN (T + TM+ aTl)

(5.124)

 

(7, +

7’м ) [ 1 + а ( 7 ’ +

7'н) ] <Т» -

 

 

Это уравнение имеет

решения

 

 

ст» -

п

2

DN (T + T„ + aTl)

(5.125)

и, сто -

+ TJ п + а (Т + Гм)]

Пренебрегая постоянной времени регулятора Тм = 0,

получаем

простую формулу

 

 

 

 

 

 

 

 

 

 

(5-126)

При Dn = 100 град2,

а

= 0,1 c~\

Т '= 10 с среднее

квадра­

тическое отклонение температуры объекта при отсутствии обратной связи через регулятор = 7,06 град.

Таким образом, при наличии определенной зоны нечувствитель­ ности среднее квадратическое отклонение температуры объекта

заключено в пределах

 

стО0)

(5.127)

где (7^ определяется формулой (5.125) или (5.126),

а ст0о— фор­

мулой (5.123) или уравнением (5.121).

 


Г л а в а 6 М Е Т О Д Ы А Н А Л И З А Д И С К Р Е Т Н Ы Х СИСТЕМ

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

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

Любая дискретная система преобразует непрерывные функции

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

вдискретные моменты времени. Это относится также к случайным процессам. Для описания случайных процессов в дискретных си­ стемах можно пользоваться теми же характеристиками, что и для непрерывных систем, если считать эти характеристики функциями дискретных аргументов. Так, математическим ожиданием случайной функции дискретного аргумента X (/г) будем называть такую функ­ цию тх (/г), которая при каждом значении аргумента А равна мате­ матическому ожиданию функции X (А) при данном А:

тх (А)

= M I X (А)].

 

 

Математическое ожидание

случайной функции

X,

наблюдаемой

в непрерывные моменты времени t = (А + е) Тп,

0 ^

е ^ 1,

1пх (А, е)

= M IX (А, е)].

 

 

Корреляционной функцией случайной функции дискретного аргу­ мента X (А) называется функция, определяемая соотношением

Кх (А, Ах) = М ([X (А) — шх (А)] [X (А2) — mx (Ах)]).

Из этого выражения при Аг = А получаем дисперсию случайной функции X (/г).

Корреляционная функция случайной функции X (А, е), наблю­

даемой в непрерывные смещенные моменты времени

t =

(/г + е) Тп,

О^

е ^

1, определяется

формулой

 

 

 

 

Кх (А, е, Ах, ex) =

М {[X (Ах е) — mx (А,

е)]

X

 

 

X [X (Ах,

8х) — тх (Ах, ех)]J •

 

 

 

 

е

1, 0 ^ 8х с 1).

 

 

При

Ах = А, ех = е эта формула определяет дисперсию случай­

ной

функции.

 

 

 

 

Аналогично могут быть обобщены определения для других мо­ ментов дискретных случайных функций. Функция плотности вероят­ ности дискретного случайного процесса также зависит от дискрет­ ного момента времени.

167


Для стационарной случайной функции дискретного аргумента, как и для непрерывной случайной функции, имеем

тх (/г) = const, Кх (/г, /гх) = 1гх (/г— /гх).

Спектральная плотность дискретной стационарной случайной функции определяется по формуле

си

£< “ >= & £

где Тп — период повторения импульсов (такт) исследуемой дискрет­ ной системы. Спектральная плотность стационарной случайной по­ следовательности равна интенсивности соответствующего дискрет­ ного белого шума Vd в интервале

СО

Интегральное каноническое представление стационарной слу­ чайной функции X (/г) выражается через Vd следующим образом [59]:

X (/г)

= тх (/г) + j

Vd (со)е'“ЛГп dco,

 

Я

 

где

 

 

М =

[Vd(со) Vd(со')] =

Sxd (со) б (со - со').

Дискретный стационарный белый шум X (h) так же, как и непре рывный, имеет постоянную интенсивность

Gd = ~ S dx = const.

* П

При этом корреляционная функция дискретного белого шума имеет вид

Gd при h1 — h2,

kx {hi lb) — О^б/,,/,.

О при /г2 =?= /г2-

Преобразование случайных функций дискретными системами может быть описано формулами, аналогичными рассмотренным в предыдущих главах для непрерывных систем. В частности для линейных дискретных систем выходная векторная функция Y (т, е) связана с входной векторной функцией X (h) линейным операто­ ром А и выражается формулой

Y (т, в) = AhX (Л),

(6.1)

где Ah — произвольный линейный оператор, применяемый к аргу­ менту h. Так как операции оператора Ah и математического ожида­ ния М перестановочны, то, применяя операцию математического ожидания к выражению (6.1), получим

ти (m, е) = Ah тх (Л).

(6.2)

168


Для корреляционной функции выходной переменной также спра­ ведлива формула, аналогичная соответствующему выражению в не­ прерывных системах:

Ку

(т,

е, Шц 8j) = AhA,nKx (h, /гД

(6.3)

Формулы (6.2)

и

(6.3) аналогичны соответственно

формулам

(2.6) и (2.8).

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

Д , (т, г, т±, в О = АкА„Гх (Л, / гД

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

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

6.2. Метод весовых функций

Для вероятностного анализа линейной дискретной системы может быть применен метод весовых функций, основанный на использова­ нии формулы (1.41), если известны весовые функции дискретной системы.

Для одномерной физически возможной нестационарной дискрет­ ной системы связь между входной и выходной случайными перемен­ ными на основании формулы (1.41) имеет вид

h

у (Л, В) = k=082 (Л, е, k) X (k). (6.4)

Применяя операцию математического ожидания к правой и левой частям выражения (6.4), получим

 

л

 

ту (/г, е) =

Е ё Д, е. Щтх (/г).

(6.5)

k=0

 

169