Turbo Paskal "Операции над матрицами"

Примечаниеот автора: Курсовой проект на Turbo Pascal, включающий в себя такие операции над матрицами, как перемножение двух матриц, транспонирование, вычисление определителя, вычисление обратной матрицы, сложение (вычитание) двух матриц
Загрузить архив:
Файл: ref-16532.zip (22kb [zip], Скачиваний: 378) скачать

Государственный Комитет Российской Федерации

по Высшему Образованию

Санкт-Петербургский

Государственный Электротехнический Университет «ЛЭТИ»

Кафедра Радиосистем

Пояснительная Записка

к Курсовой Работе по Дисциплине

«Информатика. Основы Программирования»

Вариант №9: «Операции над Матрицами»

Выполнил:студент Чапчаев В.В.

                     факультет: РТ

                     уч.группа: № 2142

Проверил:   ассистент Аникин А.П.

Санкт – Петербург

2 0 0 3

Условие задачи

Программа работает в текстовом режиме. Пользователь выбирает операцию, вводит с клавиатуры размер матрицы и значения элементов матрицы. Программа выполняет заданную операцию и выводит на экран ее результат. Для реализации в программе предлагаются следующие операции: перемножение двух матриц, транспонирование, вычисление определителя, вычисление обратной матрицы, сложение (вычитание) двух матриц.

Метод решения задачи

После запуска программы, пользователю предлагается выбрать требуемую операцию: нахождение определителя, обратной матрицы, транспонированной матрицы, сложение или вычитание двух матриц, умножение двух матриц. Это производится по средствам ввода пользователем цифры, соответствующей данной операции. Далее следует запрос о размерах данной матрицы (кол-ве строк и столбцов), после ввода которых программа проверяет, чтобы введенные значения не были отрицательными или не превышали допустимых размеров матрицы (в программе – 10х10). В случае, если перечисленные условия не выполняются, на экран выводится сообщение об ошибке, и программа завершает работу. При выполнении условий, пользователю предлагается ввести с клавиатуры исходную матрицу «А» (в программе – MAS1). Ввод производится построчно, для перехода к следующей строке, пользователю необходимо нажать клавишу  «ENTER». После того, как число строк становится равным числу, введенному пользователем, ввод матрицы прекращается. Введенная матрица выводится на экран для проверки. Далее следует проверка «пригодности» введенной матрицы к операции, выбранной пользователем, в случае «непригодности» матрицы, на экран выводится сообщение об ошибке, и программа завершает работу. В случае «пригодности» матрицы, над ней производится выбранная операция, после чего на экран выводится сообщение о результате. Если пользователем была выбрана операция, требующая вторую матрицу «В» (в программе – MAS2), пользователю предлагается ввести ее размеры, как и при вводе матрицы «А», далее программа проверяет, возможно ли проведение выбранной операции над матрицами с введенными параметрами. Если операция невозможна, то на экран выводится сообщение об ошибке, и программа завершает работу. Иначе, операция выполняется, и ее результирующая матрица «С» (в программе – MAS3) выводится на экран.

Ввод матриц возможен только с клавиатуры, непосредственно при выполнении программы. Ввод данных «из файла» не предусмотрен, т.к. этот вариант в условии задачи не оговаривается.

Исходные данные и форма их представления

Переменные

Структура

Тип

Назначение

MAS1, MAS2, MAS3

Двумерный массив

Вещественный

Матрицы «А», «В», «С»

i1, i2, i3

Переменные

Целочисленный

Счетчик строк матриц А, В и С соответственно

j1, j2, j3

Переменные

Целочисленный

Счетчик столбцов матриц А, В и С соответственно

m1, m2, m3

Переменные

Целочисленный

Кол-во строк матриц А, В и С соответственно

n1, n2, n3

Переменные

Целочисленный

Кол-во столбцов матриц А, В и С соответственно

operation

Переменные

Целочисленный

Вариант программы

det

Переменные

Вещественный

Определитель

k

Переменные

Целочисленный

Рабочая переменная

Блок – схема программы

Начало

Case operation

Operation = 1

Operation = 2

Operation = 3

Operation = 6

Ввод operation

НЕТ

Ввод m1

Ввод n1

1>n1 or n1>10

or 1>m1 or m1>10

ДА

Ввод MAS1

Вывод MAS1

Подсчет det

Вывод det

Подсчет обр.матр.

Вывод обр.матр.

Транспонир. матрицы

Вывод

транспонир. матрицы

Operation = 4,5

Ввод m2

Ввод n2

n1<>n2 or m1<>m2

