Primer proyeccion Gustavo Araujo 9no A

zt=[ 10056 18573 20195 21268 23635 23255 16675 15934 21250 24899 28636 32642 37187 41763 45789 48508]';

n=length(zt)
zt1=[0;zt(1:n-1)]
zt2=[0;0;zt(1:n-2)]
zt3=[0;0;0;zt(1:n-3)]
zt4=[0;0;0;0;zt(1:n-4)]
zt5=[0;0;0;0;0;zt(1:n-5)]
zt6=[0;0;0;0;0;0;zt(1:n-6)]
zt7=[0;0;0;0;0;0;0;0zt(1:n-7)]
x=[zt1,zt2,zt3,zt4,zt5,zt6;zt7]

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))
bet=imag(raiz(1))
r=sqrt(alfa^2+bet^2)
teta=acos(alfa/r)
t=[1:n]';
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:n]';
z1=(r.^t).*cos(teta*t)
z2=(r.^t).*sin(teta*t)
Z= A(1)*z1+A(2)*z2
plot2d(t,Z)