Autocovariance của một (2,1) quá trình ARMA - nguồn gốc của mô hình phân tích cho


13

Tôi cần rút ra các biểu thức phân tích cho hàm tự động điều khiển của quy trình ARMA (2.1) được ký hiệu là:γ ( k )γ(k)

y t = φ 1 y t - 1 + φ 2 y t - 2 + θ 1 ε t - 1 + ε tyt=ϕ1yt−1+ϕ2yt−2+θ1ϵt−1+ϵt

Vì vậy, tôi biết rằng:

γ ( k ) = E [ y t , y t - k ]γ(k)=E[yt,yt−k]

để tôi có thể viết:

γ ( k ) = φ 1 E [ y t - 1 y t - k ] + φ 2 E [ y t - 2 y t - k ] + θ 1 E [ ε t - 1 y t - k ] + E [ ε t y t - k ]γ(k)=ϕ1E[yt−1yt−k]+ϕ2E[yt−2yt−k]+θ1E[ϵt−1yt−k]+E[ϵtyt−k]

sau đó, để lấy phiên bản phân tích của hàm tự động điều khiển, tôi cần thay thế các giá trị của - 0, 1, 2 ... cho đến khi tôi nhận được đệ quy có giá trị với tất cả lớn hơn một số nguyên.k kkk

Do đó, tôi thay thế và thực hiện điều này để có được:k = 0k=0

γ ( 0 ) = E [ y t , y t ] = ϕ 1 E [ y t - 1 y t ] + ϕ 2 E [ y t - 2 y t ] + θ 1 E [ ϵ t - 1 y t ] + E [ ε t y t ]

γ(0)=E[yt,yt]=ϕ1E[yt−1yt]+ϕ2E[yt−2yt]+θ1E[ϵt−1yt]+E[ϵtyt]

bây giờ tôi có thể đơn giản hóa hai thuật ngữ đầu tiên và sau đó thay thế cho như trước:y tyt

γ(0)=ϕ1γ(1)+ϕ2γ(2)+θ1E[ϵt−1(ϕ1yt−1+ϕ2yt−2+θ1ϵt−1+ϵt)]+E[ϵt(ϕ1yt−1+ϕ2yt−2+θ1ϵt−1+ϵt)]

γ(0)=ϕ1γ(1)+ϕ2γ(2)+θ1E[ϵt−1(ϕ1yt−1+ϕ2yt−2+θ1ϵt−1+ϵt)]+E[ϵt(ϕ1yt−1+ϕ2yt−2+θ1ϵt−1+ϵt)]

sau đó tôi nhân ra tám điều khoản, đó là:

+θ1ϕ1E[ϵt−1yt−1]+θ1ϕ2E[ϵt−1yt−2]+θ21E[(ϵt−1)2]=θ21σ2ϵ+θ1E[ϵt−1ϵt]=θ1E[ϵt−1]E[ϵt]=0+ϕ1E[ϵtyt−1]+ϕ2E[ϵtyt−2]+θ1E[ϵtϵt−1]=θ1E[ϵt]E[ϵt−1]=0+E[(ϵt)2]=σ2ϵ

+θ1ϕ1E[ϵt−1yt−1]+θ1ϕ2E[ϵt−1yt−2]+θ21E[(ϵt−1)2]=θ21σ2ϵ+θ1E[ϵt−1ϵt]=θ1E[ϵt−1]E[ϵt]=0+ϕ1E[ϵtyt−1]+ϕ2E[ϵtyt−2]+θ1E[ϵtϵt−1]=θ1E[ϵt]E[ϵt−1]=0+E[(ϵt)2]=σ2ϵ

Vì vậy, tôi còn lại cần phải giải quyết bốn điều khoản còn lại. Tôi muốn sử dụng logic tương tự cho các dòng 1, 2, 5 và 6 như tôi đã sử dụng trên các dòng 4 và 7 - ví dụ cho dòng 1:

θ1ϕ1E[ϵt−1yt−1]=θ1ϕ1E[ϵt−1]E[yt−1]=0θ1ϕ1E[ϵt−1yt−1]=θ1ϕ1E[ϵt−1]E[yt−1]=0 vì .E[ϵt−1]=0E[ϵt−1]=0

Tương tự cho các dòng 2, 5 và 6. Nhưng tôi có một giải pháp mô hình cho thấy biểu thức cho đơn giản hóa thành:γ(0)γ(0)

γ(0)=ϕ1γ(1)+ϕ2γ(2)+θ1(ϕ1+θ1)σ2ϵ+σ2ϵγ(0)=ϕ1γ(1)+ϕ2γ(2)+θ1(ϕ1+θ1)σ2ϵ+σ2ϵ

Điều này cho thấy sự đơn giản hóa của tôi như được mô tả ở trên sẽ bỏ lỡ thuật ngữ với hệ số - theo logic của tôi là 0. Có phải logic của tôi bị lỗi hay là giải pháp mô hình mà tôi thấy không chính xác?ϕ1ϕ1

Giải pháp đã thực hiện cũng đề xuất rằng "tương tự" có thể được tìm thấy là:γ(1)γ(1)

γ(1)=ϕ1γ(0)+ϕ2γ(1)+θ1σ2ϵγ(1)=ϕ1γ(0)+ϕ2γ(1)+θ1σ2ϵ

và cho :k>1k>1

