Курсовая работа по численным методам

1. Методом Крылова развернуть характеристический определитель матрицы А . Исходную систему линейных уравнений решить методом Жордана-Гаусса. Решение. Метод Крылова основан на свойстве квадратной матрицы обращать в нуль свой характеристический многочлен.Согласно теореме Гамильтона-Кали, всякая квадратная матрица является корнем своего характеристического многочлена и, следовательно, обращает его в нуль. Пусть 1 характеристический многочлен. Заменяя в выражении 1 величину на , получим . 2 Возьмем произвольный ненулевой вектор . 3 Умножим обе части выражения 2 на 4 Положим , 5 т.е. 6 Учитывая 5, выражение 4 запишем в виде , 7 или в виде Решаем систему 7. Если эта система имеет единственное решение, то ее корни являются коэффициентами характеристического многочлена 1. Если известны коэффициенты и корни характеристического многочлена, то метод Крылова дает возможность найти соответствующие собственные векторы по следующей формуле 8 Здесь векторы, использованные при нахождении коэффициентов методом Крылова, а коэффициенты определяются по схеме Горнера 9 Используя все выше сказанное, развернем характеристический определитель матрицы А методом Крылова. Выберем в качестве начального следующий вектор , Вычислим Составим матричное уравнение , или Полученную систему уравнений решим методом Жордана-Гаусса. -72-61-61-110-3-3-33051-167-131-13121290 -8-619-61901190-11-889-8890-159165796519 65193100-6-5-5010-9-8-01Исходя из результатов таблицы, имеем . Таким образом характеристическое уравнение матрицы имеет вид 2. Для определения собственных чисел матрицы необходимо решить полученное характеристическое уравнение третьей степени Данное кубическое уравнение невозможно решить стандартными средствами. Воспользуемся для этой цели числовыми методами, а точнее методами приближенного вычисления. 2.1 Исследование функции. Вычислим первую и вторую производные данной функции Необходимо выбрать интервал, на котором будем находить решение.

Для отделения корней существует несколько способов.

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

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

Рискну нарушить это требование, и объединить эти два способа в один. То есть исследовать функцию аналитически и по результатам исследования построить приблизительный график функции.

Областью значений исходного уравнения является вся ось . Приравняв первую производную к нулю, мы можем получить критические точки данной функции точки минимумов и максимумов, или же точки, в которых функция не определена.Стоит отметить, что для вычисления квадратного корня, также применимы числовые методы, на которых и основаны микрокалькуляторы и программы для ЭВМ. Данные методы основаны на логарифмировании корня и последующего вычисления. вычисляется при помощи числового ряда Уравнение имеет решение Изменив знак равенства на знак неравенства или , можем найти промежутки возрастания и убывания функции.

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

Далее необходимо найти, интервалы, в которых график функции пересекает ось . Сразу можно определиться, что так при значение функции больше нуля, а при - меньше нуля, то одна из точек пересечения, будет лежать на данном интервале. Произведя не хитрые математические вычисления значения функции для , сузим интервал до . Далее рассмотрим оставшиеся два интервала.Известно, что при - значение функции отрицательно, а в первой критической точке положительно, то будем сужать этот промежуток. В данном случае применим метод половинного деления. -100-1059042-50-139492-25-19092-12-2426- 6-320-34-5-172-4 -12Таким образом получили еще один интервал . Следующий будет от и до бесконечности.

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

Первая и вторая производные функции имеют разные знаки, т.е В этом случае итерационный процесс осуществляем по формуле Для оценки точности приближение можно воспользоваться формулой , где при , точное значение корня. Итак решим наше уравнение методом хорд с точностью . 2.1 Интервал . Так как первая и вторые производные в точке, от которой мы начинаем работать имеют различные знаки, то работаем по второму варианту.

