Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Методы вычислений в Excel.pdf
Скачиваний:
187
Добавлен:
14.03.2016
Размер:
998 Кб
Скачать

ПРИМЕРЫ РЕШЕНИЯ ЗАДАЧ И РЕКОМЕНДАЦИИ К ЭКЗАМЕНУ

Задача 1.

Получить

решение

уравнения

f x x3 x 2

9x 9 0 методом деления отрезка пополам с

точностью 0,05. Интервал изоляции 4; 3,8 .

Проверим, что данных отрезок является интервалом

изоляции корня. Найдем значение

функции на

концах этого

интервала: f 4 3,

f 3,8

2,768 . Т.е.

на интервале

содержится корень уравнения. Проверим, что он единственный.

 

 

 

2x 9 0

 

 

f

x 3x

2

 

 

 

 

 

 

 

 

 

 

 

 

x

 

2

 

4 4 3 9

 

1,43

 

 

 

 

 

 

 

 

 

1

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

2

 

4 4 3 9

 

2,097

 

 

 

6

 

 

 

 

 

 

 

 

 

 

на всем интервале x x2 , а

 

Следовательно,

 

 

 

f x 0

x2

2,097 3,8 , т.е. функция

f x монотонно возрастающая

на

4; 3,8 ,

следовательно, данный интервал содержит один

корень уравнения и является интервалом изоляции.

Расчеты проведем в Excel по методу деления отрезка пополам, результаты оформим в виде таблицы.

k

a

b

c

f a

f c

 

b a

 

 

 

 

0

-4

-3.8

-3.9

-3

-0.009

0.2

1

-3.9

-3.8

-3.85

-0.009

1.405875

0.1

2

-3.9

-3.85

-3.875

-0.009

0.705078

0.05

3

-3.9

-3.875

-3.8875

-0.009

0.349705

0.025

Расчетные формулы:

88

ai

c

i 1

,

если f

a

i

1

f

c

i 1

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ai 1 ,

в противном случае

 

 

 

bi

b

 

 

,

если f a

 

 

f

c

 

 

0;

 

 

 

 

i 1

 

 

 

 

 

i 1

 

 

i 1

 

 

 

 

 

 

 

 

ci 1 ,

 

в противном случае.

 

 

 

 

ci

 

ai

bi

, i 0,1,2,...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: x=-3.8875

 

 

 

 

 

 

 

 

 

 

 

Задача 2.

Получить

 

решение

уравнения

f x x3

x 2

9x 9 0

методом

простой

итерации

с

точностью 0.001. Интервал изоляции 5, 3 .

 

 

 

Аналогично доказываем, что интервал является интервалом

изоляции.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

5 46 , f 3 18

 

 

 

 

 

 

 

 

 

на

всем

интервале

x x2

, а x2 2,097 3 ,

т.е.

f

x 0

функция

 

f x

монотонно

возрастающая

на 5, 3 ,

следовательно, данный интервал является интервалом изоляции.

Расчетные формулы:

 

 

xk 1 xk cf xk

 

 

 

x0

 

a b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

c

 

 

2

 

 

, M max f

 

 

 

 

 

 

 

 

 

m M

 

 

x a,b

x a,b

 

Найдем значения констант. Для этого вычислим значения

первой и второй производных:

 

 

 

 

 

 

 

 

2x 9

 

 

 

f x 3x

2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 0

 

 

f

x 6x

 

 

 

Экстремум производной функции находится в точке

f

x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

1 5, 3 .

 

 

 

 

 

 

 

3

 

 

 

 

 

 

89

Находим значения производной на концах отрезка

 

10 9 56

 

 

 

 

 

 

 

 

 

 

f 5 75

 

 

 

 

 

 

 

 

 

 

 

6 9 12 .

 

 

 

 

 

 

 

 

f 3 27

 

2

 

 

 

 

 

 

Таким образом, M 56, m 12, c

 

 

0,0294 .

 

12 56

 

 

 

 

 

 

 

 

 

 

 

 

xk 1 xk 0,0294 xk3 xk2 9xk 9 ,

x0

5 3

4 .

Вычисления оформляем в таблице

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

k

xk

 

f xk

 

 

xk xk 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

-4

 

-3

 

 

 

 

 

 

 

 

 

 

1

-3.9118

 

-0.35069

 

 

0.0882

 

 

 

2

-3.90149

 

-0.05198

 

 

0.01031

 

 

 

3

-3.89996

 

-0.0079

 

 

0.001528

 

 

 

4

-3.89973

 

-0.0012

 

 

0.000232

 

 

Ответ: x= -3,8997

 

 

Задача 3.

Получить

решение

уравнения

f x x3 x 2

9x 9 0 методом

Ньютона

с точностью

