Форум программистов [Powered by Invision Power Board]

Код


type matrix = array[1..100,1..100] of integer;
............
{нахождение количества 0}
function z(a:matrix;k,l:integer):integer;
var c,i,j:integer;
begin
  c:=0;
  for I := 1 to k  do
    for J := 1 to l do
      if a[i,j]=0 then c:=c+1;
end;
{умножение матрицы на число }
function scmult(a:matrix;k,l:integer;sc:integer) :matrix;
 var i,j:integer;
 b:matrix;
begin
  for I := 1 to k  do
    for J := 1 to l do
     b[i,j]:=a[i,j]*sc;
end;
 {ввод матрицы}
procedure Vvod(var mt:matrix;x,y:byte);
var i,j:byte;
begin
for i:=1 to x do
for j:=1 to y do
mt[i,j]:=random(x*y);
end;
{вывод матриц на экран}
procedure Vyvod(var mt:matrix;x,y:byte;);
var i,j:byte;
begin
for i:=1 to x do
 begin
  for j:=1 to y do
  write(mt[i,j]:4);
  writeln;
 end;
end;
{траеспонирование матриц}
procedure Trans(m1:matrix;x,y:byte;var m2:matrix);
var i,j:byte;
begin
for i:=1 to y do
for j:=1 to x do
m2[i,j]:=m1[j,i];
end;
{умножение матриц}
procedure Umn(m1,m2:matrix;x,y:byte;var m3:matrix);
var i,j,k:byte;
begin
for k:=1 to x do
for j:=1 to x do
 begin
  m3[k,j]:=0;
  for i:=1 to y do
  m3[k,j]:=m3[k,j]+m1[k,i]*m2[i,j];
 end;
end;
{возведение матриц в квадрат}
procedure Umn2(m1:matrix;x,y:byte;var m3:matrix);
var i,j,k:byte;
begin
for k:=1 to x do
for j:=1 to x do
 begin
  m3[k,j]:=0;
  for i:=1 to y do
  m3[k,j]:=m3[k,j]+m1[k,i]*m1[i,j];
 end;
end;

Код

PROCEDURE MatrixT(VAR A : TMatrix; Rows, Cols : Integer);
VAR
  i, j : Integer;
  z    : Integer;
  m    : Integer;
BEGIN
  if Rows>Cols
    then m:=Rows
    else m:=Cols;
  for i:=1 to m do begin
    for j:=Succ(i) to m do begin
      z:=A[i, j];
      A[i, j]:=A[j, i];
      A[j, i]:=z;
    end;
  end;
END;

Код


VAR
  za, zb        : Integer;
  m, n          : Integer;
  A, B, C, D, E : Matrix;
BEGIN
  n:=4;
  m:=5;
  Vvod(A,n,m);
  Vvod(B,n,m);
  WriteLn('Matrix A:');
  Vyvod(A,n,m);
  WriteLn('Matrix B:');
  Vyvod(B,n,m);
  za:=z(A, n, m);
  zb:=z(B, n, m);
  WriteLn('z(A)=', za);
  WriteLn('z(B)=', zb);
  if za>zb
  then begin
    {C:=za*(A*B^T)+zb*(A*B^T)^2}
    {D:=(A*B^T)}
    Trans(B,n,m, C);
    Umn(A,C,m,n, D);
    {C:=za*D}
    scmult(D,n,n,za, C);
    {E:=D^2}
    Umn2(D,n,n, E);
    {D:=zb*E}
    scmult(E,n,n,zb, D);
    {C:=C+D}
    SumMatr(C, D, n,n, C);   {<--- Нужно описать процедуру сложения двух матриц}
    WriteLn('Матрица C');
    Vyvod(C, n,n);
  end
  else begin
.................................
     здесь аналогичный код вычисления матрицы D
.................................
  end;
END.