Статистический анализ линейных систем автоматического управления

Московский государственный технический университет  
имени Н.Э. Баумана 
А.С. Ющенко, Д.С. Делия 
 
 
СТАТИСТИЧЕСКИЙ АНАЛИЗ  
ЛИНЕЙНЫХ СИСТЕМ 
АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ 
 
 
Методические указания к лабораторному практикуму  
по курсу «Статистическая динамика систем  
автоматического управления» 
 
 
Под редакцией А.С. Ющенко 
 
 
 
 
 
 
 
 
 
М о с к в а  
Издательство МГТУ им. Н.Э. Баумана 
2 0 0 9  
 

УДК 519.711.3 
ББК 22.161.6 
Ю985 
Рецензент д-р техн. наук, проф. В.С. Медведев 
 
Ю985 
Ющенко А.С., Делия Д.С.  
Статистический анализ линейных систем автоматического управления: Метод. указания к лабораторному практикуму по курсу «Статистическая динамика систем автоматического управления» / Под ред. А.С. Ющенко. — М.: Изд-во 
МГТУ им. Н.Э. Баумана, 2009. — 52 с.: ил. 
 
В работе приведено описание трех лабораторных работ по курсу 
«Статистическая динамика систем автоматического управления» с 
использованием ПЭВМ и пакета MATLAB 7, посвященных формированию случайных процессов с заданными характеристиками, анализу 
линейной непрерывной автоматической системы при воздействии на 
нее случайного процесса и исследованию автоматической системы 
при воздействии на нее случайного сигнала и помехи. Методические 
указания содержат необходимые математические формулы, сведения 
о пакете MATLAB 7, необходимые при проведении лабораторных работ, а также примеры проведения исследований. 
Для студентов старших курсов, обучающихся по направлениям 
подготовки 220400 «Мехатроника и робототехника» и 220200 «Управление в технических системах». 
УДК 519.711.3 
ББК 22.161.6 
Учебное издание 
Ющенко Аркадий Семенович 
Делия Данила Сергеевич 
СТАТИСТИЧЕСКИЙ АНАЛИЗ ЛИНЕЙНЫХ СИСТЕМ  
АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ 
 
Редактор С.А. Серебрякова 
Корректор  М.А. Василевская 
Компьютерная верстка С.А. Серебряковой 
Подписано в печать  16.12.2008. Формат 60×84/16. Бумага офсетная. 
Усл. печ. л. 3,02. Уч.-изд. л. 2,83.  
Изд. № 148. Тираж 200 экз. Заказ        . 
 
Издательство МГТУ им. Н.Э. Баумана 
Типография МГТУ им. Н.Э. Баумана 
105005, Москва, 2-я Бауманская ул., 5 
 
          © МГТУ им. Н.Э. Баумана, 2009 
 

ВВЕДЕНИЕ 
Данное пособие предназначено для самостоятельного проведения 
студентами цикла лабораторных работ по курсу «Статистическая динамика систем автоматического управления» на ПЭВМ с использованием 
программного пакета MATLAB. Такой курс предусмотрен в учебном 
плане ряда специальностей, связанных с автоматизацией, мехатроникой и 
робототехникой. Предполагается, что студенты знакомы с основным курсом теории автоматического управления и с теорией вероятностей. Для 
того чтобы облегчить проведение лабораторных работ, в методических 
указаниях даны основные соотношения из курса «Статистическая динамика систем автоматического управления», которые используются при 
статистическом анализе линейных стационарных автоматических систем. 
Приведено краткое описание возможностей пакета MATLAB 7 для проведения такого анализа.  
В пособии даны указания к проведению трех лабораторных работ, 
объединенных в общий цикл. Вначале нужно построить формирующий 
фильтр, позволяющий получить случайный процесс с заданными статистическими характеристиками. Затем провести исследование системы, на 
вход которой поступают детерминированный сигнал и случайная помеха. 
Наконец, предлагается исследовать замкнутую систему автоматического 
регулирования, ко входу которой приложен случайный сигнал, представляющий собой аддитивную смесь полезного случайного сигнала и случайного шума. Описание каждой лабораторной работы сопровождается 
примерами. Для того чтобы правильно интерпретировать получаемые 
результаты, мы рекомендуем предварительно выполнить статистическое 
исследование в аналитической форме с использованием приведенных в 
пособии соотношений. Более подробную информацию о методах статистического анализа автоматических систем и о возможностях пакета 
MATLAB можно найти в литературе, список которой приведен в конце 
пособия [1 – 4]. 
 
3 

 
1. ОСНОВНЫЕ МАТЕМАТИЧЕСКИЕ СООТНОШЕНИЯ, 
ИСПОЛЬЗУЕМЫЕ ПРИ СТАТИСТИЧЕСКОМ АНАЛИЗЕ 
АВТОМАТИЧЕСКИХ СИСТЕМ 
1.1. Характеристики случайных процессов 
Основными характеристиками случайных величин являются 
функция распределения вероятности и плотность распределения 
вероятности. Функция распределения случайной величины ξ, принимающей любые вещественные значения, определяется соотношением 
 
( )
{
}
F x
P
x
ξ
=
ξ <
  