0,001. Интервал изоляции 5, 3 .

Проверка значений интервала изоляции была сделана в

примере выше. Расчетные формулы метода Ньютона:

xk 1 xk

f xk

 

 

 

 

 

 

 

 

или в нашем случае

 

 

f xk

 

 

 

xk 1 xk

x3

x 2

9x

k

9

 

 

 

k

k

 

 

 

 

 

 

 

3xk2 2xk 9

 

 

 

 

 

Выбираем

 

 

нулевое

 

приближение.

f

 

 

 

x 6x 2 ,

f 5 28 , f 5 125 25 45 9 46 . Знак функции и

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

x0 5 .

Результаты представлены в таблице

90

 

 

 

 

 

 

 

 

 

k

xk

f xk

f xk

 

 

xk xk 1

 

 

 

 

 

 

 

 

 

 

 

 

 

0

-5

-46

56

 

 

 

 

 

1

-4.17857

-8.89217

35.02423

 

0.821429

2

-3.92469

-0.72721

29.36009

 

0.253886

3

-3.89992

-0.00659

28.82821

 

0.024769

4

-3.89969

-5.6E-07

28.82332

 

0.000229

Ответ: x=-3,89969

Задача 4. Решить систему линейных уравнений методом простой итерации с точностью 0,05:

9x1 x2 3x3 3 2x1 8x2 x3 1

x1 x2 6x3 6

Проверим условие диагонального преобладания:

9 1 3 4

8 2 1 3

6 1 1 2

Условия диагонального преобладания выполняются. Разрешим систему уравнений относительно xi

x1k 1

 

1

3 x2k 3x3k

 

 

9

 

x2k 1

 

1

1 2x1k x3k

 

 

8

 

x3k 1

 

1

6 x1k x2k

 

 

6

 

xi0 0

Результаты можно представить в виде таблицы

k

x1k

x2k

x3k

max

 

xik 1 xik

 

 

 

 

 

 

 

 

 

 

0

0

0

0

 

 

1

-0.33333

-0.125

-1

1

91

2

-0.01389

-0.33333

-1.07639

0.319444

3

-0.01157

-0.26302

-1.05787

0.070313

4

-0.00993

-0.26013

-1.04577

0.012105

Ответ содержится в последней строке таблицы.

Задача 5. Решить систему линейных уравнений методом Гаусса-Зейделя с точностью 0,05:

9x1 x2 3x3 3 2x1 8x2 x3 1

x1 x2 6x3 6

Аналогично проверяем условие диагонального преобладания.

Разрешим систему уравнений относительно xi

x1k 1

 

1

3 x2k 3x3k

 

 

9

 

x2k 1

 

1

1 2x1k 1 x3k

 

 

8

 

x3k 1

 

1

6 x1k 1 x2k 1

 

 

6

 

xi0 0

Результаты удобно оформить в таблицу

k

x1k

x2k

x3k

max

 

xik 1 xik

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

0

0

 

 

 

 

 

1

-0.33333

-0.20833

-1.09028

1.090278

2

0.006944

-0.25955

-1.0421

0.340278

3

-0.01481

-0.25896

-1.04563

0.02175

Ответ содержится в последней строке таблицы.

Задача 6. Для таблично заданной функции:

x

-2

1

1,5

2

f

0,1

-0,2

0,5

1,2

92

вычислить значение функции в точке z=1,2, используя формулы линейной интерполяции.

Определяем интервал, которому принадлежит z : 1;1,5 .

Расчетные формулы:

 

 

f z

fi fi 1

z xi fi .

 

 

 

 

 

 

xi xi 1

 

 

Тогда

 

 

f 1,2 0,5 0,2 1,2 1,5 0,5

0,7

0,3 0,5 0,08 .

 

1,5 1

0,5

 

Ответ: 0,08.

 

 

Задача 7. Для таблично заданной функции:

 

 

 

 

x

 

-2

 

1

 

 

 

 

 

1,5

 

2

 

 

 

 

 

 

 

 

 

 

 

f

 

0,1

 

-0,2

 

 

0,5

 

1,2

 

 

 

 

 

 

выписать базисные полиномы и вычислить значение

полинома Лагранжа в точке z=1,2; при n=3.

 

 

 

 

 

 

 

l0

z

 

z x1 z x2

z x3

 

 

 

 

1,2 1 1,2 1,5 1,2 2

 

0,00114

 

x0 x1

x0 x2 x0 x3

 

 

2 1 2 1,5 2

2

 

 

 

 

 

 

 

l1

z

 

z x0

z x2

z x3

 

 

1,2 2 1,2 1,5 1,2 2

0,512

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1 x0

x1 x2

x1 x3

 

 

1 2 1

1,5 1 2

 

 

l2

z

 

z x0 z x1

