2025-11-11T17:13:08.951657

Interpolation of functionals of stochastic sequences with stationary increments from observations with noise

Luz, Moklyachuk
The problem of optimal estimation of linear functional ${{A}_{N}}ξ=\sum\limits_{k=0}^{N}{a(k)ξ(k)}\,$ depending on the unknown values of a stochastic sequence $ξ(m)$ with stationary $n$-th increments from observations of the sequence $ξ(k)$ at points $k=-1,-2,\ldots $ and of the sequence $ξ(k)+η(k)$ at points of time $k=N+1,N+2,\ldots $ is considered. Formulas for calculating the mean square error and the spectral characteristic of the optimal linear estimate of the functional are proposed under condition of spectral certainty, where spectral densities of the sequences $ξ(m)$ and $η(m)$ are exactly known. Minimax (robust) method of estimation is applied in the case where the spectral densities are not known exactly while some sets of admissible spectral densities are given. Formulas that determine the least favorable spectral densities and the minimax spectral characteristics are proposed for some specific sets of admissible densities.
academic

استيفاء الدوال الخطية للتسلسلات العشوائية ذات الزيادات المستقرة من الملاحظات المشوبة بالضوضاء

المعلومات الأساسية

  • معرّف الورقة: 2510.22764
  • العنوان: Interpolation of functionals of stochastic sequences with stationary increments from observations with noise
  • المؤلفون: Maksym Luz, Mykhailo Moklyachuk (جامعة تاراس شيفتشينكو الوطنية في كييف)
  • التصنيف: math.ST stat.TH
  • المجلة المنشورة: Applied Statistics. Actuarial and Financial Mathematics. No. 2, 131–148, 2012
  • رابط الورقة: https://arxiv.org/abs/2510.22764

الملخص

تتناول هذه الورقة مشكلة الحصول على التقديرات المثلى للدالة الخطية ANξ=∑k=0Na(k)ξ(k)A_N\xi = \sum_{k=0}^{N} a(k)\xi(k)، وهذه الدالة تعتمد على القيم المجهولة للتسلسل العشوائي ξ(m)\xi(m) الذي يمتلك زيادات مستقرة من الرتبة n. تتضمن بيانات الملاحظة التسلسل ξ(k)\xi(k) عند النقاط k=−1,−2,…k = -1, -2, \ldots والتسلسل المشوب بالضوضاء ξ(k)+η(k)\xi(k) + \eta(k) عند النقاط k=N+1,N+2,…k = N+1, N+2, \ldots. تحت الشروط المحددة طيفياً (حيث تكون كثافات الطيف للتسلسلات ξ(m)\xi(m) و η(m)\eta(m) معروفة بدقة)، يتم اقتراح صيغ لحساب متوسط الخطأ التربيعي للتقديرات الخطية المثلى والخصائص الطيفية. عندما لا تكون كثافات الطيف معروفة بدقة بل يتم إعطاء مجموعة من كثافات الطيف المسموحة، يتم تطبيق طريقة التقدير الأصغر-الأعظم (القوي). بالنسبة لمجموعات كثافة معينة مسموحة، يتم اقتراح صيغ لتحديد كثافات الطيف الأكثر عدم ملاءمة والخصائص الطيفية الأصغر-الأعظم.

