实验四 矩阵奇异值分解.docVIP

实验四 矩阵的奇异值分解 1、原理 设A∈Cm×n，s1，s2，…，sr是A的非零奇异值，则存在m阶酉矩阵U∈Cm×n及n阶酉矩阵V，m×n矩阵D， D= = 使得A=UDVH这就是矩阵A的奇异值分解。 2、算法 第一步：求出AHA的特征值≥≥…≥＞0==…=，确定非零奇异值=，i=1,2，…，r。 第二步：分别求出矩阵AHA的对应于特征值的特征向量并将其单位正交化，得到标准正交向量组α1，α2，…，αn令V=（α1，α2，…，αn）=（V1，V2），V1=（α1，α2，…，αr），V2=（αr+1，αr+2，…，αn） 第三步：若U=（γ1，γ2，…，γr，γr+1，γr+2，…，γm）=（U1，U2），其中U1=（γ1，γ2，…，γr），U2=（γr+1，γr+2，…，γm）， 则因（Aα1，Aα2，…，Aαr）=（s1γ1，s2γ2，…，srγr） 即有U1=AV1 。其中= 第四步：解方程组AAHy=0，对基础解系单位正交化可以求得γr+1，γr+2，…，γm，令U=（γ1，γ2，…，γr，γr+1，γr+2，…，γm）。 3、程序 在matlab中有求解矩阵奇异值分解的svd函数，调用格式为[U，S，V] = svd(A)，其中U就是所求的U矩阵，S是所求的对角阵，V就是所求的酉矩阵V。 4、例子 以教材上的A=[1 0；0 1；1 0]为例来说明上述求矩阵的奇异值