pib=[2.814 3.115 3.366]'
n=length(pib)
pib1=[0;pib(1:n-1)]
pib2=[0;0;pib(1:n-2)]
x=[pib1 pib2]
b=inv(x'*x)*(x'*pib)
s=poly(0,'s')
s=s^2-b(1)*s-b(2)
raiz=roots(s)
t=[1:n]'
p1=exp(raiz(1)*t)
p2=exp(raiz(2)*t)
x=[p1 p2]
A=inv(x'*x)*(x'*pib)
z=A(1)*p1+A(2)*p2
plot(t,pib)
t=[1:0.1:n+1]'
z1=exp(raiz(1)*t)
z2=exp(raiz(2)*t)
z=A(1)*z1+A(2)*z2
plot(t,z)
No hay comentarios:
Publicar un comentario