agai.m

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

function[A]=Agai(m,n);   %matrixA函数表示一个A矩阵。
m=15;
n=7;
W=10;
B=3;%所给的仿真条件
t=linspace(pi/(60*B),pi/2,15);
%测试矩阵A
%t=[-pi/(30*B):pi/(210*B):pi/(30*B)];

A=zeros(m,n);%先行分配一个m*n大的0矩阵，有利于提高运算速度。

  y=0; 
for i=1:m 
    %y=0;
    j=7
    %y=0;
         %for j=1:n
              
        %A(i,j)=(cos(W*t(i))*sin(B*t(i))/(pi*t(i))).^(j-1)%利用for循环来构造这样一个m*n的矩阵
       
                        for L=0:1:(j-1) 
                             %y=zongjie(L)*dierbufen(6-L,t[i])*disanbufen(L,t[i])+y
                              y=zongjie(L,j)*dierbufen(j-1-L,t(i))*disanbufen(L,t(i))+y;
                        end
                        A(i,j)=y;
               
                        %end
    %y=0;
end