ДА

НЕТ

Ввод MAS2

Вывод MAS2

Operation = 4

ДА

k :=1

НЕТ

Operation = 5

ДА

k :=-1

Вывод MAS3

(сумма/разность)

1

2

3

4

5

6

7

Ввод MAS2

Вывод MAS2

m3 := m1

n3 := n2

Вывод MAS3

(произведение)

Ввод m2

Ввод n2

1>=m2 or m2>10

or 1>=n2 or n2>10

n2<>m1

ДА

НЕТ

Текст: «Ошибка»

Конец


Текст Программы:

ProgramMATRIX;{Наименование программы}

Uses CRT;

VAR i1, i2, i3 : integer;{Счетчикстрок}

    j1, j2, j3 : integer;{Счетчикстолбцов}

operation: integer;    {Вариант развития программы}

det :      real;       {Определитель}k:   integer;          {Рабочая переменная}

    {Массивы (матрицы), используемые в программе}

    MAS1,   {Матрица А}              

    MAS2,   {Матрица В}

    MAS3 : array [1..10,1..10] of real;   {МатрицаС}

BEGIN  {Начало программы}

     WriteLn ('Что Вы желаете делать с матрицами?');

     {Выбор пользователем варианта развития программы}

     WriteLn ('Если Вы желаете найти определитель матрицы, нажмите 1');

     WriteLn ('Если Вы желаете найти обратную матрицу, нажмите 2');

    WriteLn ('Если Вы желаете транспонировать матрицу, нажмите 3');

     WriteLn ('Если Вы желаете сложить матрицы, нажмите 4');

     WriteLn ('Если Вы желаете вычесть матрицы, нажмите 5');

     WriteLn ('Если Вы желаете перемножить матрицы, нажмите 6');

    ReadLn (operation); {Занесение выбранного варианта в память}

     WriteLn ('Введите кол-во строк исходной матрицы, не более 10');

     ReadLn (m1);

     WriteLn ('Введите кол-во столбцов исходной матрицы, не более 10');

     ReadLn (n1);

      If ((1>n1) or(n1>10) or (1>m1) or(m1>10)) {Условияошибки}

        then begin

               WriteLn ('ОШИБКА!!!');

               Halt;

             end

        else begin

               WriteLn ('Введите исходную матрицу');{Ввод исходной матрицы}

               for i1:=1 to m1 do

                  for j1:=1 to n1 do Read (MAS1[i1,j1]);

              end;

               for i1:=1 to m1 do   {Выводисходнойматрицы}

                      begin

                       for j1:=1 to n1 do

                        Write (MAS1[i1,j1]);

                        WriteLn;

                      end;

          Case operation of{Оператор выбора «operation»}

              1: begin{Определитель}

                    if (m1<>n1) then writeLn ('ОШИБКА!!!') {Условие ошибки}

                    else

                        begin {Формула определителя}

                         det:=(MAS1[1,1]*MAS1[2,2]*MAS1[3,3]

                         +MAS1[2,1]*MAS1[3,2]*MAS1[1,3]

                         +MAS1[1,2]*MAS1[2,3]*MAS1[3,1])

                         -(MAS1[3,1]*MAS1[2,2]*MAS1[1,3]

                         +MAS1[3,2]*MAS1[2,3]*MAS1[1,1]

                         +MAS1[2,1]*MAS1[1,2]*MAS1[3,3]);

                         WriteLn ('Opredelitel det=',det); {Выводопределителя}

                        end;

                 end;

              2: begin{Обратнаяматрица}

                   if (m1<>n1) then WriteLn ('ОШИБКА!!!') {Условиеошибки}

                   else begin

                          det:=(MAS1[1,1]*MAS1[2,2]*MAS1[3,3] {Определитель}

                          +MAS1[2,1]*MAS1[3,2]*MAS1[1,3]

                          +MAS1[1,2]*MAS1[2,3]*MAS1[3,1])

                          -(MAS1[3,1]*MAS1[2,2]*MAS1[1,3]

                          +MAS1[3,2]*MAS1[2,3]*MAS1[1,1]

                          +MAS1[2,1]*MAS1[1,2]*MAS1[3,3]);

                          if det = 0 then WriteLn ('ОШИБКА!!!') {Условиеошибки}

                          else begin {Союзная матрица}

for i1:=1 to m1 do

                                     for j1:=1 to n1 do MAS2[i1,j1]:=MAS1[j1,i1];

{Итоговаяформула}             

                                for i1:=1 to m1 do

                                     for j1:=1 to n1 do                                                                                                        MAS3[i1,j1]:=(1/det)*MAS2[i1,j1];

