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

11 чел.

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


 

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

2731. Субъекты таможенного права 250 KB
  Субъектом права в любой отрасли права принято считать носителя определенных прав и обязанностей, которым он наделен государством для реализации своих жизненных потребностей либо возложенных на него полномочий в соответствующей сфере общественных отношений.
2732. Определение коэффициентов трения с помощью наклонного маятника 135 KB
  Определение коэффициентов трения с помощью наклонного маятника Цель работы: ознакомиться со сложным механическим движением; определить коэффициенты трения различных пар материалов. МЕТОДИКА ЭКСПЕРИМЕНТА На шарик, выведенный из положения равновесия, ...
2733. Определение скорости полета тела с помощью баллистического крутильного маятника 261.5 KB
  Определение скорости полета тела с помощью баллистического крутильного маятника Цель работы: Научиться экспериментально определять скорости движения тел. МЕТОДИКА ЭКСПЕРИМЕНТА В основе экспериментального определения скорости полета тела с помощ...
2734. Моделирование электростатического поля 101.5 KB
  Моделирование электростатического поля Приборы и принадлежности: электролитическая кювета, набор электродов, понижающий трансформатор, цифровой вольтметр В7-38. Введение. Нередко различные физические поля описываются одинаковыми уравнениями. Решение...
2735. Определение удельного электрического сопротивления проволоки 111 KB
  Определение удельного электрического сопротивления проволоки Приборы: лабораторная установка ФПМ-01, мост постоянного тока Р-333. Цель работы: приобретение навыков проведения простейших измерений электрических величин, практическое применение законо...
2736. Реостат и делитель напряжения 158.5 KB
  Реостат и делитель напряжения Приборы и принадлежности: источник тока, два вольтметра, два миллиамперметра, реостат, нагрузочные резисторы. Введение. Реостат – устройство для регулирования тока или напряжения в электрических цепях путем изменен...
2737. Измерение электрического сопротивления одинарным мостом постоянного тока 148 KB
  Измерение электрического сопротивления одинарным мостом постоянного тока (мостом Уитстона) Приборы и принадлежности: реохорд, магазин сопротивлений, источник постоянного тока, гальванометр, два резистора с неизвестным сопротивлением. Введение. Для и...
2738. Гальванометр и его применение 206.5 KB
  Гальванометр и его применение Приборы и принадлежности: источник питания постоянного тока Б5-70, ампервольтметр М2018, реостат, лабораторная панель с гальванометром ЛМ, вольтметром М252, двумя магазинами сопротивлений. Введение. Гальванометрами назы...
2739. Компенсационный метод измерения электрических величин 83 KB
  Компенсационный метод измерения электрических величин Приборы и принадлежности. исследуемый гальванический элемент, нормальный элемент, источник рабочего тока, реохорд, гальванометр, двухполюсный переключатель, магазин сопротивлений...