[silas_xue]

lz 代码中矩阵乘法有误 具体位置在New=New+D(k,k)*(U(:,k)*V(:,k));
(U(:,k)*V(:,k))中U的列数和V的行数不相等 是不能做乘法的
根据我对SVD分解和阅读lz的理解 U(:,k)改为U(:,k)' 但是没有在MatlAB中调试
lz调试看看

Thx for reading.
PS:若还算满意，直接点击“Thanks”，再次登陆时亦便于查看回答是否真的帮到你了。
个人观点 仅供参考 多多交流 相互学习