خلفية البحث والدافع

  1. المشكلة البحثية: تتمثل المشكلة الأساسية التي تعالجها هذه الورقة في الحصول على تقديرات استيفاء خطية مثلى للدوال الخطية للتسلسلات العشوائية ذات الزيادات المستقرة تحت ظروف الملاحظات المشوبة بالضوضاء.
  2. أهمية المشكلة:
    • التسلسلات العشوائية ذات الزيادات المستقرة تمثل تعميماً مهماً للتسلسلات العشوائية المستقرة بشكل عام، وتشمل تسلسلات ARIMA والسلاسل الزمنية الدورية (الموسمية)
    • في التطبيقات العملية، غالباً ما تحتوي بيانات الملاحظة على ضوضاء، مما يتطلب إجراء تقديرات تحت ظروف المعلومات غير الكاملة
    • تتمتع هذه المشكلة بقيمة تطبيقية مهمة في مجالات التمويل ومعالجة الإشارات والتنبؤ الإحصائي
  3. قيود الطرق الموجودة:
    • تنطبق الطريقة الكلاسيكية لـ Kolmogorov-Wiener بشكل أساسي على الحالات التي تكون فيها كثافة الطيف معروفة
    • عندما تكون كثافة الطيف غير معروفة أو غير مؤكدة، تفتقر الطرق التقليدية إلى القوة الكافية
    • يوجد نقص في الإطار النظري المنهجي للملاحظات المشوبة بالضوضاء
  4. دافع البحث: إنشاء نظرية تقدير قوية تحت ظروف عدم التأكد من كثافة الطيف، لتوفير طرق تقدير أكثر موثوقية للتطبيقات العملية.

المساهمات الأساسية

  1. الإطار النظري: إنشاء إطار نظري شامل للاستيفاء الخطي الأمثل للتسلسلات العشوائية ذات الزيادات المستقرة من الرتبة n تحت الملاحظات المشوبة بالضوضاء
  2. صيغ التقدير: تحت الشروط المحددة طيفياً، توفير صيغ صريحة لحساب الخصائص الطيفية للتقديرات الخطية المثلى ومتوسط الخطأ التربيعي
  3. طريقة الأصغر-الأعظم: توسيع طريقة التقدير الأصغر-الأعظم (القوي) إلى مشاكل الاستيفاء مع الملاحظات المشوبة بالضوضاء
  4. كثافة الطيف الأكثر عدم ملاءمة: بالنسبة لفئات كثافة الطيف المسموحة المحددة، توفير طرق لتحديد كثافات الطيف الأكثر عدم ملاءمة والخصائص الطيفية الأصغر-الأعظم
  5. علاقات التصفية: الكشف عن الروابط الداخلية بين مشاكل الاستيفاء ومشاكل التصفية

شرح الطريقة

تعريف المهمة

المدخلات:

  • قيم التسلسل ξ(k)\xi(k) المرصودة عند النقاط k=−1,−2,…k = -1, -2, \ldots
  • قيم التسلسل المشوب بالضوضاء ξ(k)+η(k)\xi(k) + \eta(k) المرصودة عند النقاط k=N+1,N+2,…k = N+1, N+2, \ldots

المخرجات: التقدير الأمثل للدالة الخطية ANξ=∑k=0Na(k)ξ(k)A_N\xi = \sum_{k=0}^{N} a(k)\xi(k)

القيود:

  • ξ(m)\xi(m) يمتلك زيادات مستقرة من الرتبة n
  • η(k)\eta(k) غير مرتبط بـ ξ(k)\xi(k)
  • الفترة الزمنية {0,1,…,N}\{0, 1, \ldots, N\} هي فترة القيم المجهولة

معمارية النموذج

1. التمثيل الطيفي للزيادات المستقرة

التعريف: يتم تعريف الزيادة العشوائية من الرتبة n على النحو التالي: ξ(n)(m,μ)=(1−Bμ)nξ(m)=∑l=0n(−1)lCnlξ(m−lμ)\xi^{(n)}(m,\mu) = (1-B_\mu)^n\xi(m) = \sum_{l=0}^{n} (-1)^l C_n^l \xi(m-l\mu)

حيث BμB_\mu هو عامل الإزاحة بخطوة μ\mu.

التحليل الطيفي: تمتلك الزيادات المستقرة من الرتبة n التمثيل الطيفي: ξ(n)(m,μ)=∫−ππeimλ(1−e−iμλ)n1(iλ)ndZ(λ)\xi^{(n)}(m,\mu) = \int_{-\pi}^{\pi} e^{im\lambda} (1-e^{-i\mu\lambda})^n \frac{1}{(i\lambda)^n} dZ(\lambda)

حيث Z(λ)Z(\lambda) هو مقياس عشوائي متعامد.

2. تحليل الدالة

