• Введение
• Глава I. Исследование метода итераций
• 1.1 Нелинейные уравнения
• 1.2 Метод итераций
• 1.3 Геометрический смысл
• 1.4 Решение нелинейного уравнения методом итераций
• 1.5 Экономическое применение
• Глава II. Разработка программы
• 2.1 Автоматизация метода
• Заключение
• Список использованных источников и литературы
• Приложения

Введение

Если законы функционирования модели нелинейны, а моделируемые процесс или система обладают одной степенью свободы (т.е. имеют одну независимую переменную), то такая модель, как правило, описывается одним нелинейным уравнением.

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

Методы решения квадратных уравнений были известны еще древним грекам. Решение уравнений третьей и четвертой степеней были получены усилиями итальянских математиков Ш. Ферро, Н. Тартальи, Дж. Картано, Л. Феррари в эпоху Возрождения. Затем наступила пора поиска формул для нахождения корней уравнений пятой и более высоких степеней. Настойчивые, но безрезультатные попытки продолжались около 300 лет и завершились благодаря работам норвежского математика Н. Абеля. Он доказал, что общее уравне6ие пятой и более высоких степеней неразрешимы в радикалах. Решение общего уравнения n-ой степени

a₀xⁿ+a₁x^n-1+…+a_n-1x+a_n=0, a₀0

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

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

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

Объект: нелинейные уравнения.

Цели: исследовать метод итераций, автоматизировать его с применением в среде Pascal.

Достижение поставленных целей потребовало решение следующих задач:

· Изучить метод итераций;

· Исследовать его применение;

· Разработать вычислительный алгоритм метода итераций;

· Составить программы решения систем нелинейных уравнений методом итераций в среде Turbo Pascal.

Глава I. Исследование метода итераций

итерация нелинейный уравнение алгоритм

Одной из важнейших и наиболее распространённых задач прикладной математики является задача решения нелинейных уравнений, встречающихся в разных областях научных исследований. Любое уравнение в общем случае можно представить в виде

f ( x ) = 0.

Нелинейные уравнения можно разделить на два класса - алгебраические и трансцендентные.

Алгебраическими уравнениями называются уравнения, содержащие только алгебраические функции (целые, рациональные, иррациональные). Алгебраическое уравнение в общем виде можно представить многочленом n-ой степени с действительными коэффициентами:

f (x) = а₀x_n + а₁х_n-1 +... + а_n =0.

Например, х³ + х² + 2х = 0.

Трансцендентными называются уравнения, содержащие другие функции (тригонометрические, показательные, логарифмические и т.д.), например: 2x-sin x = 0.

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

Однако точное решение уравнения не всегда является необходимым. Задачу отыскания корней уравнения можно считать практически решенной, если мы сумеем найти корни уравнения с заданной степенью точности. Для этого используются приближенные (численные) методы решения.

Таким образом, уравнение типа или называется нелинейным. Решить уравнение - это значит найти такое x, при котором уравнение превращается в тождество. В общем случае уравнение может иметь 0; 1; 2;...? корней. Рассмотрим нахождение корня нелинейного уравнения с помощью метода итераций на заданном интервале [a,b].

1.2 Метод итераций

Уравнение представим в виде: .

Далее на отрезке [a,b], где функция имеет корень, выбирается произвольная точка x₀ и далее последовательно вычисляется:

(1)

Если на отрезке [a,b] выполнено условие |цґ(x)| ? q <1, то итерационный процесс сходится к корню уравнения .

Если необходимо вычислить корень с точностью е, то процесс итераций продолжается до тех пор, пока для двух последовательных приближений x_n и x_n-1 не будет выполнено:

, при этом всегда выполняется , где е задается погрешностью корня x*.

Если q ?0.5, то можно пользоваться соотношением .

Процесс определения интервала изоляции [a,b], содержащего только один из корней уравнения, называется отделением этого корня.

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

Процесс определения корней алгебраических и трансцендентных уравнений состоит из 2 этапов:

