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

14 чел.

Выполнил: Марудо А.В., 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


 

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

22151. Ионизационные преобразователи 758 KB
  приемники ионизирующих излучений Область применения ионизационных преобразователей Схемы включения ионизационных преобразователей Погрешности приборов с использованием ионизирующих излучений Заключение Контрольные вопросы Список литературы ОПРЕДЕЛЕНИЕ СУЩНОСТИ ВОПРОСА К ионизационным преобразователям обычно относят большую группу преобразователей в которых измеряемая неэлектрическая величина функционально связана с током ионной проводимости газа возникающим под действием излучения радиоактивных веществ или рентгеновских лучей....
22152. АКУСТИЧЕСКИЕ ВОЛНЫ В ТВЕРДОМ АКУСТИЧЕСКИЕ ВОЛНЫ В ТВЕРДОМ ТЕЛЕ ТЕЛЕ1 1. 487.18 KB
  АКУСТИЧЕСКИЕ ВОЛНЫ В ТВЕРДОМ ТЕЛЕ 1 1. АКУСТИЧЕСКИЕ ВОЛНЫ В ИЗОТРОПНОМ ТВЕРДОМ ТЕЛЕ. АКУСТИЧЕСКИЕ ВОЛНЫ В АНИЗОТРОПНОМ УПРУГОМ ТВЕРДОМ ТЕЛЕ 14 2.
22153. Вопросы по курсу ФОПИ 24 KB
  Приборы для измерения количества тепла.Калориметры для измерения теплотворной способности жидких и газообразных веществ. 8Погрешности измерения тепловой энергии.Методы измерения деформаций и механических напряжений тензометры.
22154. ПРЕОБРАЗОВАТЕЛИ ХОЛЛА 1.8 MB
  Эффект Холла. Параметры и характеристики датчиков Холла. Изготовление и применение датчиков Холла.
22157. ОБЛІК КОРОТКОСТРОКОВИХ ЗОБОВ’ЯЗАНЬ 62.5 KB
  Поняття, оцінка та види короткострокових зобов’язань (пасивів). Облік заборгованості по рахунках постачальників та векселів виданих (векселів до сплати). Облік інших короткострокових зобов’язань. Облік заробітної плати і зобов’язань по заробітній платі.
22158. Резистивные преобразователи перемещения 2.95 MB
  Контактные преобразователи принцип действия и основные типы контактных преобразователей 4 требования к электрической цепи область применения 6 II.2 классификация реостатных преобразователей по конструктивным особенностям...
22159. ТАХОМЕТРЫ 855.42 KB
  Приборы и датчики угловой скорости Приборы предназначенные для измерения частоты вращения называются тахометрами. Тахометры применяются для измерения частоты вращения вала двигателя и его агрегатов. Наибольшее распространение получили следующие методы измерения частоты вращения по принципу действия чувствительного элемента ЧЭ: центробежные в которых ЧЭ реагирует на центробежную силу развиваемую неуравновешенными массами при вращении вала; магнитоиндукционные основанные на зависимости наводимых в металлическом теле вихревых токов от...