VPF::подпрограммы в паскале - Форум программистов

Витас

Дата 5.3.2010, 08:38 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 5
Регистрация: 5.3.2010

Репутация: нет
Всего: нет

Даны две матрицы А размера n*m B размера m*n, состоящие из целых елементов. Если количество z(A) нулевых элементов матрицы А, больше количества z(B) нулевых элементов матрицы B, то вычислить новую матрицу С по формуле C=z(A)*A*B^T(транспонированная) +z(B)*(A*B^T)^2, иначе вычислить новую матрицу D по формуле D по формуле D=z(B)*((A^T)*B)^2+(A^T)*B*z(A)

ФедосеевПавел

Дата 5.3.2010, 17:26 (ссылка)

(нет голосов)

Загрузка ...

Опытный

Профиль
Группа: Участник
Сообщений: 291
Регистрация: 7.2.2009

Репутация: 2
Всего: 10

Мне кажется, что в условии ошибка...

Матрицы A(n,m), B(m,n). Это "сцепленные" матрицы. Их можно умножать A*B=(n,m)*(m,n)=(n,n) или
B*A=(m,n)*(n,m)=(m,m).

Теперь попробуем вычислить, например C=z(A)*A*B^T+z(B)*(A*B^T)^2
Обозначим B1=B^T. B1(n,m).

А как умножить A*B^T=A*B1=(n,m)*(n,m) ?

И тоже самое для D.

Витас

Дата 9.3.2010, 17:42 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 5
Регистрация: 5.3.2010

Репутация: нет
Всего: нет

вот процедуры. помогите дописать тело программы и описание переменных

Код


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;

Добавлено @ 17:45
и в условии матрицы А и В обе размера (n*m)

Тегами пользуйся для подсветки синтаксиса

Это сообщение отредактировал(а) volvo877 - 9.3.2010, 19:22

ФедосеевПавел

Дата 9.3.2010, 19:52 (ссылка)

(нет голосов)

Загрузка ...

Опытный

Профиль
Группа: Участник
Сообщений: 291
Регистрация: 7.2.2009

Репутация: 2
Всего: 10

Пересмотри подпрограммы scmult и Trans.
Первая даже не откомпилируется (функция должна возвращать только перечислимые типы (Integer, Real и т.д.).
Вторая будет возвращать неверный результат для неквадратных матриц. Могу предложить такой вариант

Код

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.

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Object Pascal: кроссплатформенные технологии \| Следующая тема »