يتم تحليل الدالة الهدف على النحو التالي: ANξ=BNξ−VNξA_N\xi = B_N\xi - V_N\xi

حيث:

  • BNξ=∑k=0Nbμ(k)ξ(n)(k,μ)B_N\xi = \sum_{k=0}^{N} b_\mu(k)\xi^{(n)}(k,\mu) هي دالة الزيادة
  • VNξ=∑k=−1−μnvμ(k)ξ(k)V_N\xi = \sum_{k=-1}^{-\mu n} v_\mu(k)\xi(k) هي دالة الملاحظة المعروفة

3. بناء التقدير الأمثل

باستخدام طريقة الإسقاط المتعامد في فضاء هيلبرت، يتم إسقاط BNξB_N\xi على فضاء الملاحظة: B^Nξ=∫−ππhμ(1)(λ)dZξμ(n)(λ)+∫−ππhμ(2)(λ)dZξμ(n)+ημ(n)(λ)\hat{B}_N\xi = \int_{-\pi}^{\pi} h_\mu^{(1)}(\lambda) dZ_{\xi_\mu^{(n)}}(\lambda) + \int_{-\pi}^{\pi} h_\mu^{(2)}(\lambda) dZ_{\xi_\mu^{(n)} + \eta_\mu^{(n)}}(\lambda)

نقاط الابتكار التقني

  1. تحديد الخصائص الطيفية: من خلال حل نظام معادلات خطية لتحديد الخصائص الطيفية المثلى:
    • تتكون مصفوفات (Fμ)−1(F_\mu)^{-1} و (Gμ)−1(G_\mu)^{-1} من معاملات فورييه لكثافة الطيف
    • تحقق المعاملات شروط التعامد لضمان أمثلية التقدير
  2. طريقة المعادلات المصفوفية: تحويل المشكلة الطيفية المستمرة إلى معادلات مصفوفية منفصلة: Gμccμ=GμeeμG_\mu^c c_\mu = G_\mu^e e_\mubμ,N+Gμecμ=Gμceμ+Fμceμb_{\mu,N} + G_\mu^e c_\mu = G_\mu^c e_\mu + F_\mu^c e_\mu
  3. إطار الأصغر-الأعظم: عند عدم التأكد من كثافة الطيف، البحث عن نقطة سرج (h0,f0,g0)(h^0, f^0, g^0) التي تحقق: Δ(h;f0,g0)≥Δ(h0;f0,g0)≥Δ(h0;f,g)\Delta(h; f^0, g^0) \geq \Delta(h^0; f^0, g^0) \geq \Delta(h^0; f, g)

إعداد التجارب

إطار التحقق النظري

هذه الورقة تركز بشكل أساسي على البحث النظري، وتتحقق من صحة الطريقة من خلال الاستنتاج الرياضي، بما في ذلك:

  1. التحقق من الشروط: التحقق من أن كثافة الطيف تحقق شروط التكامل (15) و (16)
  2. قابلية عكس المصفوفات: ضمان قابلية عكس المصفوفات ذات الصلة
  3. تحليل التقارب: تحليل تقارب سلاسل فورييه

تحليل الحالات الخاصة

  1. حالة الضوضاء البيضاء: عندما تكون η(n)(k,μ)\eta^{(n)}(k,\mu) ضوضاء بيضاء، المصفوفة Gμe=σ2IG_\mu^e = \sigma^2 I
  2. حالة الدعم المحدود: صيغ مبسطة عندما تكون سلسلة المعاملات ذات دعم محدود
  3. فئات كثافة الطيف المحددة: تحليل خصائص الفئات D0,μ−D_{0,\mu}^- و DM,μ−D_{M,\mu}^-

نتائج التجارب

النتائج النظرية الرئيسية

النظرية 3.1 (وجود التقدير الأمثل)

تحت الشروط (15) و (16)، يوجد التقدير الخطي الأمثل B^Nξ\hat{B}_N\xi، وتُعطى خصائصه الطيفية بالصيغ (17) و (18)، ويُحسب متوسط الخطأ التربيعي بالصيغة (19).

