- •Методическое пособие
- •Введение
- •Теория вопроса
- •Погрешности исходных данных
- •Погрешности вычислений
- •Погрешности ограничения
- •Посмотрите, сколько значащих цифр осталось?
- •Индивидуальные задания и последовательность выполнения работы
- •Исходные данные к выполнению работы Таблица 1
- •Контрольные вопросы
- •Табулирование функций
- •Теория вопроса
- •Итерационный процесс продолжается до тех пор, пока значения х1(k), х2(k), х3(k) не станут близкими с заданной погрешностью к значениям х1(k-1), х2(k-1), х3(k-1).
- •Изменим знаки у элементов с нечетной суммой индексов, тогда
- •Зададим начальное приближение
- •Исключим переменную х2 из третьего уравнения системы. Для этого умножим на 0,5 и вычтем из третьего уравнения. Получим новую систему эквивалентную двум предыдущим:
- •Вычислим первые приближения
- •Индивидуальные задания и последовательность выполнения работы
- •Варианты для лабораторных работ
- •Контрольные вопросы
- •Теория вопроса
- •Индивидуальные задания и последовательность выполнения работы
- •Теория вопроса
- •Индивидуальные задания и последовательность выполнения работы
- •Контрольные вопросы
- •Теория вопроса
- •Индивидуальные задания и последовательность выполнения работы
- •Исходные данные к выполнению работы Таблица 5
- •Контрольные вопросы:
- •Теория вопроса
- •Решением дифференциального уравнения является неизвестная функция
- •Индивидуальные задания и последовательность выполнения работы
- •Контрольные вопросы
- •Рекомендуемая литература
Индивидуальные задания и последовательность выполнения работы
Перенесите данные своего варианта в электронную таблицу и проведите интерполирование в соответствии с заданием. Постройте график функции, нанесите на него точку интерполирования, оцените погрешность интерполирования, сохраните работу в своей папке и покажите преподавателю.
Задания для проведения интерполяции. Таблица 3
№ варианта |
Y |
X |
Тип интерполяции |
1 |
2 4 3 1 |
1 2 3 4 |
линейная при Х=2,25 |
2 |
4 3 3.5 3.2 |
1 2 3 4 |
линейная при Х=3,7 |
3
|
27 32 56 100 |
18 20 22 24 |
линейная при Х=23,15 |
4 |
1 3 4 2 |
7 8 9 10 |
Квадратичная при Х= 5,5 |
5
|
2.7 5.4 3 -1 |
1 2 3 4 |
Квадратичная при Х= 1,5 |
6 |
-1 -0.87 -0.25 0 |
-1 -2 -3 -4 |
Квадратичная при Х= -3,6 |
7 |
-1 -0.87 -0.25 0 |
4 3 2 1 |
линейная при Х=1,15 |
8 |
34 28 56 4 |
34 38 42 46 |
линейная при Х=39,5 |
9 |
4 56 28 34 |
34 38 42 46 |
Квадратичная при Х= 36 |
10 |
1 3 5 7 |
4 5 6 7 |
линейная при Х=5,5 |
11 |
4 5 6 7 |
1 2 3 4 |
линейная при Х=2,25 |
12 |
1.2 4.8 4.7 4.7 |
4 5.6 6.2 8 |
Квадратичная при Х= 7 |
13 |
34 38 42 46 |
-1 -0.87 -0.25 0 |
Квадратичная при Х= -0,7 |
14 |
1.75 33.4 22.4 17.2 |
1.55 2.55 4.55 10 |
Квадратичная при Х= -0,7 |
15 |
3 10 50 200 |
10 20 30 40 |
линейная при Х=2,25 |
Лабораторная работа №4 (2ч.) РЕШЕНИЕ НЕЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ С ПОМОЩЬЮ ПЭВМ.
Цель работы: изучение численных методов решения нелинейных алгебраических уравнений.
Теория вопроса
Уравнение f(x)=0, левая часть которого представляет собой многочлен от х, степени больше единицы или содержит тригонометрические, логарифмические и другие элементарные функции называется нелинейным. Например,
Значения х, удовлетворяющие условию f(x)=0, называются корнями этого уравнения. Методы решения нелинейных уравнений делятся на прямые и итерационные. Прямые методы позволяют получить корень уравнения по некоторым конечным аналитическим формулам. Однако, такие формулы в общем случае известны только для небольшого класса нелинейных уравнений.
Наиболее простым и наглядным методом решения уравнений является графический, но его можно использовать только в том случае, когда возможно построить график функции на плоскости.
Пример 1. Решить уравнение tgx=ax (a,x>0). Обозначим y1=tgx, y2=ax. Построим на плоскости графики функций y1 и y2 и отметим, что точки, в которых y1=y2 и будут решением (или решениями) исходного уравнения. То есть, те значения переменной х, при которых y1=y2 обращают исходное уравнение в тождество и являются его корнями.
Так как функция tgx периодическая, то существует бесконечное множество тангенсоид, из которых на графике показаны только три. Причем, при и т.д., . Сказанное означает, что прямая y2=ах будет пересекать тангенсоиды бесконечное число раз, и, следовательно, исходное уравнение имеет бесконечное множество корней.
В этом примере использование графического метода сразу позволяет сделать заключение о множественности корней исходного трансцендентного уравнения. Достоинством графического метода является его наглядность, а недостатками низкая точность решения и ограниченность количества переменных.
Н айти корни уравнения можно и численными методами.
Итерационными (численными) называются методы последовательного приближения к корню уравнения. Сущность этих методов состоит в следующем. Задается точность решения е. Затем поиск корня происходит в два этапа. На первом этапе выполняется отделение корней, т.е. определяются отрезки [a,b]. На каждом из которых функция f(x) непрерывна и имеет только один корень. По сути дела, на этом этапе находят приближенные значения корней с погрешностью, задаваемой длиной отрезка [a,b]. На втором этапе с помощью некоторого алгоритма происходит последовательное уточнение значения корня до нахождения такого приближенного значения х истинной величины корня Х*, при котором выполняется условие
Каждое повторное уточнение корня называется итерацией. Число итераций заранее неизвестно и зависит от требуемой точности решения е.
Отделение корней может быть выполнено следующим образом. Вычисляется с некоторым шагом h (предположительно меньшим расстояния между корнями) ряд значений функции f(x) в интересующей нас области изменения х. Если функция f(x) непрерывна, то корни уравнения будут расположены между теми соседними координатами х, в которых функция имеет разные знаки (в частном случае одно из этих значений функции может равняться нулю) и для которых, следовательно, выполняется соотношение f(x)*f(x+h)0.
Пример 1. Отделить корень уравнения (х-1)2-е-х=0 в области [-5,5].
Просчитываем с помощью ЭТ значения функции от ха= -5 до хb=5 с шагом h=5 и получаем: корень уравнения расположен между х = -3,0 и х = -2,5.
Если функция f(x) на отрезке [a,b] непрерывна, то уточнение корня на отрезке [a,b] до величины, удовлетворяющей одному из условий обычно осуществляется одним из следующих итерационных методов: деления отрезка пополам (МДО); хорд (МХР); касательных (МК). Рассмотрим МДО.
В методе деления отрезка пополам отрезок [a,b] делим пополам и вычисляем с=(а+b)/2 и значение f (c). Если выполняется соотношение , где р – погрешность вычисления функции (например, р=0001), то с является корнем уравнения. В противном случае выбираем тот из отрезков [a,с] или [с,b], на концах которого функция f(x) имеет разные знаки. Выбранный отрезок снова обозначим через [a,b] и опять делим пополам и т.д. Процесс повторяется до тех пор, пока не будет получено точное значение корня или длина b-а очередного отрезка [a,b] не станет меньше заданной точности е. Вышеизложенное реализуется следующим алгоритмом (рис. 5). Учтите, что в алгоритме отсутствует проверка на правильность ввода границ первоначального интервала.
Преимуществом МДО является возможность его применения к любым непрерывным функциям и простота выполнения вычисленного процесса, недостатком – низкая скорость приближения к корню.