Чи повністю описує функція автокореляції стохастичний процес?


31

Чи стохастичний процес повністю описаний його функцією автокореляції?

Якщо ні, то які додаткові властивості будуть потрібні?

Відповіді:


44

Що розуміється під повним описом стохастичного процесу? Ну, математично, стохастичний процес - це сукупність випадкових змінних, по одній на кожен момент в наборі індексів , де зазвичай - вся реальна лінія або позитивна реальна лінія, і повний опис означає, що для кожного цілого числа та часових моментів ми знаємо (спільні) розподіли в випадкові величини , ,{X(t):tT}t TTn1nt1,t2,,tnTnX(t1)X(t2),X(tn). Це величезна кількість інформації: нам потрібно знати CDF для кожного моменту часу , (двовимірний) спільний CDF і для всіх варіантів моменту часу і , (тривимірні) CDF , і , і т.д. і т.д. і т.д. і т.д.X(t)tX(t1)X(t2)t1t2X(t1)X(t2)X(t3)

Таким чином, природно, люди розшукували простіші описи та більш обмежувальні моделі. Одне спрощення відбувається, коли процес інваріантний зміні часу виникнення. Що це означає, що це

  • Усі випадкові змінні в процесі мають однакові CDF: для всіх .FX(t1)(x)=FX(t2)(x)t1,t2
  • Будь-які дві випадкові величини, розділені на деяку задану кількість часу, мають той самий спільний CDF, як і будь-яка інша пара випадкових змінних, розділених однаковою кількістю часу. Наприклад, випадкові величини і розділені на секунди, як і випадкові змінні і , і таким чиномX(t1)X(t1+τ)τX(t2)X(t2+τ)FX(t1),X(t1+τ)(x,y)=FX(t2),X(t2+τ)(x,y)
  • Будь-які три випадкові величини , , розміщені і мають однаковий спільний CDF, як , , які також розташовані та один від одного,X(t1)X(t1+τ1)X(t1+τ1+τ2)τ1τ2X(t2)X(t2+τ1)X(t2+τ1+τ2)τ1τ2
  • і так далі для всіх багатовимірних CDF. Дивіться, наприклад, відповідь Пітера К. для детальної інформації про багатовимірний випадок.

Ефективно, імовірнісні описи випадкового процесу не залежать від того, що ми назвати початком по осі часу: зміщення всіх моментів часу на деяку фіксовану суму на дає однаковий імовірнісний опис випадкових змінних. Ця властивість називається стаціонарністю суворого сенсу, а випадковий процес, який користується цією властивістю, називається строго стаціонарним випадковим процесом або, простіше кажучи, стаціонарним випадковим процесом. t1,t2,,tnτt1+τ,t2+τ,,tn+τ

Зауважимо, що сувора стаціонарність сама по собі не вимагає особливої ​​форми ФДВ. Наприклад, це не говорить про те, що всі змінні є гауссовими.

Прикметник суворо підказує, що можна визначити більш слабку форму стаціонарності. Якщо -порядковий спільний CDF такий же, як -порядник CDF для всіх варіантів і , тоді, як кажуть, випадковий процес стаціонарний на замовлення і називається -постійним стаціонарним випадковим процесом. Зауважимо, що -стаціонарний стаціонарний випадковий процес також є стаціонарним, щоб замовити для кожного додатногоNthX(t1),X(t2),,X(tN)NthX(t1+τ),X(t2+τ),,X(tN+τ)t1,t2,,tNτNNthNthnn<N . (Це тому, що -спільний CDF-межі є межею -порядок CDF як підходу аргументів : узагальнення ). Те А строго стаціонарний випадковий процес є випадковим процесом , який знаходиться в нерухомому стані на всі замовлення .nthNthNnFX(x)=limyFX,Y(x,y)N

Якщо випадковий процес є стаціонарним (принаймні) порядку , то всі мають однаковий розподіл, і так, припускаючи, що середнє значення існує, однаковий для всіх . Аналогічно, однаковий для всіх і називається силою процесу. Всі фізичні процеси мають кінцеву силу, і тому прийнято вважати, що в такому випадку, і особливо в старій технічній літературі, процес називають процесом другого порядку . Вибір імені невдалий, оскільки він викликає плутанину з другим порядком 1X(t)E[X(t)]=μtE[(X(t))2]tE[(X(t))2]<стаціонарність (пор. ця моя відповідь на stats.SE ), і тому тут ми назвемо процес, для якого є кінцевим для всіх (чи є чи ні - константа) як процес з кінцевою потужністю і уникайте цієї плутанини. Але зауважте ще раз цеE[(X(t))2]tE[(X(t))2]

