Файл: Болдырев, В. С. Методы математической статистики в гидрографии и кораблевождении учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 24.10.2024

Просмотров: 100

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

рюі наиболее общий случай, когда в составе реализации находится несколько слагаемых•

y(i)=x(t) + 2(t)+z(t)

Здесь x ( t )

-

чисто случайная составляющая#

 

2 ( t )

-

периодическая (неслучайная) составляю-

 

 

 

щая #

 

 

 

 

z ( t )

- неслучайная

функция времени# она опре­

 

 

 

деляется и исключается из состава слу­

 

 

 

чайной функции сравнительно просто,

 

 

 

путем применения к^ (t)

операции нахок-

 

 

 

дения математического ожидания и по­

 

 

 

следующего

центрирования.

Следовательно,

не

нарушая

общности,

можно положить

z ( t) = 0

. Тогда

 

 

 

 

 

Пусть

 

%(t)=3c(t)+4(i)

 

 

 

 

 

 

 

 

 

 

2 ( t ) - A s i n ( o o t + f ) •

0 .2 8 )

Применим к

fp(t)

операцию нахождения

корреляционной

Функции:

 

 

 

 

 

 

 

К^{'с)-€іт -J-

s in { u jti'cP)sin(u>t'tcoT+<f>)clt ;

 

 

^

( Г )= — cos

U v ,

0 .2 9 ;

 

 

 

Таким

о б р а з о м ,

в

р е з у л ь т а т е

применения оп ерации н а ­

хож дения корреляц и он н ой функции получили н ек оторую

пер и оди ч еск ую

кривую

с ам плитудой - — .

При эт ом и н -

Т


формация и начальной фазе <f оказалась утерянной. В более общем случае, если периодическая составляющая представлена рядом,

 

 

 

П

 

 

 

 

 

? ( ^ =Yé CK5ir2(aJKt + CP«') 7

(4.30)

 

 

к=0

 

 

 

 

где

 

 

ZJT

К .

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

Первый член

этого

ряда

&0

может

быть

включен в

состав г ( і)

как

постоянное

Z

 

поэтому в даль­

слагаемое,

нейшем его

можно не

учитывать.

Применяя к функции g it)

оператор нахождения корреляционной функции, получим

*

Q2

 

/< (?) = %

- f c o s c O ^ .

(4.31)

к = І

Теперь применим оператор нахождения корреляционной функции к сумме случайной и периодической составляющих:

К^(Т)=&Л7 -J [*(0+z(£)][x(t+T)+?(t+r)]dt;

Т----’

г

П

К (T)=€im Y

c(t)x(t*T)ctt+

 

 

Т— е»о

I о

 

 

т

т

Л

 

* \ < t ) 4 ( t * T ) d t + ^ ( O y i t + T j o l t Г *

 

О

О

 

и

Первое и последнее слагаемые будут равны К (Т)

соответственно. Два других слагаемых представ­

ляют из себя взаимные корреляционные

функции мевду

чис­

230


той случайной и периодической составляющими. По условию, оба слагаемых являются независимыми. Поэтому их взаим­ ная корреляционная функция равна нулю. Итак,

K ^ Z ) = Kx i'Z ) tK y ( Z )

(4.32)

При больших Т А

О

а К?(7)

останется

незатухающей. Следовательно, при наличии в составеy ( t )

периодической составляющей K y {Z) будет

незатухающей

р2

?

 

функцией с амплитудами ‘■'Л

и частотами

<*>

2

 

к

Изложенный прием позволяет определить

наличие в со­

ставе реализации периодической составляющей.Полученная в результате обработки корреляционная функция, точнее ее периодическая часть, может быть в дальнейшем под­ вергнута гармоническому анализу. Итак, для того, чтобы убедиться в наличии или отсутствии в составе реализации периодической составляющей, достаточно найти корреля­ ционную функцию данной реализации. Если корреляционная (функция получится затухающей, то это значит, что перио­ дической составляющей практически не было. Существенное затруднение на практике представляет ограниченность длины корреляционной функции. Это ограничение возникает

за счет конструктивных особенностей большинства корреля­ торов, но его можно обойти, если обрабатывать реализа­ цию на ЭВМ. Однако при этом надо иметь в виду, что точ­ ность получения корреляционной функции падает по мере увеличения Z . Поэтому при больших Z небольшие ко­ лебания корреляционной функции свидетельствуют скорее всего о погрешностях обработки и недостатке исходных данных.

Для выделения периодической составляющей рассмотрим полученный график корреляционной функции при достаточно большом Z . Разлагая (Л.32) по Фурье,

231


= É ^ f c o s U u j Z ,

 

К= I

 

 

 

получим частоту tOk и амплитуды гармоник Ак .

Но для

получения функции

p t f )

нужно еще

найти углы <рк

, то

есть "привязать"

каждую гармонику

к оси времени.

Для

этого можно использовать либо метод последовательных

приближений, либо

метод Фуриха. Метод последовательных

приближений заключается

в подборе

положения каждой гар­

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

Метод Фуриха предполагает повторную обработку кривой (4 .32), которую запишем в нормированном виде»

 

*»с

2

К ( 7 ) - К

^

c o s c o U t .

 

k=l

 

Б результате обработки первое слагаемое дает

Ф > ~ £ ' 7 ) К ( ^ = н ------- *

Ѵк*І ^ К ( ^ Л + т с/т )Н ({< л)

Здесь с /г ~ конечный эл ем ент а р г у м ен т а 27 . Так как

! к { о ) \ > \ к к~У\ .

0*5 0

можно записать

N N

£ [ К ( Щ 2 £ К(кд+тСІТ)К(1<А)

-

>

-

 

 

А — /V

но K ( i; â ) —~ 0

при /Ѵ- ~

0 0 , следовательно, при не­

ограниченном возрастании N самое больное сдагемое стремится к нуле, а значит все остальные слагаемые и подавно стремятся к нуле. Таким образом, при достаточно

(г) , л

обращается

н нуль.

большом N

К ,жл Т )

 

х (х

 

 

 

Аналогично

 

2

 

 

N

Г

 

 

 

 

—- ccos*b

ДA іiаt o) ;

 

k*l

| j - '

2

7

 

Ar-/

L Ы

 

 

ири достаточно большом N толе обращаются в нуль.

Итак, повторная обработка на корреляторе (при пред­ варительном нормировании) дает окончательно

Очевидно, повторная обработка дает

Л

233


где

А

означает нормирование

п - і раз. Так к а к А < ^ .

те

при

п -й обработке все А^

уменьшатся. Иначе гово­

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

обработки ( Т ) практически будет представлять

собой косинусоиду частоты и) . Определив^ , построим

функцию, умножив y ( t ) на

coscot ■

90

g (t)cos соt = ЛГ(t) соs CO t+ ^

Ak s гп (Асоt + % ) соs со t

k~l

и найдем ее математическое

ожидание:

 

1 9 0

Al[g(t)coscot]=Al[f(t)coscof]*-A1 ^A^sintycot+^coscât

Lk=l

По условию, первое слагаемое будет равно нулю. Вы­ полнив указанные преобразования, получим

Следовательно,

А

м[у. ({)cosu>t\ = -rjrSinfi

и аналогично

A l [ y ( t ) s z n a ) t] = ^ c o s ff ,

откуда

(4.33)

Sin cot]

234

и

Аг-^\\P\\y(t)coscüt\ +М

В общем те случае

. ф = M[^(t)cosUu)t] ѣ

M \ g ( t ) s i n k c ü i \ ;

^=rf[M\jS(f)coskul?]J + \M{y(t)sin1<u)t]]2| * (Ч.34)

Полученные результаты ыохно применить, например, к исследованию морских течений. Если вектор скорости приливо-отливного течения описывает с постоянной угло­ вой скоростью эллипс, то в составе случайной функции скорости суммарного течения будет синусоидальная перио­ дическая составляющая. Действительно, параметрические

уравнения эллипса в системе координат/’# У , оси кото­ рой совпадают с главными осями эллипса, могут быть за­ писаны в виде:

х ' = а cos n t

^

ё sin n t

}

где

п - постоянная

угловая скорость. Перейдя к систе­

ме

координат Х О 'У »

получим:

 

x c o s ^ + <£Sincf >= a cos n t j

или

- X sin .

fy c o s

»■ é S i n n t

235