النظرية 3.2 (تقدير الدالة)

التقدير الأمثل للدالة الهدف ANξA_N\xi هو: A^Nξ=−∑k=−1−μnvμ(k)ξ(k)+∫−ππhμ(a,1)(λ)dZξμ(n)(λ)+∫−ππhμ(a,2)(λ)dZξμ(n)+ημ(n)(λ)\hat{A}_N\xi = -\sum_{k=-1}^{-\mu n} v_\mu(k)\xi(k) + \int_{-\pi}^{\pi} h_\mu^{(a,1)}(\lambda) dZ_{\xi_\mu^{(n)}}(\lambda) + \int_{-\pi}^{\pi} h_\mu^{(a,2)}(\lambda) dZ_{\xi_\mu^{(n)} + \eta_\mu^{(n)}}(\lambda)

النظرية 6.1 (كثافة الطيف الأكثر عدم ملاءمة)

في الفئة D0,μ−×D0,μ−D_{0,\mu}^- \times D_{0,\mu}^-، تأخذ كثافة الطيف الأكثر عدم ملاءمة الشكل: f0(λ)=λ2n∣1−eiλμ∣2n∑k=−∞∞fμ0(∣k∣)eiλkf^0(\lambda) = \frac{\lambda^{2n}}{|1-e^{i\lambda\mu}|^{2n} \sum_{k=-\infty}^{\infty} f_\mu^0(|k|)e^{i\lambda k}}

نتائج الحالات الخاصة

  1. تحت ظروف الضوضاء البيضاء: عندما تكون الضوضاء بيضاء وسلسلة المعاملات موجبة، تتبسط الكثافة الأكثر عدم ملاءمة إلى: f0(λ)=∣(iλ)n(1−e−iλμ)n∑k=0Nγμ(k)e−iλk∣2f^0(\lambda) = \left|\frac{(i\lambda)^n}{(1-e^{-i\lambda\mu})^n} \sum_{k=0}^{N} \gamma_\mu(k)e^{-i\lambda k}\right|^2
  2. فئات القيود المحدودة: في الفئة DM,μ−D_{M,\mu}^-، يتم الحصول على كثافة الطيف الأكثر عدم ملاءمة من خلال حل نظام معادلات خطية محدود الأبعاد.

الأعمال ذات الصلة

التطور التاريخي

  1. النظرية الكلاسيكية: الأعمال الكلاسيكية لـ Kolmogorov (1992) و Wiener (1966) في نظرية التنبؤ
  2. نظرية الزيادات المستقرة: الأعمال الرائدة لـ Yaglom (1955) و Pinsker (1955)
  3. طريقة الأصغر-الأعظم: العمل الأول لـ Grenander (1957) في طريقة الاستقراء الأصغر-الأعظم

التطور الحديث

  1. طرق التحسين المحدب: الطرق القوية لـ Franke (1985) و Kassam & Poor (1985)
  2. العمليات المتجهة: التعميم المتجه لـ Rozanov (1967) و Masyutka & Moklyachuk (2012)
  3. التسلسلات المرتبطة دورياً: دراسة التسلسلات المرتبطة دورياً لـ Dubovets'ka & Moklyachuk (2013)

الابتكار في هذه الورقة

مقارنة بالأعمال الموجودة، يتمثل الابتكار الرئيسي في هذه الورقة في:

  • معالجة منهجية لمشاكل الاستيفاء مع الملاحظات المشوبة بالضوضاء
  • إنشاء الروابط بين مشاكل الاستيفاء والتصفية
  • توفير حلول صريحة لفئات متعددة من كثافات الطيف

الخلاصة والمناقشة

الاستنتاجات الرئيسية

  1. إطار نظري شامل: إنشاء نظرية استيفاء شاملة للتسلسلات العشوائية ذات الزيادات المستقرة تحت الملاحظات المشوبة بالضوضاء
  2. صيغ صريحة: توفير صيغ حسابية صريحة للخصائص الطيفية للتقديرات المثلى ومتوسط الخطأ التربيعي
  3. القوة: توفر طريقة الأصغر-الأعظم استراتيجية تقدير قوية عند عدم التأكد من كثافة الطيف
  4. الفائدة العملية: يمكن تطبيق النتائج النظرية مباشرة على نماذج ARIMA والسلاسل الزمنية الموسمية والمشاكل العملية الأخرى

