Решение задач линейной алгебры средствами Delphi: Метод Гаусса для решения систем линейных алгебраических уравнений

Захаров Д.С., ИСТ-21д

МИНОБРНАУКИ РОССИИ

Федеральное государственное бюджетное образовательное  учреждение 
высшего профессионального образования 
«Себряковский филиал Волгоградского государственного архитектурно-строительного университета» 
(ФГБОУ ВПО «СФ ВолгГАСУ»)

Кафедра Математических и естественно-научных дисциплин

 

 

Решение задач линейной алгебры средствами Delphi: Метод Гаусса для решения систем линейных алгебраических уравнений

курсовая работа

по специальности 230400.62 «Информационные  системы и технологии (бакалавр)»

 

	Исполнитель: Захаров Дмитрий ИСТ-21д ___________________________ (подпись)
«Прошла защиту»     Оценка ________________________	Руководитель: Отрощенко Д.С., преподаватель   ___________________________  (подпись)

 

 

 

Михайловка

2012

Содержание

 

1. Введение 3

2. Содержательная часть 4

2.1. История 4

2.2. Описание метода 4

2.3. Условие совместности 5

3. Практическая часть. 6

3.1. Постановка задачи 6

3.2.     Методы решения линейных уравнений 6

3.3. Руководство пользователя 7

3.4. Листинг 9

4. Заключение 17

5. Список литературы 18

1. Введение

Для разработки курсового проекта была предоставлена задача: написать программу которая выводила бы решение системы линейных уравнений с тремя переменными методом Гаусса.

Метод Гаусса — классический метод решения системы линейных алгебраических уравнений (СЛАУ). Это метод последовательного исключения переменных, когда с помощью элементарных преобразований система уравнений приводится к равносильной системе ступенчатого (или треугольного) вида, из которой последовательно, начиная с последних (по номеру) переменных, находятся все остальные переменные.

Данная  программа писалась на языке Delphi в бесплатной среде Lazarus. Lazarus — свободная среда разработки программного обеспечения для компилятора Free Pascal (часто используется сокращение FPC— свободно распространяемый компилятор языка программирования Pascal) на языке Object Pascal. Интегрированная среда разработки предоставляет возможность кроссплатформенной разработки приложений в Delphi-подобном окружении. На данный момент является единственным инструментом быстрой разработки приложений (RAD), позволяющим Delphi-программистам создавать приложения с графическим интерфейсом для Linux (и других не-Windows) систем.

Позволяет достаточно несложно переносить Delphi-программы  с графическим интерфейсом в  различные операционные системы: Linux, FreeBSD, Mac OS X, Microsoft Windows, Android. Начиная с Delphi XE2 в самом Delphi имеется возможность компиляции программ для Mac OS X и iOS.

2. Содержательная часть

1. История

Хотя в  настоящее время данный метод повсеместно называется методом Гаусса, он был известен и до К. Ф. Гаусса. Первое известное описание данного метода — в китайском трактате «Математика в девяти книгах», составленном между I в. до н.э. и II в. н. э.

2. Описание  метода

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

Матрица называется основной матрицей системы,  — столбцом свободных членов.

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

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

Тогда переменные называются главными переменными. Все остальные называются свободными.

Если  хотя бы одно число  , где , то рассматриваемая система несовместна, т.е. у неё нет ни одного решения.

Пусть для любых .

Перенесём свободные переменные за знаки равенств и поделим каждое из уравнений системы на свой коэффициент  при самом левом        (    , где  — номер строки):

, 
где

Если  свободным переменным системы (2) придавать  все возможные значения и решать новую систему относительно главных  неизвестных снизу вверх (то есть от нижнего уравнения к верхнему), то мы получим все решения этой СЛАУ. Так как эта система получена путём элементарных преобразований над исходной системой (1), то по теореме об эквивалентности при элементарных преобразованиях системы (1) и (2) эквивалентны, то есть множества их решений совпадают.

2. Условие совместности

Упомянутое  выше условие  для всех может быть сформулировано в качестве необходимого и достаточного условия совместности:

Теорема Кронекера  — Капелли.

Система совместна тогда и только тогда, когда ранг её основной матрицы равен рангу её расширенной матрицы.

3. Практическая часть.

1. Постановка  задачи

Дана  система из n линейных уравнений  вида:

a₁₁x₁+a₁₂x₂+…+a_1nx_n=b₁      (3) 
a₂₁x₁+a₂₂x₂+…+a_2nx_n=b₂ 
a_n1x₁+a_n2x₂+…+a_nnx_n=b_n. 
Требуется получить решение системы линейных уравнения (3) при любых заданных коэффициентах a_ij  и b_i.

3.2. Методы решения линейных уравнений

Системы линейных уравнений вида (3) решаются точными и итерационными методами. Точные методы дают точное решение за конечное число операций, если все они выполнялись без погрешности. Число операций итерационных методов зависит от заданной погрешности вычислений. Метод Гаусса является точным методом решения системы линейных уравнений.

Метод Гаусса или метод последовательного исключения неизвестных основан на приведении матрицы коэффициентов А=(a_ij) к треугольному виду. При этом алгоритм решения  системы линейных уравнений следующий:

1. С помощью двух циклов с управляющими переменными  i=1,2,..,n  и  j=1,2,…,n  организуется  ввод  коэффициентов a_ij и b_i  образующих массивы A=(a_ij), B=(b_i).

2.Делается прямой ход исключения переменных путем преобразования коэффициентов по формулам:

