WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


Pages:     || 2 | 3 |
ЧИСЛЕННЫЕ МЕТОДЫ АНАЛИЗА ИЗДАТЕЛЬСТВО ГОУ ВПО ТГТУ Министерство образования и науки Российской Федерации Государственное образовательное учреждение высшего профессионального образования «Тамбовский государственный технический университет» ЧИСЛЕННЫЕ МЕТОДЫ АНАЛИЗА Методические указания для студентов специальности 240802 Тамбов Издательство ГОУ ВПО ТГТУ 2010 УДК 66.011(076.5) ББК В192.1я73-5 Ч-671 Рецензент Доктор технических наук, профессор ГОУ ВПО ТГТУ Е.Н. Туголуков Составители:

А.Н. Пахомов, Ю.В. Пахомова Ч-671 Численные методы анализа : методические указания / сост. :

А.Н. Пахомов, Ю.В. Пахомова. – Тамбов : Изд-во ГОУ ВПО ТГТУ, 2010. – 32 с. – 50 экз.

Даны методические указания к лабораторным работам по курсу «Основы математической физики и численные методы анализа»» для студентов специальности 240802.

УДК 66.011(076.5) ББК В192.1я73-5 © Государственное образовательное учреждение высшего профессионального образования «Тамбовский государственный технический университет» (ГОУ ВПО ТГТУ), 2010 Учебное издание ЧИСЛЕННЫЕ МЕТОДЫ АНАЛИЗА Методические указания Составители:

ПАХОМОВ Андрей Николаевич, ПАХОМОВА Юлия Владимировна Редактор Е.С. Кузнецова Инженер по компьютерному макетированию Т.Ю. Зотова Подписано в печать 02.10.2010 Формат 60 84/16. 1,86 усл. печ. л. Тираж 50 экз. Заказ № 529 Издательско-полиграфический центр ГОУ ВПО ТГТУ 392000, Тамбов, ул. Советская, 106, к. 14 ВВЕДЕНИЕ Современному инженеру в своей практике приходится решать задачи расчёта и моделирования технологического процесса и оборудования. Для успешного решения подобных задач необходимо уметь применять численные и аналитические методы расчёта уравнений и их систем. Выбор метода решения в значительной мере определяет быстроту и точность получаемого решения.

Выбор подходящего метода для решения уравнений зависит от характера рассматриваемой задачи. Задачи, сводящиеся к решению отдельных уравнений и их систем, можно классифицировать по числу уравнений и в зависимости от предлагаемого характера и числа решений. Одно уравнение называется линейным, алгебраическим или трансцендентным в зависимости от того, имеет оно одно решение, n решений или неопределённое число решений. Систему уравнений будем называть линейной или нелинейной в зависимости от математической природы входящих в неё уравнений [2].

Решение линейного уравнения с одним неизвестным получается достаточно просто (см. школьный курс математики).

РЕШЕНИЕ НЕЛИНЕЙНЫХ УРАВНЕНИЙ Нелинейные уравнения можно разделить на два класса – алгебраические и трансцендентные.

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

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

cosx = x, logx = x – 5, x3 = logx + x5 + 40.

Более строгое определение таково: Трансцендентное уравнение – это уравнение вида f(x) = g(x), где функции f и g являются аналитическими функциями, и по крайней мере одна из них не является алгебраической.

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

1) точные методы;

2) итерационные методы.

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

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

Пусть дано уравнение f(x) = 0, (1) где:

1) функция f(x) непрерывна на отрезке [a, b] вместе со своими производными первого и второго порядка;

2) значения f(x) на концах отрезка имеют разные знаки (f(a) f(b) < 0);

3) первая и вторая производные f (x) и f (x) сохраняют определённый знак на всём отрезке.

Условия 1 и 2 гарантируют, что на интервале [a, b] находится хотя бы один корень, а из 3 следует, что f(x) на данном интервале монотонна, и поэтому корень будет единственным.

Решить уравнение (1) итерационным методом – значит установить, имеет ли оно корни, сколько корней, и найти значения корней с нужной точностью.

Всякое значение, обращающее функцию f(x) в нуль, т.е. такое, что f() = 0, называется корнем уравнения (1) или нулём функции f(x).

