概念讲解
公式
公式里的
U就是下面运行结果的a$u,
Σ就是下面的a$d变换成的对角矩阵(diag(a$d)),
V就是下面运行结果的a$v
生成例子
x <- matrix(1:20,nrow = 4,ncol = 5)
生成了一个4✖5的矩阵,下面对这个矩阵进行奇异值分解
奇异值分解
a <- svd(x)
a
把奇异值分解得到的这三部分拿出来看一下
a$d
a$u
a$v
将这三部分运算,然后和x进行一波比较
a$u%*%diag(a$d)%*%t(a$v)
x
第一个式子的运行结果和x是一样的