WriteLn;

                                 WriteLn ('Обратная матрица:');

                                 for i1:=1 to m1 do begin {Выводобратнойматрицы}

                                                      for j1:=1 to n1 do

                                                      Write (MAS3[i1,j1]);

                                                      WriteLn;

                                                    end;

                               end;

                        end;

                 end;

              3: begin {Транспонирование матрицы}

                   for i1:=1 to m1 do

                    for j1:=1 to n1 do MAS2[i1,j1]:=MAS1[j1,i1]; {Формула}

                   WriteLn ('Транспонированная матрица: ');

                   for i1:=1 to m1 do begin {Выводтранспонированнойматрицы}

                                          for j1:=1 to n1 do

                                          Write (MAS2[i1,j1]);

                                          WriteLn;

                                      end;

                 end;

              4,5: begin {Сложение/вычитание матриц}

                         {Ввод второй матрицы}

                       WriteLn ('Введите кол-во строк второй матрицы');

                       ReadLn (m2);

                       Writeln ('Введите кол-во столбцов второй матрицы');

                       ReadLn (n2);

                       If (n2<>n1) or (m2<>m1)

                        then WriteLn ('OSHIBKA!!!') {Условие ошибки}

                        else begin

WriteLn ('Введитевторуюматрицу');

for i1:=1 to m1 do

                                   for j1:=1 to n1 do

                                    Read (MAS2[i1,j1]);

                             end;

                       for i1:=1 to m1 do   {Вывод второй матрицы}

                           begin

                             for j1:=1 to n1 do

                             Write (MAS2[i1,j1]);

                             WriteLn;

                           end;

           if operation = 4 then k := 1;

                       if operation = 5 then k := -1;

                       for i1:=1 to m1 do

                          for j1:=1 to n1 do

                            MAS3[i1,j1]:=MAS1[i1,j1]+k*MAS2[i1,j1]; {Итоговая формула}

                       writeln('Сумма/разность:');

                       for i1:=1 to m1 do

                           begin

                             for j1:=1 to n1 do Write(MAS3[i1,j1]);

                             WriteLn;

                           end;

                   end;

              6:begin {Умножение матриц}

                    {Ввод второй матрицы}

                    WriteLn ('Введите кол-во строк второй матрицы');

                    ReadLn (m2);

                   Writeln ('Введите кол-во столбцов второй матрицы');

                    ReadLn (n2);

                    If ((1>=m2) or(m2>=10) or (1>=n2) or(n2>=10) {Условиеошибки}

                                or(n2<>m1)) then WriteLn ('ОШИБКА!!!')

                       elsebegin

                              WriteLn ('Введите вторую матрицу');

                              for i2:=1 to m2 do

for j2:=1 to n2 do Read (MAS2[i2,j2]);

                            end;

        for i2:=1 to m2 do begin {Выводвторойматрицы}

                                         for j2:=1 to n2 do

                                         Write (MAS2[i2,j2]);

                                         WriteLn;

                                       end;

                    m3:=m1; n3:=n2;

                    for i3:=1 to m3 do

                        for j3:=1 to n3 do begin

                                             MAS3[i3,j3] := 0;

                                             for i2:= 1 to m2 do

                                  {Итоговая формула}

                                  MAS3[i3,j3] :=MAS3[i3,j3] + MAS1[i3,i2] * MAS2[i2,j3];

                                           end;

                    begin {Вывод произведения}

                      writeln;

                      writeln('Произведение:');

                      for i3:=1 to m1 do begin

                                           for j3:=1 to n2 do Write (MAS3[i3,j3]);

                                           WriteLn;

                                         end;

                    end;

              end;

          End; {End Case}

          ReadKey;

END.   {Конец программы}

Пример работы программы

1.   Операция нахождения определителя:

1)

2)

3)

4)                  æ1  45ö

ç5  37÷

è8  14ø

5)   det = 54.

2.   Операция нахождения обратной матрицы:

1)   При появлении на экране меню с выбором операций, выбираем 2.

2)   На предложение ввести кол-во строк матрицы, выбираем 3.

3)   Кол-во столбцов – 3.

4)   Вводим матрицу:      æ1  45ö

ç5  37÷

è8  14ø

5)   Программа выдает ответ:      æ1,851*100  9,259*1001,481*10    ö

              ç7,407*100  5,555*1001,851*100÷

              è9,259*100  1,296*10    7,407*100ø.

3.Операция транспонирования матрицы:

1)   При появлении на экране меню с выбором операций, выбираем 3.

