fullrfk

fullrfk — full rank factorization of A^k

Calling sequence

[Bk,Ck]=fullrfk(A,k)  

Parameters

A : real or complex matrix
k : integer
Bk,Ck : real or complex matrices

Description

This function computes the full rank factorization of A^k i.e. Bk*Ck=A^k where Bk is full column rank and Ck full row rank. One has range(Bk)=range(A^k) and ker(Ck)=ker(A^k).

For k=1, fullrfk is equivalent to fullrf.

Examples



A=rand(5,2)*rand(2,5);[Bk,Ck]=fullrfk(A,3);
norm(Bk*Ck-A^3,1)
 
  

See also

fullrf, range

Author

F.D (1990);