xiuzhenga.m

利用矩阵分析的方法对线性方程组进行求解分析

%矩阵A的逆niA，与W,B都有关，并且与才采点，采样范围有关
%t=linspace(pi/60,pi,15)
%y=zongjie(L)*dierbufen(6-L,t[i])*disanbufen(L,t[i])
function[A3]=xiuzhengA(A,m,n);   %matrixA函数表示一个A矩阵。
%m=20;
%n=7;
[U S V]=svd(A);
sdiag=diag(S);
A3=zeros(m,n);%先行分配一个m*n大的0矩阵，有利于提高运算速度。
epsilon=1e-10;
for i=1:(n-4)
    if sdiag(i)>epsilon
          A3=A3+U(:,i)*sdiag(i)*(V(:,i))';
   end
end