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);