65833

Решение систем линейных алгебраических уравнений

Лабораторная работа

Информатика, кибернетика и программирование

Задание 1 Условие: Решить СЛАУ методом итераций Якоби либо Зейделя: Ход работы: Итерационный метод Якоби: Достаточное условие сходимости: Чтобы выполнялось достаточное условие сходимости преобразуем систему: Блоксхема метода Якоби: Код программы: for i:=1 to 4 do x[i]:=1; e:=0.001; k:=0; while k=0 do begin dx:=0; for i:=1 to 4 do begin s:=0; for j:=1 to i1 do...

Русский

2014-08-09

741 KB

13 чел.

Выполнил: Марудо А.В., 2 курс, 3 группа

Проверил: Шапочкина Ирина Викторовна

Лабораторная работа #1(вариант #5)

Тема: Решение систем линейных алгебраических уравнений.

Цель: Найти решение СЛАУ, т.е. такую упорядоченную совокупность чисел, которая обращает все уравнения системы в верное равенство, пользуясь прямыми (точными) и итерационными методами решения систем.

Задание 1

Условие: Решить СЛАУ методом итераций (Якоби либо Зейделя):

Ход работы: 

Итерационный метод Якоби:

Достаточное условие сходимости:

Чтобы выполнялось достаточное условие сходимости, преобразуем систему:

Блок-схема метода Якоби:

Код программы:

   for i:=1 to 4 do

     x[i]:=1;

     e:=0.001;

     k:=0;

      while k=0 do

       begin

        dx:=0;

        for i:=1 to 4 do

         begin

          s:=0;

          for j:=1 to i-1 do

           s:=s+a[i,j]*x[j];

            for j:=i+1 to 4 do

             s:=s+a[i,j]*x[j];

             y:=(b[i]-s)/a[i,i];

             d:=abs(y-x[i]);

             if d>dx then

              dx:=d;

              x[i]:=y;

         end;

            if dx<e then k:=1;

       end;

Вектор невязки

for t:=1 to 4 do

  begin

   for q:=1 to 4 do

    v[t]:=v[t]+x[q]*a[t,q];

  end;

 Label12.Caption:=Floattostr(v[1]+2);

 Label13.Caption:=Floattostr(v[2]+32);

 Label14.Caption:=Floattostr(v[3]-26);

 Label15.Caption:=Floattostr(v[4]+26);

Полученные результаты:

Задание 2

Условие: Решить СЛАУ Методом Гаусса с постолбцовым выбором главного элемента:

Ход работы:

Метод Гаусса:

Идея метода: привести матрицу к нижнему треугольному виду последовательным выполнением арифметических операций построчно.

Преобразуем систему, чтобы избавиться от нулевых и близких к нулю диагональных элементов, для проведения прямого хода:

Блок-схема метода Гаусса:

Код программы:

Прямой ход

  n:=4;

  for i:=1 to n do

    begin

     l:=i;

     for z:=i+1 to n do

       begin

        if (Abs(a[z,i])<(Abs(a[l,i]))) then

        l:=z;

        if l<>i then

         begin

          for j:=i to n do

           begin

            q:=a[i,j];

            a[i,j]:=a[l,j];

            a[l,j]:=q;

           end;

          w:=b[i];

          b[i]:=b[l];

          b[l]:=w;

         end;

       end;

     for k:=i+1 to n do

       begin

        c:=a[k,i]/a[i,i];

        for j:=i+1 to n do

         begin

          a[k,j]:=a[k,j]-c*a[i,j];

         end;

        b[k]:=b[k]-c*b[i];

       end;

    end;

Обратный ход

for i:=n downto 1 do

 begin

  s:=0;

  for j:=i+1 to n do

   begin

    s:=s+a[i,j]*x[j];

   end;

  x[i]:=(b[i]-s)/a[i,i];

 end;

Вевтор невязки

for t:=1 to 4 do

  begin

   for g:=1 to 4 do

    v[t]:=v[t]+x[g]*u[t,g];

  end;

 Label12.Caption:=Floattostr(v[1]-125);

 Label13.Caption:=Floattostr(v[2]-23);

 Label14.Caption:=Floattostr(v[3]-58);

 Label15.Caption:=Floattostr(v[4]+90);

Полученные результаты:

Вывод:

В данной работе были написаны программы для нахождения решения СЛАУ прямым (Гаусса) и итерационным (Якоби) методами. Исходные системы уравнений имеют 4-й порядок. Это позволило нам использовать оба метода (для прямых методов m<=100).