a_ij =-a_ij/a_ii; a_kj=a_kj+a_ija_ik; b_i=b_i+a_ijb_i

где i=1,2,…,n-1; j=i+1,i+2,…,n; k=i+1,i+2,…,n.

В конце этих преобразований получается x_n=b_n/a_nn.

3.Организуется обратный ход (последовательное нахождение x_n-1, x_n-2 , …, x₂, x₁), проводя вычисления по формулам:

h=b_i, h=h-x_ia_ij, x_i=h/a_ii, где i=n-1, n-2,…, 2, 1; j=i+1, i+2,…,n;

В результате формируется массив X=(x_i) неизвестных x_n, x_n-1,…, x₂, x₁.

Особенностью  метода Гаусса является то, что при прямом ходе производится выбор элемента строки отличного от нуля, если элемент a_ii, стоящий на главной диагонали матрицы А, равен 0.

Этот  выбор осуществляется путем перестановки соответствующих строк матрицы  А. Это исключает деление на 0, если матрица A=(a_ij) содержит нулевые элементы на главной диагонали.

2. Руководство пользователя

Пользоваться  данной программой довольно просто: начинать работу нужно с ввода коэффициентов при переменных x₁,x₂,x₃ (для удобства в программе были изменены на переменные X, Y, Z).

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

Однако, имеется ряд правил занесения  коэффициентов в систему: в первое уравнение надлежит ставить уравнение, в котором есть первый коэффициент  при Х. В обратном случае программа  выдаст ошибку:

Затем при  нажатии на кнопку «Ответ» программа  выдаст результат:

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

4.  Листинг

unit Unit1;

 

{$mode objfpc}{$H+}

 

interface

 

uses

  Classes, SysUtils, FileUtil, Forms, Controls, Graphics, Dialogs, Grids,

  StdCtrls;

 

type

 

  { TForm1 }

 

  TForm1 = class(TForm)

    Button1: TButton;

    Button2: TButton;

    Button3: TButton;

    Label16: TLabel;

    Label17: TLabel;

    Label18: TLabel;

    Label3: TLabel;

    Label4: TLabel;

    Label8: TLabel;

    Label9: TLabel;

    Scobka3: TLabel;

    Scobka4: TLabel;

    Scobka5: TLabel;

    Uravnenie2: TLabel;

    Label1: TLabel;

    Label2: TLabel;

    Uravnenie3: TLabel;

    Uravnenie4: TLabel;

    X1: TEdit;

    X2: TEdit;

    X3: TEdit;

    O1: TEdit;

    Y1: TEdit;

    Y2: TEdit;

    O2: TEdit;

    Y3: TEdit;

    Z1: TEdit;

    Z2: TEdit;

    Z3: TEdit;

    O3: TEdit;

    Uravnenie1: TLabel;

    Reshenie: TLabel;

    Main: TLabel;

    Label5: TLabel;

    Label6: TLabel;

    Label7: TLabel;

    Label10: TLabel;

    Label11: TLabel;

    Label12: TLabel;

    Label13: TLabel;

    Label14: TLabel;

    Label15: TLabel;

    LabelY: TLabel;

    LabelX: TLabel;

    LabelZ: TLabel;

    Scobka1: TLabel;

    OtvetX: TLabel;

    OtvetY: TLabel;

    OtvetZ: TLabel;

    procedure Button1Click(Sender: TObject);

    procedure Button2Click(Sender: TObject);

    procedure Button3Click(Sender: TObject);

 

  private

    { private declarations }

  public

    { public declarations }

  end;

 

var

   px1:integer; //Переменные вводимые в ячейки

   px2:integer;

   px3:integer;

   py1:integer;

   py2:integer;

   py3:integer;

   pz1:integer;

   pz2:integer;

   pz3:integer;

   po1:integer;

   po2:integer;

   po3:integer;

   ipx2:real; //Измененные переменные(1 шаг решения системы)

   ipy2:real;

   ipz2:real;

   ipo2:real;

   ipx3:real;

   ipy3:real;

   ipz3:real;

   ipo3:real;

   kipy3:real;//Конечно измененные переменные(2 шаг решения системы)

   kipz3:real;

   kipo3:real;

   RX:real; //Окончательное решение системы

   RY:real;

   RZ:real;

   VX1:string;//Переменные для "удобного" и наглядного вывода на экран

   VY1:string;

   VZ1:string;

   VX2:string;

   VY2:string;

   VZ2:string;

   VX3:string;

   VY3:string;

   VZ3:string;

  Form1: TForm1;

 

implementation

 

{$R *.lfm}

 

{ TForm1 }

 

procedure TForm1.Button2Click(Sender: TObject);

begin

   Close;

end;

 

procedure TForm1.Button3Click(Sender: TObject);

begin

  px1:=StrtoInt(x1.text);

  px2:=StrtoInt(x2.text);

  px3:=StrtoInt(x3.text);

  py1:=StrtoInt(y1.text);

  py2:=StrtoInt(y2.text);

  py3:=StrtoInt(y3.text);

  pz1:=StrtoInt(z1.text);

  pz2:=StrtoInt(z2.text);

  pz3:=StrtoInt(z3.text);

  po1:=StrtoInt(o1.text);

  po2:=StrtoInt(o2.text);

  po3:=StrtoInt(o3.text);

  if px1=0 then showmessage ('Поставьте на первое место уравнение в котором есть коэффициент при X') else

  begin

  ipx2:=px1*(-px2/px1)+px2;

  ipy2:=py1*(-px2/px1)+py2;