Файл: Болдырев, В. С. Методы математической статистики в гидрографии и кораблевождении учеб. пособие.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