pracma (version 0.7.5)

pinv: Pseudoinverse or Generalized Inverse

Description

Compute the pseudoinverse of a matrix.

Usage

pinv(A)

Arguments

A

matrix

Value

• The pseudoinverse of matrix A.

Details

Compute the pseudoinverse B of a matrix A using the singular value decomposition svd().

The pseudoinverse $B$ solves the problem to minimize $|A x - b|$ by setting $x = B b$

s <- svd(A) D <- diag(s$d); Dinv <- diag(1/s$d) U <- s$u; V <- s$v U <- s$u; V <- s$v X = V Dinv U'

Thus B is computed as s$v %*% diag(1/s$d) %*% t(s$u).

References

Ben-Israel, A., and Th. N. E. Greville (2003). Generalized Inverses - Theory and Applications. Springer-Verlag, New York.

See Also

MASS::ginv

Examples

A <- matrix(c(7,6,4,8,10,11,12,9,3,5,1,2), 3, 4)
b <- apply(A, 1, sum)  # 32 16 20  row sum
x <- pinv(A) %*% b
A %*% x              #=> 32 16 20  as column vector