Задача нахождения корня уравнения f(x) = 0 итерационным методом состоит из двух этапов:

1) определение отрезка локализации корней – отыскание приближённого значения корня или содержащего его отрезка;

2) уточнение приближённых корней – доведение их до заданной степени точности.

Процесс определения отрезка локализации корней начинается с установления знаков функции f(x) в границах x = a и x = b области её существования.

Пример 1. Определить отрезок локализации корней уравнения:

х3 – 7х + 3 = 0. (2) Построим график этой функции (рис. 1).



Следовательно, уравнение (2) имеет три действительных корня, лежащих в интервалах [–3, –1], [0, 1] и [1, 5].

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

В инженерной практике распространён графический способ определения приближённых корней.

f(x) –––––5 x Рис. 1. Определение отрезка локализации Принимая во внимание, что действительные корни уравнения (1) – это точки пересечения графика функции f(x) с осью абсцисс, достаточно построить график функции f(x) и отметить точки пересечения f(x) с осью Ох или отметить на оси Ох отрезки, содержащие по одному корню. Построение графиков часто удаётся сильно упростить, заменив уравнение (1) равносильным ему уравнением:

f1(x) = f2(x), (3) где функции f1(x) и f2(x) более простые, чем функция f(x). Тогда, построив графики функций у = f1(x) и у = f2(x), искомые корни получим как абсциссы точек пересечения этих графиков.

Пример 2. Графически отделить корни уравнения (рис. 2):

Рис. 2. Графическое определение корня x lg x = 1. (4) Уравнение (4) удобно переписать в виде равенства:

lg x = 1/ x Отсюда ясно, что корни уравнения (4) могут быть найдены как абсциссы точек пересечения логарифмической кривой y = lg x и гиперболы y =1/x. Построив эти кривые, приближённо найдём единственный корень = 2,5 уравнения (4) или определим его содержащий отрезок [2, 3].

Итерационный процесс состоит в последовательном уточнении начального приближения х0. Каждый такой шаг называется итерацией. В результате итераций находится последовательность приближённых значений корня х1, х2,..., хn. Если эти значения с увеличением числа итераций n приближаются к истинному значению корня, то говорят, что итерационный процесс сходится.

Метод половинного деления a + b Для нахождения корня уравнения (1), принадлежащего отрезку [a, b], делим этот отрезок пополам. Если f = 0, a + b a + b a + b a + b a, то = 0 является корнем уравнения. Если f 0, то выбираем ту из половин или, b, на 2 2 2 концах которой функция f(x) имеет противоположные знаки. Новый суженный отрезок [а1, b1] снова делим пополам и производим те же самые действия.

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

Пример 3. Методом половинного деления уточнить корень уравнения f(x) x4 + 2x3 – x – 1 = 0, лежащий на отрезке [0, 1].

Последовательно имеем:

f(0) = – 1; f(1) = 1; f(0,5) = 0,06 + 0,25 – 0,5 – 1 = – 1,19;

f(0,75) = 0,32 + 0,84 – 0,75 – 1 = – 0,59;

f(0,875) = 0,59 + 1,34 – 0,88 – 1 = + 0,05;

f(0,8125) = 0,436 + 1,072 – 0,812 – 1 = – 0,304;

f(0,8438) = 0,507 + 1,202 – 0,844 – 1 = – 0,135;

f(0,8594) = 0,546 + 1,270 – 0,859 – 1 = – 0,043 и т.д.

Можно принять = (0,859 + 0,875) = 0,867.

Метод хорд В данном методе процесс итераций состоит в том, что в качестве приближений к корню уравнения (1) принимаются значения х1, х2,..., хn точек пересечения хорды АВ с осью абсцисс (рис. 3). Сначала запишем уравнение хорды AB:

y - f (a) x - a =.

f (b) - f (a) b - a Для точки пересечения хорды AB с осью абсцисс (х = х1, y = 0) получим уравнение:

f (a) x1 = (b - a).