· отделение корней, - т.е. определение интервалов изоляции [a,b], внутри которого лежит каждый корень уравнения;

· уточнение корней, - т.е. сужение интервала [a,b] до величины равной заданной степени точности е.

1.3 Геометрический смысл

Будем предполагать, что функции ц(x) и j(x) являются непрерывными. На плоскости X0Y построим графики функции Y=x и Y=ц(x). Каждый вещественный корень x* уравнения является абсциссой точки пересечения кривой Y= ц(x) с прямой Y=x. Начиная с некоторой точки A₀(x₀, ц(x₀)), строим ломаные линии A₀B₁A₁B₂A₂…(лестница), звенья которой попеременно параллельны оси 0X и оси 0Y, причем вершины A₀,A₁,A₂… лежат на кривой Y=ц(x). Общие абсциссы точек A₁ и B₁, A₂ и B₂ … представляют собой последовательные приближения x₁, х₂,…,х_k,… корня x*, которые сходятся к нему монотонно и односторонне.

Рис. 1.

На рис. 1 представлен случай, когда 0<j(x)<1, т.е. угол касательной к графику функции Y=ц(x) меньше 45⁰, т.е. a < 45⁰. Функция Y=ц(x) является возрастающей и вогнутой. Если -1<j(x)<0, то ломаная A₀B₁A₁B₂A₂… будет иметь вид спирали (рис.2). В этом случае сходимость является двусторонней, т.е. искомый корень всегда принадлежит отрезку [x_k,x_k+1]. Функция Y=ц(x) является убывающей и вогнутой.

Рис. 2.

Если |j(x)|>1, т.е. угол наклона касательной к кривой ц(x) превышает 45⁰, то в этом случае итерации сходиться не будут (рис.3).

Рис. 3.

Если же |j(x)|<1 в некоторой окрестности корня, а вдали от него это неравенство не выполняется, то итерационный процесс будет сходящимся только в том случае, если начальное приближение x₀ выбрано достаточно близко к корню (рис. 4).

Рис. 4.

При произвольном выборе начального приближения сходимости может не быть (рис. 5).

Рис. 5.

1.4 Решение нелинейного уравнения методом итераций

Рассмотрим уравнение у = ln(x) - x + 1,8. Представим его в виде x = ln(x) + 1,8.

Проверим условие сходимости, найдя производную от функции f(x) и подставив в получившееся выражения концы отрезка [2,3].

f '(x) = (ln(x) + 1,8)' = 1/x;

f '(x) = 1/2 = 0,5;

f '(x) = 1/3 = 0,3333333;

Как видим, условие сходимости итерационного процесса выполняется, т.е. f '(x) <1.

Из отрезка, на котором определен корень уравнения, выбираем произвольную точку x_o = 2. Затем последовательно вычисляем x₁, x₂, x_k, подставляя x_o в формулу (1). Получим следующее:

По полученным данным построим график функции x = ln(x) + 1,8

Рис. 6. График функции

Далее сравниваем с точностью разность между значениями F(x _k+1 ) и F(x _k ) функции до тех пор, пока разность не будет меньше точности. При выполнении условия выводим корень уравнения F(x _k+1 ).

2.1 Автоматизация метода

Программа на Pascal состоит из отдельных разделов или блоков, которые должны располагаться в следующем порядке:

- заголовок программы;

- раздел объявления переменных;

- раздел объявления процедур и функций;

- тело программы.

Заголовок программы состоит из зарезервированного слова PROGRAM и имени программы: PROGRAM MetodIteraccii;

Приложение 1

Кодовая часть.

Блок-схема метода итераций.

Начало

x _n, E

x_n+1 = ln (x_n) + 1,8

- f `(x)| < 1 +

`Не удовлетворяет-

условию сходимости' |x_n+1 - x_n| ? E x _n:= x_n+1

+

x _n+1 Конец

Приложение 3

1) y = cos(x) - x + 5

2) y = 0.7 sin(x) -x + 2.7

Размещено на Allbest.ru