Pascal ABC: Программирование матричных операций

[AnnieBond]

1) Рассчитать элементы квадратной матрицы A = (aij ), i,j = 1,2,...,n по заданной формуле;
2) Вычислить элементы вектора-столбца X = (xi), i = 1,2,...,n по заданному правилу;
3) Вычислить произведение матрицы А на столбец Х;
5) вычисления значения Y по заданной формуле.
Размерность задачи: n=5.

[ByAlex]

В первом примере индексы правильные или нет? Там Aij или Aji?

[AnnieBond]

Все 100% правильно.Я сделала только первый пункт, дальше не знаю

Program Lab05;
Const
     n=5;
Var
   i,j:Integer;
   A:Array [1..n,1..n] of Real;
  Begin
//Вычисление элементов матрицы
   For i:=1 to n do
      For j:=1 to n do
         Begin
             A[j,i]:=(n*j)/((n-j+1)*(n+i));
         End;
//Выводим матрицу
For i:=1 to n do
   Begin
        For j:=1 to n do
            Write(A[j,i]:12:5);
           Writeln;
   End;
   End.

[Stilet]

//Второе
   For i:=1 to n do
      For j:=1 to n do begin s:=0;
        For k:=1 to i do s:=s+a[i,k];
        x[i]:=x[i]+a[i,j]-s;
      end;

[ByAlex]

//Выводим матрицу
For i:=1 to n do
Begin
   For j:=1 to n do
      Write(A[j,i]:12:5);
   Writeln;
End;

Неправильно выводишь элементы матрицы!!! Ты выводишь элементы столбца в строку. Надо write(A[i,j]:4:2,' ')! Да и заполняешь тоже не так, как нужно. Следи за индексами.

[AnnieBond]

ByAlex, я же говорю, здесь строки и столбцы так и должны быть названы, все правильно. Размерность сказал преподаватель.

[BDA]

program Lab05;

const
  n = 5;

var
  i, j: Integer;
  A: Array [1..n, 1..n] of Real;

begin
  //Вычисление элементов матрицы
  for j := 1 to n do
    for i := 1 to n do
    begin
      A[j, i] := (n * j) / ((n - j + 1) * (n + i));
    end;
  //Выводим матрицу
  for j := 1 to n do
  begin
    for i := 1 to n do
      Write(A[j, i]:12:5);
    Writeln;
  end;
end.

При заполнении, циклы должны располагаться так, чтобы программа работала быстрее.
При выводе нужно циклы тоже поменять местами, иначе матрица выводится транспонированной.

Если не верите, что вывод неправильный, то возьмите любой элемент, лежащий не на главной диагонали, например, A23.
Элемент находится во 2 строке, в 3 столбце.
Рассчитываем по формуле значение в этой ячейке: 5*2/((5-2+1)(5+3))=10/(4*8)=0.3125, а у вас другое

Чтобы не плодить сообщения, буду добавлять код в это сообщение.
Оказалось, что подсчет вектора от Stilet'а работает неверно.

program Lab05;

const
  n = 5;

var
  i, j, k: Integer;
  A: Array [1..n, 1..n] of Real;
  X, Z: Array [1..n] of Real;
  s, c, y: real;
  t: boolean;

begin
  //Вычисление элементов матрицы
  for j := 1 to n do
    for i := 1 to n do
    begin
      A[j, i] := (n * j) / ((n - j + 1) * (n + i));
    end;
  //Выводим матрицу
  writeln('Matrix A');
  for j := 1 to n do
  begin
    for i := 1 to n do
      Write(A[j, i]:12:5);
    Writeln;
  end;
  //Рассчитываем вектор Х
  for i := 1 to n do
  begin
    x[i] := 0;
    for j := 1 to n do 
      x[i] := x[i] + a[i, j];     
    s := 0;
    for k := 1 to i do 
      s := s + a[i, k];
    x[i] := x[i] - n * s;
  end;
  //Выводим вектор Х
  writeln('Vector X');
  for i := 1 to n do
    Write(x[i]:12:5);
  writeln;
  //Рассчитываем вектор Z
  for i := 1 to n do
  begin
    z[i] := 0;
    for j := 1 to n do
      z[i] := z[i] + a[i, j] * x[j];
  end;
  //Выводим вектор Z
  writeln('Vector Z=AX');
  for i := 1 to n do
    Write(z[i]:12:5);
  writeln;
  //Рассчитываем c
  c := 0;
  for i := 1 to n do
    c := c + x[i];
  c := c / n;
  //Рассчитываем y
  y := 1;
  t := false;
  for i := 1 to n do
    if x[i] > c then
    begin
      t := true;
      y := y * x[i];
    end;
  //Выводим y
  if t then
    writeln('y=', y:0:5)
  else
    writeln('y can not calculate');   
  readln;
end.

[AnnieBond]

спасибо а как выполнить 3 и 4 пункты?

[AnnieBond]

Спасибо большое!!!

[ByAlex]

Сообщение от AnnieBond

ByAlex, я же говорю, здесь строки и столбцы так и должны быть названы, все правильно. Размерность сказал преподаватель.

Это понятно! Заполнять можно и так, но вывод у вас не правильный!!! У вас результирующая матрица получается транспонированной.