Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Автоматизация научных исследований..pdf
Скачиваний:
17
Добавлен:
15.11.2022
Размер:
6.27 Mб
Скачать
(2 .2 1 )

Учитывая соотношения (2.14), (2.15), получим

\lSx(f), 0 < / < оо

Gx(f) = - О, при других / ,

где

(2.17)

Из соотношений (2.16), (2.13) определим оценку односторонней спектральной плотности Gx( fk).

Имеем

-

2

^

а

N

(2.18)

 

rijNM ы

 

2

 

 

 

Определим А",(У^.)|

из соотношения (2.12), получим

 

 

^ ( Л ) = а д ) + УЙ(Л).

(2.19)

где

 

п . ..

. ^

 

 

 

 

( 2пкпл

 

 

Ш к ) = & Ц Хшcos

У

 

 

 

/ 1=0

 

 

 

 

W-1

 

2пкп

(2.20)

 

Q(.A) = - A'X x/„sin

N

 

 

/1=0

 

 

N

/=!,«.; к = 0,1,...,—.

Из соотношения (2.19) имеем

У

|*(Л )Г = ? ( / , ) + £</*); i = i,"j;* = o,i,...,y.

Таким образом, оценка Gx(fk) определяется с использованием соотношений (2.18), (2.21), (2.20).

2.4. Получение вторым способом сглаженной оценки спектральной плотности

Рассмотрим произвольную функцию v(r), такую, что Т не явля­

ется ее периодом, и пусть (рис. 2.5)

х(0 =

(2.22)

Функция v(/) задана на бесконечном интервале времени. Функция u(t) есть косинусоидальное сглаживающее окно Ханна, определяемое формулой

1 - cos'

, 0 < t < Т;

u(t) = <

\ Т )

О, в остальных случаях.

Преобразование Фурье функции м(/) имеет вид

где f x= Y J , а (/, (/) задана формулой вида

sin я/Г

-МГ

£/,(/) = Л

,

.. я/Г

(2.23)

(2.24)

(2.25)

UlW - f i ) = -T

s \ n n ( f - f x)T

e-Mr

 

 

n( f ~ f \ ) T

J

 

V ,(f + f x) =- T

sin7t(/ +f\)T

e-Mr

 

 

n ( f + f t)T

J

 

График функции |t/(/)| показан на рис. 2.6.

 

 

Найдем преобразование Фурье от функции x(t). Имеем

 

г

00

 

 

X ( f ) = \x{t)e-JU"dt = \U(a)V(f - a)da,

(2.27)

О

-оо

 

 

где V(f) - преобразование Фурье от функции v(/); а - частота, Гц.

X {fk) = - V{fk) - - F (/M) - - V(/M ),

(2.28)

где

 

V(fk)=\v(t)e-J2'k'ndt.

(2.29)

Предположим теперь, что в пределах каждой полосы частот шириной

Д/ = 1/7’ функция

v(t) ведет себя как ограниченный по частоте бе­

лый шум. Тогда

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

V"(f)

и V(g), вычисленное для любых частот /

и g из набора кА/ - к/Т,

имеет вид

 

 

 

 

M[F‘(/)F(g)] = {0,{ * g;

(2.30)

 

1Л /

= g,

 

где М - символ математического ожидания; V‘(f) - величина, ком­ плексно-сопряженная к величине V ( f ) .

С учетом (2.30) из (2.28) получаем

при любых f k =к/Т (к = 0,1,2,...,N / 2). Эта величина характеризует потери, вызванные использованием спекгрального окна Ханна (2.23) для оценки спектральной плотности. Поэтому при оценивании спек­ тральной плотности по формулам (2.13), (2.18) выражение (2.12) нуж­

но умножить на масштабный множитель В частности, при вы­

числении спектральной плотности по формулам (2.13), (2.18) и ис­ пользовании сглаживающего окна Ханна следует принимать

~

1я N~] ( !

 

2 7Ш

exp

- jln k n

(2.32)

 

V-> н=0

1-COS

----

 

N

 

 

 

N

 

 

 

 

где f k =k/(NAt)

(к = 0,1,...,N/2).

 

 

 

 

 

 

 

Сглаженная оценка Gx{fk) односторонней спектральной плотно­

сти запишется в виде

 

 

 

 

 

 

 

 

~

9

"</

~

•,

 

 

 

Д/

(2.33)

<?,</*) =—J - Z

 

 

к = 0,1,...,-.

 

ndNAt

,=)

 

 

 

 

 

2

 

Определим |x /(/t )| Из соотношения (2. 32) получим

 

 

X,(fk) = Pi(fk) + jQ,(fkl

(2.34)

где [3]

 

 

 

 

 

 

 

 

 

 

8 ЛМ

 

 

 

 

 

 

^2пкп^

 

^ ( / ,) = A ' ) / l ? . x" l 1" cos^

J

 

cosV N У

 

 

8 (V-1

 

 

2

----

^

. ( 2пкп\

(2.35)

Q,(fi ) = - A L - ' £ X,I

1 - COS

sinl------

 

VJ н=0

 

 

 

N

 

V N

 

 

N

 

 

 

 

 

 

 

i = \,nd;k = 0X...,

 

 

 

 

 

 

 

 

Из соотношения (2.34) имеем

 

 

 

 

 

 

 

 

X ,(Л)Г

- Л2(Л ) + й ’(Л); /

 

 

 

=

у.

(2.36)

Таким образом, оценка Gx{fk) определяется с использованием соотношений (2.33), (2.36), (2.35).