γ(k)=ϕ1γ(k−1)+ϕ2(k−2)γ(k)=ϕ1γ(k−1)+ϕ2(k−2)

Tôi hy vọng câu hỏi là rõ ràng. Bất kỳ trợ giúp sẽ được nhiều đánh giá cao. Cảm ơn bạn trước.

Đây là một câu hỏi liên quan đến nghiên cứu của tôi, và không chuẩn bị cho bất kỳ kỳ thi hoặc khóa học.

Câu trả lời:


8

Nếu quy trình ARMA là nhân quả, có một công thức chung cung cấp các hệ số tự động hóa.

Hãy xem xét các nguyên nhân quá trình trong đó là tiếng ồn trắng có giá trị trung bình bằng 0 và phương sai . Theo thuộc tính nhân quả, quy trình có thể được viết là trong đó biểu thị -weights.ARMA(p,q)ARMA(p,q)yt=p∑i=1ϕiyt−1+q∑j=1θjϵt−j+ϵt,

yt=∑i=1pϕiyt−1+∑j=1qθjϵt−j+ϵt,
ϵtϵtσ2ϵσ2ϵyt=∞∑j=0ψjϵt−j,
yt=∑j=0∞ψjϵt−j,
ψjψjψψ

Phương trình đồng nhất chung cho các hệ số tự động tương quan của một quá trình nguyên nhân là với các điều kiện ban đầu ARMA(p,q)ARMA(p,q)γ(k)−ϕ1γ(k−1)−⋯−ϕpγ(k−p)=0,k≥max(p,q+1),

γ(k)−ϕ1γ(k−1)−⋯−ϕpγ(k−p)=0,k≥max(p,q+1),
γ(k)−p∑j=1ϕjγ(k−j)=σ2ϵq∑j=kθjψj−k,0≤k<max(p,q+1).
γ(k)−∑j=1pϕjγ(k−j)=σ2ϵ∑j=kqθjψj−k,0≤k<max(p,q+1).

2

Lỗi tính toán của bạn trong câu hỏi ban đầu của bạn nằm ở

θ1ϕ1E[ϵt−1yt−1]=θ1ϕ1E[ϵt−1]E[yt−1]=0(mistaken)

θ1ϕ1E[ϵt−1yt−1]=θ1ϕ1E[ϵt−1]E[yt−1]=0(mistaken)

Bạn không thể tách kỳ vọng - và không độc lập.E[ϵt−1yt−1]E[ϵt−1yt−1]ϵt−1ϵt−1yt−1yt−1


Như bạn có thể thấy từ bản cập nhật của tôi (bên dưới) tôi đã nhận ra điều này ngay sau khi hoàn thành bài đăng - nhưng cảm ơn rất nhiều vì sự giúp đỡ của bạn!
— thủy văn

1

ĐỒNG Ý. Vì vậy, quá trình viết bài thực sự chỉ cho tôi giải pháp.

Hãy xem xét các thuật ngữ Kỳ vọng 1, 2, 5 và 6 từ phía trên mà tôi nghĩ nên là 0.

Ngay lập tức cho các điều khoản 5 - - và 6 - : các thuật ngữ này chắc chắn bằng không, vì và độc lập với và .E[ϵtyt−1]E[ϵtyt−1]E[ϵtyt−2]E[ϵtyt−2]yt−1yt−1yt−2yt−2ϵtϵtE[ϵt]=0E[ϵt]=0

Tuy nhiên, các điều khoản 1 và 2 trông như thể Kỳ vọng có hai biến tương quan. Vì vậy, hãy xem xét các biểu thức cho và do đó:yt−1yt−1yt−2yt−2

yt−1=ϕ1yt−2+ϕ2yt−3+θ1ϵt−2+ϵt−1yt−2=ϕ1yt−3+ϕ2yt−4+θ1ϵt−3+ϵt−2

yt−1=ϕ1yt−2+ϕ2yt−3+θ1ϵt−2+ϵt−1yt−2=ϕ1yt−3+ϕ2yt−4+θ1ϵt−3+ϵt−2

Và nhớ lại thuật ngữ 1 - . Nếu chúng ta nhân cả hai mặt của biểu thức cho với và sau đó lấy Kỳ vọng, thì rõ ràng tất cả các thuật ngữ ở phía bên phải ngoại trừ giá trị cuối cùng trở thành số không (vì các giá trị của , và độc lập với và ) để cung cấp:ϕ1θ1E[ϵt−1yt−1]ϕ1θ1E[ϵt−1yt−1]yt−1yt−1ϵt−1ϵt−1yt−2yt−2yt−3yt−3ϵt−2ϵt−2ϵt−1ϵt−1E[ϵt−1]=0E[ϵt−1]=0

E[ϵt−1yt−1]=E[(ϵt−1)2]=σ2ϵ

E[ϵt−1yt−1]=E[(ϵt−1)2]=σ2ϵ

Vì vậy, thuật ngữ 1 trở thành . Đối với thuật ngữ 2, cần phải rõ ràng rằng, theo cùng một logic, tất cả các thuật ngữ đều bằng không.+ϕ1θ1σ2ϵ+ϕ1θ1σ2ϵ

Do đó, câu trả lời mô hình ban đầu là chính xác.

Tuy nhiên, nếu bất cứ ai có thể đề xuất một cách khác để có được một giải pháp chung (ngay cả khi lộn xộn), tôi sẽ rất vui khi nghe nó!

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookie và Chính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.