
книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов
.pdf68 Г лапа 2
процесс на входе системы есть изменение разности потенциалов,
как показано на рис. 2.6, |
где С — емкость, |
R — сопротивление, |
|||||||
L — индуктивность, |
e(t) — приложенное напряжение, |
і(() |
— ре-' |
||||||
|
|
зультирующий |
процесс — сила |
то |
|||||
L |
Ці) |
ка. Следует напомнить, |
что |
і{() = |
|||||
|
|
= d q ( i ) / d t , |
где |
q ( t ) — заряд. |
|
|
|||
|
|
Положим, что в качестве вход |
|||||||
|
|
ного процесса задается приложенное |
|||||||
|
|
напряжение; в качестве процесса на |
|||||||
|
|
выходе системы выбран результи |
|||||||
|
|
рующий заряд. |
Как и в случае рас |
||||||
с |
|
смотренных в |
подразд. 2.4.1 меха |
||||||
Р и с . 2.6. Электрическая сн- |
нических |
систем, |
для |
того |
чтобы ' |
||||
найти соответствующую |
этому |
слу |
|||||||
стема с колебаниями |
напря |
чаю частотную |
характеристику, |
не |
|||||
жения на входе. |
|
обходимо |
прежде |
всего |
получить |
||||
|
|
дифференциальное |
уравнение, |
|
опи |
сывающее данную систему. Как следует из основного закона тео рии электрических цепей, сумма всех значений падения напря жения в элементах цепи равна нулю, т. е.
e(i) + ec(t) + eR(l) + eL(t)=Q, |
(2.39) |
|
где |
|
|
ес { і ) = — ^-<?(0— падение напряжения на емкости, |
(2.39а) |
|
eR(t)= —Rq(t) —падение напряжения |
|
> |
на сопротивлении, |
||
|
|
(2.396) |
eL(t) = —Lq (t)—падение напряжения |
на индуктивности. |
|
|
|
(2.39в) |
Отсюда находим дифференциальное уравнение, описывающее систему
Lq(t) + Rq (t) + ~ q ( t ) = e(t). |
(2.40) |
Отметим, что между уравнением (2.40) и уравнением движения (2.19) механической системы с вынуждающей силой на входе су ществует аналогия. Применяя описанную в подразд. 2.4.1 мето дику анализа, сразу получаем частотную характеристику даннойщ
простой электрической |
системы |
^ |
|
Н (/) - ,= |
Л |
(2n/)2L + 2я//Я |
(2.41) |
с |
Частотные характеристики физических систем |
69 |
І£де индекс е — q обозначает, что частотная характеристика H[f) связывает напряжение на входе системы с зарядом на выходе. Заметим, что величина H(J)e_q имеет размерность кулои/вольт.
График функции H{j)e_Qимеет тот же вид, что и график частот ной характеристики на рис. 2.3, описывающей механиче скую систему. При этом коэффициент затухания С и собственная частота незатухающих колебаний /„ в электрической цепи опреде ляются равенствами
(2.42а)
(2.426)
fn 2л V ТС"
Можно показать, что между механической и электрической систе мами существует полная аналогия, о чем свидетельствуют данные табл. 2.2.
Таблица 2.2
Понятия-аналоги механических и электрических систем
В х о д н о й п р о ц е с с
Выходной процесс
Постоянные параметры
Э л е к т р и ч е с к а я с и с т е м а :
в ход н о й |
п р о ц е с с — н а п р я ж е н и е |
||
|
Н а п р я ж е н и е е U) |
||
Заряд q |
(/) |
|
|
Сила |
|
|
|
|
тока і 1) = dqjdl |
||
Индуктипность( |
L |
Сопротивление R
Емкость С
М е х а н и ч е с к а я с и с т е м а : в хо д н о й п р о ц е с с — в ы н у ж
д а ю щ а я |
с и л а |
Г и л а |
F (О |
Смещение у (/) Скорость V (і) = dyjdt
Масса m
Коэффициент торможения
Модуль упругости 1 Jk
Чаще используются частотные характеристики электрических систем, связывающие напряжение как входной процесс с силой тока на выходе
Я ( / ) „ , = R + i |
1 |
(2.43) ' |
2л[С |
где H(f)e_i имеет размерность ампер/вольт. Функция, обратная величине (2.43), которую можно обозначить через #(/)£_е, называется_импедансом. Импеданс
tf(fli-«=Ä + / ( 2 n / L ~ 2 ^ с ). |
(2.44) |
Заметим,"что, пользуясь данными табл. 2.2, можно записать ана лог уравнения (2.44) для механической системы
і
. |
Я (/)^ f= с + / ( 2я fm----- |
~ ~ у |
(2.45) |
о
Таблица 2.3
Параметры-аналоги некоторых физических систем
С и с т е м а |
П р о ц е с с па в х о д е |
П р о ц е с с на в ы хо д е |
П о ст о я н н ы е п ар ам ет р ы |
Электрическая
Механическая (воз- вратно-поступатель ное движение)
Механическая (враща тельное движение)
Акустическая
Термическая
Магнитная
Напряжение |
Сила тока |
Сила |
Скорость |
Вращающий мо |
Угловая скорость |
мент |
|
Давление |
Скорость частиц |
Температура |
Поток тепла |
Магнитодвижущая Магнитный поток сила
Индуктивность |
Сопротивление |
Емкость |
Масса |
Коэффициент тормо |
Модуль упругости |
|
жения |
|
Момент инерции |
Угловой коэффициент |
Угловой коэффициент |
|
торможения |
упругой деформа |
|
|
ции |
Коэффициент инер |
Акустический импе |
Акустическая емкость |
ции (акустичес |
данс |
|
кая масса) |
|
|
|
Температурное сопро |
Термоемкость |
|
тивление |
|
|
Магнитное сопротив |
|
|
ление |
|
------- Л /----------------------------------------
Частотные характеристики физических систем |
71 |
^Функцию (2.45) часто называют механическим импедансом ввиду того, что она аналогична более употребительному термину «элек трический импеданс».
2.4.3. Другие системы
Аналитический метод решения, использованный в подразд. 2.4.1, позволяет, по крайней мере теоретически, получить частотную характеристику любой устойчивой строго определенной физической системы (т. е. любой системы, динамический режим работы которой можно описать с помощью точного уравнения). ^Кроме того, различные физические системы могут обладать ана логичными параметрами, как показано на примере механической
и электрической систем в подразд. 2.4.2. Перечень параметрованалогов некоторых обычных физических систем приведен
втабл. 2.3.
2.5.Практические соображения
Вразд. 2.4 был описан аналитический метод определения ча стотных характеристик физических систем. Для того чтобы упро стить выкладки и дать общее представление об идее метода рассматривались лишь самые простые механические системы. На основании этих примеров не следует делать вывод, что аналитиче ский путь определения частотных характеристик для физических
^систем всегда столь прост.
Рассмотрим, например, такую механическую систему, как сплошная упругая конструкция, параметры которой (масса, коэф фициент торможения, коэффициент жесткости) не сосредоточены в точке, как предполагалось при рассмотрении примеров в под разд. 2.4.1, а распределены некоторым образом в пространстве. Можно назвать различные входные и выходные процессы такой механической системы, представляющие интерес для исследова теля. Частотная характеристика для каждой комбинации про цессов на входе и выходе обычно имеет много пиков, соответст вующих различным значениям резонансной частоты, в отличие от рассмотренной в подразд. 2.4.1 системы с единственным значением этой частоты. Частотные характеристики сравнительно несложных сплошных конструкций типа балок и плит можно найти с доста точной степенью точности аналитическим путем. Однако в случае более сложных конструкций необходимо затратить много усилий, ^тобы получить хотя бы грубое представление об искомой частот ной характеристике.
По изложенным выше причинам для определения или по край ней мере для проверки формы частотных характеристик реальных физических систем используются эмпирические методы. Наиболее
72 Глапа 2
простой эмпирический метод состоит в том, что к входу системѣ? прилагается гармонический процесс н измеряется амплитуда и фаза процесса на выходе при изменении частоты гармоники. Как показано в разд. 2.3, отношение амплитуд процессов на выходе и входе системы дает амплитудную частотную характеристику, а сдвиг фаз процесса на выходе относительно входного процесса определяет фазовую частотную характеристику.
К сожалению, встречаются и такие случаи, когда трудно моде лировать сигнал на входе исследуемой системы гармонической функцией. Рассмотрим, например, задачу определения динамиче ской реакции корабля на создаваемые волнами нагрузки, само лета — на турбулентные порывы ветра или строительной конст* рукции на нагрузки, создаваемые избыточным акустическим давле нием. Ввиду больших размеров исследуемых систем и сложного характера распределения приложенных к ним нагрузок весьма трудно найти частотные характеристики эмпирическим путем, подавая на вход системы калиброванный гармонический про цесс. В таких случаях можно получить ценные сведения, прила гая моделируемую нагрузку в виде гармонического процесса к моделям, уменьшенным до приемлемого размера. С другой сто роны, иногда более целесообразно вычислять искомые частотные характеристики по результатам измерений процессов на входе и выходе физических систем, находящихся в естественных для них условиях. В естественных же условиях перечисленные выше системы подвергаются воздействиям, носящим характер случай ных или переходных процессов, но никак не гармонически^#, Однако даже и в таких случаях можно определить частотные ха рактеристики путем анализа процессов на выходе системы при использовании в качестве входов переходных или случайных процессов. Эти вопросы детально рассмотрены в гл. 5 и 10.
Упражнения
1.Запишите уравнение, выражающее условие линейности физической системы.
2.Найдите, каким из перечисленных ниже требований должна
удовлетворять физическая система, чтобы ее можно было описать
однозначной весовой |
функцией h(т): |
|
||
а) |
постоянство |
параметров; б) линейность; в) физическ |
||
осуществимость; |
г) |
устойчивость. |
|
|
3. Докажите, что физическая система, обладающая весовой |
||||
функцией |
h(x) = |
Ае?х, неустойчива при а > 0. |
#4 |
4.Найдите весовую функцию изображенной на рис. 2.2 ме ханической системы с вынуждающей силой на входе.
5.Найдите частотную характеристику физической системы, обладающей весовой функцией h(x) = Ае~ах, где а > 0.
Частотные характеристики физических систем |
73 |
6. Пусть изображенная на рис. 2.1 механическая |
система |
'характеризуется следующими параметрами: коэффициент жест
кости |
пружины |
/е = |
2000 |
кг/с2, |
коэффициент торможения |
с = 40 кг/с, масса in = |
2 кг (в единицах СИ). |
||||
Вычислите |
частоту |
незатухающих колебаний /я; |
|||
а) |
собственную |
||||
б) коэффициент затухания £; |
|
||||
в) |
резонансную частоту /г; |
частотной характеристики |
|||
г) |
пиковое значение |
амплитудной |
\ r n \ -
7. Пусть масса изображенной на рис. 2.1 механической си стемы выведена из состояния равновесия, а затем освобождена. ^Докажите, что при возникающих колебаниях интервалы времени между прохождениями положения равновесия составляют
где £2< 1,0.
8. Докажите, что резонансная частота изображенной на рис. 2.2 механической системы с вынуждающей силой на входе составляет, согласно уравнению (2.25),
fr = fnV 1 -2 £ 2 , где С2 < 0 ,5 .
9. Полоса пропускания по уровню половинной энергии ре зонансной физической системы определяется соотношением Вг =
= f t — /і. где
і |
|
|
I Н (/l) І 2=1 Я (/a) ! 2= 4 “ I Я (/r) I 2- |
|
|
Докажите, |
что для изображенной на рис. |
2.2 механической систе |
|||
мы с вынуждающей силой на входе при малых значениях |
I ве |
||||
личина |
Вг ~ |
2£/,. |
аналога изображенной |
||
10. |
Составьте схему электрического |
||||
на рис. |
2.4 |
механической системы со смещением основания на |
|||
входе. Выразите значения параметров-аналогов цепи через |
пара |
||||
метры механической системы k, с и ш. |
|
|
МАТЕМАТИЧЕСКИЕ ОСНОВЫ АНАЛИЗА
СТАЦИОНАРНЫХ СЛУЧАЙНЫХ ПРОЦЕССОВ
В этой главе, которая представляет собой развитие гл. 1, кратко излагаются математические основы методов анализа стационарных случайных процессов. Приводятся элементарные и более сложные положения теории случайных процессов, начиная от исходных принципов, которые составляют необходимую осно ву для практических приложений теории к задачам анализа. Описаны важнейшие вероятностные свойства случайных величин, стационарных и эргодических случайных процессов, гауссовских случайных процессов, приведены сведения о линейных преобра зованиях и теоремы о дискретном представлении процессов.
3.1.Вероятностные основы анализа случайных величин
3.1.1. Случайная величина
В основе теории вероятностей лежит понятие множества/ определяемого как совокупность объектов (называемых также точками или элементами), относительно каждого из которых можно судить, принадлежит он данному множеству или не при надлежит. В частности, возможные исходы опыта (или результа ты измерения) образуют множество точек, называемое выбороч ным пространством. Эти точки можно различными способами группировать в множества, называемые событиями, и при со ответствующих условиях возможно поставить в соответствие каж дому событию определенную вероятность. Эти вероятности всегда принадлежат отрезку 10,1], причем вероятность невозможного события равна нулю, вероятность достоверного события — еди нице. Выборочные пространства могут быть конечными или бесконечными.
Рассмотрим выборочное пространство точек, представляющих возможные исходы некоторого опыта (или результаты измерения). Случайной величиной x(k) называется функция множества, опре-? деленная для точек k выборочного пространства, т. е. случайная величина x(k) есть действительное число, принимающее значения от — оо до + оо, соответствующие каждой возможной выбороч-
Математические основы анализа случайных процессов |
75 |
^ной точке /г. Иначе говоря, случайный исход опыта с индексом k '■.может быть представлен действительным числом х(/г), называе мым случайной величиной. Пусть задано множество всех возмож ных исходов опыта. Из подмножеств этого множества можно образовать вполне аддитивный класс множеств, и каждому собы
тию, т. |
е. подмножеству из этого класса, может |
быть |
поставлена |
|||
в соответствие некоторая |
вероятность. |
|
|
|
||
Пусть x{k) — некоторая случайная величина. Тогда для любо |
||||||
го фиксированного числа х случайное событие x{k) < |
х опреде |
|||||
ляется |
как множество возможных исходов |
k, |
для |
которых |
||
x(k) |
< |
X. Пользуясь теоретико-множественной терминологией тео |
||||
рии |
вероятностей, можно |
определить функцию |
распределения |
Р(х) как вероятность, которую можно поставить в соответствие множеству точек /г, удовлетворяющих заданному неравенству x(k) < X. Заметим, что множество точек k, для которых x(k) ^ х, представляет собой подмножество совокупности всех точек k,
удовлетворяющих неравенству x(k) ^ оо. Можно |
записать |
|
Р (х )= P[x(k)^x]. |
(3.1) |
|
Очевидно, |
|
|
Р(а) < Р (Ъ), |
а < Ъ, |
(3.2) |
Р( — оо)=0, |
Р(оо)= 1. |
(3.3) |
Если случайная величина может принимать любые значения в некотором интервале (что и будет предполагаться в дальнейшем), то можно определить одномерную плотность распределения соот
ношением |
|
|
|
р ( х )= і:т |
1-РГ*<*1*Х* + Л*]-_ |
(3.4) |
|
Заметим, что |
р (х)^г 0, |
(3.5) |
|
|
|||
|
j" р (x)dx = |
1, |
(3.6) |
— СО |
|
|
|
Р(х)= |
л: |
|
(3.7) |
Jp(E)d6, |
^ J ^ = p(x). |
Для обобщения на случай дискретной случайной величины, как в примере 3.1, в плотность распределения вводятся дельта-функ- дии.
П р и м е р 3.1. Дискретное распределение. Пусть производится бросание монеты; очевидно, здесь возможны два исхода — появ ление герба или решетки, причем вероятность этих событий оди накова и- равна 1/2. Случайная величина x(k) может принимать.
76 |
Глава 3 |
только два значения — х (герб) и х (решетка), которым можно* поставить в соответствие произвольные действительные числа) Пусть, например, х (герб) = а и х (решетка) = Ь, где а и Ь — действительные числа и b > а. Тогда функция распределения определяется следующим образом:
О при X < а,
Р(х) = у- при а < X < Ь,
1 Г.ГЛИ
Плотность распределения в этом случае имеет вид
|
р (х) = |
б (х—а) + - у б (х— Ь), |
где б (х — а) |
и б (х — Ь) — дельта-функции Дирака. |
|
П р и м е р |
3.2. Равномерное (прямоугольное) распределение. |
Пусть опыт состоит в случайном выборе точки в замкнутом интер вале [а, Ь]. Непрерывную случайную величину x(k) можно опре делить как числовое значение выбранной точки. Соответствующая функция распределения имеет вид
|
' 0 |
при X < а, |
|
Р (х)= | |
° при а <; X < |
Ь, |
|
|
, 1 |
при X > Ь. |
.> |
Следовательно, плотность |
распределения |
|
|
. . |
( (Ь ~ а )~1 при а < х< Ь , |
||
Р(х) = { п |
|
|
|
|
(0 при других X. |
|
|
На этом пример 3.2 заканчивается. |
принимать значения |
||
Пусть случайная |
величина x(k) может |
в интервале (—• оо, оо). Среднее значение (математическое ожи дание) случайной величины х(/г) определяется предельным пере ходом от произведения всех возможных значений x(k) на соот
ветствующие им вероятности |
к интегралу |
|
|
|
|
СО |
|
|
М [х (/г)] = |
Гхр (x)dx —[А*, |
(3.8) |
где М [ |
] есть символ математического ожидания |
или усредне |
|
ния по k. |
Аналогичным образом математическое ожидание любой |
Математические основы анализа случайных процессов |
77 |
■\
P-действительной однозначной непрерывной функции g(x) случай ной величины хЦі) определяется в виде
|
|
СО |
|
М[£(*(£))]= j g(x)p(x)dx, |
(3.9) |
||
|
|
— СО |
|
где р(х) — плотность |
распределения случайной величины x(k)1. |
||
В частности, положив |
g(x) = |
х2, найдем среднее значение |
квад |
рата случайной величины x(k) |
|
||
|
|
СО |
|
M[x2(k)]= |
Г x*p(x)dx = 4f%. |
(3.10) |
Дисперсия случайной величины x{k) определяется как сред ний квадрат отклонения случайной величины x{k) от ее математи ческого ожидания. Здесь g(x) = (х — рД2 и
СО
М [(X (k) - р,.)2] = J ( Х - Р х)*р (x)dx = ¥ ] - р .2= а ; . (3.11)
— СО
По определению стандартное отклонение ох случайной ве личины x(k) есть положительное значение квадратного корня из дисперсии. Стандартное отклонение измеряется в тех же едини-
^ цах, что и среднее значение.
Плотность распределения функции случайной величи ны. Пусть x{k) — случайная величина с плотностью распределения р{х), g{x) — некоторая действительная однозначная непрерыв ная функция X. Если производная dg/dx существует и отлийиа от нуля, то можно определить плотность распределения p(g), соот ветствующую случайной величине g(x(k)) = g(k). Для обобщения на случай многозначной обратной функции допустим, что обрат-
*) Следует |
сделать такое замечание |
Пусть g4Ul |
= |
g{x\ |
при £(•*) |
|||
> 0; |
g ‘(x) — 0 |
пои g(x) $ 0: g ' {х) = 0 |
при |
gix) |
0; g |
(л) |
= р(д) пj>и |
|
g{x) |
< 0. |
Математическое ожидание неличины |
g{x(k)) |
определяется толь |
||||
ко тогда, |
когда |
хотя бы один из двух интегралов |
|
|
|
л>
со |
со |
f £+ W p (x)dx, |
j V (*)р (x)dx |
— CO |
— on |
конечен. Если оба они конечны, то математическое ожидание конечно; если первый бесконечен, а второй конечен, M|g(jr(ft))| = <х ; если первый конечен, а второй бесконечен, М [£{а(А))] = —оо.— Прим. ред.