Программа для вычисления следа матриц (с промежуточными вычислениями)

@motomaniak · Регистрация: 03.06.2013

Author24 — интернет-сервис помощи студентам

Приветствую всех!

Пишу небольшую курсовую, никак не могу разобраться в ряде мелочей. Из меня не б-г весть какой программист, поэтому прошу сразу камнями не закидывать, наверняка в приведённом коде есть глупые и грубые ошибки (оттого-то он и не работает).

В частности, выдаётся ошибка вида "Can't evaluate constant expression" в каждой процедуре. Если с английского на русский я это могу перевести, то с русского на человеческий - уже нет. Какая constant, если оно var?

Здесь нет параметров процедур, поскольку я не очень понимаю, надо ли их описывать вообще, ведь в каждой процедуре и так присутствует описание типов переменных. Сразу скажу, что с настолько серьёзным (для меня серьёзным) уровнем программирования сталкиваюсь впервые, а чтение манов в Сети просветления как-то не дало. Возможно, мозг просто плавится на питерской жаре

Подскажите ещё по этому вопросу, если не сложно.

Итак, задание:
1. Сформировать произвольную квадратную матрицу A;
2. Вычесть из матрицы A единичную матрицу;
3. Сформировать диагональную матрицу C из элементов, которые определяются как суммы элементов столбцов полученной матрицы;
4. Найти след матрицы C.

Обязательное требование - максимальное использование подпрограмм (поэтому сделано так, как сделано... я бы их количество урезал, но раз такое требование, лучше перебдеть...), ввод данных из файла и вывод результатов в файл. До работы с файлами пока не дошёл: это, на мой взгляд, последнее, что надо дорабатывать в программе.

И вот как я это попытался реализовать с высоты своих скромных познаний:

Pascal

program kursovaja;
 
uses crt;
 
var x,i,j:integer;
    s:real;
    M:array[1..x,1..x] of real;
    A:array[1..x,1..x] of integer;
 
(*Ввод размерности матриц*)
procedure EnterSize;
var x:integer;
begin
  repeat
  writeln('Введите размерность квадратной матрицы n=',x);
  readln(x);
  if x<1 then writeln('Ошибка ввода! Размерность квадратной матрицы может быть только целым числом, большим или равным 1. Повторите ввод');
  until x>=1;
end;
 
(*Ввод матрицы с клавиатуры*)
procedure EnterMatrix;
var x,i,j:integer;
    M:array[1..x,1..x] of real;
begin
writeln('Введите элементы квадратной матрицы');
for i:=1 to x do
for j:=1 to x do begin
                 writeln('M[',i,',',j,']=');
                 readln(M[i,j]);
                 end;
end;
 
(*Вывод матрицы на экран в "человекопонятном" виде*)
procedure ShowMatrix;
var x,i,j:integer;
    M:array[1..x,1..x] of real;
begin
for i:=1 to x do begin
                 for j:=1 to x do
                 write(M[i,j]:6:2);
                 writeln;
                 end;
end;
 
(*Вычитание из матрицы единичной матрицы*)
procedure SubtractMatrix;
var x,i:integer;
    M:array[1..x,1..x] of real;
    A:array[1..x,1..x] of integer;
begin
for i:=1 to x do
              A[i,i]:=1;
for i:=1 to x do begin
for j:=1 to x do
              M[i,i]:=M[i,i]-A[i,i];
end;
end;
 
(*Формирование диагональной матрицы из сумм элементов столбцов исх. матрицы*)
procedure FormMatrix;
var x,i,j:integer;
    M:array[1..x,1..x] of real;
begin
for i:=1 to x do
for j:=1 to x do begin
                 M[i,i]:=M[i,i]+M[i,j];
                 end;
end;
 
(*Нахождение следа полученной матрицы*)
procedure FindTrace;
var x,i:integer;
    M:array[1..x,1..x] of real;
begin
for i:=1 to x do begin
                 s:=0;
                 s:=s+M[i,i];
end;
writeln('След матрицы Tr(M)= ',s);
end;
 
 
(*Основной модуль программы*)
 
begin
 
clrscr;
 
EnterSize;
EnterMatrix;
ShowMatrix;
SubtractMatrix;
FormMatrix;
FindTrace;
 
readkey;
end.

Заранее весьма признателен за любую помощь. Работаю исключительно сам, пытаюсь дойти своим умом, но пока не вполне выходит.

@motomaniak · 25.12.2013, 14:42 **[ТС]**

Несмотря на то, что за говнокод в первом сообщении мне безумно стыдно, проблему я решил (давненько уже) и решил выложить рабочую программу, которая у меня сохранилась.

Может, пригодится кому.

Вместо byte лучше использовать integer, это я так извращался, другим не рекомендую.

Pascal

program kursov10;
 
uses crt;
 
const maxcount=255;
 
type matrix=array[1..maxcount,1..maxcount] of real;
 
var i,j,n:byte;
    A:matrix;
 
{Ввод размерности матрицы}
procedure MatrixSize (var x:byte);
begin
repeat
  writeln('Введите размерность матрицы N=');
  readln(x);
  if (x<1) then writeln('Ошибка! Размерность матрицы должна выражаться натуральным числом. Повторите ввод.');
until (x>=1);
end;
 
{Ввод матрицы с клавиатуры}
procedure EnterMatrix (x:byte; var M:matrix);
var i,j:byte;
begin
  writeln('Введите элементы квадратной матрицы:');
  for i:=1 to x do
  for j:=1 to x do begin
  writeln('M[',i,',',j,']=');
  readln(M[i,j]);
  end;
end;
 
{Вывод матрицы на экран}
procedure PrintMatrix (x:byte; var M:matrix);
var i,j:byte;
begin
  for i:=1 to x do begin
  for j:=1 to x do
  writeln('M[',i,',',j,']=',M[i,j]);
  end;
end;
 
{Вычитание единичной матрицы из исходной}
procedure SubtractMatrix (x:byte; var M:matrix);
var i,j:byte;
begin
  for i:=1 to x do
  for j:=1 to x do begin
                   if i=j then M[i,j]:=(M[i,j]-1);
                   end;
end;
 
{Формирование диагональной матрицы C по заданию}
procedure FormNewMatrix (x:byte; var M:matrix);
var i,j:byte;
    s:real;
begin
  for j:=1 to x do begin
                   s:=0;
  for i:=1 to x do begin
                   s:=s+M[i,j];
                   end;
  for i:=1 to x do begin
                   if j=i then M[i,j]:=S else M[i,j]:=0;
                   end;
                   end;
end;
 
{Нахождение следа матрицы}
procedure TraceMatrix (x:byte; var M:matrix);
var i,j:byte;
    tr:real;
begin
  tr:=0;
  for i:=1 to x do
  for j:=1 to x do begin
                   if i=j then Tr:=Tr+M[i,j];
                   end;
  writeln ('След матрицы Tr=',Tr);
end;
 
 
begin
clrscr;
 
MatrixSize(n);
EnterMatrix(n,A);
PrintMatrix(n,A);
SubtractMatrix(n,A);
writeln('После вычитания единичной матрицы исходная матрица А принимает следующий вид:');
PrintMatrix(n,A);
FormNewMatrix(n,A);
writeln('Так выглядит новая матрица С:');
PrintMatrix(n,A);
TraceMatrix(n,A);
 
readkey;
end.

	25.12.2013, 14:42

Опции темы