f (b) - f (a) Пусть для определённости f (x) > 0 при а х b (случай f (x) < 0 водится к нашему, если записать уравнение в виде f(x) = 0). Тогда кривая = f(x) будет выпукла вниз и, следовательно, расположена ниже своей хорды АВ. Возможны два случая: 1) f(а) > 0 (рис. 3, а) и 2) f(b) < 0 (рис. 3, б).

а) б) а) б) Рис. 3. Метод хорд В первом случае конец а неподвижен и последовательные приближения:

x0 = b;

f (xi ) xi+1 = xi i f (xi )- f (a)(x - a), (i = 0, 1, 2,...) (5) образуют ограниченную монотонно убывающую последовательность, причём a < <... < xi +1 < xi <... < x1 < x0.

Во втором случае неподвижен конец b, а последовательные приближения:

x0 = а;

f (xi ) xi+1 = xi f (b)- f (xi )(b - xi ) (6) образуют ограниченную монотонно возрастающую последовательность, причём x0 < x1 <... < xi < xi+1 <... < < b.

Обобщая эти результаты, заключаем:

1) неподвижен тот конец, для которого знак функции f (х) совпадает со знаком её второй производной f (х);

2) последовательные приближения xn лежат по ту сторону корня, где функция f (х) имеет знак, противоположный знаку её второй производной f (х).

Итерационный процесс продолжается до тех пор, пока не будет обнаружено, что xi – xi – 1 <, где – заданная предельная абсолютная погрешность.

Пример 4. Найти положительный корень уравнения f (x) x3 – 0,2x2 – 0,2х – 1,2 = с точностью = 0,01.

Прежде всего определяем отрезок локализации корня.

Так как f (1) = –0,6 < 0 и f (2) = 5,6 > 0, то искомый корень лежит в интервале [1, 2]. Полученный интервал велик, поэтому разделим его пополам.

Так как f (1,5) = 1,425 > 0, то 1< < 1,5.

Так как f (x) = 6x – 0,4 > 0 при 1 < х < 1,5 и f (1,5) > 0, то воспользуемся формулой (5) для решения поставленной задачи:

0,x1 = 1+ (1,5 -1) = 1,15 ; x1 – x0 = 0,15 >, 1,425 + 0,следовательно, продолжаем вычисления;

f (х1) = –0,173;

0,x2 = 1,15 + (1,5 -1,15)= 1,190 ; x2 – x1 = 0,04 >, 1,425 + 0,f (х2) = –0,036;

0,x3 = 1,190 + (1,5 -1,190)= 1,198 ; x3 – x2 = 0,008 <.

1,425 + 0,Таким образом, можно принять = 1,198 с точностью = 0,01.





Точный корень уравнения = 1,2.

Метод Ньютона Отличие этого итерационного метода от предыдущего состоит в том, что вместо хорды на каждом шаге проводится касательная к кривой y = f(x) при = хi и ищется точка пересечения касательной с осью абсцисс (рис. 4). При этом не обязательно задавать отрезок [а, b], содержащий корень уравнения (1), достаточно найти лишь некоторое начальное приближение корня x = х0.

Рис. 4. Метод Ньютона Применяя метод Ньютона, следует руководствоваться следующим правилом: в качестве исходной точки х0 выбирается тот конец интервала [а, b], которому отвечает ордината того же знака, что и знак f (х).

Уравнение касательной, проведённой к кривой y = f(x) через точку В0 с координатами х0 и f(х0), имеет вид y - f (x0) = f (x0)(x - x0).

Отсюда найдём следующее приближение корня х1 как абсциссу точки пересечения касательной с осью Ох (y = 0) f (x0) x1 = x0 -.

f (x0) Аналогично могут быть найдены и следующие приближения как точки пересечения с осью абсцисс касательных, проведённых в точках В1, В2 т.д. Формула для i + 1 приближения имеет вид f (xi ) xi+1 = xi -. (7) f (xi ) Для окончания итерационного процесса может быть использовано или условие f(xi) <, или условие близости двух последовательных приближений xi – xi – 1 <.

Итерационный процесс сходится, если f(х0) f (х0) > 0.

Метод последовательных приближений Для использования метода последовательных приближений исходное нелинейное уравнение f(х) = 0 заменяется равносильным уравнением x = (x). (8) Пусть известно начальное приближение корня х = х0. Подставляя это значение в правую часть уравнения (8), получим новое приближение:

