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

  • «КАЗАНСКИЙ ГОСУДАРСТВЕННЫЙ ЭНЕРЕГЕТИЧЕСКИЙ УНИВЕРСИТЕТ» Кафедра: «Инженерная графика»Практическая работа №1

  • Работу выполнил: Студент группы ЗПТ-1-19Северьянова Л.Н. Руководитель: Гимадиев Р. Ш. Казань 2021 Теоретические положения

  • Рабочее задание. 24 Вариант

  • мми. ММИ1 24В. Практическая работа 1 Решение систем линейных алгебраических уравнений (слау) методом итераций


    Скачать 84.88 Kb.
    НазваниеПрактическая работа 1 Решение систем линейных алгебраических уравнений (слау) методом итераций
    Дата16.01.2022
    Размер84.88 Kb.
    Формат файлаdocx
    Имя файлаММИ1 24В.docx
    ТипПрактическая работа
    #332466


    МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ

    РОССИЙСКОЙ ФЕДЕРАЦИИ
    Государственное образовательное учреждение

    высшего профессионального образования

    «КАЗАНСКИЙ ГОСУДАРСТВЕННЫЙ

    ЭНЕРЕГЕТИЧЕСКИЙ УНИВЕРСИТЕТ»

    Кафедра: «Инженерная графика»


    Практическая работа №1
    «Решение систем линейных алгебраических уравнений (СЛАУ)

    методом итераций»


    Дисциплина: «Методы моделирования и исследования»

    Работу выполнил:

    Студент группы ЗПТ-1-19
    Северьянова Л.Н.

    Руководитель: Гимадиев Р. Ш.

    Казань 2021

    Теоретические положения

    Пусть дана система уравнений


    (1)

    учитывая (1) преобразуется к виду

    , (2)

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

    (3)

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

    Доказательство. Пусть - решение, тогда

    , (4)

    тогда для нормы имеем оценку

    ,

    , поскольку , то

    (5)

    Из (2) однородной системы при вытекает единственное решение, а следовательно, существование и единственное решение систем (2) и (1) при любом .

    Пусть , где -решение системы (2). Вычитая

    (4)-(3) получаем

    (6)

    т.е. , следовательно , где - - ая степень матрицы ,

    - погрешность - ой итерации.

    , (7)

    т.е. норма погрешности не медленнее геометрической прогрессии со знаменателем

    Из выражения (7) определяется число необходимых итераций при решении для достижения заданной точности
    . (8)

    Первая и вторая норма матрицы ( ) вычисляются по формулам

    , . (9)

    Если хотя бы одна из норм меньше единицы, то итерационный процесс сходится.
    Рабочее задание. 24 Вариант Методом простой итерации решить систему линейных уравнений с точностью до 0,001, предварительно оценив число необходимых для этого шагов



    Система уравнений имеет вид



    Пусть - точное решение, - -ое приближение



    - погрешность нулевого приближения. За погрешность принимается норма столбца . А первая норма основной матрицы составляет

    ,

    итерационный процесс сходится.



    k<

    При k> 90 достигается точность и итерации производятся по формуле в таблице 1

    При примем



    При k=1 и далее

    Таблица


    Метод итераций

    k

    x1

    x2

    x3

    x4

    Погрешность

    0

    0,43

    -1,8

    -0,8

    1,7




    1

    -0,2041

    -2,5124

    -0,9131

    1,7574

    0,7124

    2

    -0,2313

    -3,0307

    -1,2088

    1,5717

    0,5183

    3

    -0,2100

    -3,3110

    -1,2522

    1,3344

    0,2804

    4

    -0,1266

    -3,4737

    -1,2663

    1,1975

    0,1627

    5

    -0,0739

    -3,5532

    -1,2519

    1,1273

    0,0796

    6

    -0,0420

    -3,5915

    -1,2425

    1,0990

    0,0382

    7

    -0,0277

    -3,6082

    -1,2363

    1,0884

    0,0167

    8

    -0,0216

    -3,6155

    -1,2337

    1,0850

    0,0073

    9

    -0,0194

    -3,6186

    -1,2327

    1,0839

    0,0031

    10

    -0,0185

    -3,6200

    -1,2324

    1,0835

    0,0014

    11

    -0,0182

    -3,6207

    -1,2323

    1,0833

    0,0006


    Сходимость в тысячных долях имеет место уже при k=11

    Результаты 11 итерации округляем до 0,001

    Ответ:









    -0,018

    -3,621

    -1,232

    1,083


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