¿Cómo puedo modificar el algoritmo de Iteración de Arnoldi para que lea una LTI realization?
Publicado por Alexis (1 intervención) el 05/05/2018 19:37:29
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:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
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
Valora esta pregunta
0