svd

用迭代优化做奇异值分解。 之前翻了无数本书都没找到怎么对标准正交矩阵做约束，只好自己想想，然后想起了标准正交矩阵左乘自己的转置是单位矩阵。