
книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов
.pdf118 |
Глава 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) |
—в
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 |
|
2Ѣ |
6. Случайная величина х имеет плотность распределения |
||
р ( * ) = 4 - 6 (X + 3) + 4 - б ( * - 2 ) + — 4 = . е - * ѵ . |
||
*• |
* |
8 у 2л |
Определите
ка) среднее значение и дисперсию;
б) |
вероятность события, состоящего в том, |
что х > 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) — 0» |
если |
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) = J2аР (z ) d z = Р [г < |
Za] = 1 —a, |
(4.15а) |
|
|
— СО |
|
|
|
|
|
оо |
|
|
* |
1 — Р (г а) = |
Г р ( z ) d z = Р [Z > Za] = a. |
(4.156) |
|
■Щ, |
|
J |
|
|
|
|
*a |
|
|
величину |
za, удовлетворяющую формулам (4.15), называют |
|||
100а %-ной |
точкой. |
|
|
|