nia.m

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

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