стаціонарний процес першого порядку не повинен бути процесом обмеженої потужності.

Розглянемо випадковий процес, який є стаціонарним для замовлення . Тепер, оскільки спільний розподіл і збігається з функцією спільного розподілу і , і значення залежить лише від . Ці очікування є скінченними для процесу з кінцевою потужністю, і їх значення називається функцією автокореляції процесу: - це функція , часу поділ випадкових величин і , і не залежить від2X(t1)X(t1+τ)X(t2)X(t2+τ)E[X(t1)X(t1+τ)]=E[X(t2)X(t2+τ)]τRX(τ)=E[X(t)X(t+τ)]τX(t)X(t+τ)tзовсім. Зауважте також, що і тому функція автокореляції є рівною функцією її аргументу.

E[X(t)X(t+τ)]=E[X(t+τ)X(t)]=E[X(t+τ)X(t+ττ)]=RX(τ),

Стаціонарний випадковий процес кінцевої потужності другого порядку має такі властивості

  1. Його середнє значення - константаE[X(t)]
  2. Його функція автокореляції є функцією , часовим поділом випадкових величин і , і робить взагалі не залежать від .RX(τ)=E[X(t)X(t+τ)]τX(t)X(t+τ)t

Припущення про стаціонарність певною мірою спрощує опис випадкового процесу, але для інженерів та статистиків, зацікавлених у побудові моделей за експериментальними даними, оцінка всіх цих CDF є нетривіальним завданням, особливо коли є лише відрізок одного вибіркового шляху (або реалізація) на якій можна проводити вимірювання. Два вимірювання, які порівняно легко зробити (оскільки інженер вже має необхідні інструменти на своєму робочому столі (або програми в MATLAB / Python / Octave / C ++ у своїй бібліотеці програм), є значення постійного струму з і автокореляційної функціїx(t) 1T0Tx(t)dtx(t)Rx(τ)=1T0Tx(t)x(t+τ)dt(або його перетворення Фур'є, спектр потужності ). Вважаючи ці вимірювання як оцінку середнього значення та функції автокореляції процесу з кінцевою потужністю, веде до дуже корисної моделі, про яку ми поговоримо далі.x(t)


Випадковим процесом з кінцевою потужністю називають процес стаціонарного широкого сенсу (WSS) (також слабко стаціонарний випадковий процес, який, на щастя, також має той самий ініціалізм WSS), якщо він має постійне середнє значення і функцію його автокореляції залежить лише від різниці в часі (або ).RX(t1,t2)=E[X(t1)X(t2)]t1t2t2t1

Зауважте, що визначення нічого не говорить про CDF випадкових змінних, що містять процес; це повністю обмеження моментів першого та другого порядку випадкових змінних. Звичайно, стаціонарний (або -стаціонарний (для ) або стаціонарний) випадковий процес є випадковим процесом WSS, але навпаки це не повинно бути правдою.NthN>2

Процес WSS не повинен бути стаціонарним для будь-якого замовлення.

Розглянемо, наприклад, випадковий процес де приймає чотири однаково вірогідні значення і . (Не лякайтеся: чотири можливі вибіркові шляхи цього випадкового процесу - це лише чотири сигнали сигналу сигналу QPSK). Зауважимо, що кожен - дискретна випадкова величина, яка, як правило, приймає чотири однаково вірогідні значення і , Легко побачити, що взагалі і{X(t):X(t)=cos(t+Θ),<t<}Θ0,π/2,π3π/2X(t)cos(t),cos(t+π/2)=sin(t),cos(t+π)=cos(t)cos(t+3π/2)=sin(t)X(t)X(s)мають різні розподіли, і тому процес навіть не є стаціонарним. З іншого боку, для кожного той час, як Коротше кажучи, процес має нульове середнє значення, а його функція автокореляції залежить лише від різниці в часі , і тому процес є широким сенсом стаціонарним. Але він не є стаціонарним першого порядку, і тому не може бути стаціонарним і для вищих порядків.

E[X(t)]=14cos(t)+14(sin(t))+14(cos(t))+14sin(t)=0
t
E[X(t)X(s)]=14[cos(t)cos(s)+(cos(t))(cos(s))+sin(t)sin(s)+(sin(t))(sin(s))]=12[cos(t)cos(s)+sin(t)sin(s)]=12cos(ts).
ts

Навіть для WSS-процесів, які є стаціонарними (або строго стаціонарними) випадковими процесами, мало що можна сказати про конкретні форми розподілів випадкових змінних. Коротко,

Процес WSS не обов'язково є стаціонарним (у будь-якому порядку), а середня і автокореляційна функція процесу WSS недостатня для того, щоб дати повний статистичний опис процесу.

Нарешті, припустимо , що стохастичний процес передбачається бути гаусовим процесом ( «довести» це з достатнім ступенем впевненості не є тривіальним завданням). Це означає, що для кожного , - гауссова випадкова величина і для всіх натуральних чисел та вибір моментів , , , випадкових величин , , це спільно випадкові величини Гаусса . Зараз спільна функція щільності Гаусса повністюtX(t)n2nt1t2,tnNX(t1)X(t2),X(tn)визначається за допомогою засобів, дисперсій та коваріацій випадкових змінних, і в цьому випадку, знаючи середню функцію (вона не повинна бути константою, як це потрібно для широкого сенсу -стаціонарність) і функція автокореляції для всіх (вона не повинна залежати лише від як це потрібно для широкого стаціонарності) достатньо, щоб повністю визначити статистику процесу.μX(t)=E[X(t)]RX(t1,t2)=E[X(t1)X(t2)]t1,t2t1t2

Якщо гауссовский процес є процесом WSS, то це також строго стаціонарний гауссовский процес. На щастя для інженерів та процесорів сигналів, багато процесів фізичного шуму можуть бути добре змодельовані як WSS-гауссові процеси (і, отже, суворо стаціонарні процеси), так що експериментальне спостереження за функцією автокореляції легко забезпечує всі спільні розподіли. Крім того, оскільки Гауссові процеси зберігають свій гауссовий характер при проходженні через лінійні системи, а функція автокореляції на виході пов'язана з функцією автокореляції на вході як

Ry=hh~RX
щоб статистику виходу можна було також легко визначити, процес WSS в цілому і WSS Гаусса процес, зокрема, мають велике значення в інженерних програмах.


Чи не могли б ви прокоментувати "Білий шум" у цьому сенсі? За визначенням автокореляція при є дисперсією випадкових змінних. Чи означає це, що AWGN (Additive White Gaussian Noise) має нескінченну дисперсію? Я запитую це, тому що зазвичай люди пишуть , це неправильно? Чи слід писати ? Спасибі. τ=0n(t) N(0,N0/2)n(t) N(0,δ(0)N0/2)
Рой

1
@Drazick Будь ласка, задайте окреме запитання.
Діліп Сарват

1
Це фантастичний міні-курс у визначенні стаціонарних процесів. Я ніколи не бачив нічого подібного - викладено так методично і чітко. Вікі спільноти?
abalter

@Dilip Sarwate Вибачте за моє незнання. У прикладі. Чому E [X (t)] = 0 для всіх t? Ви припускали ергодичність? Як ви отримали функцію щільності ймовірності X (t) з функції щільності ймовірності тети для обчислення очікуваного значення? E [X (t) X (s)] = E [cos (t + theta) * cos (s + theta)] так? Які кроки ви зробили, щоб спростити цей вираз і дійти до того, що ви написали? Спасибі
VMMF

1
@VMMF Там немає NO ергодичності використовується. - дискретна випадкова величина, оскільки - дискретна випадкова величина, і вона з однаковою ймовірністю приймає значення і . Ерго, . приймає значення , , і з однаковою ймовірністю . Отже,X(t)=cos(t+Θ)Θ±cos(t)±sin(t)14E[X(t)]=0X(t)X(s)cos(t)cos(s)(cos(t))(cos(s))=cos(t)cos(s)sin(t)sin(s)(sin(t))(sin(s))=sin(t)sin(s)14E[X(t)(X(s)]=12(cos(t)cos(s)+sin(t)sin(s))=12cos(ts) . Отже,
Діліп Сарват
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.