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

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

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

68 Г лапа 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 V ТС"

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

Таблица 2.2

Понятия-аналоги механических и электрических систем

В х о д н о й п р о ц е с с

Выходной процесс

Постоянные параметры

Э л е к т р и ч е с к а я с и с т е м а :

в ход н о й

п р о ц е с с — н а п р я ж е н и е

 

Н а п р я ж е н и е е U)

Заряд q

(/)

 

Сила

 

 

 

тока і 1) = dqjdl

Индуктипность(

L

Сопротивление R

Емкость С

М е х а н и ч е с к а я с и с т е м а : в хо д н о й п р о ц е с с — в ы н у ж ­

д а ю щ а я

с и л а

Г и л а

F

Смещение у (/) Скорость V (і) = dyjdt

Масса m

Коэффициент торможения

Модуль упругости 1 Jk

Чаще используются частотные характеристики электрических систем, связывающие напряжение как входной процесс с силой тока на выходе

Я ( / ) „ , = R + i

1

(2.43) '

где 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))| = <х ; если первый конечен, а второй бесконечен, М [£{а(А))] = —оо.— Прим. ред.

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