Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов

.pdf
Скачиваний:
71
Добавлен:
25.10.2023
Размер:
16.48 Mб
Скачать

118

Глава 3

где <j>(() — произвольная

заданная функция, для которой эт

интеграл существует. При любых заданных ф (t) и пределах ин­ тегрирования (а, Ь) характеристика I есть случайная величина, зависящая от конкретной выборочной функции х(і). Для иссле­ дования статистических свойств случайной величины / область интегрирования (а, Ь) обычно разбивают на интервалы At и рас-'

сматривают приближение в виде линейной комбинации

.-

N

 

 

 

 

 

Іы« 2 * (ш )ф (

*

'

(

3

• 139)

і=1

 

 

 

 

 

Сходимость IN к / может быть определена различными крите^ риями. Говорят, что последовательность [IN\ сходится к /

1) в среднем квадратическом, если

lim М ПlN—/ | а]= 0 ;

(3,140а)

N—>-оо

2)по вероятности, если для любого е > 0

lim Р [ 11N—/ 1> е]= 0 .

(3.1406)

N *~со

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

смысле сходимости в среднем квадратическом.

3.5.2.Теоремы о дискретном представлении случайных процессов

Реализации случайного процесса с~непрерывным временем часто представляются и анализируются в дискретной форме. Связанные с этим проблемы будут рассмотрены в гл. 7 и 9. Здесь же приведены две важные теоремы о дискретном представлении случайных процессов, необходимые для понимания излагаемого далее материала.

Пусть реализация x(t) случайного процесса (х(/)} задана в интервале времени от 0 до Г секунд и равна нулю вне этого ин­ тервала. Преобразование Фурье этой реализации

т

X (/)= j* X

dt.

(3. Щ

b

 

 

Для того чтобы получить периодическую функцию с периодом Т секунд, предположим, что реализация x(t) непрерывно повторяется.

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

119

^Основное приращение частоты / = 1 IT. Разлагая

функцию в

ряд Фурье, находим

 

 

 

 

 

00

 

х(0 =

2 Апе™мт,

(3.142)

где

 

— СО

 

 

т

 

 

 

 

Ап=

~ ^ x(t)e~2*WTdt.

(3.143)

 

 

о

 

Из формулы (3.141) следует, что

 

і

 

Г

 

X ( - £ - j

=

j1X (t)e-™WTd t = A J .

(3.144)

о

Таким образом, величина Х(пІТ) определяет значения коэффи­ циентов Ап и, следователььо, ординаты х(і) при всех t. Вид функ­ ции x(t) в свою очередь определяет величины X(f) при всех зна­ чениях /. Этот вывод составляет содержание теоремы о дискрет­ ном представлении процесса в частотной области. Основное при­ ращение частоты / = 1 IT называется коинтервалом Найквиста.

Пусть преобразование Фурье X(J) некоторой реализации задано в интервале частот от —В до В Гц и равно нулю вне этого интервала. Интервал физически осуществимых частот состав­ ляет 0 — В Гц. Обратное преобразование Фурье имеет вид

Ч

в

(3.145)

x(t)=

j X (f)en-*ll‘df.

—в

Для того чтобы получить периодическую функцию частоты с периодом 2В Гц, положим, что функция X(f) непрерывно повто­ ряется. Основное приращение времени составляет і = 1/2В. Теперь

X (/) =

NT Спе-*МВ ,

(3.146)

где

в

 

 

(3.147)

СП = ^ Г

1 * ( № ІПиВѴ-

 

—в

 

щз формулы (3.145) следует, что

 

в

 

 

X { - й - ) = J

X{f)t*M 4f=ZBCa.

(3.148)

—в

х ( п І 2 В )

120 Глава 3

Таким образом, величина определяет значения коэффи#*' шіентов С„ и, следовательно, функцию X(f) при всех значениях /. Вид этой функции в свою очередь определяет ординаты х{1) при

всех значениях

t. Этот

вывод

составляет

содержание т е о р е м ы

о д и с к р е т н о м

п р е д с т а в л е н и и

п р о ц е с с а

во

в р е м е н н о й

о б л а с т и 1'*.

Основное приращение

времени 1 /2 В

называется

и н т е р в а л о м

Н а й к в и с т а .

Предположим теперь, что реализация х ( І ) задана только в

интервале

времени от'0 до Т секунд, а ее

преобразование Фурье

X ( f ) — в

интервале частот от —В до В

Гц. Это двойственное

предположение теоретически невозможно в силу

п р и н ц и п а

н е о п ­

р е д е л е н н о с т и [1]. В действительности, однако,

оно может

быть/;

приближенно справедливо для конечных интервалов времени и для полосовых фильтров. Полагая, что на функции х ( і) и Х(/) нало­ жены такие ограничения, касающиеся интервалов времени и частот, можно показать, что для определения функции х(1) при всех значениях t необходимо знать лишь конечное число дис­ кретных значений х(І) или Х(/). Согласно формуле (3.144), снимая дискретные значения функции X ( f ) в точках, разделен­ ных по шкале частот коннтервалом Найквиста М Т в промежутке от —В до В , можно найти число дискретных значений, которое необходимо для описания функции х ( і ) . Это число равно

N = - — = 2BT.

(3.149)

Согласно формуле (3.148), снимая дискретные значения функціи^ х(1) в точках, разделенных по шкале времени интервалом Найк­ виста 1/2Ö в промежутке от 0 до Т , можно найти, что

N = T j W = 2ßT-

(з.ібо)

Таким образом, требуется одинаковое число дискретных значений при выборке их через коинтервал Найквиста по шкале частот и при выборке через интервал Найквиста по шкале времени.

Упражнения

1. Рассмотрите, являются ли плотность распределения р ( х ) , автокорреляционная функция Rx(t) и спектральная плотность

Gx(f)

стационарного случайного

процесса

а)

четными;

 

б)

неотрицательными;

Д

в)

ограниченными в нуле.

 

:) Доказательство этой теоремы принадлежит акад. В. А. Котельнико­ ву.— Прим. ред.

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

