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

Пензенский приборостроительный колледж на тему Метод касательных решения нелинейных уравнений Выполнил Ст-т 22п группы ЛЯПИН Р.Н. Проверила Ковылкино 1999 г. ЗАДАНИЕ НА КУРСОВУЮ РАБОТУ студент Ляпин Р.Н. группа 22п 1. Тема Метод касательных решения нелинейных уравнений. 2. Изучить теоретический материал по заданной теме.

3. Составить блок схему алгоритма решения задачи . 4. Написать программу на языке Турбо-Паскаль для решения задачи в общем виде. 5. Выполнить программу с конкретными значениями исходных данных. 6. Определить корни уравнения х3 0,1 х2 0,4 х 1,2 0 аналитически и уточнить один из них с точностью до 0,01 методом касательных 7. Срок представления работы к защите 10 мая 1999 г.

8. Исходные данные для исследования научная и техническая литература. Руководитель курсовой работы Кривозубова С.А. Задание принял к исполнению Ляпин Р.Н. РЕФЕРАТ Курсовая работа содержит страниц, 1 график, 5 источников. Перечень ключевых понятий производная, метод касательных, программирование, нелинейное уравнение. Объект исследования Корни нелинейного уравнения.

Цель работы Определение корней нелинейного уравнения. Методы исследования изучение работ отечественных и зарубежных авторов по данной теме. Полученные результаты изучен метод касательных решения нелинейных уравнений рассмотрена возможность составления программы на языке программирования Турбо-Паскаль 7.0 Область применения в работе инженера. СОДЕРЖАНИЕ стр.

ВВЕДЕНИЕ 1. Краткое описание сущности метода касательных метода секущих Ньютона 2. Решение нелинейного уравнения аналитически 3. Блок схема программы 4. Программа на языке PASCAL 5. Результаты выполнения программы 13 СПИСОК ИСПОЛЬЗОВАННИХ ИСТОЧНИКОВ 14 ВВЕДЕНИЕ Процедура подготовки и решения задачи на

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

Подготовка исходных данных . 6. Ввод программы и исходных данных в ЭВМ. 7. Отладка программы. 8. Тестирование программы. 9. Решение задачи на ЭВМ и обработка результатов. В настоящей курсовой работе условие задачи дано в математической формулировке, поэтому необходимость в выполнении этапов 1 и 2 отпадает и сразу можно приступить к разработке алгоритма решения задачи на

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

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

между ними. Внутри блоков указывается информация, характеризующая выполняемые ими функции. Блоки схемы имеют сквозную нумерацию. Конфигурация и размеры блоков, а также порядок построения схем определяются ГОСТ 19.002-80 и ГОСТ 19.003-80. На этапе 4 составляется программа на языке Турбо-Паскаль. При описании программы необходимо использовать характерные приемы программирования и учитывать специфику языка. В качестве языка программирования выбран язык

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

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

Задание при обработке на ЭВМ проходит ряд шагов компиляцию, редактирование компоновку и выполнение. Обработка результатов решения задачи осуществляется с помощью ЭВМ. Выводимые результаты оформлены в виде, удобном для восприятия. 1. Краткое описание сущности метода касательных метода секущих Ньютона Пусть на отрезке a b отделен корень с уравнения f x 0 и f -функция непрерывна на отрезке a b,

а на интервале a b существуют отличные от нуля производные f и f . Так как f x 0 , то запишем уравнение f x 0 в виде x x f x f x 1 Решая его методом итераций можем записать xn1 x n f x n f x n 2 Если на отрезке ab f x f x 0, то нул евое приближение выбираем x0a. Рассмотрим геометрический смысл метода . Рассмотрим график функции yfx.

Пусть для определенности f x 0 и f x 0 рис. 1. Проведем касательную к графику функции в точке B b, f b. Ее уравнение будет иметь вид y f b f b x b Полагая в уравнении y 0 и учитывая что f x 0, решаем его относительно x. Получим x b f b f b Нашли абсциссу x1 точки c1 пересечения касательной с осью ox x1 b f b f b Проведем касательную к графику функции в точке b1 x1 f x1.Найдем абсциссу x2 точки с2 пересечения касательной