z x3

 

 

 

 

1,2 2 1,2

1 1,2

2

0,585

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 x0

x2 x1 x2 x3

 

 

1,5 2 1,5 1 1,5

2

l3

z

 

z x0 z x1

z x2

 

 

 

1,2 2 1,2 1 1,2 1,5

0,096

 

x3 x0

x3 x1 x3 x2

 

 

 

 

2 2 2 1 2 1,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L3 z f 0 l0 z f1l1 z f 2 l2 z f3l3 z

 

 

 

 

 

 

 

 

 

0,00114*0,1 0,512*0,2 0,585*0,5 0,096*1,2 0,0748

Ответ: 0,0748.

Примеры решения задач на интегрирование можно найти выше.

93

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

таблично:

 

x

-1

-0,5

0

1

 

 

f

2

3

4

4,5

 

В

данном

задании

переменная

x меняется с постоянным

шагом 0,5. При использовании формулы трапеций решение ищется в виде:

f x0

f xn

n 1

 

2

4,5

 

 

 

I h

 

 

 

f xi

0,5

 

 

3 4

 

5,125 .

 

2

 

 

2

 

 

i 1

 

 

 

 

 

В случае же, если шаг не постоянный, например:

 

 

 

x

-1

 

 

-0,6

0

 

 

0,8

 

 

 

 

f

2

 

 

3

4

 

 

4,5

 

 

необходимо

пользоваться общей формулой трапеций:

 

n

xi f xi 1

 

 

 

 

 

 

 

 

 

I

f

xi xi 1

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

2

 

 

 

 

 

 

 

 

 

Т.е. в данном случае решение ищется как

 

 

I 3 2

0,6 1

4 3

0 0,6

4,5 4

0,8 0 6,5

2

2

2

 

 

 

 

 

 

 

 

Аналогично для формул левых и правых прямоугольников:

n

xi xi 1

,

n

.

I лев f xi 1

Iправ f xi xi xi 1

i 1

 

 

i 1

 

94

ВВЕДЕНИЕ .............................................................................

3

ПРИБЛИЖЕННОЕ РЕШЕНИЕ НЕЛИНЕЙНЫХ

 

АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ....................................

5

Постановка задачи..............................................................

5

Приближенные (итерационные) методы решения НАУ10

Метод деления отрезка пополам (дихотомии)...........

10

Метод простой итерации..............................................

12

Метод релаксации.........................................................

14

Метод Ньютона (касательных) ....................................

16

Метод хорд ....................................................................

20

МЕТОДЫ РЕШЕНИЯ СИСТЕМ ЛИНЕЙНЫХ

 

АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ..................................

23

Постановка задачи............................................................

23

Прямые методы решения СЛАУ.....................................

24

Метод Крамера..............................................................

24

Метод обратной матрицы.............................................

26

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

27

Метод прогонки ................................................................

31

Итерационные методы решения линейных

 

алгебраических систем.....................................................

33

Метод простой итерации..............................................

33

Метод Якоби..................................................................

33

Метод Гаусса-Зейделя..................................................

38

АППРОКСИМАЦИЯ ФУНКЦИЙ ......................................

41

Постановка задачи интерполяции...................................

41

Локальная интерполяция..............................................

42

Кусочно-постоянная интерполяция ............................

42

Кусочно-линейная интерполяция................................

43

Кубический интерполяционный сплайн.....................

45

Глобальная интерполяция................................................

49

Полином Лагранжа.......................................................

50

Подбор эмпирических формул....................................

53

 

95

Метод наименьших квадратов.....................................

55

ЧИСЛЕННОЕ ИНТЕГРИРОВАНИЕ..................................

60

Постановка задачи........................................................

60

Формулы прямоугольников.........................................

61

Формула трапеций........................................................

64

Формула Симпсона.......................................................

65

ЧИСЛЕННОЕ РЕШЕНИЕ ОБЫКНОВЕННЫХ

 

ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ..........................

69

Постановка задачи............................................................

69

Приближенные методы решения задачи Коши для ОДУ

первого порядка ................................................................

70

Метод Эйлера....................................................................

71

Модифицированный метод Эйлера.................................

72

Методы Рунге-Кутты........................................................

73

Численные методы решения систем ОДУ первого

 

порядка...............................................................................

75

МЕТОД КОНЕЧНЫХ РАЗНОСТЕЙ РЕШЕНИЯ

 

КРАЕВЫХ ЗАДАЧ ДЛЯ ОДУ.............................................

80

Постановка задачи........................................................

80

Аппроксимация производных......................................

81

ПРИМЕРЫ РЕШЕНИЯ ЗАДАЧ И РЕКОМЕНДАЦИИ К

 

ЭКЗАМЕНУ...........................................................................

88

96