х1 = (х0).

Далее, подставляя каждый раз новое значение корня в (8), получаем последовательность значений:

xi + 1 = (xi), (i = 0, 1, …, n). (9) Геометрически метод итерации может быть пояснён следующим образом. Построим на плоскости хОу графики функций у = х и у = (х). Каждый действительный корень уравнения (8) является абсциссой точки пересечения М кривой у = (х) с прямой у = х (рис. 5, а).

а) б) Рис. 5. Сходящиеся итерационные процессы Отправляясь от некоторой точки А0 [x0, (x0)], строим ломаную А0В1А1В2А2 … («лестница»), звенья которой попеременно параллельны оси Ох и оси Оу, вершины А0, А1, А2, … лежат на кривой у = (х), а вершины В1, В2, В3, …, – на прямой у = х.

Общие абсциссы точек А1 и В1, А2 и В2, …, очевидно, представляют собой соответственно последовательные приближения х1, х2, … корня.

Возможен также другой вид ломаной А0В1А1В2А2 … – «спираль» рис. 5, б). Решение в виде «лестницы» получается, если производная (х) положительна, а решение в виде «спирали», если (х) отрицательна.

На рисунке 5, а, б кривая у = (х) в окрестности корня пологая, т.е. (x) < 1, и процесс итерации сходится. Однако, если рассмотреть случай, где (x) > 1, то процесс итерации может быть расходящимся (рис. 6). оэтому для практического применения метода итерации нужно выполнение достаточного условия сходимости итерационного процесса.

Пусть функция (х) определена и дифференцируема на отрезке [a, b], причём все её значения (х) [a, b].

Рис. 6. Расходящийся итерационный процесс Тогда, если существует правильная дробь q такая, что (x) q < 1 при a < x < b, то:

1) процесс итерации xi + 1 = (xi), (i = 0, 1, …, n) сходится независимо от начального значения х0 [a, b];

2) предельное значение = lim xn является единственным корнем уравнения х = (х) на отрезке [a, b].

n Пример 5. Уравнение f(x) x3 – x – 1 = 0 (10) имеет корень [1, 2], так как f(1) = – 1 < 0 и f(2) = 5 > 0.

Уравнение (10) можно записать в виде х = х3 – 1. (11) Здесь (х) = х3 – 1 и (х) = 3х2;

поэтому (х) 3 при 1 х и, следовательно, условия сходимости процесса итерации не выполнены.

Если записать уравнение (10) в виде x = x +1, (12) то будем иметь:

(x)= x +1 и (x)=.

33 (x +1) 1 Отсюда 0 < (x)< < при 1 х 2 и, значит, процесс итерации для уравнения (12) быстро сойдётся.

33 Найдём корень уравнения (10) с точностью до 10–2. Вычисляем последовательные приближения хn с одним запасным знаком по формуле xi = xi +1, (i = 0, 1, 2,..., n).

Найденные значения помещены в табл. 1.

1. Значения последовательных приближений xi i 0 1 2 3 xi 1 1,260 1,312 1,322 1,С точностью до 10–2 можно положить = 1,324.

РЕШЕНИЕ СИСТЕМ ЛИНЕЙНЫХ УРАВНЕНИЙ Способы решения систем линейных уравнений делятся на две группы:

1) точные методы, представляющие собой конечные алгоритмы для вычисления корней системы (решение систем с помощью обратной матрицы, правило Крамера, метод Гаусса и др.);

2) итерационные методы, позволяющие получить решение системы с заданной точностью путём сходящихся итерационных процессов (метод итерации, метод Зейделя и др.).

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

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

Рассмотрим систему n линейных алгебраических уравнений относительно n неизвестных х1, х2, …, хn:

a11x1 + a12x2 +... + a1nxn = b1 ;

a x1 + a22x2 +... + a2nxn = b2 ;

(13) L an1x1 + an2x2 +... + annxn = bn.

В соответствии с правилом умножения матриц рассмотренная система линейных уравнений может быть записана в матричном виде Ах = b, (14) где:

Pages:     || 2 | 3 |










© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.