121

'Определяют ли эти функции г) среднее значение процесса; д) дисперсию процесса; е) плотность распределения;

ж) автокорреляционную функцию; з) спектральную плотность?

Какие из этих утверждений справедливы для функций р(х), Rx(т) и Gx(f) произвольного случайного процесса?

2. Докажите уравнения (3.48) и (3.49).

3. Докажите уравнения (3.62) и (3.63).

4.Докажите теорему Парсеваля [формула (3.105)] и объяс­ ните, каким образом получено уравнение (3.108).

5.Подставив выражения (3.146) и (3.148) в формулу (3.145),. покажите, что функция х(/) может быть выражена через дискрет­ ную функцию х(пі2В) следующим образом:

 

 

п

СО

*(я/25)

/ ~2В

X (0 = 2

п

—00

 

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

р ( * ) = 4 - 6 (X + 3) + 4 - б ( * - 2 ) + — 4 = . е - * ѵ .

*•

*

8 у

Определите

ка) среднее значение и дисперсию;

б)

вероятность события, состоящего в том,

что х > 1.

7.

Функция Rx(т) = 25 е-4 1т 1 cos 2зт/т -f

16. Определите

а) среднее значение и дисперсию;

 

б)

соответствующую одностороннюю спектральную плотность.

8. Какая из нижеследующих функций может быть односто­

ронней спектральной

плотностью?

 

 

 

п m _

f2 + 9

I)2

Л /п

 

/2+ 1

l U )

(Р +

4)(f +

 

 

/4 + 5р + 6>

 

f *

+

4

 

 

g4(/)=

exp ( — j R )

g8(/)= f * — 4/2 + 3 '

 

Г- + 2

9. Дана двусторонняя

спектральная

плотность

 

 

83 (/) +

2 о ( і

10

),

/ < І 0 .

 

 

Ш

 

 

S A D = .0,

 

 

 

 

1/1 > 10.

Определите а) среднее значение и дисперсию;

б) соответствующую автокорреляционную функцию.

122

 

 

Глава 3

10.

Какое из нижеследующих свойств всегда справедливо-

для пары эргодических случайных процессов?

а) Rxy{оо) = р,р,.

 

б)

Если Rxy(0) = 0, то рл = 0 или ри = 0.

в)

Rxy(x) —

если

Rx(x) = 0 или Ry(x) = 0-

Г)

\ R xy(x ) \ 2 <

R x (r)

R y ( r ) .

д)

Gxy(0) = 0, если р* = 0 или р„ = 0.

е)

\GXy ( f ) \ ^ G x(0)Gy(0).

ж)

Gxyif) = 0,

если Gx{f) = 0 или Gy(f) = 0.

з)

