1
<?xml version="1.0" encoding="ISO-8859-1" standalone="no"?>
2
<!DOCTYPE MAN SYSTEM "../../manrev.dtd">
4
<LANGUAGE>eng</LANGUAGE>
6
<TYPE>Scilab Function</TYPE>
7
<DATE>April 1993</DATE>
8
<SHORT_DESCRIPTION name="fullrf"> full rank factorization</SHORT_DESCRIPTION>
10
<CALLING_SEQUENCE_ITEM>[Q,M,rk]=fullrf(A,[tol]) </CALLING_SEQUENCE_ITEM>
15
<PARAM_NAME>A</PARAM_NAME>
17
<SP>: real or complex matrix</SP>
21
<PARAM_NAME>tol</PARAM_NAME>
23
<SP>: real number (threshold for rank determination)</SP>
27
<PARAM_NAME>Q,M</PARAM_NAME>
29
<SP>: real or complex matrix</SP>
33
<PARAM_NAME>rk</PARAM_NAME>
35
<SP>: integer (rank of <VERB>A</VERB>)</SP>
42
Full rank factorization : <VERB>fullrf</VERB> returns <VERB>Q</VERB> and <VERB>M</VERB> such
43
that <VERB>A = Q*M</VERB>
44
with range(<VERB>Q</VERB>)=range(<VERB>A</VERB>) and ker(<VERB>M</VERB>)=ker(<VERB>A</VERB>),
45
<VERB>Q</VERB> full column rank , <VERB>M</VERB> full row rank,
46
<VERB>rk</VERB> = rank(<VERB>A</VERB>) = #columns(<VERB>Q</VERB>) = #rows(<VERB>M</VERB>).</P>
47
<P><VERB>tol</VERB> is an optional real parameter (default value is <VERB>sqrt(%eps)</VERB>).
48
The rank <VERB>rk</VERB> of <VERB>A</VERB> is defined as the number of singular values
49
larger than <VERB>norm(A)*tol</VERB>.</P>
51
If A is symmetric, <VERB>fullrf</VERB> returns <VERB>M=Q'</VERB>.</P>
55
A=rand(5,2)*rand(2,5);
58
[X,d]=rowcomp(A);Y=X';
59
svd([A,Y(:,1:d),Q]) //span(Q) = span(A) = span(Y(:,1:2))
79
<AUTHOR>F.D.; </AUTHOR>