Scilab Website | Contribute with GitLab | Mailing list archives | ATOMS toolboxes
Scilab Online Help
2025.0.0 - English


range

range (span) of A^k

Syntax

[X,dim]=range(A,k)

Arguments

A

real square matrix

k

integer

X

orthonormal real matrix

dim

integer (dimension of subspace)

Description

Computation of Range A^k ; the first dim rows of X span the range of A^k. The last rows of X span the orthogonal complement of the range. X*X' is the Identity matrix

Examples

A=rand(4,2)*rand(2,4);   // 4 column vectors, 2 independent.
[X,dim]=range(A,1);dim   // compute the range

y1=A*rand(4,1);          //a vector which is in the range of A
y2=rand(4,1);            //a vector which is not  in the range of A
norm(X(dim+1:$,:)*y1)    //the last entries are zeros, y1 is in the range of A
norm(X(dim+1:$,:)*y2)    //the last entries are not zeros

I=X(1:dim,:)'            //I is a basis of the range
coeffs=X(1:dim,:)*y1     // components of y1 relative to the I basis

norm(I*coeffs-y1)        //check

See also

  • fullrfk — full rank factorization of A^k
  • rowcomp — row compression, range

Used Functions

The range function is based on the rowcomp function which uses the svd decomposition.

Report an issue
<< kernel Kernel rowcomp >>

Copyright (c) 2022-2024 (Dassault Systèmes)
Copyright (c) 2017-2022 (ESI Group)
Copyright (c) 2011-2017 (Scilab Enterprises)
Copyright (c) 1989-2012 (INRIA)
Copyright (c) 1989-2007 (ENPC)
with contributors
Last updated:
Thu Oct 24 11:13:08 CEST 2024