с осью Ox x2 x1 f x1 f x1 Вообще xk1 x k f x k f x k 3 Таким образом, формула 3 дает последовательные приближения xk корня, получаемые из уравнения касательной , проведенной к графику функции в точке b k x k f x k0 метод уточнения корня c ab уравнения f x 0 с помощью формулы 3 называется методом касательной или методом Ньютона. Геометрический смысл метода касательных состоит в замене дуги y f x касательной, одной к одной

из крайних точек . Начальное приближение x 0 a или x0 b брать таким, чтобы вся последовательность приближения х k принадлежала интервалу ab . В случае существования производных f , f , сохраняющих свои знаки в интервале, за х0 берется тот конец отрезка ab, для которого выполняется условие f х0 f х0 0. Для оценки приближения используется общая формула c-x k-1 f x k1m , где m min f x на отрезке ab . На практике проще пользоваться другим правилом Если на отрезке ab выполняется условие 0 m f x и e -

заданная точность решения, то неравенство x k1-x k e влечет выполнение неравенства c-x k-1 e . В этом случае процесс последовательного приближения продолжают до тех пор, пока не выполнится неравенство c-x k-1 e . 2. Решение нелинейного уравнения аналитически Определим корни уравнения х3 0,1х2 0,4х 1,2 0 аналитически. Находим f x х3 0,1х2 0,4х 1,2 f x 3х2 0,1х 0,4 f 1 2,5 0 f 0 1,2 0 f 1 0,3 0 x- -101 sign f x

Следовательно, уравнение имеет действительный корень, лежащий в промежутке 0 1 . Приведем уравнение к виду x j x , так , чтобы j x 1 при 0 x 1. Так как max f x f 1 3 0,1 0,4 3,5 то можно взять R 2. Тогда j x x f x R x 0,5 х3 0,05 х2 0,2 х 0,6 0,5 х3 0,05 х2 0,8 х 0,6. Пусть х0 0 , тогда х n1 j х n. Вычисления расположим в таблице. nхnх2nх3nj хn.f x11110,85-0,1736320,850,72250,6141250,93 681250,0846530,93681250,877617660,822163 1940,89448752-0,0465140,894487520,800107 9230,7156865520,9177413440,02428850,9177 413440,8422491740,7729668890,905597172-0

,0130660,9055971720,8201062380,742685890 ,9121294810,00692370,9121294810,83198019 0,7588736590,908667746-0,003780,90866774 60,8256770720,7502661240,9105172810,0019 6890,9105172810,8290417190,7548568120,90 953-0,00105100,909530,8272508840,7524122 530,9100579950,000559110,9100579950,8282 055550,7537150870,909778575-0,0003120,90 97785750,8276970550,7530210480,909927483 0,000159130,9099274830,8279680250,753390 8610,909848155-8,5E-05140,9098481550,827 8236650,7531938340,9098904244,5E-05150,9 098904240,8279005830,7532988120,90986790 4-2,4E-05160,9098679040,8278596020,75324 28810,9098799021,28E-05170,9098799020,82 78814370,7532726810,90987351-6,8E-06180, 909873510,8278698030,7532568040,90987691 63,63E-06190,9098769160,8278760020,75326 52630,909875101-1,9E-06200,9098751010,82 78726990,7532607560,9098760681,03E-06 График функции y х3 0,1х2 0,4х 1,2 3. Блок схема программы 4. Программа на языке PASCAL 7.0 program metodkasatelНазвание программы uses Crt Модуль дисплейных функций var Блок описаний переменных xn,xn1,a,b,c,mx,y0,x0 real function f1x1Real Real Основная функция begin f1 x1x1x1-0.5-0.05x1x10.8x10.6 end function f2x4Real

Real Производная от основной функции begin f2 x4x4x40.5x4x40.1x4x40.4x4 1.2 end begin Начало основного тела программы Clrscr Очистка экрана перед выполнением программы a0b1c0.01 Writeln От A,a, до B,b Вывод на экран Writeln Погрешность с,c Readln Ожидание нажатия клавиши Enter xnb xn1 f1xn y0f2b while ABSy0 c do Проверка по точности вычисления корня begin

