Контрольная работа по "Математике"

Министерство образования Российской Федерации Чайковский технологический институт Ижевского государственного технического университета                                       Типовой расчет по дисциплине: Вычислительная математика         Вариант 30                   Выполнил: студент гр. АСОИиУ1-02                                                                                          Мищихин Е.А.                                                                                             Проверил: доцент                                                                                                                Морозова А.Р.               г. Чайковский 2004 г.

Содержание:

1. Приближенное решение уравнения f(x)=0 методом деления пополам………………3
2. Метод простых итераций решения уравнения f(x)=0…………………………………5
3. Приближенное решение уравнения f(x)=0 методом Ньютона………………………..7
4. Вычисление корней многочлена………………………………………………………..8
5. Решение линейной системы методом простых итераций…………………………….9
6. Решение системы нелинейных уравнений методом Ньютона………………………11
7. Решение системы линейных уравнений методом Гаусса…………………………....13
8. Решение системы линейных уравнений методом простых итераций………………15

14. Численное интегрирование…………………………………………………………….17

24. Приближенное решение задачи  Коши методом Эйлера……………………………..19

25. Приближенное решение задачи  Коши методом Рунге-Кутта……………………….20

27. Приближенное решение задачи  Коши методом Адамса-Башфорта………………...22

50. Аппроксимация функции по  методу наименьших квадратов……………………….24

Приложения…………………………………………………………………………………26

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1. Приближенное  решение уравнения f(x)=0 методом деления пополам

(метод  биссекции)

Теория:

               Пусть требуется решить уравнение f(x)=0. Предположим, что корень уравнения отделен и лежит на интервале [a; b].

             Выбираем точку с, лежащую внутри интервала: середину интервала

Вычисляем значение f(с). В качестве нового интервала изоляции мы примем ту из двух половинок интервала [a; b], на концах которой функция имеет разные знаки. Таким путем можно как угодно сузить участок, на котором находится корень.

 

Решение:

График функции f(x) (рис.1)

 

i	a	b	c	f(a)	f(b)	f(c)
1	-1,00000	1,00000	0,00000	1,00000	-1,00000	1,57080
2	0,00000	1,00000	0,50000	1,57080	-1,00000	1,19312
3	0,50000	1,00000	0,75000	1,19312	-1,00000	0,55151
4	0,75000	1,00000	0,87500	0,55151	-1,00000	0,02886
5	0,87500	1,00000	0,93750	0,02886	-1,00000	-0,32674
6	0,87500	0,93750	0,90625	0,02886	-0,32674	-0,13716
7	0,87500	0,90625	0,89063	0,02886	-0,13716	-0,05173
8	0,87500	0,89063	0,88281	0,02886	-0,05173	-0,01087
9	0,87500	0,88281	0,87891	0,02886	-0,01087	0,00913
10	0,87891	0,88281	0,88086	0,00913	-0,01087	-0,00084
11	0,87891	0,88086	0,87988	0,00913	-0,00084	0,00416
12	0,87988	0,88086	0,88037	0,00416	-0,00084	0,00166
13	0,88037	0,88086	0,88062	0,00166	-0,00084	0,00041
14	0,88062	0,88086	0,88074	0,00041	-0,00084	-0,00021
15	0,88062	0,88074	0,88068	0,00041	-0,00021	0,00010
16	0,88068	0,88074	0,88071	0,00010	-0,00021	-0,00005
17	0,88068	0,88071	0,88069	0,00010	-0,00005	0,00002
18	0,88069	0,88071	0,88070	0,00002	-0,00005	-0,00002
19	0,88069	0,88070	0,88070	0,00002	-0,00002	0,00000
20	0,88070	0,88070	0,88070	0,00000	-0,00002	-0,00001

 

 

 

 

 

 

 

 

 

 

 

 

 

F (-1) =0+1=1;

F (1) =0-1=-1;

 

2. Метод  простых итераций решения уравнения f(x)=0.

 

Теория:

Пусть требуется  решить уравнение f(x)=0. Предположим, что корень отделен и лежит на интервале [a;b].

Каким-то образом  выразим  .

Пусть - нулевое приближение. Для получения следующего приближения в правую часть уравнения вместо x подставим :

Следующие приближения  получаются аналогично:

,

,

………………

………………

Если последовательность имеет предел, то этот предел будет являться корнем уравнения.

Одно из значений с достаточно большим номером можно принять за приближенное значение корня. Однако может случиться, что последовательность не имеет предела, и тогда метод итераций не приводит к цели. Если последовательность этих приближений сходиться, то есть имеет некоторое предельное значение, то этот предел и будет являться корнем данного уравнения.

Если последовательность расходиться, то данное уравнение не имеет решения.

Также должно выполняться условие  .

Решение:

График функции f(x)(рис.1)

x	F(x)
0,00000	1,00000	-1,00000
1,00000	0,73505	0,26495
0,73505	0,96030	-0,22524
0,96030	0,79553	0,16477
0,79553	0,93591	-0,14038
0,93591	0,82607	0,10984
0,82607	0,91940	-0,09333
0,91940	0,84434	0,07505
0,84434	0,90788	-0,06353
0,90788	0,85606	0,05181
0,85606	0,89977	-0,04370
0,89977	0,86384	0,03593
0,86384	0,89405	-0,03021
0,89405	0,86909	0,02496
0,86909	0,89004	-0,02094

 

 

 