Результаты вычисления приведены в таблице. -4,0-3,0-66,04,0,0740741-4,0-3,1142857-6 6,0-2,36883970,0438674-4,0-3,0440850-66, 01,59017360,0294477-4,0-3,0901012-66,0-0 ,98796930,0182957-4,0-3,0610770-66,0,645 65780,0119566-4,0-3,0798611-66,0-0,40867 780,0075681-4,0-3,0678974-66,0,26407720, 0048903-4,0-3,0755972-66,0-0,16840770,00 31187-4,0-3,0706743-66,0,10831070,002005 8-4,0-3,0738353-66,0-0,06928330,0012830- 4,0-3,0718112-66,0,04447290,0008236-4,0- 3,0731096-66,0-0,02848360,0005275-4,0-3, 0722776-66,0,01826900,0003383-4,0-3,0728 111-66,0-0,01170680,0002168-4,0-3,072469 2-66,0,00750610,0001390-4,0-3,0726884-66 ,0-0,00481090,0000891-4,0-3,0725479-66,0 ,00308430,0000571-4,0-3,0726380-66,0-0,0 0197700, 2.2.2 Интервал . Так как первая и вторые производные в точке, от которой мы начинаем работать имеют различные знаки, то работаем по второму варианту.

Результаты вычисления приведены в таблице. ,04,04,0-10,0-0,23,03,28571434,0-0,87463 56-0,04859093,03,23444984,0-0,0423087-0, 00235053,03,23199594,0-0,0019734-0,00010 963,03,23188154,0-0,0000919-0, 2.2.3 Интервал . Так как первая и вторые производные в точке, от которой мы начинаем работать имеют одинаковые знаки, то работаем по первому варианту.

Результаты вычисления приведены в таблице. ,06,0-12,04,0,675,7506,0-2,01562504,0,33 593755,83376626,0-0,16130144,0,02688365, 84020986,0-0,01201984,0,00200335,8406885 6,0-0,00089094,0,00014855,84072406,0-0,0 0006604,0, Итак, корнями уравнения будут 2.3 Метод касательных метод Ньютона.

В век повальной компьютеризации не есть хорошо считать при помощи логарифмической линейки. Поэтому, разработаем алгоритм и прикладную программу для решения кубических уравнений методом Ньютона.Ниже приведена блок-схема алгоритма и листинг программы, реализующей данный алгоритм на языке С. Также привожу текст, которая выдает данная программа при решении исходного уравнения. метод Ньютона длЯ решениЯ кубических уравнений include math.h include iostream.h double a40, b30, c20, prec0.0 double minim0, maxim0 void Hellovoid void Input void Derivative void Calculation double CalcFundouble double CalcFirstdouble double CalcSeconddouble mainvoid Hello Input Derivative Calculation return 0 void Hellovoid cout Программа длЯ решениЯ кубических уравнений методом касательных метод Ньютона.nn void Input cout Кубическое уравнение имеет вид endl a1x3a2x2a3xa40 endl endl for int i0i 4i cout Введите значение коэффициента a i1 cin ai cout endl Необходимо указать интервал поиска решениЯ. endl Введите нижнюю границу поиска cin minim cout Введите верхнюю границу поиска cin maxim whileminimmaximminim maxim cout nНижнЯЯ граница должна быть меньше верхней и не может быть ей равна. endl Повторите ввод нижней границы cin minim cout Повторите ввод верхней границы cin maxim cout Введите допустимую погрешность cin prec void Derivative b0a03 b1a12 b2a2 c0b02 c1b1 cout nnn Исходное уравнение имеет вид nn a0 x3 a1 x2 a2 x a3 0nn ПерваЯ производнаЯ имеет вид nn fx b0 x2 b1 x b2 nn ВтораЯ производнаЯ имеет вид nn fx c0 x c1 nn void Calculation double x0, m0 cout endl Xn fXn fXnm endl endl if absCalcFunminimabsCalcSecondminim 0 xminim else xmaxim if CalcFirstminim CalcFirstmaxim mabsCalcFirstmaxim else mabsCalcFirstminim cout cout.width15cout.precision10 cout x cout cout.width15cout.precision10 cout CalcFunx cout cout.width15cout.precision10 cout fabsCalcFunxm cout n whilefabsCalcFunxm prec xx-CalcFunxCalcFirstx cout cout.width15cout.precision10 cout x cout cout.width15cout.precision10 cout CalcFunx cout cout.width15cout.precision10 cout fabsCalcFunxm cout n cout double CalcFundouble x return a0xxxa1xxa2xa3 double CalcFirstdouble x return b0xxb1xb2 double CalcSeconddouble x return c0xc1 Программа длЯ решениЯ кубических уравнений методом касательных метод Ньютона.

