Курсовые / curstau / CANON

                    { Љ ­®­ЁзҐбЄЁҐ д®а¬л }

 unit CANON;
   interface
       uses ALG,GB,ENEX;
       procedure CONT (n    :integer;A,B      :matr;var Sco        :matr);
       procedure OBS  (n    :integer;A,C      :matr;var Sob        :matr);
       procedure CFCON(n,m,r:integer;A,B,C,Sco:matr;var Au, Bu, Cu :matr);
       procedure CFOBS(n,m,r:integer;A,B,C,Sob:matr;var Aob,Bob,Cob:matr);


   implementation

                    {ўлзЁб«Ґ­ЁҐ ¬ ваЁжл гЇа ў«пҐ¬®бвЁ }

     procedure CONT(n:integer;A,B:matr;var Sco:matr);

              var i,j,k :integer;           {                 n-1   }
                  E,T,B1:matr;              { Sco=[B: A*B:..:A  *B] }
     begin
       for i:=1 to n do for j:=1 to n do
       begin
         E[i,j]:=0; E[i,i]:=1;
         B1[i,j]:=B[i,j];
       end;
       k:=1;
      repeat
       MUL(n,n,1,E,B1,T);
       for i:=1 to n do  Sco[i,k]:=T[i,1];
       for i:=1 to n do for j:=1 to n do
       begin
         B1[i,j]:=T[i,j];
         E[i,j]:=A[i,j];
       end;
       k:=k+1;
      until k=n+1;
     end;
     {------------------------------------------}


                  {ўлзЁб«Ґ­ЁҐ ¬ ваЁжл ­ Ў«о¤ Ґ¬®бвЁ  }

     procedure OBS(n:integer;A,C:matr;var Sob:matr);
              var i,j,k      :integer;
                  E,T,C1     :matr;
     begin
       for i:=1 to n do for j:=1 to n do
       begin
         E[i,j]:=0; E[i,i]:=1;
         C1[i,j]:=C[i,j];
       end;
       k:=1;
      repeat
       MUL(1,n,n,C1,E,T);
       for j:=1 to n do  Sob[k,j]:=T[1,j];
       for i:=1 to n do for j:=1 to n do
       begin
         C1[i,j]:=T[i,j];
         E[i,j]:=A[i,j];
       end;
       k:=k+1;
      until k=n+1;
     end;
     {------------------------------------------}


            {ЇаҐ®Ўа §®ў ­ЁҐ Є Є ­®­. д®а¬Ґ гЇа ў«пҐ¬®бвЁ}

     procedure CFCON (n,m,r:integer;A,B,C,Sco:matr;var Au, Bu, Cu :matr);
              var i,j,k          :integer;
                  E,O,T,Tk,ISco  :matr;
     begin
       INV(n,Sco,ISco);
       for j:=1 to n do O[i,j]:=0;
       MUL(1,n,n,O,ISco,T);
       for i:=1 to n do for j:=1 to n do
       if i=j then E[i,j]:=1 else E[i,j]:=0;
       k:=1;
      repeat
       MUL(1,n,n,T,E,Tk);
       for j:=1 to n do  Sco[k,j]:=Tk[k,j];
       for i:=1 to n do for j:=1 to n do E[i,j]:=A[i,j];
       k:=k+1;
      until k=n+1;
     end;
     {------------------------------------------}


               {ЇаҐ®Ўа §®ў ­ЁҐ Є Є ­®­. д®а¬Ґ ­ Ў«о¤ Ґ¬®бвЁ}

     procedure CFOBS (n,m,r:integer;A,B,C,Sob:matr;var Aob,Bob,Cob:matr);
              var ISob,Sa  :matr;
     begin
       INV(n,Sob,ISob);
       MUL(n,n,n,Sob,A,Sa);
       MUL(n,n,n,Sa,ISob,Aob);
       MUL(n,n,1,Sob,B,Bob);
       MUL(1,n,n,C,ISob,Cob);
     end;
     {------------------------------------------}
   end.