(1) 
и представляет собой вероятность того, что случайная величина ξ 
принимает значения, меньшие заданного значения x. 
Плотность распределения вероятности случайной величины ξ 
может быть определена по функции распределения вероятности  
с использованием формулы  
 
( )
( )
.
dF x
p x
ξ
ξ
=
 
(2) 
dx
Случайный процесс определяется множеством случайных реализаций ξ(t), 0 ≤ t ≤ T. Фиксируя произвольным образом моменты 
времени ti, i = 1, 2, …, N, можно получить N-мерную случайную величину ξ = [ξ(t1), ξ(t2),…,ξ(tN)], т. е. случайный вектор, компонентами которого являются случайные величины, представляющие собой 
значения реализаций ξ(t) в дискретные моменты времени. Таким 
образом, случайный процесс характеризуется множеством функций 
распределения вероятности, определяющих векторную случайную 
величину ξ: 
 
4 

2
1
2
1
1
2
2
( ,
,...,
, , ,..., )
{ (
,
( )
, ..., ( )
},
k
k
k
k
k
F
x x
x t t
t
P
t
x
t
x
t
x
ξ
=
ξ
) <
ξ
<
ξ
<
 
(3) 
где k = 1, 2, …, N, или соответствующим множеством плотностей 
распределения вероятности: 
(
,
,...,
,
,
,...,
)
p
x
x
x
t
t
t
=
ξ
k
k
k
1
2
1
2
k
 
 
(4)
 
ξ
                    
(
,
,...,
,
,
,...,
),
...
F
x
x
x
t
t
t
x
x
x
∂
= ∂
∂
∂
k
k
k
k
1
2
1
2
1
2
где k = 1, 2, …, N. 
Используя плотности распределения вероятности, можно определить моменты различного порядка для случайного процесса 
ξ(t). Наиболее часто применяют начальный момент первого порядка (математическое ожидание): 
∞
 
1
( )
{ ( )}
( , )
m t
M
t
xp
x t dx
ξ
ξ
−∞
=
ξ
= ∫
  
(5) 
и центральный момент второго порядка (корреляционную (автокорреляционную) функцию): 


( ,
)
{ ( ) ( )}
K
t t
M
t
t
=
ξ
ξ
=
ξ
1
2
1
2
∞
∞
 
(6)
 
                    
(
( ))(
( ))
(
,
, ,
)
.
x
m t
x
m t
p
x x t t
dx dx
1
1
2
2
2
1
2
1
2
1
2
=
−
−
∫∫
ξ
ξ
ξ
−∞−∞
Из последнего выражения можно также найти дисперсию случайного процесса: 
 
( )
{ ( ) ( )}.
D t
M
t
t
ξ
=
ξ
ξ


  
(7) 
Взаимнокорреляционную функцию двух случайных процессов 
ξ(t) и η(t) определяют по формуле 


K
t t
M
t
t
( , )
{ ( ) ( )}
=
ξ
η
=
ξη
1
2
1
2
∞∞
 
(
( ))(
( ))
( ,
, , )
.
x
m t
x
m t
p
x x t t dx dx
  
(8)
 
1
1
2
2
2
1
2
1
2
1
2
=
−
−
∫∫
ξ
ξ
ξη
−∞−∞
 
5 

Напомним, что случайный процесс ξ(t) называется стационарным (в широком смысле), если его математическое ожидание постоянно, а корреляционная функция зависит только от разности 
аргументов: 
 
( )
const
m t
m
ξ
ξ
=
=
,  
(9) 
 
1
2
1
2
( , )
(
)
( ).
K t t
K t
t
K
ξ
ξ
ξ
=
−
=
τ  
(10) 
Для стационарных процессов можно определить спектральную 
плотность случайного процесса как преобразование Фурье корреляционной функции: 
 
( )
( )
.
j
S
K
e
d
∞
−ωτ
ξ
ξ
−∞
ω =
τ
τ
∫
  
(11) 
По заданной спектральной плотности можно определить корреляционную функцию: 
j
K
S
e
d
 
1
( )
( )
.
2
∞
ωτ
ξ
ξ
−∞
τ =
ω
ω
π ∫
  
(12) 
Из последней формулы следует и выражение для вычисления дисперсии стационарного случайного процесса по его спектральной 
плотности: 
∞
 
1
(0)
( )
.
2
D
K
S
d
ξ
ξ
ξ
−∞
=
=
ω
ω
π ∫
  
(13) 
При статистической обработке случайных сигналов используются методы и формулы математической статистики. 
Пусть на интервале времени [0, T] экспериментально получено 
n реализаций случайного процесса ξ(t), которые мы обозначим 
xi(t), i = 1, 2, …, n. Тогда оценка математического ожидания случайного процесса может быть определена по формуле 
n
1
ˆ ( )
( ).
 
1
i
i
m t
x t
n
ξ
=
= ∑
  
(14) 
Оценка (14) является несмещенной, т. е. 
ˆ
{
}
,
M m
m
ξ
ξ
=
 и состоятельной, поскольку выполняется условие 
ˆ
lim
{
}
0.
n
D mξ
→∞
=
 
 
6