Кубическое уравнение имеет вид a1x3a2x2a3xa40 Введите значение коэффициента a1 1 Введите значение коэффициента a2 -6 Введите значение коэффициента a3 -9 Введите значение коэффициента a4 58 Необходимо указать интервал поиска решениЯ. Введите нижнюю границу поиска -4 Введите верхнюю границу поиска -3 Введите допустимую погрешность 0.00005 Исходное уравнение имеет вид 1x3-6x2-9x580 ПерваЯ производнаЯ имеет вид fx3x2-12x-9 ВтораЯ производнаЯ имеет вид fx6x-12 Xn fXn fXnm -4 -66 1.2 -3.24137931 -9.922506048 0.183750112 -3.079817529 -0.40621762 0.007522548518 -3.07261683-0.0007897932301.462580056e-0 5 Программа длЯ решениЯ кубических уравнений методом касательных метод Ньютона.

Кубическое уравнение имеет вид a1x3a2x2a3xa40 Введите значение коэффициента a1 1 Введите значение коэффициента a2 -6 Введите значение коэффициента a3 -9 Введите значение коэффициента a4 58 Необходимо указать интервал поиска решениЯ. Введите нижнюю границу поиска 3 Введите верхнюю границу поиска 4 Введите допустимую погрешность 0.00005 Исходное уравнение имеет вид 1x3-6x2-9x580 ПерваЯ производнаЯ имеет вид fx3x2-12x-9 ВтораЯ производнаЯ имеет вид fx6x-12 Xn fXn fXnm 3 4 0.4 3.2 0.159122085 0.01768023167 3.231855174 0.0003411376333.790418145e-05 Программа длЯ решениЯ кубических уравнений методом касательных метод Ньютона.

Кубическое уравнение имеет вид a1x3a2x2a3xa40 Введите значение коэффициента a1 1 Введите значение коэффициента a2 -6 Введите значение коэффициента a3 -9 Введите значение коэффициента a4 58 Необходимо указать интервал поиска решениЯ. Введите нижнюю границу поиска 5 Введите верхнюю границу поиска 6 Введите допустимую погрешность 0.00005 Исходное уравнение имеет вид 1x3-6x2-9x580 ПерваЯ производнаЯ имеет вид fx3x2-12x-9 ВтораЯ производнаЯ имеет вид fx6x-12 Xn fXn fXnm 6 4 0.67 5.851851852 0.2601229487 0.04335382479 5.840787634 0.001413241032 0.000235540172 5.8407268624.255405933e-087.092343222e-0 9 2.4 Метод итераций.

