Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Вычислительная математика.pdf
Скачиваний:
3
Добавлен:
24.11.2025
Размер:
3.89 Mб
Скачать

ЭЛЕМЕНТЫ ЧИСЛЕННЫХ МЕТОДОВ

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

1.1 Дана система линейных алгебраических уравнений

a11x1 + a12 x2 + a13x3 = b1,a21x1 + a22 x2 + a23x3 = b2 ,a31x1 + a32 x2 + a33x3 = b3.

Составьте программу, которая реализует алгоритм одного из прямых методов для решения системы линейных алгебраических уравнений порядка n и вычисляет одновременно обратную матрицу для матрицы системы. Примените составленную программу к данной системе.

Номер

i

ai1

ai2

ai3

bi

варианта

1

1

0,21

– 0,45

– 0,20

1,91

 

2

0,30

0,25

0,43

0,32

 

3

0,60

– 0,35

– 0,25

1,83

2

1

– 3

0,5

0,5

– 56,5

 

2

0,5

– 6,0

0,5

– 100

 

3

0,5

0,5

– 3

– 210

3

1

0,45

– 0,94

– 0,15

– 0,15

 

2

– 0,01

0,34

0,06

0,31

 

3

– 0,35

0,05

0,63

0,37

4

1

0,63

0,05

0,15

0,34

 

2

0,15

0,10

0,71

0,42

 

3

0,03

0,34

0,10

0,32

5

1

– 0,2

1,60

– 0,10

0,30

 

2

– 0,30

0,10

– 1,50

0,40

 

3

1,20

– 0,20

0,30

– 0,60

6

1

0,30

1,20

– 0,20

– 0,60

 

2

– 0,10

– 0,20

1,60

0,30

 

3

0,05

0,34

0,10

0,32

7

1

0,20

0,44

0,81

0,74

 

2

0,58

– 0,29

0,05

0,02

 

3

0,05

0,34

0,10

0,32

8

1

6,36

11,75

10

– 41,4

 

2

7,42

19,03

11,75

– 49,49

 

3

5,77

7,48

6,36

– 27,67

9

1

– 9,11

1,02

– 0,73

– 1,25

 

2

7,61

6,25

– 2,32

2,33

 

3

– 4,64

1,13

– 8,88

– 3,75

10

1

– 9,11

– 1,06

– 0,67

– 1,56

 

2

7,61

6,35

– 2,42

2,33

 

3

– 4,64

1,23

– 8,88

– 3,57

11

1

1,02

– 0,73

– 9,11

– 1,25

 

2

6,25

– 2,32

7,62

2,33

 

3

1,13

– 8,88

4,64

– 3,75

44

12

1

0,06

0,92

0,03

– 0,82

 

2

0,99

0,01

0,07

0,66

 

3

1,01

0,02

0,99

– 0,98

13

1

0,10

– 0,07

– 0,96

– 2,04

 

2

0,04

– 0,99

– 0,85

– 3,73

 

3

0,91

1,04

0,19

– 1,67

14

1

0,62

0,81

0,77

– 8,18

 

2

0,03

– 1,11

– 1,08

0,08

 

3

0,97

0,02

– 1,08

0,06

15

1

0,63

– 0,37

1,76

–9,29

 

2

0,90

0,99

0,05

0,12

 

3

0,13

– 0,95

0,69

0,69

1.2 Составьте программу, которая реализует алгоритм метода простой итерации и метода Зейделя решения линейной системы порядка n. Примените составленную программу к данной системе и сравните полученные результаты. Вычисления производите до достижения заданной точности ε =103 или до тех пор, пока число итераций не превысит 104 .

2 Интерполирование алгебраическими многочленами

2.1 По заданной таблице значений функции найдите формулу интерполяционного многочлена Лагранжа. Составьте программу, реализующую данную задачу. Постройте график интерполяционного многочлена Лагранжа и отметьте на нем узловые точки

Mi (xi , yi ), i = 0, 1, 2, 3 .

Вариант

x0

x1

x2

x3

y0

y1

y2

y3

1

– 1

0

3

4

3

5

2

– 6

2

2

3

5

6

4

1

7

2

3

0

2

3

5

– 1

– 4

2

– 8

4

7

9

13

15

2

– 2

3

– 4

5

– 3

– 1

3

5

7

– 1

4

– 6

6

1

2

4

7

– 3

– 7

2

8

7

– 1

1

2

4

4

9

1

6

8

2

4

5

7

9

– 3

6

– 2

9

– 4

– 2

0

3

2

8

5

10

10

– 1

1,5

3

5

4

– 7

1

– 8

11

2

4

7

8

– 1

– 6

3

12

12

– 9

– 7