2)   На предложение ввести кол-во строк матрицы, выбираем 3.

3)   Кол-во столбцов – 3.

4)   Вводим матрицу:      æ1  45ö

ç5  37÷

è8  14ø

5)   Программа выдает ответ:      æ1  58ö

              ç4  31÷

              è5  74ø.

4.Операция сложения матриц:

1)   При появлении на экране меню с выбором операций, выбираем 3.

2)   На предложение ввести кол-во строк матрицы, выбираем 3.

3)   Кол-во столбцов – 3.

4)   Вводим матрицу:      æ1  45ö

ç5  37÷

è8  14ø

5)   На предложение ввести кол-во строк матрицы, выбираем 3.

6)   Кол-во столбцов – 3.

7)   Вводим матрицу:      æ    2  810ö

ç10  614÷

                                          è16  2    8ø.

8)   Программа выдает ответ:      æ    3  1215ö

              ç15    921÷

                                                        è24    312ø.

5.Операция вычитания матриц:

1)   При появлении на экране меню с выбором операций, выбираем 3.

2)   На предложение ввести кол-во строк матрицы, выбираем 3.

3)   Кол-во столбцов – 3.

4)   Вводим матрицу:      æ1  45ö

ç5  37÷

è8  14ø

5)   На предложение ввести кол-во строк матрицы, выбираем 3.

6)   Кол-во столбцов – 3.

7)   Вводим матрицу:      æ    2  810ö

ç10  614÷

                                          è16  2    8ø.

8)   Программа выдает ответ:     æ-1  -4-5ö

             ç-5  -3-7÷

             è-8  -1-4ø.

6.Операция перемножения матриц:

1)   При появлении на экране меню с выбором операций, выбираем 3.

2)   На предложение ввести кол-во строк матрицы, выбираем 3.

3)   Кол-во столбцов – 3.

4)   Вводим матрицу:      æ1  45ö

ç5  37÷

è8  14ø.

5)   На предложение ввести кол-во строк матрицы, выбираем 3.

6)   Кол-во столбцов – 3.

7)   Вводим матрицу:      æ    2  810ö

ç10  614÷

                                          è16  2    8ø.

8)   Программа выдает ответ:     æ122  42106ö

             ç152  72148÷

             è   90   78126ø.

7.Первый исключительный случай:

1)   При появлении на экране меню с выбором операций, выбираем 1.

2)   На предложение ввести кол-во строк матрицы, выбираем 0.

3)   Кол-во столбцов – 3.

4)

8.Второй исключительный случай:

1)   При появлении на экране меню с выбором операций, выбираем 1 (2, 3).

2)   На предложение ввести кол-во строк матрицы, выбираем 3.

3)   Кол-во столбцов – 2.

4)   Вводим матрицу:      æ1  4ö

ç5  3÷

è8  1ø.

5)

9.Третий исключительный случай:

1)   При появлении на экране меню с выбором операций, выбираем 4.

2)   На предложение ввести кол-во строк матрицы, выбираем 3.

3)

4)   Вводим матрицу:      æ1  45ö

ç5  37÷

è8  14ø.

5)   На предложение ввести кол-во строк матрицы, выбираем 2.

6)   Кол-во столбцов – 2.

6)

10.Четвертый исключительный случай:

1)   При появлении на экране меню с выбором операций, выбираем 6.

2)   На предложение ввести кол-во строк матрицы, выбираем 3.

3)   Кол-во столбцов – 3.

4)   Вводим матрицу:     æ1  45ö

ç5  37÷

è8  14ø.

5)   На предложение ввести кол-во строк матрицы, выбираем 2.

6)   Кол-во столбцов – 2.

7)

Выводы по работе

Для проверки результатов, выдаваемых программой, использовался математический пакет MathCAD 2001 Professional.

Существенный недостаток программы в том, что при нахождении обратной матрицы и определителя, программа может работать только с матрицами, размером не превышающие 3х3. Это связано со значительным усложнением программы, в случаи необходимости выполнения вышеперечисленных операций над матрицами большего размера.

При написании программы использовалась следующая литература:

1) Фаронов В.В. «TurboPascal 7.0, начальный курс. Учебное пособие».

2) Мизрохи С.В. «TurboPascal и объектно-ориентированное програмирование».

3) Pascal».

4) Боревич Э.И. «Определители и матрицы».

5) Гутер Р.С., Резниковский П.Т. «Программирование и вычислительная математика».

6)

7) Горстко А.Б., Ерусалимский Я.М. «Математика. Общий курс».