
книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов
.pdf48 |
Глава 1 |
Фурье автокорреляционной функции, так и взаимная спектралі^ ная плотность двух реализаций представляет преобразование Фурье взаимной корреляционной функции. Поскольку взаимная корреляционная функция не обладает свойством четности, взаим ная спектральная плотность есть обычно комплексная величина
GxyU )~C xy( n - jQ xy(f), |
(1.47) |
где действительная часть Cxy(f) называется синфазной составляю щей, а мнимая часть Qxy(f) — квадратурной составляющей взаим ной спектральной плотности.
В функции частоты синфазную составляющую можно предста вить как отношение среднего произведения функций x(t) и y(t^ в узком интервале частот от / до / + А/ к ширине этого интервала.' Такое же определение можно дать квадратурной составляющей взаимной спектральной плотности, за тем исключением, что либо x(t), либо у(і), но не обе функции сразу сдвинуты во времени та ким образом, что составляющие с частотой / будут сдвинуты по
фазе на |
90°. |
|
Удобно выразить взаимную спектральную плотность в показа |
||
тельной |
форме: |
(1.48) |
|
Оху(П=\Оху(П\е-іѲ*Уи)\ |
где модуль I Gxy(f) | и аргумент Qxy(f) связаны с величинами Cxy(f) и ОхѵФ формулами
I Охв if)\ = Y C 2xy( f ) + Q \ y(n, |
(1 -49) |
|
(/)= arctg Qxyif) |
(1.50f |
|
|
Cxy (f) |
|
Приведем еще одно полезное |
соотношение |
(1.51) |
і О х Л я |
і к м / а д - |
При использовании взаимной спектральной плотности для ре шения физических задач целесообразно ввести в рассмотрение дей ствительную величину
„2 |
f f \ __ 1G x y (f ) 12 |
(1.52) |
||
Х У |
(f)-- |
G x |
( f ) G y ( f) |
|
|
|
|
которая называется функцией когерентности. В том случае, когда при некотором значении частоты y%y(f) = 0, говорят, что функции х(і) и у(() на данной частоте некогерентны, или, другими словами, некоррелированы. Если функции х{і) и y(t) статистически независимы, то при всех значениях частоты yly(f) = 0. В том слу^ чае, когда при всех значениях частоты ylyif) = 1, говорят, чта функции x(t) и y{t) полностью когерентны.
Типичный график взаимной спектральной плотности как функдин частоты (зависимость Gxy от /) для двух реализаций случай-
Основные характеристики физических процессов |
49 |
s ных процессов приведен на рис. 1.21. Этот график |
называется |
' взаимным спектром. Заметим, что график состоит из двух ча стей, определяющих модуль и фазу.
Применение. Взаимная спектральная плотность, как и взаимная корреляционная функция, используется во многих
•случаях. Некоторые примеры ее применения приведены ниже.
Измерение частотной характеристики. В этом случае исполь зуется важное соотношение, связывающее взаимную спектраль.
. ную плотность с осмовными^характеристиками исследуемых фи.
зических систем. Рассмотрим, например, электрическую цепь, обладающую частотной характеристикой #(/). Предположим, что напряжение на входе цепи является стационарным случайным процессом со спектральной плотностью Gx(f). Как отмечалось в подразд. 1.3.4, сигнал на выходе электрической цепи является стационарным .случайным процессом со спектром, определяемым формулой (1.39). Однако выполняется и другое важное соотно шение, а именно: взаимная спектральная_плотность сигналов на входе и выходе имеет вид
Gxy{f) = H if)Gx{f)- |
(1.53 |
. Эта формула позволяет определить полную частотную характери- - стику линейной системы по данным измерений взаимной спек тральной плотности из соотношения H(f) — Gxy(f)/Gx(f). Довери тельные интервалы, определяющие точность измерения #(/), можно найти, вычислив соответствующую функцию когерент-
, ности. Эти вспрссы будут детально рассмотрены в гл. 5 и 6.
4— * 2 4 4
50 Глава 1
Измерение времени задержки. В подразд. 1.4.2 эта задача при~^
водилась в качестве, примера использования взаимной корреля-' ционной функции. В данном случае фазовый угол Qxy(f), величина которого определяется взаимной спектральной плотностью, свя зывающей сигналы на входе и выходе некоторой системы, дает сдвиг фазы сигнала с частотой / при прохождении его через си стему. Следовательно, время задержки в системе сигнала с неко торой частотой / описывается формулой т = QxU(f)!2nf. Отметим, что по данным измерений взаимной спектральной плотности мож но вычислить время задержки как функцию частоты, тогда как взаимная корреляционная функция не дает такой возможности.
Теория линейного прогнозирования и фильтрации. Спектраль ная и взаимная спектральная плотности имеют много различных' приложений в теории линейного прогнозирования и фильтрации. Решение этих задач требует определения оптимального по некото рому критерию линейного фильтра, который позволяет передать и прогнозировать желаемый сигнал, исключив посторонний шум. Характеристики такого фильтра определяются различными спек трами и взаимными спектрами, связывающими сигнал на входе плюс шум с сигналом на выходе плюс шум.
Рассмотрим следующий пример из этой области. Допустим,, что необходимо получить оценку сигнала у(() путем линейного преобразования входного сигнала x(t). Положим далее, что сиг нал у(і) искажается посторонним шумом п((), который не зависит от x(t). Если пока не затрагивать вопроса о физической осущест вимости искомого фильтра, то оптимальная весовая функция,, отыскиваемая по критерию минимума среднего квадрата ошибки,.' имеет вид
Я у / х ( / ) = - % ( 1 . 5 4 >
где вертикальная черта в индексе у\х означает, что y(t) опреде ляется при заданном сигнале х(().
Минимальное значение среднего квадрата ошибки для опти мальной системы
СО
= |
(1.55) |
о |
|
где y\yif) — функция когерентности, |
связывающая х(() и y(t} |
и вычисляемая по формуле (1.52). Нетрудно видеть, что в том слу чае, когда когерентность на всех частотах равна единице, функ ция у(1) получается путем линейного преобразования функции x(J) с нулевым средним квадратом ошибки. При помощи соотно шения (1.55) на практике определяют, насколько близка некото рая реальная система к оптимальной.
Основные характеристики физических процессов |
61 |
^Упражнения
Определите, какие утверждения в упражнениях 1—4 верны. 1. Стационарный случайный процесс:
а) процесс с дискретным временем;
б) |
процесс с непрерывным временем; |
в) |
эргодический процесс; |
г) |
имеет не зависящие от времени осредненные по ансамблю |
характеристики; |
|
д) |
обладает тем свойством, что осредненные по времени и по |
.ансамблю характеристики одинаковы. 2. Эргодический случайный процесс:
а) процесс |
с дискретным |
временем; |
|
б) |
процесс |
с непрерывным временем; |
|
в) |
стационарен; |
времени осредненные по ансамблю |
|
г) |
имеет не зависящие от |
характеристики; д) обладает тем свойством, что осредненные по времени и по
ансамблю характеристики одинаковы.
3. Одну выборочную функцию можно использовать для полу чения всех статистических характеристик случайного процесса, ■если этот процесс
а) детерминированный; б) эргодический; в) стационарный;
г) обладает всеми этими свойствами.
4. Автокорреляционная функция стационарного случайного процесса
а) |
убывает с ростом |т |; |
б) |
зависит только от сдвига т; |
в) |
стремится к постоянной величине с ростом |т |; |
г) |
всегда не отрицательна. |
5.Как изменяются ответы на вопросы упражнения 4, если ■стационарный случайный процесс не содержит периодических ■составляющих?
6.Задан стационарный случайный процесс с нулевой спек тральной плотностью при / = 0. Докажите, что автокорреляцион
ная функция такого процесса удовлетворяет соотношению
|
|
|
с ° 1 |
|
|
|
|
|
I R (r)dr = |
0. |
|
|
|
|
—со і |
|
|
7. |
Задан стационарный случайный процесс с автокорреля |
||||
ционной функцией R(т) |
= А е-а Іт і, |
где |
0. Найдите |
||
а) |
среднее |
значение |
р; |
|
|
іб) |
среднее |
значение |
квадрата ¥ 2; |
|
|
4*
52 |
Глава 1 |
в) спектральную плотность G(/). |
|
8. Полагая, |
что случайный процесс, рассматриваемый в уп |
ражнении 7, является нормальным, запишите выражение для. плотности распределения процесса.
9. |
Запишите выражение для автокорреляционной функции |
а) |
гармонического процесса; |
б) |
периодического прямоугольного волнового процесса; |
в) |
периодического треугольного волнового процесса. |
10. |
Совместная плотность распределения двух случайных про |
цессов |
удовлетворяет соотношению р(х, у) = р(х) р(у). |
Укажите утверждения из перечисленных ниже, которые всегда верны:
а) взаимная корреляционная функция этих процессов равна нулю при всех значениях сдвига т;
б) взаимная спектральная плотность этих процессов равна нулю при всех значениях частоты /;
в) функция когерентности, связывающая процессы, равна нулю при всех значениях частоты /.
S 3
ГЛАВА 2
ЧАСТОТНЫЕ ХАРАКТЕРИСТИКИ
ФИЗИЧЕСКИХ СИСТЕМ
ч
Прежде чем подробно рассматривать вопросы измерения и анализа случайных физических процессов, следует уточнить не которые необходимые понятия и основные соотношения, характе ризующие динамический режим физических систем. В этой главе приводятся теоретические формулы, которые могут быть исполь зованы для описания частотных характеристик идеальных систем,, и на примерах простых физических систем даются пояснения ос новных идей.
2.1. Линейные системы с постоянными параметрами
Идеальной называется система с постоянными параметрами,. обладающая свойством линейности сигналов в двух определенных точках — на входе, или в точке приложения воздействия, и на выходе, или в точке определения реакции системы. Говорят, что система,имеет постоянные параметры, если все основные свойства ее инвариантны во времени. Например, простая пассивная элек трическая цепь является системой с постоянными параметрами, если сопротивления, емкости и индуктивности всех ее элементов не меняются во времени.' Система линейна, если частотные ха рактеристики обладают свойствами аддитивности и однородности.
.Понятие аддитивности означает, что реакция системы на сумму входных сигналов равна сумме реакций на каждый отдельно взя тый входной сигнал. Понятие однородное™ означает, что реак ция системы на любой сигнал, умноженный на некоторую по стоянную, равна этой постоянной, умноженной на реакцию систе мы на входной сигнал.
Запишем эти утверждения в аналитической форме. Пусть сим вол f(x) обозначает реакцию системы на входной сигнал х. Систе ма обладает свойством линейности, если для двух любых входных
сигналов Хх и хг и постоянной с справедливы соотношения |
|
f(xi + x2)= f(x 1) + f(x 2), |
(2.1а). |
f{cx) = cf{x). |
(2.16) |
Формула (2.1а) выражает свойство аддитивности, а |
формула |
(2.16) — свойство однородности. |
|
■54 |
Глава 2 |
Предположение о постоянстве параметров вполне приемлемо для многих физических систем, с которыми приходится иметь дело на практике. Например, обычно не наблюдается заметных измене ний основных характеристик электрических цепей или механиче ских устройств в пределах любых представляющих практический интерес интервалов времени. Такие системы встречаются, конеч но, далеко не всегда. Сопротивление электрического резистора может меняться вследствие сильного нагрева, а прочность соору жения также может меняться при повреждении, вызванном уста лостью металла под воздействием непрерывной вибрации. Кроме того, некоторые физические системы конструируют именно как
•системы с переменными параметрами. В качестве яркого примера ' можно привести системы связи, в которых используются электрон ные приборы. Однако примеры такого рода обычно представляют -собой весьма специфические частные случаи.
Большие ограничения накладывает предположение о линей ности реальных систем. При наличии экстремальных условий на входе частотные характеристики всех реальных физических си стем нелинейны. Например, при увеличении приложенного к электрическому конденсатору напряжения наблюдается резкая криволинейность характеристики и, следовательно, сила тока уже не будет прямо пропорциональна напряжению. Точно так же металлический трос при увеличении нагрузки в конце концов разрывается, и поэтому напряжение в тросе не будет пропорцио нально приложенной нагрузке. Еще более осложняет дело то об стоятельство, что нелинейность обычно возникает постепенно, а.' не обнаруживается в какой-либо момент. Например, соотношение между напряжением в тросе и приложенной к нему нагрузкой
•обнаруживает отклонение от линейности задолго до того, как происходит окончательный обрыв троса. Тем не менее, не рискуя допустить больших ошибок, частотные характеристики многих физических систем можно считать, по крайней мере в некотором ограниченном диапазоне приложенного на входе воздействия, линейными.
-2.2. Основные динамические свойства физических систем
Динамические свойства линейной системы с постоянными пара- "метрами можно охарактеризовать при помощи весовой функции' h{т), которая определяется как реакция системы в некоторый мо мент t на единичную импульсную функцию, поданную на вход системы в момент / — т. Полезность весовой функции как характе ристики системы связана со следующим. При любом входном
Частотные характеристики физических систем |
5S- |
g процессе x(t) процесс на выходе системы y(t) определяется ишпегра1лом свертки
СО |
|
у (і) = jh (x )x (t—x)dx. |
(2.2)’ |
— СО
Таким образом, выходной процесс y(t) представляет собой взве шенную (на бесконечном интервале времени) сумму всех значе ний входного процесса
Для того чтобы линейная система с постоянными параметрами была физически осуществимой, необходимо, чтобы система реаги
ровала только на прошлые |
значения сигнала. Это означает, |
что |
/г (х) |
=-0 при т < 0 . |
(2.3) |
Следовательно, для реальных физических систем нижний предел интегрирования в формуле (2.1) в действительности равен нулю,, а не минус бесконечности.
Говорят, что линейная система с постоянными параметрами- устойчива, если при произвольной допустимой ограниченной функ
ции на входе системы функция на выходе также |
является ограни |
||
ченной. На |
основании соотношения (2.2) можно записать |
||
|
|
СО |
|
|
Ы О І Н |
j* А СО* (*—*)* К |
|
|
|
—со |
|
|
со |
|
|
? |
< j |
|А(т) || х ( і—т) \dx. |
(2.4> |
—СО
Втом случае, когда функция х(і) на входе системы ограничена,, существует конечная константа А, такая, что
| jc( / ) K j4 при всех t. |
(2.5).- |
|
Из формулы (2.4) следует, что |
-ч |
|
00 |
|
|
J |
I h (т) I dx. |
(2.6> |
— СО
Поэтому если весовая функция h(x) линейной системы с постоян
ными параметрами |
абсолютно интегрируема, т. е. |
|
|
СО |
|
■f |
j jh(x) JdxC oo, |
(2.7)- |
—oo
то функция на выходе будет ограничена и система, следователь но, обладает свойством устойчивости.
'56 |
Глава 2 |
Линейную систему с постоянными параметрами можно охарак-с* теризовать также передаточной функцией Н(р), которая опреде ляется как преобразование Лапласа весовой функции /г(т):
w |
|
Я (р )= J h (x)e~ptdx, р=а-\- jb. |
(2.8) |
Критерий устойчивости линейной системы с постоянными пара-
.метрами в случае ее физической осуществимости принимает ин тересную форму, если связать его с передаточной функцией Н{р). В этом случае система является устойчивой, если функция Н(р) не имеет полюсов в правой полуплоскости комплексной пере--з
■менной р или |
на |
мнимой |
оси (т. е. при а > 0 полюсы отсутст |
вуют). Верно |
и |
обратное |
утверждение — система неустойчива, |
если функция Я(р) имеет хотя бы один полюс в правой полупло скости комплексной переменной р или на мнимой оси.
Важным свойством линейных систем с постоянными парамет рами является сохранение частотной структуры сигналов. Рас смотрим линейную систему с постоянными параметрами, обла дающую весовой функцией h{т). Из формулы (2.2) следует, что при произвольной форме входного процесса х(і) п-я производная по времени процесса на выходе y(t) имеет вид
dny (t) |
■= |
: (2-9) |
|
dt" |
|||
|
|
Предположим далее, что на вход системы подается гармоническое' колебание
X (t)~ X sin (2л// + Ѳ). |
(2.10) |
|
•Бторая производная функции x(t) |
|
|
d2x (0 |
= —4л2f*x(t). |
(2. 11) |
di2 |
||
К,ак видно из соотношения |
(2.9), вторая производная выходного |
|
.процесса y(t) |
|
|
d 2dty [3t ) |
|
(2.12) |
Таким образом, процесс на выходе y(t) также имеет форму сину(~ч соиды с той же частотой, что и частота процесса ’faa входе х(/): Отсюда следует, что при помощи линейной системы с постоянны ми параметрами нельзя преобразовывать частоту, а можно ме нять лишь амплитуду и фазу приложенного на входе процесса.
Частотные характеристики физических систем |
5Г |
ft.3. Частотные характеристики
Если линейная система с постоянными параметрами физически, осуществима и устойчива, то ее динамические свойства можно опи сать при помощи частотной характеристики H{f), которая оп ределяется как преобразование Фурье весовой функции h1*), т. е.
со |
|
ГЯ (/) = hJ(x)e~2ni^di. |
(2.13) |
L |
|
о
Отметим, что нижний предел интегрирования равен нулю, а не минус бесконечности, так как Ңт) = 0 при т < 0. Частотная характеристика представляет собой просто частный случай пере даточной функции, когда в показателе экспоненты р = а + jb,. а = 0 и b = 2л/. Для физически осуществимых устойчивых си стем передаточную функцию можно заменить частотной характе ристикой без потери информации.
Применяя преобразование Фурье к левой и правой частям фор мулы (2.2), получают важное соотношение для частотной характе ристики линейной системы с постоянными параметрами. Обозна чив преобразование Фурье процесса на входе системы х(і) через
X(J), а процесса на выходе у(і) — через_У(/), |
из формулы (2.2) |
можно найти, что |
|
Y ( / ) = # (f)X (/). |
(2.14) |
Таким образом, использовав частотную характеристику систе м ы и выполнив преобразование Фурье процессов на ее входе и 'выходе, можно свести интеграл свертки (2.2) к простому алгебраи ческому выражению (2.14).
Частотная характеристика в общем случае является комплекс ной величиной, которую удобно представить через ее модуль и аргумент. Для^этого следует переписать H(f) в показательной
форме |
|
[ Н( П=\ Н( П\ е~ІФІП. |
(2.15) |
шодуль \Н(Г)\ называется амплитудной частотной характери стикой системы, а аргумент ф{[) — фазовой частотной характе ристикой. Используя эти понятия, частотной характеристике можно дать простую физическую интерпретацию. Пусть на вход системы поступает гармоническое колебание с частотой /, сущест вующее, как мы предполагаем, на бесконечном интервале времени. Процесс на выходе системы, как показано в разд. 2.2, также будет Синусоидальным той же частоты. Отношение амплитуд процессов на выходе и входе системы определяет амплитудную частотную характеристику системы |Я (/)|, а сдвиг фаз процесса на выходе системы относительно процесса на входе определяет фазовую частотную характеристику 0(/).