svd_equation.m

来自「利用矩阵分析的方法对线性方程组进行求解分析」· M 代码 · 共 23 行

function x=SVD_equation(m,n,A,b)
%用奇异值分解求解病态线性方程组
%A为系数矩阵,b为右向量
%m=length(A);%返回行数
epsilon=1e-10;
x=zeros(n,1);
[U,S,V]=svd(A);
sdiag=diag(S);
for i=1:(n-0)
    if sdiag(i)>epsilon
           x=x+(U(:,i)'*b)/sdiag(i)*V(:,i);
   end
end
%function x=svd_equation(A,B)
%eps=1e-10;n=length(A);x=zeros(n,1);
%[u,s,v]=svd(A);sigma=diag(s);
%for i=1:n
   % if abs(sigma(i))>=eps
      %  x=x+(u(:,i)'*B)/sigma(i)*v(:,i);
      %end 
      %end