Файл: Железнов, Ю. Д. Статистические исследования точности тонколистовой прокатки.pdf

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

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

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

Добавлен: 21.10.2024

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

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

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

p ( T )

Рис. 20. Кривая нормированной автокорреляционной функции

автокорреляционной функции обычно аппроксимируют аналити­ ческой функцией вида

R (т) = R (0) е~%1т I

(72)

при отсутствии скрытых периодичностей и

R (т) = R (0) е~а 1т 1cos сот

(73)

при наличии скрытых периодичностей. Для сложных процессов используют сумму функций вида (72 и 73). На рис. 20 приведена аппроксимируемая по формуле (73) кривая. Более точная аппрок­ симация может быть получена в виде суммы двух функций вида (73).

Т а б л и ц а 5. Расчетная таблица автокорреляционной функции

­

 

 

вы

 

1

 

Номер борки

ц*

ч~

 

° Ч

0

 

 

II

1

*1

О

о о

X

 

2

*2

X

*2*2

3

*3

X

* А

 

4

*4

X

X4T4

38

 

о

О о

*30

X

*38*38

39

*39

X

*39*39

40

*40

X

*40*40

 

 

V

 

30,9

 

2 j

 

10"4

 

р(т) =

 

1

= R )/R (0)

 

п

1 2 3 . . . 38 39 40

О о

о о

о о

О О

О о

о о

 

*1*3

*1*4

*1*38

*1*30

*1*40

*2*3

*2*4

*2*5

*2*39

*2*40

 

*3*4

*3*5

*3*0

*3*40

 

 

*4*5

*4*0

*4*7

 

 

 

О о

О о

 

 

 

 

*38*39

*38*40

 

 

 

 

*39*40

 

 

 

 

 

16 X

2,21 X

5,03х

—4,3 X

5 X

20 X

X Ю-4 х ю -*

X Ю~4

х ю - 4 х ю - 4 X 10-4

0,51

0,071

—0,16

—0,14

0,16

0,64

40


10. Спектральная плотность случайного процесса

Спектральная плотность не содержит дополнительной инфор­ мации по сравнению с автокорреляционной функцией, так как они связаны взаимным преобразованием Фурье. Однако опреде­ ление спектральной плотности в некоторых задачах становится необходимым. Спектральная плотность S (со) позволяет непосред­ ственно судить о частотных свойствах случайного процесса. Она характеризует интенсивность случайного процесса при различных частотах, отличающихся друг от друга на бесконечно малую ве­ личину.

По физическому смыслу спектральная плотность — это вели­ чина, которая пропорциональна средней мощности процесса в ин­ тервале частот от со до со, + Асо.

Спектральная плотность функции х (t) в общем виде опреде­ ляется по формуле [8]

 

 

 

S (со) = lim -^=- |F (/со) |2,

(74)

 

 

 

Г->оо л

 

где

Т

длина

реализации, с;

 

 

F (/со) —

преобразование Фурье функции х (t).

случай­

 

По формуле (74)

вычислить спектральную плотность

ного процесса почти невозможно, так как для определения изо­ бражения Фурье F (/со) необходимо знать аналитическое выраже­ ние случайного процесса х (/). Поэтому спектральную плот­ ность случайного процесса определяют в соответствии с теоремой Винера— Хинчина. Автокорреляционная функция и спектральная плотность представляют собой взаимные преобразования Фурье:

 

 

СО

 

S(C 0)=4 r

 

1

(75)

 

— СО

 

Я (т )=

Jсо

5 (со) е;МТ dco.

(76)

Спектральная плотность S (со) и автокорреляционная функция являются четной и вещественной функцией, поэтому формулы (75) и (76) удобно представить в следующем виде:

00

S (со) = —

f R (т) cos сот dv,

(77)

К

J

 

 

О

 

аналогично

 

 

00

 

 

R (т) = 2 |

S (со) cos сот dco.

(78)

О

 

 

41


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

оо

 

2 | S (со) dco = х2 (t).

(79)

о

 

Кроме этого, имеется связь между спектральной

плотностью

и видом функции х (/), которая описывает данный случайный про­ цесс. На рис. 21 показаны кривые спектральной плотности, соот­ ветствующие случайным процессам, приведенным на рис. 17.

Как видно из рис. 21, чем медленнее изменяется величина х (/), тем уже график спектральной плотности, и чем шире график спек­ тральной плотности, тем быстрее изменяется величина х (t). Пики на рис. 21 указывают на то, что в случайном процессе преобладают частоты, соответствующие этим пикам.

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