القيود

  1. التعقيد الحسابي: ينطوي على عكس مصفوفات لا نهائية الأبعاد، مما يتطلب تقريب بالقطع في الحسابات العملية
  2. قيود الشروط: يتطلب أن تحقق كثافة الطيف شروط تكامل محددة
  3. قيود الخطية: يقتصر على التقديرات الخطية فقط، وقد تحقق التقديرات غير الخطية أداءً أفضل
  4. افتراضات الضوضاء: يفترض أن الضوضاء غير مرتبطة بالإشارة، بينما قد توجد ارتباطات في الواقع

الاتجاهات المستقبلية

  1. الخوارزميات الحسابية: تطوير خوارزميات حسابية فعالة
  2. التعميم غير الخطي: توسيع النطاق إلى تقدير الدوال غير الخطية
  3. الضوضاء المرتبطة: النظر في حالات الضوضاء المرتبطة بالإشارة
  4. التعميم متعدد الأبعاد: توسيع النطاق إلى حالات الحقول العشوائية متعددة الأبعاد

التقييم المتعمق

المزايا

  1. الصرامة النظرية: الاستنتاج الرياضي دقيق وإثبات النظريات شامل
  2. ابتكار الطريقة: تطبيق مبتكر لطريقة الأصغر-الأعظم على مشاكل الاستيفاء مع الضوضاء
  3. اكتمال النتائج: توفير نظام نظري شامل من الوجود إلى البناء
  4. القيمة العملية: توفر النتائج النظرية إرشادات لها قيمة عملية

أوجه القصور

  1. نقص التجارب العددية: عدم توفير أمثلة عددية للتحقق من النتائج النظرية
  2. تحليل التعقيد الحسابي: تحليل غير كافٍ لتعقيد الحسابات العملية
  3. التقارب: تحليل محدود لتقارب التقريب بالقطع
  4. الدراسات المقارنة: نقص المقارنة مع الطرق الأخرى

التأثير

  1. المساهمة الأكاديمية: توفير أدوات نظرية جديدة لنظرية تقدير العمليات العشوائية
  2. الإمكانات التطبيقية: آفاق تطبيقية في السلاسل الزمنية المالية ومعالجة الإشارات
  3. الأساس النظري: توفير أساس نظري للبحث ذي الصلة اللاحق
  4. الأهمية المنهجية: الأهمية المنهجية لتطبيق طريقة الأصغر-الأعظم

السيناريوهات المعمول بها

  1. النمذجة المالية: تقدير المعاملات والتنبؤ بنماذج ARIMA
  2. معالجة الإشارات: الاستيفاء والتصفية للإشارات المشوبة بالضوضاء
  3. التنبؤ الاقتصادي: تحليل السلاسل الزمنية الاقتصادية ذات الاتجاه والموسمية
  4. التطبيقات الهندسية: مشاكل تقدير الحالة في أنظمة التحكم

المراجع

تستشهد هذه الورقة بالأدبيات المهمة في هذا المجال، بما في ذلك:

  • Kolmogorov (1992): العمل الأساسي في نظرية العمليات العشوائية
  • Yaglom (1955, 1987): النظرية الكلاسيكية للعمليات ذات الزيادات المستقرة
  • Grenander (1957): العمل الرائد في طريقة الأصغر-الأعظم
  • سلسلة أعمال Moklyachuk: التطور المنهجي لنظرية التقدير القوي

التقييم الشامل: هذه ورقة ممتازة ذات طابع نظري قوي، وقد قدمت مساهمات مهمة في نظرية تقدير العمليات العشوائية. على الرغم من افتقارها للتحقق العددي، فإن الإطار النظري شامل والطريقة مبتكرة، مما يعطيها تأثيراً مهماً على تطور هذا المجال.