Как и для предыдущего метода, привожу блок-схему алгоритма решения и листинг программы, реализующей этот алгоритм на языке программирования С. метод итераций длЯ решениЯ кубических уравнений include math.h include iostream.h double a40, b30, prec0.0 double minim0, maxim0 void Hellovoid void Input void Derivative void Calculation double CalcFundouble double CalcFirstdouble mainvoid Hello Input Derivative Calculation return 0 void Hellovoid cout Программа длЯ решениЯ кубических уравнений методом итераций.nn void Input cout Кубическое уравнение имеет вид endl a1x3a2x2a3xa40 endl endl for int i0i 4i cout Введите значение коэффициента a i1 cin ai cout endl Необходимо указать интервал поиска решениЯ. endl Введите нижнюю границу поиска cin minim cout Введите верхнюю границу поиска cin maxim whileminimmaximminim maxim cout nНижнЯЯ граница должна быть меньше верхней и не может быть ей равна. endl Повторите ввод нижней границы cin minim cout Повторите ввод верхней границы cin maxim cout Введите допустимую погрешность cin prec void Derivative b0a03 b1a12 b2a2 void Calculation double x0, xold0, m0 cout endl Xn fXn Xn1-Xn endl endl iffabsCalcFirstminim fabsCalcFirstmaxim mxxoldminim else mxxoldmaxim mfabs1CalcFirstm cout cout.width15cout.precision10 cout x cout cout.width15cout.precision10 cout CalcFunx cout n ifCalcFirstx 0 do xoldx xxold-mCalcFunxold cout cout.width15cout.precision10 cout x cout cout.width15cout.precision10 cout CalcFunx cout cout.width15cout.precision10 cout fabs CalcFunx - CalcFunxold cout n while fabs CalcFunx - CalcFunxold prec else do xoldx xxoldmCalcFunxold cout cout.width15cout.precision10 cout x cout cout.width15cout.precision10 cout CalcFunx cout cout.width15cout.precision10 cout fabs CalcFunx - CalcFunxold cout n while fabs CalcFunx - CalcFunxold prec cout double CalcFundouble x return a0xxxa1xxa2xa3 double CalcFirstdouble x return b0xxb1xb2 Программа длЯ решениЯ кубических уравнений методом итераций. Кубическое уравнение имеет вид a1x3a2x2a3xa40 Введите значение коэффициента a1 1 Введите значение коэффициента a2 -6 Введите значение коэффициента a3 -9 Введите значение коэффициента a4 58 Необходимо указать интервал поиска решениЯ. Введите нижнюю границу поиска -4 Введите верхнюю границу поиска -3 Введите допустимую погрешность 0.00005 Xn fXn Xn1-Xn -4 -66 -3.24137931 -9.922506048 56.07749395 -3.127327517 -3.12093462 6.801571427 -3.091454705 -1.064778438 2.056156183 -3.079215872 -0.372281515 0.6924969227 -3.074936774 -0.131239433 0.241042082 -3.073428275 -0.04639844126 0.08484099175 -3.07289496 -0.01642029825 0.02997814301 -3.072706221-0.005813178631 0.01060711962 -3.072639403-0.002058264249 0.003754914382 -3.072615744-0.000728799396 0.001329464852 -3.072607367-0.0002580606280.00047073876 78 -3.072604401-9.137721784e-00.00016668341 08 -3.072603351-3.235601088e-05.902120696e- 05 -3.072602979-1.145703711e-02.089897377e- 05 Программа длЯ решениЯ кубических уравнений методом итераций. Кубическое уравнение имеет вид a1x3a2x2a3xa40 Введите значение коэффициента a1 1 Введите значение коэффициента a2 -6 Введите значение коэффициента a3 -9 Введите значение коэффициента a4 58 Необходимо указать интервал поиска решениЯ. Введите нижнюю границу поиска 3 Введите верхнюю границу поиска 4 Введите допустимую погрешность 0.00005 Xn fXn Xn1-Xn 3 4 3.2 0.159122085 3.840877915 3.231062338 0.01338370012 0.1457383849 3.231805877 0.001151957391 0.01223174272 3.2318698759.934183961e-05 0.001052615552 3.2318753948.568402322e-069.077343728e-0 5 3.231875877.390497921e-07 7.82935253e-06 Программа длЯ решениЯ кубических уравнений методом итераций.

Кубическое уравнение имеет вид a1x3a2x2a3xa40 Введите значение коэффициента a1 1 Введите значение коэффициента a2 -6 Введите значение коэффициента a3 -9 Введите значение коэффициента a4 58 Необходимо указать интервал поиска решениЯ. Введите нижнюю границу поиска 5 Введите верхнюю границу поиска 6 Введите допустимую погрешность 0.00005 Xn fXn Xn1-Xn 6 4 5.851851852 0.2601229487 3.739877051 5.842217669 0.0346921878 0.2254307609 5.840932773 0.004788677115 0.02990351069 5.8407554140.0006639855431 0.004124691572 5.8407308229.212373716e-050.000571861805 9 5.840727411.278267885e-057.934105832e-05 5.8407269371.773688694e-061.100899016e-0 5 Решив уравнение , получили корень МетодКорень 1Корень 2Корень 3Хорд-3,0726383,2318815,840724Касательны х Ньютона-3,0726163,2318555,840726Итераций -3,0726023,2318755,840726 Для дальнейших расчетов будем использовать среднее арифметическое значение полученных корней Используя полученные значения, определим собственные значения исходной матрицы. Собственные вектора матрицы А определим по формуле Для нашей матрицы, данная формула примет следующий вид Коэффициенты определяются по схеме Горнера Для имеем Для имеем Для имеем Далее можем найти собственные векторы . Для контроля полученных значений, развернем исходную матрицу А , и определим ее собственные векторы методом непосредственного развертывания.

Характеристический многочлен для данной матрицы имеет вид . Находим . Число диагональных миноров второго порядка у матрицы второго порядка . Выписываем эти миноры и складываем их . И, в заключение, находим Таким образом, характеристическое уравнение имеет вид Данное уравнение идентично уравнению, полученному при помощи метода Крылова.

Нет смысла заново его решать.

Воспользуемся уже вычисленными корнями их средним значением.

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

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

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