Gx{oo) = Gy(оо) = GXy{оо) = 0.

и)

Если р(х, у) = р(х)р(у), то Rxy{r) = 0.

к)

Если р{X, у) = 0, то GXy(f) = 0.

4

Л

ОСНОВНЫЕ ПОЛОЖЕНИЯ МАТЕМАТИЧЕСКОЙ СТАТИСТИКИ

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

4.1. Выборочные величины и оценивание параметров

Рассмотрим случайную величину х, определение которой приведено в подразд. 3.1.1. Индекс k, обозначающий точку вы­ борочного пространства, для простоты опустим. Введем теперь в рассмотрение два основных параметра, которые характеризуют среднее положение и степень изменчивости случайной величины X — среднее значение этой величины и ее дисперсию. Согласно формулам (3.8) и (3.11), среднее значение и дисперсия опре­ деляютсяследующим образом:

СО

 

(4.1)

— СО

 

со

 

=м [(*— ^ ) 2]== J ( Х— |А,)ар (x)dx,

(4.2)

—00

 

где р(х) — плотность распределения случайной величины х. На практике точно определить эти параметры, разумеется, нель­ зя, так как плотность распределения, как правило, не известна.

124 Глава 4

Следовательно, остается довольствоваться лишь оценками сред­ него значения и дисперсии, основанными на конечном числе на­

блюденных величин.

 

получе­

Один из существующих способов, используемых для

ния оценок среднего значения

и дисперсии величины х на осью-,

не N независимых измеренных значений этой величины,

заклю­

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

 

_ .

X N

(4.3)

 

= -дГ 2 хі’

 

1=1

 

/V

(4.4)

і=Ч

Здесь X и si представляют собой выборочное среднее и выборочную

дисперсию соответственно. Знак (л) у р* и а? означает, что эти выборочные параметры используются в качестве оценок среднего значения и дисперсии величины х. Нижний индекс у si означает, что этот параметр представляет собой смещенную оценку диспер­ сии (этот вопрос рассмотрен ниже). Число наблюденных значений случайной величины, используемых для получения оценок (вы­ борочных величин), называется объемом выборки.

Формулы (4.3) и (4.4)

не единственные,

которые можно ис­

пользовать для

вычисления

оценок среднего

значения и диспер­

сии случайной

величины х.

Подходящие оценки среднего значе-/

ння и дисперсии можно также

получить,

например, путем деле­

ния сумм, входящих

в соотношения (4.3) и (4.4), на (N — 1), а

не на N. Сами по себе

оценки

не являются

правильными или

неправильными,

так как

они

определены

в

некоторой степени

произвольно. Тем не менее некоторые оценки можно считать «хорошими» или «лучшими» по сравнению с другими.

Для того чтобы определить качество или «доброкачествен­ ность» оценки, необходимо принять во внимание три обстоятель­ ства. Во-первых, желательно, чтобы математическое ожидание

оценки было равно определяемому параметру, т.

е.

М [Ф ]=Ф ,

(4.5)

где Ф — оценка параметра Ф. Если это условие выполняется,, оценка называется несмещенной. Во-вторых, желательно, чтобь*^ средний квадрат ошибки данной оценки был не больше среднего квадрата ошибки при любой другой оценке, т. е.

М К ® !— Ф )2] < м [(<&, — ф )21,

(4.6>

Основные положения математической статистики

125

'где Фг — рассматриваемая оценка, Фг ■— любая другая

оцен­

ка. При выполнении этого условия оценка называется эффек­

тивной.

В-третьих, желательно, чтобы по мере увеличения объема

выборки

оценка приближалась к оцениваемому

параметру с

ве­

роятностью, стремящейся к единице. Это значит, что при

лю­

бом е >

.0

 

 

 

Нгп Р [| Ф—Ф| > е ] = 0 .

(4.7а)

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

lim М [(Ф—Ф)2]= 0 .

(4.76)

N

 

Заметим, что условия (4.7) выражают в сущности

определяемые

соотношениями (3.140) требования сходимости: а) по вероятности и б) в среднем квадратическом.

Рассмотрим оценку среднего значения (4.3). Математическое

ожидание выборочного

среднего значения х

составляет

 

М [х]=М

1 N

'

N

 

 

(4.8)

 

= ~ff

=

F S * «

2 * «

 

 

/«=і

_

г=і

 

 

 

Следовательно,

как

это

вытекает ттз