Для каждой из систем была найдена своя совокупность чисел-решений. Если подставить в исходные системы найденные переменные, уравнения системы обращаются в верные тождества.

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

БГУ

Физический факультет

2011/2012 учебный год

Минск

PAGE   \* MERGEFORMAT 1


 

А также другие работы, которые могут Вас заинтересовать

41148. Лексические и фразеологические нормы русского литературного языка 142 KB
  Лексическое значение слова. Многозначные слова омонимы синонимы антонимы паронимы.Многозначные слова нужно отличать от омонимов. Если связь между значениями слов утрачена, значит это омонимы. Омонимы – слова одинаковые по звучанию, но разные по значению (жать руку- жать серпом)
41149. Способы задания плоскости на эпюре 592.5 KB
  Способы задания плоскости на эпюре Из курса элементарной геометрии известно что через три точки не лежащие на одной прямой можно провести плоскость и при том только одну. Таким образом положение плоскости в пространстве логично определить задать тремя точками точки А В С табл. Кроме этого положение плоскости в пространстве определяют: прямая АВ и точка С не лежащая на прямой табл.
41150. Изобарная и изохорная теплоемкости 80 KB
  3 где индекс х обозначает условия протекания процесса подвода теплоты индекс указывает при каких условиях подводится теплота: при постоянном давлении или при постоянном объеме.2 то можно предположить что поскольку при постоянном объеме внешняя теплота равна изменению внутренней энергии а при постоянном давлении изменению теплосодержания рабочего тела то в общем случае это должны быть разные количества теплоты. Другими словами количество теплоты для увеличения температуры на 1 градус при постоянном давлении будет отличаться от...
41151. Преобразование треугольника в эквивалентную звезду. Метод контурных токов 134 KB
  Преобразование треугольника в эквивалентную звезду. Преобразованием треугольника в эквивалентную звезду называется такая замена части цепи соединенной по схеме треугольником цепью соединенной по схеме звезды при которой токи и напряжения в остальной части цепи сохраняются неизменными. под эквивалентностью треугольника и звезды понимается то что при одинаковых напряжениях между одноименными зажимами токи входящие в одноименные выводы одинаковы. Преобразование треугольника в звезду.
41152. Программные средства обеспечения безопасности передачи данных в компьютерных сетях 320.5 KB
  Введение Интенсивное развитие глобальных компьютерных сетей появление новых технологий поиска информации привлекают все больше внимания к сети Internet со стороны частных лиц и различных организаций. Развитие глобальных сетей привело к многократному увеличению количества не только пользователей но и атак на компьютеры подключенные к сети Internet. При подключении к Internet локальной или корпоративной сети необходимо позаботиться об обеспечении информационной безопасности в этой сети. В сфере компьютерных сетей межсетевой экран представляет...
41153. Применение Первого начала термодинамики для процессов идеального газа 218 KB
  Т =const Если Т =const то после подстановки в уравнение состояния 4.4 получим pV = const или pv = const 4. В соотвествии с 1 внутрення энергия пропорциональна температуре U  T или dU  dT другими словами dU = CVdT где CV коэффициент пропорциональности CV теплоемкость газа при постоянном объеме Дж кг К Если Т =const то dT =0 и U = 0.7: pv = const следовательно pv=...
41154. Формирование уставного капитала и имущества организации, их состав и назначение. Определение потребности в основных и оборотных средствах, показатели использования 135.5 KB
  Определение потребности в основных и оборотных средствах показатели использования.Основные средства предприятия. Состав основных средств структура и воспроизводственная характеристика. Классификация основных средств основных фондов и ее экономическое значение.
41155. ТОВАР В СИСТЕМЕ МАРКЕТИНГА 1.06 MB
  Основные журналы по изучению товара в системе маркетинга Содержание товара как инструмента комплекса маркетинга. Услуги как особый вид товара в системе маркетинга. Концепция товарного ассортимента решения относительно товара продуктовой линии и упаковки.
41156. УЧЕТ ДЕНЕЖНЫХ СРЕДСТВ И ИХ ЭКВИВАЛЕНТОВ 138 KB
  Управление денежными средствами становится все более важным изза огромной сложности финансовых рынков. Правильное раскрытие и классификация денежных средств и их эквивалентов необходимы для точной оценки ликвидности компании. В МСБУ 7 Отчеты о движении денежных средств даны следующие определения денежным средствам их эквивалентам и потокам денежных средств: Денежные средства включают наличные деньги и вклады до востребования. Эквивалент денежных средств краткосрочные высоколиквидные вложения...