Pronostico del PIB

-Pronostico del 2008
zt=[21250 24899 28636 32642 37187 41763 45789 48508]';
n=length(zt)
zt1=[0;zt(1:n-1)]
zt2=[0;0;zt(1:n-2)]
x=[zt1,zt2]
b=inv(x'*x)*(x'*zt)
s=poly(0,'s')
s=s^2-b(1)*s-b(2)
raiz=roots(s)
alfa=real(raiz(1))
beta=imag(raiz(1))
r=sqrt(alfa^2+beta^2)
teta=acos(alfa/r)
t=[1:n]'
z1=(raiz(1).^t)
z2=(raiz(2).^t)
x=[z1 z2]
a=inv(x'*x)*(x'*zt)
z=a(1)*z1+a(2)*z2
plot(t,[zt z])
plot(t,zt)
e=zt-z
sqrt(e'*e)
plot(t,e)
t=[1:n+1]';
z1=(raiz(1).^t)
z2=(raiz(2).^t)
z=a(1)*z1+a(2)*z2
plot(t,z)


-Pronostico de PIB 1998
vector=[18573 20195 21268 23635 23255]';
n=length(V)
v1=[0;V(1:n-1)]
v2=[0;0;V(1:n-2)]
x=[v1,v2]b=inv(x'*x)*(x'*V)
s=poly(0,'s')
s=s^2-b(1)*s-b(2)
raiz=roots(s)t=[1:n]'
v1=(raiz(1).^t)
v2=(raiz(2).^t)
x=[v1 v2]
a=inv(x'*x)*(x'* V)
z=a(1)*v1+a(2)*v2
plot(t, V)
e=V-z
sqrt(e'*e)
plot(e)
t=[1:n+1]';
v1=(raiz(1).^t)
v2=(raiz(2).^t)
z=a(1)*v1+a(2)*v2
plot(t,z)