– 4

– 1

3

– 3

4

– 9

13

0

1

4

6

7

– 1

8

2

14

– 8

– 5

0

2

9

– 2

4

6

15

– 7

– 5

– 4

– 1

4

– 4

5

10

2.2 Вычислите одно значение заданной функции для промежуточного значения аргумента с помощью интерполяционного многочлена Лагранжа и оцените погрешность интерполяции. Составьте программу, реализующую данную задачу.

45

Вариант

Таблица

x

1

1

3,8

2

2

3,5

3

3

0.5

4

4

4,8

5

1

4,1

6

2

3,9

7

3

3,3

8

4

4,0

9

1

2,9

10

2

5,3

11

3

4,1

12

4

7,6

13

1

4,4

14

2

2,5

15

3

5,2

Таблица 1

Таблица 2

 

 

 

 

 

x

f (x) = (1/ x) lg x + x2

 

x

f (x) = ln2,3x 0,8/ x

 

 

 

1,2

0,3486

1,3

1,7777

 

 

 

 

1,9

1,0537

2,1

4,5634

 

 

 

 

3,3

1,7844

3,7

13,8436

 

 

 

 

4,7

2,2103

4,5

20,3952

 

 

 

 

5,4

2,3712

6,1

37,3387

 

 

 

 

6,8

2,6322

7,7

59,4051

 

 

 

 

7,5

2,7411

8,5

72,3593

 

 

 

 

 

 

 

Таблица 3

 

Таблица 4

 

 

 

 

 

 

 

 

 

 

 

 

 

x

f (x) = 2,1 sin0,37x

 

 

x

f (x) =1,7 3

x

cos(0,4 0,7x)

 

– 3,2

– 1,9449

 

 

 

 

 

 

 

 

 

 

2,6

 

2,1874

 

– 0,8

– 0.6126

 

 

 

 

 

 

 

 

 

 

3,3

 

2,8637

 

0,4

0,3097

 

 

 

 

 

 

 

 

4,7

 

3,8161

 

2,8

1,8068

 

 

 

 

 

 

 

 

6,1

 

3,8524

 

4,0

2,0913

 

 

 

 

 

 

 

 

7,5

 

3,1905

 

6,4

1,4673

 

 

 

 

 

 

 

 

8,2

 

2,8409

 

7,6

0,6797

 

 

 

 

 

 

 

 

9,6

 

2,6137

 

 

2.3 Составьте интерполяционный

 

 

 

 

 

 

многочлен Ньютона для

функции из задания 1

с помощью программы для компьютера.

46

2.4 Дана таблица значений функции y = sin x .

x

sin x

x

sin x

x

sin x

1,1

0,89121

1,6

0,99957

2,1

0,86321

1,2

0,93204

1,7

0,99166

2,2

0,80850

1,3

0,96356

1,8

0,97385

2,3

0,74571

1,4

0,98545

1,9

0,94630

2,4

0,67546

1,5

0,99749

2,0

0,90930

2,5

0,59847

Пользуясь первой или второй интерполяционными формулами Ньютона при n = 2 , вычислите sin x для следующих значений аргумента x и укажите оценку остаточного члена

R2 .

1) 1,151; 2) 1,218; 3) 1,345; 4) 1,421; 5) 1,538; 6) 1,609; 7) 1,732; 8) 1,849; 9) 1,929; 10) 2,031;

11)2,173; 12) 2,218; 13) 2,313; 14) 2,437; 15) 2,478.

2.5С помощью программы для компьютера уплотните часть таблицы заданной функции, пользуясь первой или второй интерполяционными формулами Ньютона.

Для выполнения задания 5 по заданной таблице значений функции с равноотстоящими значениями аргумента составьте таблицу конечных разностей и определите порядок интерполяционного полинома Ньютона. В зависимости от расположения участка [a,b] уплот-

нения таблицы с шагом H относительно исходной таблицы выберите первую или вторую интерполяционную формулу Ньютона. В программе необходимо совершить подсчет погрешности метода по выбранной формуле.

Вариант

Таблица

a

b

H

1

1

0,65

0,75

0,01

2

2

0,30

0,45

0,025

3

3

1,45

1,55

0,01

4

4

1,20

1,40

0,02

5

2

0,10

0,20

0,01

6

3

1,10

1,30

0,02

7

4

1,05

1,25

0,025

8

1

0,70

0,90

0,02

9

3

1,25

1,50

0,025

10

4

1,00

1,10

0,01

11

1

0,60

0,70

0,01

12

2

0,15

0,35

0,025

13

3

1,15

1,25

0,01

14

1

0,65

0,85

0,025

15

2

0,20

0,40

0,02

47

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]