Dugaan saya adalah bahwa persamaan karakteristik yang Anda gunakan berbeda dengan persamaan saya. Biarkan saya melanjutkan dalam beberapa langkah untuk melihat apakah kami setuju.
Pertimbangkan persamaan
λ2−ϕ1λ−ϕ2=0
Jika adalah akar dari persamaan karakteristik "standar" dan pengaturan , tampilan diperoleh dari penulisan ulang yang standar sebagai berikut:
Oleh karena itu, kondisi alternatif untuk stabilitas adalah bahwa semua akar tampilan pertama berada di dalam lingkaran unit, .z1−ϕ1z−ϕ2z2=0z−1=λ1−ϕ1z−ϕ2z2⇒z−2−ϕ1z−1−ϕ2⇒λ2−ϕ1λ−ϕ2===000
AR(2)|z|>1⇔|λ|=|z−1|<1
Kami menggunakan representasi ini untuk memperoleh segitiga stasioneritas dari proses , yaitu stabil jika tiga kondisi berikut terpenuhi: AR(2)AR(2)
- ϕ2<1+ϕ1
- ϕ2<1−ϕ1
- ϕ2>−1
Ingatlah bahwa Anda dapat menulis akar tampilan pertama (jika nyata) sebagai
untuk menemukan dua kondisi pertama.λ1,2=ϕ1±ϕ21+4ϕ2−−−−−−−√2
Kemudian, adalah stasioner iff , maka (jika nyata):
Semakin besar keduanya dibatasi oleh , atau:
Secara analog, .AR(2)|λ|<1λi−1<ϕ1±ϕ21+4ϕ2−−−−−−−√2⇒−2<ϕ1±ϕ21+4ϕ2−−−−−−−√<<12
λiϕ1+ϕ21+4ϕ2−−−−−−−√<2ϕ1+ϕ21+4ϕ2−−−−−−−√⇒ϕ21+4ϕ2−−−−−−−√⇒ϕ21+4ϕ2⇒ϕ21+4ϕ2⇒ϕ2<<<<<22−ϕ1(2−ϕ1)24−4ϕ1+ϕ211−ϕ1
ϕ2<1+ϕ1
Jika kompleks, maka dan begituModulus kuadrat dari bilangan kompleks adalah kuadrat dari real ditambah kuadrat dari bagian imajiner. Karenanya,
Ini stabil jika , maka jika atau , seperti yang akan ditampilkan. (Pembatasan dihasilkan dari adalah redundan mengingat dan .)λiϕ21<−4ϕ2λ1,2=ϕ1/2±i−(ϕ21+4ϕ2)−−−−−−−−−−√/2.
λ2=(ϕ1/2)2+(−(ϕ21+4ϕ2)−−−−−−−−−−√/2)2=ϕ21/4−(ϕ21+4ϕ2)/4=−ϕ2.
|λ|<1−ϕ2<1ϕ2>−1ϕ2<1ϕ22<1ϕ2<1+ϕ1ϕ2<1−ϕ1
Merencanakan segitiga stasioneritas, juga menunjukkan garis yang memisahkan kompleks dari akar nyata, kita dapatkan
Diproduksi dalam R menggunakan
phi1 <- seq(from = -2.5, to = 2.5, length = 51)
plot(phi1,1+phi1,lty="dashed",type="l",xlab="",ylab="",cex.axis=.8,ylim=c(-1.5,1.5))
abline(a = -1, b = 0, lty="dashed")
abline(a = 1, b = -1, lty="dashed")
title(ylab=expression(phi[2]),xlab=expression(phi[1]),cex.lab=.8)
polygon(x = phi1[6:46], y = 1-abs(phi1[6:46]), col="gray")
lines(phi1,-phi1^2/4)
text(0,-.5,expression(phi[2]<phi[1]^2/4),cex=.7)
text(1.2,.5,expression(phi[2]>1-phi[1]),cex=.7)
text(-1.75,.5,expression(phi[2]>1+phi[1]),cex=.7)