matlab警告出现奇异精度

最新推荐文章于 2023-10-29 21:16:11 发布

找不到活命的氧气

最新推荐文章于 2023-10-29 21:16:11 发布

阅读量1.2w

点赞数 3

文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/u011059424/article/details/127377338

版权

该博客介绍了如何处理矩阵A在使用A运算时出现奇异精度警告的问题。通过引入Matlab代码，利用奇异值分解(SVD)和伪逆阵(pinv)来消除这种警告。当矩阵A不满秩时，伪逆阵能够确保计算结果的准确性。提供的svd_solve函数实现了一个自定义的解法，其核心是迭代地应用SVD更新解向量，直到达到预设的精度标准。

摘要由CSDN通过智能技术生成

针对矩阵A和矩阵b，使用A\b后出现警告奇异精度，可使用如下matlab代码消除。具体原因查看其他博主写的关于伪逆阵pinv，写的很详细。

使用伪逆阵，A为满秩时，跟上面A\b出现的结果一样，但A的秩不为满秩时，可以防止出现矩阵奇异精度。

function varargout = svd_solve(A,b)
N = length(b);
if nargin < 3
ep = 1e-8;
end
[U,S,V] = svd(A);
sigma = diag(S);
x = zeros(N,1);
k = 1;
while 1
if abs(sigma(k))<ep
A=A-sigma(k)*U(:,k)*V(:,k)';
else
x=x+(U(:,k)'*b)/sigma(k)*V(:,k);
end
k=k+1;
if k>N
break
end
end
[varargout{1:2}]=deal(x,A);