Методы решения нелинейных уравнений

МИНИСТЕРСТВО  ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ 

Международная «Лига развития науки  и образования» (Россия)

Международная ассоциация развития науки, образования  и культуры России (Италия) 

Международный «ИНСТИТУТ УПРАВЛЕНИЯ»

(г. Архангельск) 
 

КУРСОВАЯ  РАБОТА

ПО ДИСЦИПЛИНЕ

«Информатика  и программирование»

по теме: «Методы решения нелинейных уравнений» 
 
 
 

 

Архангельск

2005

 

Аннотация

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

Содержание

Введение              4

1. Метод  решения нелинейного уравнения  методом секущих                7

1.1. Общая характеристика методов решения нелинейных уравнений     7

1.2. Метод  секущих                    10

1.3. Тестовый  пример                             12

1.4. Разработка  алгоритма решения нелинейных  уравнений    15

2. Вычисление  значения функции при помощи  интерполяционной формулы  16

2.1. Общая  характеристика методов интерполяционной  функции    16

2.2. Интерполяционная  формула Ньютона       24

2.3. Тестовый  пример           25

2.4. Разработка  алгоритма и программы вычисления функции    26

Заключение            27

Список  литературы           28

Приложение 1.1            29

Приложение 2.1            31

  Введение

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

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

1. Численные методы линейной алгебры

  Метод Гаусса для решения системы линейных алгебраических уравнений. Устойчивость метода Гаусса. Использование метода Гаусса для вычисление обратной матрицы. Метод квадратного корня.

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

  Одношаговые итерационные методы решения систем линейных алгебраических уравнений. Каноническая форма записи. Примеры одношаговых итерационных методов. Достаточное условие сходимости.

  Необходимое и достаточное условие сходимости одношаговых стационарных итерационных методов. Теорема о сходимости одношаговых  стационарных итерационных методов. Оценка скорости сходимости. Неявный итерационный метод с чебышевским набором параметров. Оценка скорости сходимости.

  Численная устойчивость итерационного метода с чебышевским набором параметров. Упорядоченный набор итерационных параметров (пример).

  Одношаговые итерационные методы вариационного  типа. Формула для вычисления итерационного  параметра.

  Примеры итерационных методов вариационного  типа (метод скорейшего спуска; метол  минимальных невязок; метод минимальных  поправок; метод минимальных погрешностей) . Каноническая форма записи двухшаговых итерационных методов вариационного типа.

  Примеры двухшаговых итерационных методов (метод сопряженных градиентов, сопряженных  невязок, сопряженных поправок, сопряженных  погрешностей).

  Полная  и частичная проблема собственных значений.

  Степенной метод решения частичной проблемы собственных значений. Решение полной проблемы собственных значений методом  вращений. Метод обратной итерации.

2. Решение нелинейных уравнений и систем уравнений.

  Решение нелинейных уравнений. Методы разделения корней. Примеры численных методов решения нелинейных уравнений (метод простой итерации, метод Ньютона, модифицированный метод Ньютона, метод секущих).

  Сходимость  метода простой итерации.

  Метод Эйткена ускорения сходимости.

  Сходимость  метода Ньютона.

  Решение систем нелинейных уравнений.

  Примеры (применение метода простой итерации; сравнение скорости сходимости метода простой итерации и метода Ньютона; применение метода Ньютона для решения  системы двух нелинейных уравнений). 
 

  

3. Интерполяция  и приближение  функций.

  Постановка  задачи интерполирования алгебраическими  многочленами. Интерполяционная формула Лагранжа. Интерполяционная формула Ньютона (разделенные разности, схема Горнера).

  Интерполирование  с кратными узлами (существование и единственность многочлена Эрмита, погрешность интерполирования с кратными узлами). Пример (многочлен Эрмита третьей степени). Сходимость интерполяционного процесса.

  Интерполирование  сплайнами. Кубический сплайн. Наилучшее  приближение функции, заданной таблично (пример). Наилучшее приближение в гильбертовом пространстве. 
 
 
 
 
 
 
 
 
 
 

  Глава 1. Метод решения нелинейного уравнения методом секущих

  1.1. Общая характеристика методов

  1. Условие задания

  При заданных пяти вариантах допустимой ошибки e заданным численным методом вычислить приближенное значение корня функционального уравнения вида f (x) = 0, если известно, что это уравнение имеет единственный корень на отрезке [a, b].

  В проекте должно быть предусмотрено:

  - построение графика функции  f (x) на отрезке [a, b],

  - проверка корректности введенных  значений исходных данных (выполнение  условия a < b, выполнение условия e > 0),

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

  2. Содержание пояснительной записки

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

    условие задачи;     

    условия заданного варианта задания;      

    описание заданного численного  метода;     

    блок-схему алгоритма подзадачи  вычисления корня;     

    программу процедуры вычисления  корня;     

    главную программу;     

    результаты вычислений значения корня для заданных пяти вариантов допустимой ошибки

  Варианты  численного метода:

  1) метод простых итераций,

  2) метод Ньютона,

  3) метод проб,

  4) метод секущих,

  5) метод хорд.

  Упрощенный  метод Ньютона:   , n=0,1,…

  Метод ложного положения:   ,   n=0,1,…;

  c - фиксированная точка из окрестности  корня

  Метод секущих:  , n=0,1,…

  Метод Стеффенсена:  , n=0,1,… 
 
 

3. Описание Метод секущих

  Метод секущих, так же, как и метод  проб, использует не одно, а два начальных  приближения, которые мы обозначим  соответственно x_n1 и x_n2. Перед выполнением первой итерации воспользуемся правилом (6) для определения значений этих приближений.

  При выполнении каждой очередной итерации для вычисления следующего приближения по методу хорд проведем прямую линию (секущую) MN через точки с координатами (x_n1, f(x_n1)) и (x_n2, f (x_n2)), а абсциссу точки пересечения секущей MN с осью х возьмем в качестве значения следующего приближения x_s к корню (рис. 3).

  Принятое  правило нахождения следующего приближения  приводит к расчетной формуле:

  Из  трех приближений к корню оставим  два последних (отбрасываем самое старое x_n1). В методе секущих это делается по следующему правилу:

  x_n1 = x_n2;  x_n2 = x_s.                           

  Выполнение  итераций можно прекратить при выполнении условия

  |x_n2 - x_n1|< e,                                   

  а полученное значение приближения  x_s взять в качестве искомого значения корня x_w.

  Расчетные формулы должны быть применены в  алгоритме вычисления корня по методу секущих.

  Обратим внимание на то, что формула  имеет  много общего с формулой Ньютона. Знаменатель в формуле  есть не что иное, как среднее значение производной f^‘(x) на отрезке [x_n1, x_n2].  

1.2. Метод секущих

  Пусть на отрезке [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)

  Решая его методом итераций можем записать :

  xn+1 = x n– ( f(x n)/f ’(x n)) (2)

  Если  на отрезке [a;b] f ’(x) * f “(x) > 0, то нул  – евое приближение выбираем x0=a. Рассмотрим геометрический смысл метода . Рассмотрим график функции y=f(x). Пусть для определенности 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))