[kristxpher]

Este es el problema que tengo que resolver. Hasta el momento tengo este algoritmo pero no sé como modificarlo para LTI. Ayuda por favor.

Algoritmo de Arnoldi:

Código: [Seleccionar]

function [Q,H]=arnoldi(A,q,m)
    n=size(A,1);
    Q=zeros(n,m+1);
    H=zeros(m+1,m);
    Q(:,1)=q/norm(q);
    for j=1:m
        z=A*Q(:,j);
    end
    for i=1:j
        H(i,j)=dot(z,Q(:,i));
        z=z-H(i,j)*Q(:,i);
    end
    H(j+1,j)=norm(z);
    if H(j+1,j)==0, break, end
    Q(:,j+1)=z/H(j+1,j);
end