مدل خودهمبسته

در آمار و پردازش سیگنال، مدل خود همبسته ، نوعی از فرایند تصادفی است که غالبا جهت مدلسازی و پیش‌بینی انواع مختلفی از پدیده‌های طبیعی و اجتماعی به کار می‌رود.
تعریف
عبارت (AR(p به مدل خودهمبسته مرتبه p اشاره دارد و به صورت زیر تعریف می‌شود:
که پارامترهای مدل، عدد ثابت و نوفه سفید می‌باشد. عددثابت در بسیاری از مواقع جهت سادگی حذف می‌شود.
یک مدل خودهمبسته می‌تواند به صورت خروجی یک فیلتر پاسخ ضربه با قطب‌های نامحدود که ورودی آن نوفه سفید است، در نظرگرفته شود.
اعمال برخی محدودیت‌ها بر پارامترهای این مدل، برای فرایند مانا بودن آن الزامی است. برای مثال فرایند (AR(۱ با  ۱ | ≥ ۱| مانا نخواهد بود. به صورت کلی‌تر، برای اینکه مدل (AR(p مانا باشد، ریشه‌های چندجمله‌ای باید درون دایره واحد قرار گیرند، یعنی برای هر ریشه باید داشته باشیم:
مثال، یک فرایند (AR(۱
یک فرایند (AR(۱ به صورت زیر داریم:
که نوفه سفید با میانگین صفر و واریانس می‌باشد. (توجه کنید که زیرنویس کنارگذاشته شده‌است) این فرایند ماناست در صورتی‌که ، زیرا در آنصورت مشابه خروجی یک فیلتر پایدار که وردی آن نوفه سفید است، می‌باشد. (اگر در آنصورت واریانس نامحدود خواهد داشت، و بنابراین مانا نخواهد بود) در نتیجه فرض ، میانگین ( برای تمام مقادیر t یکسان خواهد بود، با قراردادن میانگین برابر خواهیم داشت:
و در نتیجه:
به صورت خاص اگر ، در آنصورت میانگین برابر صفر خواهد بود.
نشان داده می‌شود که واریانس برابر خواهد بود با:
که واریانس می‌باشد.
اتوکواریانس برابر خواهد بود با:
تابع اتوکواریانس با ثابت زمانی تنزیل می‌یابد. (برای اثبات کافی است به صورت نوشته شود. توجه کنید که در این صورت خواهیم داشت: )
تابع چگالی طیفی، تبدیل فوریه تابع اتوکواریانس خواهد بود. در شرایط گسسته، تبدیل فوریه گسسته زمان تابع اتوکواریانس برابر خواهد بود با:
\Phi(\omega)=
\frac{1{\sqrt{2\pi\,\sum_{n=-\infty^\infty B_n e^{-i\omega n
\frac{1{\sqrt{2\pi\,\left(\frac{\sigma_\varepsilon^2{1+\varphi^2-2\varphi\cos(\omega)\right)
این عبارت به علت ساختار گسسته ، متناوب است، که در عبارت کسینوسی مخرج خود را نشان می‌دهد.
اگر فرض کنیم که زمان نمونه‌برداری ( ) نسبت به ثابت زمانی( ) خیلی کوچکتر باشد، می‌توانیم از تقریب پیوسته استفاده کنیم:
که فرم لورنتزین تابع چگالی طیفی را نتیجه می‌دهد:
\Phi(\omega)=
\frac{1{\sqrt{2\pi\,\frac{\sigma_\varepsilon^2{1-\varphi^2\,\frac{\gamma{\pi(\gamma^2+\omega^2)
که فرکانس زاویه‌ای متناسب با ثابت زمانی می‌باشد.
نمایش دیگر برای از طریق جایگزین کردن با در فرمول اصلی بدست می‌آید. با N بار تکرار این جایگزینی خواهیم داشت:
با میل کردن N به سمت بینهایت، به صفر میل کرده و خواهیم داشت:
درآنصورت به صورت مجموعه‌ای از چملات اخلال به علاوه میانگین ثابت درمی‌آید. در واقع این نمایش معادل یک میانگین متحرک از مرتبه بینهایت، (∞)MA خواهد بود. اگر یک فرایند گوسی باشد، نیز فرایند گوسی خواهد بود. در حالات دیگر، قضیه حد مرکزی بیان می‌کند که ، زمانی‌که مقداری نزدیک به یک داشته باشد، به صورت تقریبی توزیع نرمال خواهد داشت.
محاسبه پارامترهای AR
مدل (AR(p با معادله زیر داریم:
یک تناظر مستقیم میان پارامترهای و تابع کواریانس فرایند وجود دارد، این تناظر می‌تواند به گونه‌ای معکوس شود که پارامترها از روی تابع خودهمبستگی (که از کواریانس‌ها بدست می‌آید) تعیین و محاسبه شوند. این محاسبه از طریق حل معادلات یول-واکر انجام می‌شود:
\gamma_m = \sum_{k=1^p \varphi_k \gamma_{m-k + \sigma_\varepsilon^2\delta_{m,0
که با p + 1، m = 0, … , p معادله خواهیم داشت. تابع خودهمبستگی Xو انحراف معیار فرایند جمله اخلال ورودی و تابع دلتای کرونکر، ضربه کرونکر، است.
از آنجا که قسمت آخر معادله، ، تنها در زمانی‌که m = ۰ باشد، غیر صفر است، معادله غالبا از طریق نمایش ماتریس برای m> 0 حل می‌شود، در این صورت خواهیم داشت:
\begin{bmatrix
\gamma_1 \\
\gamma_2 \\
\gamma_3 \\
\vdots \\
\end{bmatrix
\begin{bmatrix
\gamma_0 - \gamma_{-1 - \gamma_{-2 - \dots \\
\gamma_1 - \gamma_0 - \gamma_{-1 - \dots \\
\gamma_2 - \gamma_{1 - \gamma_{0 - \dots \\
\vdots - \vdots - \vdots - \ddots \\
\end{bmatrix
\begin{bmatrix
\varphi_{1 \\
\varphi_{2 \\
\varphi_{3 \\
\vdots \\
\end{bmatrix
برای m = ۰ داریم:
\gamma_0 = \sum_{k=1^p \varphi_k \gamma_{-k + \sigma_\varepsilon^2
که به اما اجازه می‌دهد را حل کنیم.
معادلات بالا، معادلات یول-واکر، از طریق جایگزین کردن کواریانس‌های نظری با مقادیر تخمین زده شده، یک مسیر برای تخمین پارامترهای مدل (AR(p فراهم می‌کنند. یک راه برای محاسبه کواریانس‌های تخمین زده شده استفاده از برازش حداقل مربعات، رگرسیون خطی، مقادیرX t بر p مقدار قبلی خود می‌باشد.
استخراج
معادله معرف فرایند (AR(p:
با ضرب هر دو طرف معادله در X t − m و گرفتن امید انتظاری خواهیم داشت:
بنابر تعریف تابع خودهمبستگی داریم: . مقادیر جمله اخلال از یکدیگر مستقل بوده و X t − m به ازای m> 0، مستقل از  t است، یعنی برای m> 0 داریم:، E = ۰. برای m = ۰:
E\varepsilon_t X_{t
E\left\varepsilon_t \left(\sum_{i=1^p \varphi_i\,X_{t-i+ \varepsilon_t\right)\right
\sum_{i=1^p \varphi_i\, E\varepsilon_t\,X_{t-i + E\varepsilon_t^2
0 + \sigma_\varepsilon^2,
برای m ≥ ۰ داریم:
بعلاوه داریم:
E\left\sum_{i=1^p \varphi_i\,X_{t-i X_{t-m\right
\sum_{i=1^p \varphi_i\,EX_{t X_{t-m+i
\sum_{i=1^p \varphi_i\,\gamma_{m-i,
از جایگذاری معادله بالا در معادله ، معادلات یول-واکر تولید می‌شود، برای m ≥ ۰:
برای m <0 داریم