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

  • Эксперимент 1.

  • n кол-во итераций eps

  • Эксперимент 4.

  • N

  • Методы решения систем линейных алгебраических уравнений


    Скачать 151.49 Kb.
    НазваниеМетоды решения систем линейных алгебраических уравнений
    Дата22.11.2018
    Размер151.49 Kb.
    Формат файлаdocx
    Имя файлаChislennye_metody_reshenie_SLAU_Bagramov_Robert_09-405_Avtosokhr.docx
    ТипРешение
    #57324
    страница7 из 10
    1   2   3   4   5   6   7   8   9   10



    Более точное решение, которое практически накладывается на точное.

    Другие эксперименты.

    Экспериментальная часть для метода Якоби.


    Эксперимент 1.

    1.1 Зафиксируем значение eps=


    n

    кол-во итераций

    eps

    погрешность

    h*h

    5

    10

    0,008

    0,038836019

    0,04

    6

    14

    0,00462963

    0,040036738

    0,027777778

    7

    21

    0,002915452

    0,033625109

    0,020408163

    8

    29

    0,001953125

    0,029840483

    0,015625

    9

    38

    0,001371742

    0,027583694

    0,012345679

    10

    49

    0,001

    0,024749754

    0,01

    11

    61

    0,000751315

    0,023142206

    0,008264463


    При n≥11 заданная точность не достигается.

    Исходя из того, что метод Якоби чувствителен к значениям вектора начального приближения, проведем этот эксперимент с другим начальным приближением. Что при определенных значениях даст нам уменьшение погрешности при n≥11.

    1.2 Зафиксируем значение eps=


    n

    кол-во итераций

    eps

    погрешность

    h*h

    5

    12

    0,008

    0,009662575

    0,04

    6

    23

    0,00462963

    0,002809168

    0,027777778

    7

    27

    0,002915452

    0,007521441

    0,020408163

    8

    44

    0,001953125

    0,001536698

    0,015625

    9

    48

    0,001371742

    0,006712947

    0,012345679

    10

    73

    0,001

    0,000949578

    0,01

    11

    75

    0,000751315

    0,006375368

    0,008264463
















    15

    153

    0,000296296

    0,004656866

    0,004444444

    20

    345

    0,000125

    0,00045461

    0,0025

    25

    484

    0,000064

    0,002845788

    0,0016

    30

    845

    3,7037E-05

    0,000308867

    0,001111111

    31

    783

    3,35672E-05

    0,00228595

    0,001040583

    32

    972

    3,05176E-05

    0,000277248

    0,000976563
















    33

    900

    2,78265E-05

    0,00214597

    0,000918274


    В данном случае заданная точность не достигается при n≥33.

    Вывод: Получено экспериментальное подтверждение зависимости метода Якоби от вектора начального приближения. Таким образом, мы смогли достичь заданной точности при n до 33.



    Эксперимент 2

    Зафиксируем значение n = 30 . Будем изменять порядок eps.

    n

    кол-во итераций

    eps

    погрешность

    h*h

    30

    275

    0,001111111

    0,071470535

    0,001111111

    30

    845

    3,7037E-05

    0,000308867

    0,001111111

    30

    1413

    1,23457E-06

    0,002486984

    0,001111111

    30

    1981

    4,11523E-08

    0,002565232

    0,001111111


    Мы видим, что при значении eps= и n=30 заданная точность уже не достигается. Исследуем более подробно случай, когда eps=.

    n

    кол-во итераций

    eps

    погрешность

    h*h

    5

    5

    0,04

    0,099535082

    0,04

    10

    31

    0,01

    0,064396689

    0,01

    11

    22

    0,008264463

    0,134577096

    0,008264463


    Мы видим, что при eps= n целесообразно брать до 10 включительно. Далее при увеличении n заданная точность уже не достигается.
    Эксперимент 3

    Рассмотрим, как зависит количество шагов итерации, погрешность от уменьшения порядка eps.

    Будем менять n при значении

    n

    кол-во итераций

    eps

    погрешность

    h*h

    10

    92

    0,0001

    0,009117493

    0,01

    15

    239

    1,97531E-05

    0,005883609

    0,004444444

    20

    470

    0,00000625

    0,004282809

    0,0025

    25

    786

    0,00000256

    0,003367734

    0,0016

    30

    1193

    1,23E-06

    0,002771514

    0,001111111


    Будем менять n при значении

    n

    кол-во итераций

    eps

    погрешность

    h*h

    10

    134

    0,00001

    0,007646425

    0,01

    15

    353

    1,31687E-06

    0,005121005

    0,004444444

    20

    692

    3,125E-07

    0,003853261

    0,0025

    25

    1159

    1,02E-07

    0,003087237

    0,0016

    30

    1761

    4,11523E-08

    0,002574763

    0,001111111


    Вывод: При достаточно больших n и eps= и точность достигается. (Если точность достигается при eps= то при тем более точность достигается при тех же n)
    Эксперимент 4.

    Рассмотрим, как зависит метод от начального приближения.

    Пусть начальный вектор =(0,0, …,0) и eps=. Тогда:

    n

    кол-во итераций

    eps

    погрешность

    h*h

    10

    91

    0,0001

    0,009166465

    0,01

    15

    239

    1,98E-05

    0,00586313

    0,004444444

    20

    469

    0,00000625

    0,004277529

    0,0025

    25

    783

    0,00000256

    0,00336791

    0,0016

    30

    1192

    1,23E-06

    0,002767627

    0,001111111


    Теперь возьмем вектор =(0.05 ,…, 0.05).

    n

    кол-во итераций

    eps

    погрешность

    h*h

    10

    90

    0,0001

    0,009003361

    0,01

    15

    238

    1,97531E-05

    0,005841149

    0,004444444

    20

    468

    0,00000625

    0,004275496

    0,0025

    25

    784

    0,00000256

    0,003364847

    0,0016

    30

    1192

    1,23E-06

    0,002768912

    0,001111111



    По сравнению с предыдущими результатами, при =(0.05 ,…, 0.05) особых изменений не происходит. Количество итераций меняется от -1 до 1 при n от 10 до 30.

    Теперь возьмем вектор = (0.1 ,…, 0.1).


    N

    итер

    eps

    погрешность



    5

    8

    0.008

    0.013

    0.04

    10

    58

    0.001

    0,0026

    0.01

    15

    119

    0.00029

    0.00029

    0.0044




    1   2   3   4   5   6   7   8   9   10


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