Ниже приведем расчет спектральной плотности для случай­ ного процесса, показанного на рис. 16.

Для расчета- S (со) по формуле (77) аппроксимируем кривую

автокорреляционной функции (рис. 22)

выражением вида

[8]:

R (т) = R (0) е—“ 1т 1cos |)т.

(80)

Коэффициенты аппроксимации: R (0)

= 30,9 10~4, а =

0,39;

Р - 3,2.

 

 

Таким образом, в первом приближении

R{ т) = 30,9-10 4е 0,39 1т 1cos 3,2т.

Рис. 21. Зависимость спектральной плотности от вида случайного процесса

42

Рис. 22. Спектральная плотность слуРис. 23. Спектральная плотность случайного чайного процесса, полученная с попроцесса (см. рис. 16), вычисленная на ЦВМ мощью аппроксимации (см. рис. 16)

Тогда спектральная плотность определится с помощью интег­ рала:

S(co) = - ^ J

30,9 -10'

4е 0,39 I х I cos 3,2т cos сот =

= 30,9-10-4 ■0,39 —

L 0,392 +

(3,2 — со)2 ^ 0.392 + (3,2 + со)2 . • (81)

я

По формуле (81) рассчитана спектральная плотность данной реализации и построена кривая S (со) (см. рис. 22).

Недостатком этого метода является то, что в этом случае мы не можем выделить все преобладающие частоты. Из кривой автокорреляционной функции видно, что в процессе преобладают как минимум две частоты. Для получения пиков в графике спек­ тральной плотности, соответствующих частотам, необходимо ав­ токорреляционную функцию интегрировать по всему спектру. Это достигается путем использования в расчетах цифровых вы­ числительных машин (ЦВМ).

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

11. Преобразование случайного процесса в динамической системе. Структура динамической системы; стационарность и линейность

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

43


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

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

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

Данная задача легко решается для линейных систем. В этом случае приходим к следующей задаче: зная математическое ожи­

дание х (t) и автокорреляционную функцию Rx (т) входного слу­ чайного сигнала и оператор динамической системы (передаточную

функцию), определить математическое ожидание у (t), дисперсию Dy (t) и автокорреляционную функцию Ry (т) выходного пара­ метра у (/).

Что же касается закона распределения, то нужно отметить, что он может изменяться при прохождении через линейную си­ стему.

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

На практике мы имеем дело с системами, состояние которых

вполне определено в момент времени

а в последующие моменты

времени параметры, характеризующие

эту систему, изменяются

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

динамики непрерывной

системы

(рис.

24)

можно

представить

в виде

 

 

 

 

 

 

 

- % - =

F (xi ■ ■ ■

хп, у1 ..

. Уп,

}),

(82)

где у х. .

.уп — выходные величины;

 

 

 

х х. .

.хп — входные величины.

 

 

 

44


fcuc. 24. Структурная схема динамической системы

Часто на практике представляют интерес не все выходные величины, а одна из них; тогда уравнение (82) можно записать в следующем виде:

 

М 0 4 ^

+ 6п- 1 (/)5 Й +

+ M 0 ffi = * (0 .

. (83)

где

х (t) —

внешнее возмущение;

 

bn (t).

. .b0 (t) —

коэффициенты, характеризующие параметры

 

 

системы.

 

 

Если в процессе работы все без исключения параметры системы остаются без изменения, т. е. коэффициенты bn (t) = ■■•= b0 (t) = = const, то такая система называется внутренне стационарной. В противном случае система нестационарна.

Пусть на вход системы подан сигнал x(t), тогда реакция си­ стемы на этот сигнал будет у (/).

Если при сдвиге входного сигнала той же формы выходной сиг­ нал не изменится, а только сдвинется по времени на т, то такая система называется внешне стационарной [9] (рис. 25). Внешняя стационарность системы относительна. Система может быть ста­ ционарной относительно какого-нибудь входа и выхода, а также при всех комбинациях входов и выходов. Система может обладать внутренней стационарностью и не быть внешне стационарной, и наоборот. Абсолютно стационарной называется система, которая стационарна и внутренне, и внешне.

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

А ( Д *<•(')) = £ ^ - ( 0 :

(84)

в частном случае для двух входов и выходов можно записать

■4 [*i (О ~Ь х2 (^)] = Axi (i) -f~ Ах2 (t)•

Формула (84) показывает справедливость принципа суперпо­ зиции (наложения) для линейных систем.

На практике линейные системы встречаются реже, чем нели­ нейные. Чаще имеем дело с непрерывными системами, линейными

45