0,89004	0,87268	0,01735
0,87268	0,88722	-0,01453
0,88722	0,87515	0,01207
0,87515	0,88524	-0,01009
0,88524	0,87685	0,00839
0,87685	0,88387	-0,00701
0,88387	0,87803	0,00584
0,87803	0,88290	-0,00487
0,88290	0,87884	0,00406
0,87884	0,88223	-0,00339
0,88223	0,87941	0,00282
0,87941	0,88176	-0,00236
0,88176	0,87980	0,00196
0,87980	0,88144	-0,00164
0,88144	0,88007	0,00136
0,88007	0,88121	-0,00114
0,88121	0,88026	0,00095
0,88026	0,88106	-0,00079
0,88106	0,88040	0,00066
0,88040	0,88095	-0,00055
0,88095	0,88049	0,00046
0,88049	0,88087	-0,00038
0,88087	0,88055	0,00032
0,88055	0,88082	-0,00027
0,88082	0,88060	0,00022
0,88060	0,88078	-0,00018
0,88078	0,88063	0,00015
0,88063	0,88075	-0,00013
0,88075	0,88065	0,00011
0,88065	0,88074	-0,00009
0,88074	0,88066	0,00007
0,88066	0,88072	-0,00006
0,88072	0,88067	0,00005
0,88067	0,88072	-0,00004
0,88072	0,88068	0,00004
0,88068	0,88071	-0,00003
0,88071	0,88068	0,00003
0,88068	0,88071	-0,00002
0,88071	0,88069	0,00002
0,88069	0,88070	-0,00001
0,88070	0,88069	0,00001
0,88069	0,88070	-0,00001
0,88070	0,88069	0,00001
0,88069	0,88070	-0,00001
0,88070	0,88069	0,00001
0,88069	0,88070	0,00000

 

 

Ответ:

 

3.Приближенное  решение уравнения f(x)=0 методом Ньютона.

Теория:

 

Пусть требуется  решить уравнение f(x)=0. Допустим, что корень отделен и находиться на интервале [a;b].

Чтобы найти  корень уравнения проведем касательную в т. B(b; f(b)) и точка пересечения Ох и касательной будет являться первым приближением.

Чтобы найти  составим уравнение касательной:

В качестве приближенного  корня уравнения примем абсциссу точки пересечения с осью Ох. Полагая  в уравнении касательной у=0, находим  абсциссу точки пересечения

За исходную точку берем ту точку, в которой  знак функции совпадает со знаком второй производной.

 

Решение:

График функции f(x) (рис.1).

x	F`(x)	F(x)
0,99	-12,9149	-0,77047
0,9303427	-6,31171	-0,28067
0,8858741	-5,21298	-0,02674
0,880744	-5,11831	-0,00024
0,8806961	-5,11745	0,00000
0,8806961	-5,11745	0,00000

 

Ответ: х=0,8807.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.Вычисление  корней многочлена.

Решение:

Вычислим корни  многочлена методом касательных (методом  Ньютона)

График функции f(x) (рис.2)

Из графика  видны интервалы изоляции корней:

[-3;-2], [-1; 0]

 

1)[-3;-2]

 

x	f(x)	f’(x)
-3	198,00000	-727
-2,727648	51,86437	-373,021
-2,588609	9,06134	-248,108
-2,552087	0,51102	-220,457
-2,549769	0,00196	-218,769
-2,54976	0,00000	-218,762

 

2) [-1; 0]

 

x	f(x)	f’(x)
-1	-8,00000	29
-0,724138	-1,46222	18,7304
-0,646071	-0,09731	16,2784
-0,640093	-0,00053	16,10092
-0,64006	0,00000	16,09995

 

Ответ: x=-2,54976; x=-0,64006.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.Решение линейной  системы уравнений методом простых  итераций.

Теория:

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

 

Рассмотрим систему  n линейных уравнений с n неизвестными:

Решив первое уравнение относительно второе – относительно , и т.д., систему приведем к виду:

Запишем ее в  более общем виде:

Зададим какие-либо начальные значения неизвестных

.

Подставляя  их в правые части системы, получим первые приближения:

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

 

 

 

 

 

Решение:

График функций X(y)=0,5sin(y)-0,8 b Y(x)=1-cos(x-1)(рис.3).

x	Y
0	0
1,220735	0,459698
1,057198	0,024263
1,214058	0,001635
1,220293	0,022823
1,214461	0,024167
1,214085	0,022909
1,214437	0,022829
1,214459	0,022904
1,214438	0,022908
1,214437	0,024904

 

Ответ: x=1.2144; y=0.0229

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6.Решение системы нелинейных уравнений методом Ньютона.

 

 

Теория:

 

Пусть требуется решить систему нелинейных уравнений:

Допустим найдены нулевые  приближения

Чтобы найти  решение системы необходимо найти  якобиан:

Следующие приближения  находим по формулам:

 

Решение:

 

 

 

 

     

Ответ:      

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.Решение системы  линейных уравнений методом Гаусса. Решить систему линейных уравнений  вида Ax=B.

Теория:

Метод Гаусса один из наиболее распространенных способов решения систем линейных уравнений. Он является точным: если точно выполнить  все действия, то мы получим точное решение системы. Но практически решение получается приближенным из-за округления при делении. В основе метода лежит идея последовательного исключения неизвестных.

 

Рассмотрим  систему линейных уравнений четвертого порядка:

Разделим новое  уравнение на (пусть ; в противном случае переставим уравнение так, чтобы это условие было выполнено). Находя из первого уравнения, получим: