weinia.m

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

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