Тело цикла xnxn1 xn1f1xn y0 f2xn1 Печать промежуточного результата Writelnxn,xn, xn1,xn1, fxn1,y0 Readln Ожидание нажатия клавиши Enter end Конец тела цикла WritelnКонечные значения Печать полученного результата Writeln xn1,xn1, fxn1,y0 Readln Ожидание нажатия клавиши Enter end. Конец основного тела программы5.

Результаты выполнения программы От A 0.0E00 до B 1.0E00 Погрешность с 1.0E-08 От A 0.0E00 до B 1.0E00 Погрешность с 1.0E-08 xn 8.50E-01 xn1 9.3681250000E-01 fxn1 8.4649960270E-02 xn 9.3681250000E-01 xn1 8.9448751986E-01 fxn1-4.6507647892E-02 xn 8.9448751986E-01 xn1 9.1774134381E-01 fxn1 2.4288343840E-02 xn 9.1774134381E-01 xn1 9.0559717189E-01 fxn1-1.3064617920E-02 xn 9.0559717189E-01 xn1 9.1212948085E-01 fxn1 6.9234699658E-03 xn 9.1212948085E-01 xn1 9.0866774587E-01

fxn1-3.6990702320E-03 xn 9.0866774587E-01 xn1 9.1051728099E-01 fxn1 1.9678960780E-03 xn 9.1051728099E-01 xn1 9.0953333295E-01 fxn1-1.0493249720E-03 xn 9.0953333295E-01 xn1 9.1005799543E-01 fxn1 5.5884091853E-04 xn 9.1005799543E-01 xn1 9.0977857497E-01 fxn1-2.9781681224E-04 xn 9.0977857497E-01 xn1 9.0992748338E-01 fxn1 1.5865717614E-04 xn 9.0992748338E-01 xn1 9.0984815480E-01 fxn1-8.4537703515E-05 xn 9.0984815480E-01 xn1 9.0989042365E-01 fxn1 4.5040009354E-05 xn 9.0989042365E-01 xn1 9.0986790364E-01 fxn1-2.3997676180E-05

xn 9.0986790364E-01 xn1 9.0987990248E-01 fxn1 1.2785800209E-05 xn 9.0987990248E-01 xn1 9.0987350958E-01 fxn1-6.8122881203E-06 xn 9.0987350958E-01 xn1 9.0987691573E-01 fxn1 3.6295678001E-06 xn 9.0987691573E-01 xn1 9.0987510095E-01 fxn1-1.9338276616E-06 xn 9.0987510095E-01 xn1 9.0987606786E-01 fxn1 1.0303429008E-06 xn 9.0987606786E-01 xn1 9.0987555269E-01 fxn1-5.4896190704E-07 xn 9.0987555269E-01 xn1 9.0987582717E-01 fxn1 2.9248803912E-07 xn 9.0987582717E-01 xn1 9.0987568093E-01 fxn1-1.5583464119E-07 xn 9.0987568093E-01

xn1 9.0987575885E-01 fxn1 8.3031409304E-08 xn 9.0987575885E-01 xn1 9.0987571733E-01 fxn1-4.4236003305E-08 xn 9.0987571733E-01 xn1 9.0987573945E-01 fxn1 2.3572283681E-08 xn 9.0987573945E-01 xn1 9.0987572766E-01 fxn1-1.2558302842E-08 xn 9.0987572766E-01 xn1 9.0987573394E-01 fxn1 6.6920620156E-09 Конечные значения xn1 9.0987573394E-01 fxn1 6.6920620156E-09 СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ 1. Алексеев В. Е Ваулин

А.С Петрова Г. Б. Вычислительная техника и программирование. Практикум по программированию Практ .пособие М. Высш. шк 1991. 400 с. 2. Абрамов С.А Зима Е.В. Начала программирования на языке Паскаль. М. Наука, 1987. 112 с. 3. Вычислительная техника и программирование Учеб. для техн. вузов А.В. Петров, В.Е. Алексеев, А.

С. Ваулин и др. М. Высш. шк 1990 479 с. 4. Гусев В.А Мордкович А.Г. Математика Справ. материалы Кн. для учащихся. 2-е изд. М. Просвещение, 1990. 416 с. 5. Марченко А.И Марченко Л.А. Программирование в среде Turbo Pascal 7.0 К. ВЕК, М. Бином Универсал, 1998. 496 с.