Santiago Feijó (Deber PIB)

Pronóstico del PIB desde el año 1993 hasta 2007

--Años desde 1993 hasta 2007
zt=[18573 20195 21268 23635 23255 16675 15934 21250 24899 28636 32642 37187 41763 45789 48508]';
q=length(zt)
zt1=[0;zt(1:q-1)]
zt2=[0;0;zt(1:q-2)]
zt3=[0;0;0;zt(1:q-3)]
x=[zt1,zt2,zt3,zt4]
c=inv(x'*x)*(x'*zt)
R=poly(0,'R')
R=R^2-c(1)*R-c(2)
rz=roots(R)
Alpha=real(rz(1))
Beta=imag(rz(1))
r=sqrt(Alpha^2+Beta^2)
Teta=acos(Alpha/r)
t=[1:q]';
z1=(r.^t).*cos(Teta*t)
z2=(r.^t).*sin(Teta*t)
X=[z1,z2]
a=inv(X'*X)*(X'*zt)
Z=a(1)*z1+a(2)*z2
plot2d(t,[zt Z])
t=[1:0.1:q]';
z1=(r.^t).*cos(Teta*t)
z2=(r.^t).*sin(Teta*t)
Z= a(1)*z1+a(2)*z2
plot2d(t,Z)