This paper focuses on the problem of the mean square optimal estimation of linear functionals which depend on the unknown values of a multidimensional stationary stochastic sequence.
Estimates are based on observations of the sequence with an additive stationary noise sequence.
The aim of the paper is to develop methods of finding the optimal estimates of the functionals in the case of missing observations.
The problem is investigated in the case of spectral certainty where the spectral densities of the sequences are exactly known.
Formulas for calculating the mean-square errors and the spectral characteristics of the optimal linear estimates of functionals are derived under the condition of spectral certainty.
The minimax (robust) method of estimation is applied in the case of spectral uncertainty, where spectral densities of the sequences are not known exactly while sets of admissible spectral densities are given. Formulas that determine the least favorable spectral densities and the minimax spectral characteristics of the optimal estimates of functionals are proposed for some special sets of admissible densities.
पेपर ID : 2511.07228शीर्षक : Extrapolation Problem for Multidimensional Stationary Sequences with Missing Observationsलेखक : Oleksandr Masyutka, Mikhail Moklyachuk, Maria Sideiसंस्थान : Taras Shevchenko National University of Kyivवर्गीकरण : math.ST (सांख्यिकी सिद्धांत), stat.THप्रकाशित पत्रिका : Statistics, Optimization and Information Computing, Vol. 7, मार्च 2019, pp 97-117पेपर लिंक : https://arxiv.org/abs/2511.07228 यह पेपर अनुपलब्ध अवलोकनों की स्थिति में बहुआयामी स्थिर यादृच्छिक अनुक्रमों की माध्य-वर्ग इष्टतम बहिर्वेशन समस्या का अध्ययन करता है। अनुमान योजक स्थिर शोर के साथ अनुक्रम अवलोकनों पर आधारित है। अध्ययन वर्णक्रमीय निश्चितता (spectral certainty) और वर्णक्रमीय अनिश्चितता (spectral uncertainty) दोनों स्थितियों में किया गया है। वर्णक्रमीय निश्चितता की स्थिति में, इष्टतम रैखिक अनुमान की माध्य-वर्ग त्रुटि और वर्णक्रमीय विशेषताओं की गणना के लिए सूत्र प्राप्त किए गए हैं; वर्णक्रमीय अनिश्चितता की स्थिति में, अत्यंत प्रतिकूल (minimax-robust) विधि लागू की गई है, और सबसे प्रतिकूल वर्णक्रमीय घनत्व और अत्यंत प्रतिकूल वर्णक्रमीय विशेषताओं को निर्धारित करने के लिए सूत्र प्रस्तावित किए गए हैं।
इस पेपर द्वारा हल की जाने वाली मूल समस्या है: अनुपलब्ध अवलोकनों की स्थिति में, बहुआयामी स्थिर यादृच्छिक अनुक्रमों के रैखिक कार्यात्मकों का इष्टतम अनुमान कैसे लगाया जाए? विशेष रूप से:
अवलोकन मॉडल : अवलोकन अनुक्रम ξ ( j ) + η ( j ) \xi(j) + \eta(j) ξ ( j ) + η ( j ) है, जहाँ ξ ( j ) \xi(j) ξ ( j ) संकेत अनुक्रम है और η ( j ) \eta(j) η ( j ) शोर अनुक्रम हैअनुपलब्ध पैटर्न : अवलोकन बिंदु j ∈ Z − ∖ S j \in \mathbb{Z}^- \setminus S j ∈ Z − ∖ S हैं, जहाँ S = ⋃ l = 1 s { − M l − N l , … , − M l } S = \bigcup_{l=1}^{s}\{-M_l-N_l, \ldots, -M_l\} S = ⋃ l = 1 s { − M l − N l , … , − M l } अनुपलब्ध अवलोकन खंडों को दर्शाता हैअनुमान लक्ष्य : रैखिक कार्यात्मक A ξ = ∑ j = 0 ∞ a ( j ) ⊤ ξ ( j ) A\xi = \sum_{j=0}^{\infty} a(j)^\top \xi(j) A ξ = ∑ j = 0 ∞ a ( j ) ⊤ ξ ( j ) सैद्धांतिक मूल्य : Kolmogorov-Wiener शास्त्रीय पूर्वानुमान सिद्धांत को अनुपलब्ध अवलोकन परिदृश्य तक विस्तारित करता हैव्यावहारिक महत्व : वास्तविक अनुप्रयोगों में, सेंसर विफलता, डेटा संचरण व्यवधान आदि अक्सर अवलोकन अनुपलब्धता का कारण बनते हैंमजबूतता आवश्यकता : व्यावहारिक रूप से वर्णक्रमीय घनत्व अक्सर अज्ञात या अनुमानित होता है, मजबूत अनुमान विधि की आवश्यकता होती हैपूर्ण अवलोकन धारणा : पारंपरिक विधियाँ (Wiener, Yaglom, Rozanov आदि) पूर्ण अवलोकन मानती हैंवर्णक्रमीय निश्चितता धारणा : अधिकांश विधियों को सटीक रूप से ज्ञात वर्णक्रमीय घनत्व की आवश्यकता होती है, जो व्यावहारिक रूप से पूरी नहीं होतीएकल-चर सीमा : बहुआयामी स्थिति में सिद्धांत और विधियाँ अपेक्षाकृत अविकसित हैंइस पेपर की नवीनता इसमें निहित है:
Hilbert अंतरिक्ष प्रक्षेपण विधि को अनुपलब्ध अवलोकन परिदृश्य तक विस्तारित करना वर्णक्रमीय अनिश्चितता के तहत अत्यंत प्रतिकूल मजबूत अनुमान सिद्धांत विकसित करना बहुआयामी स्थिति के लिए एक संपूर्ण सैद्धांतिक ढांचा प्रदान करना सैद्धांतिक ढांचा : अनुपलब्ध अवलोकनों के तहत बहुआयामी स्थिर अनुक्रमों की बहिर्वेशन समस्या का एक संपूर्ण सैद्धांतिक प्रणाली स्थापित करता हैवर्णक्रमीय निश्चितता स्थिति :
इष्टतम रैखिक अनुमान की स्पष्ट वर्णक्रमीय विशेषता सूत्र प्राप्त करता है (सूत्र 10) माध्य-वर्ग त्रुटि के लिए सटीक गणना सूत्र देता है (सूत्र 11) वर्णक्रमीय अनिश्चितता स्थिति :
अत्यंत प्रतिकूल मजबूत अनुमान विधि विकसित करता है सबसे प्रतिकूल वर्णक्रमीय घनत्व के लक्षण वर्णन समीकरण प्रस्तावित करता है कई विशेष स्वीकार्य वर्णक्रमीय घनत्व वर्गों के लिए ठोस समाधान देता है विशेष स्थितियाँ : शोर-मुक्त अवलोकन, असंबंधित शोर आदि विशेष स्थितियों के लिए निष्कर्ष प्रदान करता हैगणना विधि : ऑपरेटर समीकरणों और फूरियर गुणांकों के माध्यम से एक गणनीय ढांचा स्थापित करता हैइनपुट :
अवलोकन अनुक्रम: { ξ ( j ) + η ( j ) , j ∈ Z − ∖ S } \{\xi(j) + \eta(j), j \in \mathbb{Z}^- \setminus S\} { ξ ( j ) + η ( j ) , j ∈ Z − ∖ S } अनुपलब्ध समुच्चय: S = ⋃ l = 1 s { − M l − N l , … , − M l } S = \bigcup_{l=1}^{s}\{-M_l-N_l, \ldots, -M_l\} S = ⋃ l = 1 s { − M l − N l , … , − M l } कार्यात्मक गुणांक: { a ( j ) , j = 0 , 1 , … } \{a(j), j=0,1,\ldots\} { a ( j ) , j = 0 , 1 , … } संतुष्ट करते हैं ∑ j = 0 ∞ ∑ k = 1 T ∣ a k ( j ) ∣ < ∞ \sum_{j=0}^{\infty}\sum_{k=1}^{T}|a_k(j)| < \infty ∑ j = 0 ∞ ∑ k = 1 T ∣ a k ( j ) ∣ < ∞ आउटपुट :
इष्टतम अनुमान: A ^ ξ = ∫ − π π h ( e i λ ) ⊤ ( Z ξ ( d λ ) + Z η ( d λ ) ) \hat{A}\xi = \int_{-\pi}^{\pi} h(e^{i\lambda})^\top (Z_\xi(d\lambda) + Z_\eta(d\lambda)) A ^ ξ = ∫ − π π h ( e iλ ) ⊤ ( Z ξ ( d λ ) + Z η ( d λ )) माध्य-वर्ग त्रुटि: Δ ( h ; F , G ) = E ∣ A ξ − A ^ ξ ∣ 2 \Delta(h; F, G) = E|A\xi - \hat{A}\xi|^2 Δ ( h ; F , G ) = E ∣ A ξ − A ^ ξ ∣ 2 बाधा शर्तें :
न्यूनतमता शर्त: ∫ − π π Tr ( F ( λ ) + G ( λ ) ) − 1 d λ < ∞ \int_{-\pi}^{\pi} \text{Tr}(F(\lambda) + G(\lambda))^{-1}d\lambda < \infty ∫ − π π Tr ( F ( λ ) + G ( λ ) ) − 1 d λ < ∞ इस पेपर की मूल विधि Kolmogorov के Hilbert अंतरिक्ष प्रक्षेपण सिद्धांत पर आधारित है:
Hilbert अंतरिक्ष निर्माण :H = L 2 ( Ω , F , P ) H = L_2(\Omega, \mathcal{F}, P) H = L 2 ( Ω , F , P ) : शून्य माध्य, परिमित विचरण यादृच्छिक चर द्वारा उत्पन्नH s ( ξ + η ) H_s(\xi + \eta) H s ( ξ + η ) : अवलोकन मानों { ξ k ( j ) + η k ( j ) : j ∈ Z − ∖ S , k = 1 , … , T } \{\xi_k(j) + \eta_k(j): j \in \mathbb{Z}^- \setminus S, k=1,\ldots,T\} { ξ k ( j ) + η k ( j ) : j ∈ Z − ∖ S , k = 1 , … , T } द्वारा उत्पन्न बंद रैखिक उप-अंतरिक्षइष्टतम अनुमान लक्षण वर्णन : इष्टतम अनुमान A ^ ξ \hat{A}\xi A ^ ξ A ξ A\xi A ξ का H s ( ξ + η ) H_s(\xi+\eta) H s ( ξ + η ) पर लंबकोणीय प्रक्षेपण है, संतुष्ट करता है:A ^ ξ ∈ H s ( ξ + η ) \hat{A}\xi \in H_s(\xi + \eta) A ^ ξ ∈ H s ( ξ + η ) A ξ − A ^ ξ ⊥ H s ( ξ + η ) A\xi - \hat{A}\xi \perp H_s(\xi + \eta) A ξ − A ^ ξ ⊥ H s ( ξ + η ) वर्णक्रमीय अपघटन का उपयोग करते हुए:
ξ ( j ) = ∫ − π π e i j λ Z ξ ( d λ ) , A ξ = ∫ − π π A ( e i λ ) ⊤ Z ξ ( d λ ) \xi(j) = \int_{-\pi}^{\pi} e^{ij\lambda} Z_\xi(d\lambda), \quad A\xi = \int_{-\pi}^{\pi} A(e^{i\lambda})^\top Z_\xi(d\lambda) ξ ( j ) = ∫ − π π e ijλ Z ξ ( d λ ) , A ξ = ∫ − π π A ( e iλ ) ⊤ Z ξ ( d λ )
जहाँ A ( e i λ ) = ∑ j = 0 ∞ a ( j ) e i j λ A(e^{i\lambda}) = \sum_{j=0}^{\infty} a(j)e^{ij\lambda} A ( e iλ ) = ∑ j = 0 ∞ a ( j ) e ijλ
लंबकोणीयता शर्त के माध्यम से, वर्णक्रमीय विशेषता h ( e i λ ) h(e^{i\lambda}) h ( e iλ ) संतुष्ट करती है:
( A ( e i λ ) ) ⊤ ( F ( λ ) + F ξ η ( λ ) ) − ( h ( e i λ ) ) ⊤ F ζ ( λ ) = ( C ( e i λ ) ) ⊤ (A(e^{i\lambda}))^\top(F(\lambda) + F_{\xi\eta}(\lambda)) - (h(e^{i\lambda}))^\top F_\zeta(\lambda) = (C(e^{i\lambda}))^\top ( A ( e iλ ) ) ⊤ ( F ( λ ) + F ξ η ( λ )) − ( h ( e iλ ) ) ⊤ F ζ ( λ ) = ( C ( e iλ ) ) ⊤
जहाँ F ζ ( λ ) = F ( λ ) + F ξ η ( λ ) + F η ξ ( λ ) + G ( λ ) F_\zeta(\lambda) = F(\lambda) + F_{\xi\eta}(\lambda) + F_{\eta\xi}(\lambda) + G(\lambda) F ζ ( λ ) = F ( λ ) + F ξ η ( λ ) + F η ξ ( λ ) + G ( λ ) , C ( e i λ ) = ∑ j ∈ U c ( j ) e i j λ C(e^{i\lambda}) = \sum_{j \in U} c(j)e^{ij\lambda} C ( e iλ ) = ∑ j ∈ U c ( j ) e ijλ , U = S ∪ { 0 , 1 , … } U = S \cup \{0,1,\ldots\} U = S ∪ { 0 , 1 , … }
फूरियर गुणांक प्रस्तुत करते हुए:
B ( k − j ) = 1 2 π ∫ − π π ( F ζ ( λ ) ) − 1 e − i ( k − j ) λ d λ B(k-j) = \frac{1}{2\pi}\int_{-\pi}^{\pi} (F_\zeta(\lambda))^{-1}e^{-i(k-j)\lambda}d\lambda B ( k − j ) = 2 π 1 ∫ − π π ( F ζ ( λ ) ) − 1 e − i ( k − j ) λ d λ
R ( k − j ) = 1 2 π ∫ − π π ( F ( λ ) + F ξ η ( λ ) ) ( F ζ ( λ ) ) − 1 e − i ( k − j ) λ d λ R(k-j) = \frac{1}{2\pi}\int_{-\pi}^{\pi} (F(\lambda) + F_{\xi\eta}(\lambda))(F_\zeta(\lambda))^{-1}e^{-i(k-j)\lambda}d\lambda R ( k − j ) = 2 π 1 ∫ − π π ( F ( λ ) + F ξ η ( λ )) ( F ζ ( λ ) ) − 1 e − i ( k − j ) λ d λ
अज्ञात गुणांक c ( k ) , k ∈ U c(k), k \in U c ( k ) , k ∈ U ऑपरेटर समीकरण द्वारा निर्धारित होते हैं:
R a = B c Ra = Bc R a = B c
जहाँ ऑपरेटर B , R B, R B , R संबंधित ब्लॉक मैट्रिक्स द्वारा परिभाषित होते हैं, अनुपलब्ध अवलोकनों की संरचना पर विचार करते हुए।
( h ( e i λ ) ) ⊤ = ( A ( e i λ ) ) ⊤ ( F ( λ ) + F ξ η ( λ ) ) ( F ζ ( λ ) ) − 1 − ( ∑ k ∈ U ( B − 1 R a ) ( k ) e i k λ ) ⊤ ( F ζ ( λ ) ) − 1 (h(e^{i\lambda}))^\top = (A(e^{i\lambda}))^\top(F(\lambda) + F_{\xi\eta}(\lambda))(F_\zeta(\lambda))^{-1} - \left(\sum_{k \in U}(B^{-1}Ra)(k)e^{ik\lambda}\right)^\top(F_\zeta(\lambda))^{-1} ( h ( e iλ ) ) ⊤ = ( A ( e iλ ) ) ⊤ ( F ( λ ) + F ξ η ( λ )) ( F ζ ( λ ) ) − 1 − ( ∑ k ∈ U ( B − 1 R a ) ( k ) e ikλ ) ⊤ ( F ζ ( λ ) ) − 1
Δ ( h ; F , G ) = ⟨ R a , B − 1 R a ⟩ + ⟨ Q a , a ⟩ \Delta(h; F, G) = \langle Ra, B^{-1}Ra \rangle + \langle Qa, a \rangle Δ ( h ; F , G ) = ⟨ R a , B − 1 R a ⟩ + ⟨ Q a , a ⟩
जहाँ Q Q Q फूरियर गुणांकों द्वारा परिभाषित रैखिक ऑपरेटर है।
सबसे प्रतिकूल वर्णक्रमीय घनत्व (परिभाषा 3.1): ( F 0 , G 0 ) ∈ D (F^0, G^0) \in \mathcal{D} ( F 0 , G 0 ) ∈ D को सबसे प्रतिकूल कहा जाता है, यदि
Δ ( h ( F 0 , G 0 ) ; F 0 , G 0 ) = max ( F , G ) ∈ D Δ ( h ( F , G ) ; F , G ) \Delta(h(F^0, G^0); F^0, G^0) = \max_{(F,G) \in \mathcal{D}} \Delta(h(F,G); F, G) Δ ( h ( F 0 , G 0 ) ; F 0 , G 0 ) = max ( F , G ) ∈ D Δ ( h ( F , G ) ; F , G )
अत्यंत प्रतिकूल वर्णक्रमीय विशेषता (परिभाषा 3.2): h 0 ∈ H D h^0 \in H_{\mathcal{D}} h 0 ∈ H D को अत्यंत प्रतिकूल कहा जाता है, यदि
min h ∈ H D max ( F , G ) ∈ D Δ ( h ; F , G ) = max ( F , G ) ∈ D Δ ( h 0 ; F , G ) \min_{h \in H_{\mathcal{D}}} \max_{(F,G) \in \mathcal{D}} \Delta(h; F, G) = \max_{(F,G) \in \mathcal{D}} \Delta(h^0; F, G) min h ∈ H D max ( F , G ) ∈ D Δ ( h ; F , G ) = max ( F , G ) ∈ D Δ ( h 0 ; F , G )
अत्यंत प्रतिकूल समस्या बाधित अनुकूलन के समतुल्य है:
max ( F , G ) ∈ D ( ⟨ R a , B − 1 R a ⟩ + ⟨ Q a , a ⟩ ) \max_{(F,G) \in \mathcal{D}} (\langle Ra, B^{-1}Ra \rangle + \langle Qa, a \rangle) max ( F , G ) ∈ D (⟨ R a , B − 1 R a ⟩ + ⟨ Q a , a ⟩)
अबाधित अनुकूलन में परिवर्तित:
Δ D ( F , G ) = − Δ ( h ( F 0 , G 0 ) ; F , G ) + δ ( ( F , G ) ∣ D ) → inf \Delta_{\mathcal{D}}(F,G) = -\Delta(h(F^0, G^0); F, G) + \delta((F,G)|\mathcal{D}) \to \inf Δ D ( F , G ) = − Δ ( h ( F 0 , G 0 ) ; F , G ) + δ (( F , G ) ∣ D ) → inf
जहाँ δ \delta δ सूचक फलन है।
सबसे प्रतिकूल वर्णक्रमीय घनत्व उप-अवकल शर्त द्वारा निर्धारित होता है:
0 ∈ ∂ Δ D ( F 0 , G 0 ) 0 \in \partial \Delta_{\mathcal{D}}(F^0, G^0) 0 ∈ ∂ Δ D ( F 0 , G 0 )
Lagrange गुणक विधि और उप-अवकल के रूप का उपयोग करते हुए, विशिष्ट लक्षण वर्णन समीकरण प्राप्त किए जा सकते हैं।
पेपर कई विशेष वर्गों पर विचार करता है, उदाहरण के लिए:
D 0 1 = { F ( λ ) ∣ 1 2 π ∫ − π π Tr F ( λ ) d λ = p } \mathcal{D}^1_0 = \left\{F(\lambda) \left| \frac{1}{2\pi}\int_{-\pi}^{\pi} \text{Tr}F(\lambda)d\lambda = p\right.\right\} D 0 1 = { F ( λ ) 2 π 1 ∫ − π π Tr F ( λ ) d λ = p }
D 1 U V = { G ( λ ) ∣ Tr V ( λ ) ≤ Tr G ( λ ) ≤ Tr U ( λ ) , 1 2 π ∫ − π π Tr G ( λ ) d λ = q } \mathcal{D}^{UV}_1 = \left\{G(\lambda) \left| \text{Tr}V(\lambda) \leq \text{Tr}G(\lambda) \leq \text{Tr}U(\lambda), \frac{1}{2\pi}\int_{-\pi}^{\pi}\text{Tr}G(\lambda)d\lambda = q\right.\right\} D 1 U V = { G ( λ ) Tr V ( λ ) ≤ Tr G ( λ ) ≤ Tr U ( λ ) , 2 π 1 ∫ − π π Tr G ( λ ) d λ = q }
सबसे प्रतिकूल वर्णक्रमीय घनत्व समीकरण (प्रमेय 4.1):
( r G 0 ( λ ) ) ∗ ( r G 0 ( λ ) ) ⊤ = α 2 ( F 0 ( λ ) + G 0 ( λ ) ) 2 (r^0_G(\lambda))^*(r^0_G(\lambda))^\top = \alpha^2(F^0(\lambda) + G^0(\lambda))^2 ( r G 0 ( λ ) ) ∗ ( r G 0 ( λ ) ) ⊤ = α 2 ( F 0 ( λ ) + G 0 ( λ ) ) 2
( r F 0 ( λ ) ) ∗ ( r F 0 ( λ ) ) ⊤ = ( β 2 + γ 1 ( λ ) + γ 2 ( λ ) ) ( F 0 ( λ ) + G 0 ( λ ) ) 2 (r^0_F(\lambda))^*(r^0_F(\lambda))^\top = (\beta^2 + \gamma_1(\lambda) + \gamma_2(\lambda))(F^0(\lambda) + G^0(\lambda))^2 ( r F 0 ( λ ) ) ∗ ( r F 0 ( λ ) ) ⊤ = ( β 2 + γ 1 ( λ ) + γ 2 ( λ )) ( F 0 ( λ ) + G 0 ( λ ) ) 2
जहाँ α 2 , β 2 \alpha^2, \beta^2 α 2 , β 2 Lagrange गुणक हैं, γ 1 ( λ ) ≤ 0 \gamma_1(\lambda) \leq 0 γ 1 ( λ ) ≤ 0 (जब Tr G 0 ( λ ) > Tr V ( λ ) \text{Tr}G^0(\lambda) > \text{Tr}V(\lambda) Tr G 0 ( λ ) > Tr V ( λ ) हो तो 0), γ 2 ( λ ) ≥ 0 \gamma_2(\lambda) \geq 0 γ 2 ( λ ) ≥ 0 (जब Tr G 0 ( λ ) < Tr U ( λ ) \text{Tr}G^0(\lambda) < \text{Tr}U(\lambda) Tr G 0 ( λ ) < Tr U ( λ ) हो तो 0)।
पेपर निम्नलिखित पर भी विचार करता है:
D 0 2 × D 2 U V \mathcal{D}^2_0 \times \mathcal{D}^{UV}_2 D 0 2 × D 2 U V : विकर्ण तत्व बाधाएंD 0 3 × D 3 U V \mathcal{D}^3_0 \times \mathcal{D}^{UV}_3 D 0 3 × D 3 U V : भारित ट्रेस बाधाएंD 0 4 × D 4 U V \mathcal{D}^4_0 \times \mathcal{D}^{UV}_4 D 0 4 × D 4 U V : मैट्रिक्स असमानता बाधाएंD ϵ × D δ 1 \mathcal{D}_\epsilon \times \mathcal{D}^1_\delta D ϵ × D δ 1 : ϵ \epsilon ϵ -प्रदूषण और δ \delta δ -पड़ोस मॉडलप्रत्येक वर्ग के लिए संबंधित लक्षण वर्णन समीकरण दिए गए हैं।
पेपर एक ठोस द्विआयामी अनुक्रम बहिर्वेशन उदाहरण प्रदान करता है:
समस्या सेटअप :
कार्यात्मक: A 1 ξ = a ( 0 ) ⊤ ξ ( 0 ) + a ( 1 ) ⊤ ξ ( 1 ) A_1\xi = a(0)^\top\xi(0) + a(1)^\top\xi(1) A 1 ξ = a ( 0 ) ⊤ ξ ( 0 ) + a ( 1 ) ⊤ ξ ( 1 ) , जहाँ a ( 0 ) = a ( 1 ) = ( 1 , 1 ) ⊤ a(0) = a(1) = (1,1)^\top a ( 0 ) = a ( 1 ) = ( 1 , 1 ) ⊤ अनुक्रम: ξ 1 ( n ) = ξ ( n ) \xi_1(n) = \xi(n) ξ 1 ( n ) = ξ ( n ) , ξ 2 ( n ) = ξ ( n ) + η ( n ) \xi_2(n) = \xi(n) + \eta(n) ξ 2 ( n ) = ξ ( n ) + η ( n ) अनुपलब्ध समुच्चय: S = { − 3 , − 2 } S = \{-3, -2\} S = { − 3 , − 2 } वर्णक्रमीय घनत्व:
f ( λ ) = 1 ∣ 1 − b 1 e i λ ∣ 2 , g ( λ ) = 1 ∣ 1 − b 2 e i λ ∣ 2 f(\lambda) = \frac{1}{|1-b_1e^{i\lambda}|^2}, \quad g(\lambda) = \frac{1}{|1-b_2e^{i\lambda}|^2} f ( λ ) = ∣1 − b 1 e iλ ∣ 2 1 , g ( λ ) = ∣1 − b 2 e iλ ∣ 2 1 वर्णक्रमीय घनत्व मैट्रिक्स:
F ( λ ) = ( f ( λ ) f ( λ ) f ( λ ) f ( λ ) + g ( λ ) ) F(\lambda) = \begin{pmatrix} f(\lambda) & f(\lambda) \\ f(\lambda) & f(\lambda) + g(\lambda) \end{pmatrix} F ( λ ) = ( f ( λ ) f ( λ ) f ( λ ) f ( λ ) + g ( λ ) ) व्युत्क्रम वर्णक्रमीय घनत्व मैट्रिक्स :
( F ( λ ) ) − 1 = ( 1 f ( λ ) + 1 g ( λ ) − 1 g ( λ ) − 1 g ( λ ) 1 g ( λ ) ) = B ( − 1 ) e − i λ + B ( 0 ) + B ( 1 ) e i λ (F(\lambda))^{-1} = \begin{pmatrix} \frac{1}{f(\lambda)} + \frac{1}{g(\lambda)} & -\frac{1}{g(\lambda)} \\ -\frac{1}{g(\lambda)} & \frac{1}{g(\lambda)} \end{pmatrix} = B(-1)e^{-i\lambda} + B(0) + B(1)e^{i\lambda} ( F ( λ ) ) − 1 = ( f ( λ ) 1 + g ( λ ) 1 − g ( λ ) 1 − g ( λ ) 1 g ( λ ) 1 ) = B ( − 1 ) e − iλ + B ( 0 ) + B ( 1 ) e iλ फूरियर गुणांक :
B ( 0 ) = ( 2 + b 1 2 + b 2 2 − 1 − b 2 2 − 1 − b 2 2 1 + b 2 2 ) , B ( 1 ) = B ( − 1 ) = ( − b 1 − b 2 b 2 b 2 − b 2 ) B(0) = \begin{pmatrix} 2+b_1^2+b_2^2 & -1-b_2^2 \\ -1-b_2^2 & 1+b_2^2 \end{pmatrix}, \quad B(1) = B(-1) = \begin{pmatrix} -b_1-b_2 & b_2 \\ b_2 & -b_2 \end{pmatrix} B ( 0 ) = ( 2 + b 1 2 + b 2 2 − 1 − b 2 2 − 1 − b 2 2 1 + b 2 2 ) , B ( 1 ) = B ( − 1 ) = ( − b 1 − b 2 b 2 b 2 − b 2 ) ऑपरेटर मैट्रिक्स : ब्लॉक मैट्रिक्स B B B का निर्माण, अनुपलब्ध स्थान { − 3 , − 2 } \{-3, -2\} { − 3 , − 2 } और भविष्य की स्थिति { 0 , 1 , 2 , … } \{0, 1, 2, \ldots\} { 0 , 1 , 2 , … } पर विचार करते हुएवर्णक्रमीय अपघटन : गुणनखंडन का उपयोग करते हुए
( F ( λ ) ) − 1 = ( ∑ j = 0 ∞ ψ ( j ) e − i j λ ) ⋅ ( ∑ j = 0 ∞ ψ ( j ) e − i j λ ) ∗ (F(\lambda))^{-1} = \left(\sum_{j=0}^{\infty}\psi(j)e^{-ij\lambda}\right) \cdot \left(\sum_{j=0}^{\infty}\psi(j)e^{-ij\lambda}\right)^* ( F ( λ ) ) − 1 = ( ∑ j = 0 ∞ ψ ( j ) e − ijλ ) ⋅ ( ∑ j = 0 ∞ ψ ( j ) e − ijλ ) ∗ जहाँ ψ ( 0 ) = ( 1 1 0 − 1 ) \psi(0) = \begin{pmatrix} 1 & 1 \\ 0 & -1 \end{pmatrix} ψ ( 0 ) = ( 1 0 1 − 1 ) , ψ ( 1 ) = ( − b 1 − b 2 0 b 2 ) \psi(1) = \begin{pmatrix} -b_1 & -b_2 \\ 0 & b_2 \end{pmatrix} ψ ( 1 ) = ( − b 1 0 − b 2 b 2 ) व्युत्क्रम ऑपरेटर गणना : B 11 − 1 ( i , j ) = ( Θ ∗ Θ ) ( i , j ) = ∑ l = 0 min ( i , j ) ( θ ( i − l ) ) ∗ θ ( j − l ) B^{-1}_{11}(i,j) = (\Theta^*\Theta)(i,j) = \sum_{l=0}^{\min(i,j)}(\theta(i-l))^*\theta(j-l) B 11 − 1 ( i , j ) = ( Θ ∗ Θ ) ( i , j ) = ∑ l = 0 m i n ( i , j ) ( θ ( i − l ) ) ∗ θ ( j − l ) वर्णक्रमीय विशेषता :
( h 1 ( e i λ ) ) ⊤ = − ( b 2 + b 2 2 − 2 ( b 1 + b 1 2 ) , − b 2 − b 2 2 ) e − i λ (h_1(e^{i\lambda}))^\top = -(b_2 + b_2^2 - 2(b_1 + b_1^2), -b_2 - b_2^2)e^{-i\lambda} ( h 1 ( e iλ ) ) ⊤ = − ( b 2 + b 2 2 − 2 ( b 1 + b 1 2 ) , − b 2 − b 2 2 ) e − iλ
माध्य-वर्ग त्रुटि :
Δ ( h 1 ; F ) = 10 + 8 b 1 + 4 b 1 2 + 2 b 2 + b 2 2 \Delta(h_1; F) = 10 + 8b_1 + 4b_1^2 + 2b_2 + b_2^2 Δ ( h 1 ; F ) = 10 + 8 b 1 + 4 b 1 2 + 2 b 2 + b 2 2
यह उदाहरण प्रदर्शित करता है:
अनुपलब्ध अवलोकनों की ब्लॉक संरचना को कैसे संभालें वर्णक्रमीय अपघटन का उपयोग करके गणना को कैसे सरल बनाएं इष्टतम वर्णक्रमीय विशेषता का स्पष्ट रूप पेपर उदाहरण 2.1 के माध्यम से सैद्धांतिक ढांचे की व्यवहार्यता को सत्यापित करता है:
वर्णक्रमीय विशेषता की सरलता : इष्टतम वर्णक्रमीय विशेषता में परिमित समर्थन होता है (केवल λ − 1 \lambda^{-1} λ − 1 पद में गैर-शून्य), अनुपलब्ध अवलोकनों के प्रभाव को स्थानीय रूप से प्रतिबिंबित करता हैत्रुटि की गणनीयता : माध्य-वर्ग त्रुटि अभिव्यक्ति पैरामीटर b 1 , b 2 b_1, b_2 b 1 , b 2 का एक सरल बहुपद है, विश्लेषण और अनुकूलन के लिए सुविधाजनक हैपैरामीटर प्रभाव :b 1 , b 2 b_1, b_2 b 1 , b 2 जितना बड़ा, त्रुटि उतनी बड़ी (संकेत और शोर की स्वसंबंध बढ़ता है)त्रुटि b 1 b_1 b 1 के प्रति संवेदनशीलता अधिक है (संकेत स्वसंबंध का प्रभाव अधिक महत्वपूर्ण है) मौजूदा विधियों की तुलना में:
पूर्णता : समस्या मॉडलिंग से लेकर ठोस गणना तक एक संपूर्ण ढांचा प्रदान करता हैसामान्यता : बहुआयामी अनुक्रमों, किसी भी अनुपलब्ध पैटर्न पर लागू होता हैमजबूतता : अत्यंत प्रतिकूल विधि वर्णक्रमीय अनिश्चितता को संभालती हैगणनीयता : ऑपरेटर समीकरणों और फूरियर गुणांकों के माध्यम से कार्यान्वयनपेपर कई प्रमेयों द्वारा गारंटियाँ प्रदान करता है:
प्रमेय 2.1 : वर्णक्रमीय निश्चितता स्थिति में इष्टतम समाधान का अस्तित्व और अद्वितीयताप्रमेय 4.1, 5.1 : विभिन्न स्वीकार्य वर्गों के तहत सबसे प्रतिकूल वर्णक्रमीय घनत्व का लक्षण वर्णननिष्कर्ष 2.1-2.4, 4.1-4.2, 5.1-5.2 : विशेष स्थितियों के सरलीकृत परिणामKolmogorov (1941) : पहली बार स्थिर अनुक्रम पूर्वानुमान समस्या के वर्णक्रमीय विधि का प्रस्तावWiener (1949) : निरंतर समय फिल्टरिंग सिद्धांत विकसित कियाYaglom (1955, 1987) : स्थिर प्रक्रिया के संबंधित सिद्धांत का व्यवस्थित अध्ययनRozanov (1967) : बहुआयामी स्थिर प्रक्रिया सिद्धांतHannan (1970) : बहु-चर समय श्रृंखला विश्लेषणBondon (2002, 2005) : अधूरे अतीत की पूर्वानुमान समस्याCheng & Pourahmadi (1996, 1998) : L p ( w ) L^p(w) L p ( w ) अंतरिक्ष में चरम समस्याएं और प्रक्षेपKasahara, Pourahmadi & Inoue (2009) : अनुपलब्ध मान पूर्वानुमान की द्वैत विधिPelagatti (2015) : अप्रेक्षणीय घटकों का समय श्रृंखला मॉडलिंगGrenander (1957) : पहली बार स्थिर प्रक्रिया बहिर्वेशन की अत्यंत प्रतिकूल विधि का प्रस्तावKassam & Poor (1985) : संकेत प्रसंस्करण में मजबूत तकनीकों का सर्वेक्षणFranke (1984, 1985) : समय श्रृंखला की मजबूत पूर्वानुमान और प्रक्षेपFranke & Poor (1984) : अत्यंत प्रतिकूल मजबूत फिल्टरिंगVastola & Poor (1983) : वर्णक्रमीय अनिश्चितता का Wiener फिल्टरिंग पर प्रभाव विश्लेषणMoklyachuk (2008, 2015) : स्थिर अनुक्रम कार्यात्मकों की मजबूत अनुमानMoklyachuk & Masyutka (2008-2012) : बहुआयामी स्थिर प्रक्रिया की अत्यंत प्रतिकूल पूर्वानुमानMoklyachuk & Sidei (2015-2017) : अनुपलब्ध अवलोकनों के तहत प्रक्षेप, बहिर्वेशन और फिल्टरिंगLuz & Moklyachuk (2015-2016) : स्थिर वृद्धि प्रक्रिया की अनुमानमौजूदा कार्य की तुलना में:
व्यवस्थितता : अनुपलब्ध अवलोकन बहिर्वेशन सिद्धांत का पहली बार व्यवस्थित अध्ययनपूर्णता : वर्णक्रमीय निश्चितता और अनिश्चितता दोनों स्थितियों को एक साथ संभालता हैसामान्यता : कई अनुपलब्ध पैटर्न और स्वीकार्य वर्णक्रमीय घनत्व वर्गों पर विचार करता हैव्यावहारिकता : स्पष्ट गणना सूत्र और ऑपरेटर समीकरण प्रदान करता हैसैद्धांतिक ढांचा : अनुपलब्ध अवलोकनों के तहत बहुआयामी स्थिर अनुक्रमों की बहिर्वेशन का एक संपूर्ण सैद्धांतिक प्रणाली सफलतापूर्वक स्थापित करता हैवर्णक्रमीय निश्चितता परिणाम :इष्टतम वर्णक्रमीय विशेषता ऑपरेटर समीकरण R a = B c Ra = Bc R a = B c और सूत्र (10) द्वारा निर्धारित होती है माध्य-वर्ग त्रुटि सूत्र (11) के माध्यम से सटीक रूप से गणना की जा सकती है विधि संबंधित और असंबंधित शोर दोनों पर लागू होती है वर्णक्रमीय अनिश्चितता परिणाम :सबसे प्रतिकूल वर्णक्रमीय घनत्व उप-अवकल शर्त 0 ∈ ∂ Δ D ( F 0 , G 0 ) 0 \in \partial\Delta_{\mathcal{D}}(F^0, G^0) 0 ∈ ∂ Δ D ( F 0 , G 0 ) द्वारा लक्षणित होता है कई विशेष स्वीकार्य वर्गों के लिए स्पष्ट Lagrange समीकरण दिए गए हैं अत्यंत प्रतिकूल अनुमान में सैडल बिंदु गुण होते हैं गणना विधि : फूरियर गुणांकों और ऑपरेटर मैट्रिक्स के माध्यम से एक गणनीय ढांचा कार्यान्वित किया गया हैगणना जटिलता :अनंत-आयामी ऑपरेटर समीकरण को हल करने की आवश्यकता (व्यावहारिक रूप से काटने की आवश्यकता) व्युत्क्रम ऑपरेटर B − 1 B^{-1} B − 1 की गणना कठिन हो सकती है अनुपलब्ध खंड जितने अधिक, मैट्रिक्स आयाम उतना अधिक सैद्धांतिक धारणाएं :न्यूनतमता शर्त (1) या (12) के पूरा होने की आवश्यकता ऑपरेटर B B B की व्युत्क्रमणीयता मानता है (Salehi 1979 देखें) कार्यात्मक गुणांकों को निरपेक्ष योग शर्त (3) को संतुष्ट करना चाहिए वर्णक्रमीय अनिश्चितता :केवल विशिष्ट स्वीकार्य वर्णक्रमीय घनत्व वर्गों पर विचार किया गया है सबसे प्रतिकूल वर्णक्रमीय घनत्व का संख्यात्मक समाधान जटिल हो सकता है डेटा से स्वीकार्य वर्ग का अनुमान कैसे लगाएं इस पर चर्चा नहीं की गई व्यावहारिकता :बड़े पैमाने पर संख्यात्मक प्रयोगों की कमी वास्तविक डेटा अनुप्रयोग के साथ एकीकृत नहीं अन्य विधियों के साथ संख्यात्मक तुलना की कमी पेपर द्वारा संकेत दी गई अनुसंधान दिशाएं:
एल्गोरिदम विकास :ऑपरेटर समीकरणों को हल करने के लिए कुशल संख्यात्मक एल्गोरिदम बड़े पैमाने की समस्याओं के लिए अनुमानित विधियाँ काटने के आयाम को स्वचालित रूप से चुनना सैद्धांतिक विस्तार :गैर-स्थिर अनुक्रमों का सामान्यीकरण आवधिक संबंधित अनुक्रम (आंशिक कार्य पहले से है) स्थिर वृद्धि अनुक्रम (आंशिक कार्य पहले से है) अनुप्रयोग अनुसंधान :संकेत प्रसंस्करण में वास्तविक समस्याएं वित्तीय समय श्रृंखला विश्लेषण सेंसर नेटवर्क डेटा संलयन सांख्यिकीय अनुमान :डेटा से वर्णक्रमीय घनत्व का अनुमान स्वीकार्य वर्ग चयन विधि विश्वास अंतराल और परिकल्पना परीक्षण गणितीय आधार दृढ़ : Hilbert अंतरिक्ष सिद्धांत और उत्तल अनुकूलन सिद्धांत पर आधारितप्रमाण पूर्ण : प्रमेय, निष्कर्ष तार्किक रूप से स्पष्ट, शर्तें स्पष्टप्रतीक नियमित : गणितीय प्रतीकों का उपयोग नियमित, समझने में आसानअनुपलब्ध अवलोकन प्रबंधन : अनुपलब्ध संरचना को ऑपरेटर मैट्रिक्स में चतुराई से एम्बेड करता हैअत्यंत प्रतिकूल ढांचा : वर्णक्रमीय अनिश्चितता के तहत मजबूत अनुमान को व्यवस्थित रूप से विकसित करता हैबहुआयामी सामान्यीकरण : बहुआयामी स्थिति की जटिलता को सफलतापूर्वक संभालता हैकई स्थितियाँ : संबंधित/असंबंधित शोर, शोर के साथ/बिना अवलोकन को कवर करता हैकई वर्णक्रमीय वर्ग : 8 विभिन्न स्वीकार्य वर्णक्रमीय घनत्व वर्गों पर विचार करता हैस्पष्ट सूत्र : गणनीय स्पष्ट अभिव्यक्तियाँ प्रदान करता हैऐतिहासिक संदर्भ स्पष्ट : Kolmogorov से नवीनतम कार्य तकसंदर्भ व्यापक : 41 संदर्भ पत्र शामिल हैंस्थिति सटीक : मौजूदा कार्य के साथ संबंध स्पष्ट रूप से बताता हैकेवल एक उदाहरण : उदाहरण 2.1 बहुत सरल है (द्विआयामी, सरल अनुपलब्ध पैटर्न)संख्यात्मक तुलना की कमी : अन्य विधियों के साथ संख्यात्मक तुलना नहीं की गईवास्तविक डेटा नहीं : वास्तविक डेटा सेट पर सत्यापन नहीं किया गयाप्रतीक भारी : बड़ी संख्या में मैट्रिक्स, ऑपरेटर प्रतीक, पढ़ने की दहलीज अधिकसंरचना जटिल : ऑपरेटर मैट्रिक्स की ब्लॉक संरचना का विवरण पर्याप्त रूप से सहज नहीं हैचित्र की कमी : समझने में सहायता के लिए कोई चार्ट नहींगणना लागत : एल्गोरिदम जटिलता और गणना दक्षता पर चर्चा नहीं की गईपैरामीटर चयन : स्वीकार्य वर्ग पैरामीटर के व्यावहारिक चयन के लिए कोई मार्गदर्शन नहींसॉफ्टवेयर कार्यान्वयन : कोड या सॉफ्टवेयर पैकेज नहींव्युत्क्रमणीयता धारणा : ऑपरेटर B B B की व्युत्क्रमणीयता शर्त पर्याप्त रूप से स्पष्ट नहीं हैअभिसरण : अनंत-आयामी समस्या के काटने की त्रुटि विश्लेषण की कमीस्थिरता : संख्यात्मक स्थिरता पर चर्चा नहीं की गईसैद्धांतिक योगदान : ★★★★☆अनुपलब्ध अवलोकन बहिर्वेशन सिद्धांत में खाली स्थान भरता है भविष्य के अनुसंधान के लिए व्यवस्थित ढांचा प्रदान करता है विधि नवीनता : ★★★★☆अनुपलब्ध अवलोकनों को संभालने के लिए ऑपरेटर समीकरण विधि नवीन है अत्यंत प्रतिकूल ढांचे का व्यवस्थित विकास मूल्यवान है अनुप्रयोग संभावना : ★★★☆☆सिद्धांत पूर्ण लेकिन व्यावहारिकता सत्यापन की प्रतीक्षा में है अधिक वास्तविक अनुप्रयोग मामलों की आवश्यकता है पुनरुत्पादनीयता : ★★☆☆☆सैद्धांतिक सूत्र पूर्ण लेकिन एल्गोरिदम विवरण अपर्याप्त है कोड और संख्यात्मक प्रयोगों की कमी है समय श्रृंखला विश्लेषण : अनुपलब्ध डेटा प्रबंधन के लिए सैद्धांतिक उपकरण प्रदान करता हैसंकेत प्रसंस्करण : सेंसर डेटा संलयन में लागू किया जा सकता हैवित्तीय इंजीनियरिंग : उच्च आवृत्ति व्यापार डेटा की अनुपलब्धता प्रबंधनसांख्यिकी : मजबूत अनुमान सिद्धांत का विकाससेंसर नेटवर्क : सेंसर विफलता के कारण डेटा अनुपलब्धतासंचार प्रणाली : डेटा पैकेट हानि का संकेत पुनर्निर्माणवित्तीय समय श्रृंखला : अनियमित व्यापार समय की पूर्वानुमानपर्यावरण निगरानी : मौसम स्टेशन डेटा अनुपलब्धता का प्रक्षेपगैर-स्थिर प्रक्रिया : विधि स्थिरता मानती हैगैर-रैखिक प्रणाली : केवल रैखिक कार्यात्मकों पर विचार करता हैउच्च-आयामी बड़े पैमाने : गणना जटिलता बहुत अधिक हो सकती हैपूरी तरह अज्ञात वर्णक्रमीय : कुछ पूर्व जानकारी की आवश्यकता हैसमय श्रृंखला सिद्धांत शोधकर्ता : ★★★★★व्यवस्थित सैद्धांतिक ढांचा प्रदान करता है संकेत प्रसंस्करण इंजीनियर : ★★★☆☆सिद्धांत-केंद्रित, गणित पृष्ठभूमि की आवश्यकता है सांख्यिकी शोधकर्ता : ★★★★☆मजबूत अनुमान विधि संदर्भ मूल्य है अनुप्रयुक्त डेटा वैज्ञानिक : ★★☆☆☆व्यावहारिक एल्गोरिदम और कोड की कमी है पेपर इस क्षेत्र के शास्त्रीय और अग्रणी कार्यों का हवाला देता है:
आधारभूत कार्य :Kolmogorov (1992): यादृच्छिक प्रक्रिया पूर्वानुमान सिद्धांत Wiener (1966): फिल्टरिंग और पूर्वानुमान सिद्धांत Yaglom (1987): संबंधित सिद्धांत पद्धति :Grenander (1957): अत्यंत प्रतिकूल विधि Franke (1984, 1985): मजबूत पूर्वानुमान Pshenichnyj (1971): उत्तल अनुकूलन अनुपलब्ध अवलोकन :Bondon (2002, 2005) Pourahmadi et al. (2007, 2009) लेखकों की श्रृंखला कार्य : अनुसंधान की निरंतरता और गहराई प्रदर्शित करता हैयह एक सैद्धांतिक रूप से कठोर, विधि रूप से व्यवस्थित उच्च गुणवत्ता का शैक्षणिक पेपर है। मुख्य लाभ:
अनुपलब्ध अवलोकन बहिर्वेशन समस्या का संपूर्ण सैद्धांतिक ढांचा स्थापित करता है वर्णक्रमीय निश्चितता और अनिश्चितता दोनों स्थितियों को एक साथ संभालता है कई विशेष स्थितियों के लिए स्पष्ट समाधान प्रदान करता है मुख्य कमियाँ:
प्रायोगिक सत्यापन कमजोर, केवल एक सरल उदाहरण है व्यावहारिक विचार अपर्याप्त, एल्गोरिदम और कोड की कमी है पठनीयता में सुधार की आवश्यकता है, प्रतीक भारी हैं अनुशंसा सूचकांक : ★★★★☆ (सैद्धांतिक शोधकर्ता) / ★★★☆☆ (अनुप्रयुक्त शोधकर्ता)
पेपर समय श्रृंखला विश्लेषण और मजबूत अनुमान क्षेत्र में महत्वपूर्ण सैद्धांतिक योगदान करता है, लेकिन एल्गोरिदम कार्यान्वयन और व्यावहारिक अनुप्रयोग पहलुओं में पूरक और सत्यापन के लिए बाद के कार्य की आवश्यकता है।