уравнения (4.5),

оценка

=

X является несмещенной. Средний квадрат ошибки выбороч­

ного

среднего

х

 

 

 

 

 

м[(*-нд8]= м

 

 

N*1 м

S ( * i — р*)

 

 

 

 

 

 

- \і=і

 

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

N

■М [(*-р,)8] = ж М Х !(ч —іч)а — дга -{Noiy- N

(4.9)

j= i

126

Глава 4

Следовательно, согласно уравнению (4.76), оценка цх = х со­ стоятельна. Можно показать, что эта оценка также и эффективна.

Рассмотрим теперь оценку выборочной дисперсии (4.4). Мате­ матическое ожидание выборочной дисперсии si есть

М (s|] = М 1

У ( Хі—х)*

N

 

N

і=і

_1=1

Далее,

 

 

 

2

“ *)2=2 (*i—h.+ V-x— x)2 =

 

 

i = i

i = i

 

 

=2 (xi— p*)a— 2 (x ~-^x) 2 (x‘~ ^ +

N (x ~

 

i = i

/ = i

 

 

=

2 ^ - ^

2- 2 (X ~ ^ W (х— Цх) + N

( x — ixx) 2=

 

1=1

 

 

 

 

 

= 2 ^ - ^ 2- Л' ^ - ^ а-

(4Л0>

 

 

i = i

 

 

Так как

М[(Хі — p,z)2] = а2 и Щ(х — р,л)2] = а2/М,

то

 

М[sl\=-jj- (N ol-ol) =

а2.

(4.11 /

Следовательно, оценка = s* смещена. Хотя выборочная ди­ сперсия si представляет собой смещенную оценку дисперсии а2, она является состоятельной и эффективной.

Из формулы (4.11) следует, что несмещенную оценку диспер­ сии а2 можно получить, вычисляя несколько отличную выбо­ рочную дисперсию

s2= a2= - ^ - r 2 і~ х)2.

,(4.12)

/=1

Статистика (4.12) есть несмешанная оценка дисперсии а2. Кроме того, она удовлетворяет условиям эффективности и состоятель­ ности. Поэтому выборочную дисперсию (4.12) обычно считаютъ «лучшей» оценкой по сравнению с выборочной дисперсией (4Л}Т В дальнейшем в качестве оценки истинной дисперсии случайной величины будет использоваться выборочная дисперсия, опреде­ ляемая формулой (4.12).

Основные положения математической статистики

127

^|.2. Наиболее важные функции распределения

Несколько примеров теоретических функций распределения приведено в разд. 3.1. Наиболее важная из них с точки зрения прикладной статистики — это гауссовское {нормальное) распреде­ ление. В статистике широко используются еще три функции рас­ пределения, связанные с нормально распределенными случайными величинами. К ним относятся распределение х2, /-распределе­ ние и ^-распределение. Все эти распределения, включая нор­ мальное, будут описаны и рассмотрены ниже. Примеры их ис­ пользования при анализе приведены в последующих разделах.

V 2.I. Нормальное распределение

Нормальные плотность и функция распределения случайной величины X определяются уравнениями (3.35) в подразд. 3.1.3. Можно представить нормальное распределение в более удобной форме, если воспользоваться нормированной величиной

* = -* 7 ^ ■

(4.13)

Подставив равенство (4.13) в формулы (3.35), получим нормаль­ ные плотность и функцию распределения нормированной случай­ ной величины с нулевым математическим ожиданием и единичной дисперсией (р,г = 0, о%= 1):

 

Р (z )= C / 2я )~ге 2/2,

(4.14а)

 

 

Z

 

 

 

Р (г)=

(/2 іГ ) _1 Г

d t

(4.146)

 

 

— с о

 

 

Для последующего рассмотрения удобно обозначить

величину z,

которая соответствует заданной вероятности*/3^ ) ^

1 — а, сим­

волом Za.

Тогда

 

 

 

 

Р (Za) = JР (z ) d z = Р <

Za] = 1 —a,

(4.15а)

 

— СО

 

 

 

 

 

оо

 

 

*

1 — Р (г а) =

Г р ( z ) d z = Р [Z > Za] = a.

(4.156)

■Щ,

 

J

 

 

 

 

*a

 

 

величину

za, удовлетворяющую формулам (4.15), называют

100а %-ной

точкой.

 

 

 

Соседние файлы в папке книги из ГПНТБ