Главная страница

красиво. Дана система линейных алгебраических уравнений (слау)


Скачать 154.93 Kb.
НазваниеДана система линейных алгебраических уравнений (слау)
Анкоркрасиво
Дата11.04.2023
Размер154.93 Kb.
Формат файлаdocx
Имя файлаBIE21-01_-_LAB (1) (1).docx
ТипДокументы
#1054993

Дана система линейных алгебраических уравнений (СЛАУ)



Au f , (1)


a11

a12

a13

a14

u1

f1

     

где


A a21 a22

a23 a24 ,

u u2 , f f2 .

a a

a a

u  f

31 32

33 34

3   3

a a

a a

u  f

41 42

43 44

4   4



  1. Решить СЛАУ (1) методом Гаусса с выбором главного элемента по столбцам.

  2. Привести СЛАУ (1) к виду


u Bu F ,

приемлемому для применения метода простой итерации, и проверить выполнение достаточного условия сходимости этого метода, т. е.

B q1, (2)




где



B E  A ,
F f ,

1 0 0




0
E 0 1 0

0 1

0 0 0

0



0
0 , значение  выбрать таким, чтобы




1

выполнялось неравенство (2).



Решить эту СЛАУ методом простой итерации
uk1 Buk F

с точностью 0,01 , предварительно преобразовав ее к виду, приемлемому

для применения этого метода. В качестве начального приближения выбрать

u0 0 , где 0 0, 0, 0, 0 .

T


В качестве условия остановки итерационного процесса (для достижения заданной точности ) использовать неравенство

uk1 uk

1 q , (3)



где q константа из неравенства (2). Имеет место неравенство






εk

uk1 uk ,

где


εk uk U , U вектор точного решения СЛАУ,

uk1 и uk
  • векторы


приближенного решения СЛАУ, полученные на k1 и k итерациях соответственно. Следовательно, если выполняется неравенство


uk1 uk

,


равносильное неравенству (3), то справедливо и неравенство



εk ,

обеспечивающее нахождение вектора uk


заданной точностью .

приближенного решения СЛАУ с



В качестве нормы

uk1 uk

вектора


uk1 uk

взять октаэдрическую



норму, определяемую равенством




4

u
2 uii1

, где u u1, u2 , u3, u4

T ,



а в качестве нормы B матрицы B норму, согласованную с октаэдрической

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



4
B  max

b, где

B b4 .

2 1 j4

ij

i1

ij i, j1













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