84330 (Курсовая работа по численным методам)
Описание файла
Документ из архива "Курсовая работа по численным методам", который расположен в категории "". Всё это находится в предмете "математика" из , которые можно найти в файловом архиве . Не смотря на прямую связь этого архива с , его также можно найти и в других разделах. Архив можно найти в разделе "рефераты, доклады и презентации", в предмете "математика" в общих файлах.
Онлайн просмотр документа "84330"
Текст из документа "84330"
1. Методом Крылова развернуть характеристический определитель матрицы А=
. Исходную систему линейных уравнений решить методом Жордана-Гаусса.
Решение. Метод Крылова основан на свойстве квадратной матрицы обращать в нуль свой характеристический многочлен.
Согласно теореме Гамильтона-Кали, всякая квадратная матрица является корнем своего характеристического многочлена и, следовательно, обращает его в нуль.
Пусть
– (1)
характеристический многочлен.
Заменяя в выражении (1) величину 
на 
, получим
. (2)
Возьмем произвольный ненулевой вектор
. (3)
Умножим обе части выражения (2) на 
:
(4)
Положим
, (5)
т.е.
(6)
Учитывая (5), выражение (4) запишем в виде
, (7)
или в виде
Решаем систему (7). Если эта система имеет единственное решение, то ее корни 
являются коэффициентами характеристического многочлена (1).
Если известны коэффициенты и корни 
характеристического многочлена, то метод Крылова дает возможность найти соответствующие собственные векторы по следующей формуле:
(8)
Здесь 
– векторы, использованные при нахождении коэффициентов методом Крылова, а коэффициенты 
определяются по схеме Горнера
(9)
Используя все выше сказанное, развернем характеристический определитель матрицы А= методом Крылова.
Выберем в качестве начального следующий вектор:
, 
Вычислим
Составим матричное уравнение
, или 
Полученную систему уравнений решим методом Жордана-Гаусса.
|
|
|
|
|
|
| |
| 1 | 9 | 2 | 0 | -72 | -61 | -61 |
| -1 | 1 | 0 | -3 | -3 | -3 | |
| 30 | 5 | 1 | -167 | -131 | -131 | |
| 2 | 1 | 2/9 | 0 | -8 | -61/9 | -61/9 |
| 0 | 11/9 | 0 | -11 | -88/9 | -88/9 | |
| 0 | -15/9 | 1 | 657/9 | 651/9 | 651/9 | |
| 3 | 1 | 0 | 0 | -6 | -5 | -5 |
| 0 | 1 | 0 | -9 | -8 | -8 | |
| 0 | 0 | 1 | 58 | 59 | 59 | |
| 4 | 1 | 0 | 0 | |||
| 0 | 1 | 0 | ||||
| 0 | 0 | 1 |
Исходя из результатов таблицы, имеем 
.
Таким образом характеристическое уравнение матрицы имеет вид
2. Для определения собственных чисел матрицы необходимо решить полученное характеристическое уравнение третьей степени
Данное кубическое уравнение невозможно решить стандартными средствами. Воспользуемся для этой цели числовыми методами, а точнее методами приближенного вычисления.
2.1 Исследование функции.
Вычислим первую и вторую производные данной функции
Необходимо выбрать интервал, на котором будем находить решение.
Для отделения корней существует несколько способов. Наиболее популярные из них – графический и аналитический.
В литературе рассматриваются эти способы по отдельности. По заданию курсовой работы требуется отделить корни каждым из этих способов. Рискну нарушить это требование, и объединить эти два способа в один. То есть исследовать функцию аналитически и по результатам исследования построить приблизительный график функции.
Областью значений исходного уравнения является вся ось 
.
Приравняв первую производную к нулю, мы можем получить критические точки данной функции (точки минимумов и максимумов, или же точки, в которых функция не определена).
Стоит отметить, что для вычисления квадратного корня, также применимы числовые методы, на которых и основаны микрокалькуляторы и программы для ЭВМ. Данные методы основаны на логарифмировании корня и последующего вычисления.
вычисляется при помощи числового ряда
Уравнение имеет решение 
, 
. Изменив знак равенства на знак неравенства (), можем найти промежутки возрастания и убывания функции.
Функция возрастает на промежутке 
и убывает на промежутке 
. Подставив в исходное уравнение значения критических точек, имеем в результате для 
и для 
.
Приравняв вторую производную к нулю, мы можем найти точку перегиба и, соответственно, найти интервал, на котором функция выпуклая и вогнутая.
Далее необходимо найти, интервалы, в которых график функции пересекает ось 
.
Сразу можно определиться, что так при значение функции больше нуля, а при 
- меньше нуля, то одна из точек пересечения, будет лежать на данном интервале. Произведя не хитрые математические вычисления значения функции для 
, сузим интервал до 
.
Далее рассмотрим оставшиеся два интервала.
Известно, что при 
- значение функции отрицательно, а в первой критической точке положительно, то будем сужать этот промежуток. В данном случае применим метод половинного деления.
|
|
|
| 0 | 58 |
| -100 | -1059042 |
| -50 | -139492 |
| -25 | -19092 |
| -12 | -2426 |
| -6 | -320 |
| -3 | 4 |
| -5 | -172 |
| -4 | -66 |
|
|
|
| 4 | -10 |
| 100 | 939158 |
| 50 | 109608 |
| 25 | 11708 |
| 12 | 814 |
| 6 | 4 |
| 5 | -12 |
. Следующий будет от 
и до бесконечности.
Произведем аналогичные вычисления и получим промежуток 
На основании произведенного анализа построим график исходной функции.
2.2 Метод хорд.
Сразу необходимо заметить, что существуют два случая (варианта) при решении методом хорд.
Случай первый. Первая и вторая производные функции имеют одинаковые знаки, т.е. 
.
В этом случае итерационный процесс осуществляем по формуле
Случай второй. Первая и вторая производные функции имеют разные знаки, т.е. 
.
В этом случае итерационный процесс осуществляем по формуле
Для оценки точности приближение можно воспользоваться формулой
,
где 
при 
, 
– точное значение корня.
Итак решим наше уравнение методом хорд с точностью 
.
2.2.1 Интервал .
Так как первая и вторые производные в точке, от которой мы начинаем работать имеют различные знаки, то работаем по второму варианту.
Результаты вычисления приведены в таблице.
|
|
|
|
|
|
| -4,0000000 | -3,0000000 | -66,0000000 | 4,0000000 | 0,0740741 |
| -4,0000000 | -3,1142857 | -66,0000000 | -2,3688397 | 0,0438674 |
| -4,0000000 | -3,0440850 | -66,0000000 | 1,5901736 | 0,0294477 |
| -4,0000000 | -3,0901012 | -66,0000000 | -0,9879693 | 0,0182957 |
| -4,0000000 | -3,0610770 | -66,0000000 | 0,6456578 | 0,0119566 |
| -4,0000000 | -3,0798611 | -66,0000000 | -0,4086778 | 0,0075681 |
| -4,0000000 | -3,0678974 | -66,0000000 | 0,2640772 | 0,0048903 |
| -4,0000000 | -3,0755972 | -66,0000000 | -0,1684077 | 0,0031187 |
| -4,0000000 | -3,0706743 | -66,0000000 | 0,1083107 | 0,0020058 |
| -4,0000000 | -3,0738353 | -66,0000000 | -0,0692833 | 0,0012830 |
| -4,0000000 | -3,0718112 | -66,0000000 | 0,0444729 | 0,0008236 |
| -4,0000000 | -3,0731096 | -66,0000000 | -0,0284836 | 0,0005275 |
| -4,0000000 | -3,0722776 | -66,0000000 | 0,0182690 | 0,0003383 |
| -4,0000000 | -3,0728111 | -66,0000000 | -0,0117068 | 0,0002168 |
| -4,0000000 | -3,0724692 | -66,0000000 | 0,0075061 | 0,0001390 |
| -4,0000000 | -3,0726884 | -66,0000000 | -0,0048109 | 0,0000891 |
| -4,0000000 | -3,0725479 | -66,0000000 | 0,0030843 | 0,0000571 |
| -4,0000000 | -3,0726380 | -66,0000000 | -0,0019770 | 0,0000366 |
