Главная страница
Навигация по странице:

  • Методом Зейделя

  • -0,00395 0,153616 -0,1371

  • Методом Гаусса

  • -2,15725 0,29576

  • лабораторная работа 3. Лабораторная работа № 3. Лабораторная работа 3 Решение систем линейных уравнений Цель работы Решение систем линейных уравнений с помощью различных методов


    Скачать 20.41 Kb.
    НазваниеЛабораторная работа 3 Решение систем линейных уравнений Цель работы Решение систем линейных уравнений с помощью различных методов
    Анкорлабораторная работа 3
    Дата19.11.2020
    Размер20.41 Kb.
    Формат файлаdocx
    Имя файлаЛабораторная работа № 3.docx
    ТипЛабораторная работа
    #151903

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

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

    1.Цель работы

    Решение систем линейных уравнений с помощью различных методов.

    2.Задание: вариант 2

    Методами Зейделя и Гаусса решить систему линейных уравнений с точностью ε=10-3.

    (1)

    3.Расчётная часть

    Методом Зейделя

    А) Краткое теоретическое введение

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

    (2)

    Если теперь задать для неизвестных их начальные приближенные значения x(0)n,n=1,2,3 то система (2) позволяет вычислить более точные значения неизвестных на первом шаге (на первой итерации):



    Используя найденные значения неизвестных, можно еще более уточнить их на второй итерации:



    и так далее.

    Общую формулу определения i-го неизвестного на k-й итерации для системы n уравнений можно записать так:



    i=1,2,3,…,n; k=1,2,3,…

    Итерационный процесс продолжается до тех пор, пока все значения xi(k), не станут достаточно близкими к xi(k-1).



    Б) Ручной расчёт





    Первый шаг:



    продолжать расчёт

    Второй шаг:



    продолжать расчёт

    Третий шаг:



    продолжать расчёт

    И так далее

    Следующие шаги приведены в таблице:

    i

    x1

    x2

    x3

    d1

    d2

    d3

    d

    комментарий

    0

    0

    0

    0

     

     

     

     

     

    1

    0,050633

    0,131765

    -0,13605

    0,050633

    0,131765

    0,136047

    0,136047

    продолжать

    2

    0,001803

    0,152483

    -0,13663

    0,04883

    0,020718

    0,000588

    0,04883

    продолжать

    3

    -0,00364

    0,15352

    -0,13708

    0,00544

    0,001037

    0,000448

    0,00544

    продолжать

    4

    -0,00395

    0,153616

    -0,1371

    0,000317

    9,66E-05

    1,53E-05

    0,000317

    закончить

    В) Результат

    ; ;


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

    А) краткое теоретическое введение

    В общем виде система трех линейных уравнений с тремя неизвестными имеет вид:

    (1)

    Предположим, что хотя бы один из коэффициентов при неизвестном x1 отличен от нуля; пусть таким является коэффициент a11.

    Разделим все члены первого уравнения системы (1) на a11, затем все члены второго уравнения – на a21, все члены третьего уравнения – на a31. Если a21 или a31 равен нулю, то соответствующее уравнение входит в новую систему без изменений.

    (2)

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

    (3)

    Число, равное суме всех коэффициентов данного уравнения, включая его свободный член, будем называть строчной суммой. Контрольной суммой будем называть то число, которое на первом этапе совпадает с соответствующей строчной суммой, а на последующих этапах получается из контрольных сумм предыдущего этапа в результате того же преобразования, которое проводится над соответствующим уравнением.

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

    При замене правых частей системы уравнений соответствующими строчными суммами получается система уравнений такая, что если x1, x2, x3 есть решение первоначальной системы, то x1+1, x2+1, x3+1 будет решение новой системы. Числа x’1=x1+1; x’2=x2+1; x’3=x3+1 удовлетворяют первому уравнению новой системы (1).

    Наряду с системами (1), (2), (3), будем иметь в виду еще три системы (1’), (2’), (3’). Левые части параллельных пар систем (1) и (1’), (2) и (2’), (3) и (3’) соответственно совпадают, а правые части систем (1’), (2’), (3’) получаются заменой свободных членов системы (1), (2), (3) соответствующими строчными суммами. Системы (1), (2), (3) эквивалентны между собой; системы (1’), (2’), (3’) – также эквивалентны между собой. Причем система (3’) имеет вид:



    Используя «обратный ход», из системы (3) последовательно находим x3, x2, x1; параллельно из системы (3’) последовательно находим x’3, x’2, x’1. Числа x3 и x’3, x2 и x’2 , x1 и x’1 должны отличаться на 1 (расхождение в запасных знаках во внимание не принимается).



    Прямой ход

    Система

    Строчные суммы

    Контрольные суммы

    Этапы







    I







    II







    III







    IV







    V

    , где

    .

    Обратный ход




    Контроль

    Этапы

    x3=b”3/a”33

    x2=b”2-a”23.x3

    x1=b’1-a’12.x2-a’13.x3

    x’3=S*3/a”3

    x’2=S’’2-a”23.x’3

    x’1=S’1-a’12.x’2-a’13.x’3

    x’3=x3+1

    x’2=x2+1

    x’1=x1+1

    VI

    VII

    VIII


    Б) Ручной расчёт

    коэффициенты трёх неизвестных

    свободные

    члены

    строчные

    суммы

    контрольные

    суммы

    X1

    X2

    X3

    2,7

    14,9

    1,3

    2,1

    21

    21

    1,3

    2,8

    -9,3

    1,7

    -3,5

    -3,5

    15,8

    4,1

    -1,7

    0,8

    19

    19

    1

    5,51852

    0,48148

    0,77778

    7,77778

    7,77778

    1

    2,15385

    -7,15385

    1,30769

    -2,69231

    -2,69231

    1

    0,25949

    -0,10759

    0,05063

    1,20253

    1,20253




    -3,36467

    -7,63533

    0,52991

    -10,47009

    -10,47009




    -5,25902

    -0,58908

    -0,72714

    -6,57525

    -6,57525




    1

    2,26926

    -0,15749

    3,11177

    3,11177




    1

    0,11201

    0,13827

    1,25028

    1,25028







    -2,15725

    0,29576

    -1,86149

    -1,86149





    Решение данной системы

    Решение контрольной системы

    Контроль

    x3= 0,29576/(–2,15725)= –0,1371

    x'3= –1,86149/(–2,15725)= 0,86290

    x'3= x3+1

    x2= –0,15749 – 2,26926∙(–0,1371) = 0,15363

    x'2=3,11177 – 2,26926∙0,86290= 1,15363

    x'2= x2+1

    x1= 0,77778 – 0,48148∙(–0,1371)

    – 5,51852∙0,15363= –0,00398

    x'1= 7,77778 – 0,48148∙1,15363

    – 5,51852∙1,15363= 0,99602

    x'1= x1+1

    В) Результат

    ; ;

    Работу выполнил студент группы 1-31 Баринов Владимир


    написать администратору сайта