МИНОБРНАУКИ РОССИИ
Санкт-Петербургский государственный
электротехнический университет
«ЛЭТИ» им. В.И. Ульянова (Ленина)
Кафедра математического обеспечения и применения ЭВМ
отчет
по практической работе №2
по дисциплине «Вычислительная математика»
Тема: Изучение понятия обусловленности вычислительной задачи
|
Студент гр. 8383 |
|
Ларин А. |
|
Преподаватель |
|
Сучков А.И. |
2019
Цель работы.
Исследовать обусловленность нахождения корня уравнения для линейной функции.
Основные теоретические положения.
Под обусловленностью вычислительной задачи понимают чувствительность ее решения к малым погрешностям входных данных.
Задачу называют хорошо обусловленной, если малым погрешностям входных данных отвечают малые погрешности решения, и плохо обусловленной, если возможны сильные изменения решения. Количественной мерой степени обусловленности вычислительной задачи является число обусловленности, которое можно интерпретировать как коэффициент возможного возрастания погрешностей в решении по отношению к вызвавшим их погрешностям входных данных. Пусть между абсолютными погрешностями вводных данных Х и уравнения Y установлено неравенство:
,
(1)
где x* и y* - приближенные входные данные и приближенное решение. Тогда величина называется абсолютным числом обусловленности.
Если же установлено неравенство
,
(2)
между относительными ошибками данных и решения, то величину называют относительным числом обусловленности. Для плохо обусловленной задачи >> 1. Грубо говоря, если = 10N, где относительное число обусловленности, то порядок N показывает число верных цифр, которое может быть утеряно в результате по сравнению с числом верных цифр входных данных.
Ответ на вопрос о том, при каком значении задачу следует признать плохо обусловленной, зависит, с одной стороны, от предъявляемых требований к точности решения и, с другой, от уровня обеспечиваемой точности исходных данных. Например, если требуется найти решение с точностью 0.1%, а входная информация задается с точностью 0.02%, то уже значение = 10 сигнализирует о плохой обусловленности. Однако, при тех же требованиях к точности результата, гарантия, что исходные данные задаются с точностью не ниже 0.0001%, означает, что при = 103 задача хорошо обусловлена.
Если рассматривать задачу вычисления корня уравнения Y = f(X), то роль числа обусловленности будет играть величина

где
- корень уравнения.
Исследовать обусловленность задачи нахождения корня уравнения f(x)=0 для линейной функции f(x) = c(x-d). Значения функции f(x) следует вычислять приближенно с точностью Delta (точность входных данных), варьируемой в пределах от 0.1 до 0.000001. Корень уравнения будем находить методом бисекции с задаваемой точностью Eps. Порядок выполнения работы следующий:
Графически
или аналитически отделить корень
уравнения
,
т.е. найти отрезки
,
на которых функция удовлетворяет
условиям применимости метода бисекции.
Составить
подпрограмму вычисления функции
для параметров
и
,
вводимых с клавиатуры. Предусмотреть
округление вычисленных значений
функции
с использованием программы-функции
Round
с точностью Delta,
также вводимой с клавиатуры.
Составить головную программу, вычисляющую корень уравнения с заданной точностью Eps и содержащую обращение к подпрограмме F, программам-функциям BISECT, Round и представление результатов.
Провести вычисления по программе, варьируя значения параметров c (тангенс угла наклона прямой), Eps (точность вычисления корня) и Delta (точность задания сходных данных).
Проанализировать полученные результаты и обосновать выбор точности Eps вычисления корня. Сопоставить полученные теоретические результаты с экспериментальными данными.
Метод бисекции состоит в построении последовательности вложенных друг в друга отрезков, на концах которых функция имеет разные знаки. Каждый последующий отрезок получается делением пополам предыдущего. Этот процесс позволяет найти нуль функции f(x) с любой заданной точностью.
Если рассматривать
задачу вычисления корня уравнения Y
= f(X),
то роль числа обусловленности будет
играть величина
.
Определим абсолютное число обусловленности задачи вычисления корня:

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

Обозначим правую
часть неравенства
как
Пусть d
= 11. Найдем отрезок
,
на котором функция будет удовлетворять
условиям применимости метода бисекции,
то есть на концах данного отрезка
значение функции должны быть разными
по знаку и в его пределах должен быть
лишь один корень.
В табл. 1 рассмотрим случай для с = 5, то есть f(x) = 5(x-11).
Проанализировав результаты табл.1, был получен отрезок [10;15].
Таблица 1 – Поиск отрезка
|
Значение x |
Значение f(x) |
|
5 |
-30 |
|
10 |
-5 |
|
15 |
20 |
|
20 |
45 |
Для построения табл.2-5 была использована программа, которая представлена в приложении А.
В табл. 2 Eps и Delta постоянные, изменяется c. d = 11, a = 10, b = 15.
Таблица 2 – Изменение параметра с
|
Значение c |
Значение Eps |
Значение Delta |
Значение X |
Значение _max |
Значение |
Значение _max |
|
200 |
0.01 |
0.1 |
10.996094 |
0.1 |
0.005 |
хор. |
|
20 |
0.01 |
0.1 |
10.996094 |
0.1 |
0.05 |
хор. |
|
0.2 |
0.01 |
0.1 |
11.015625 |
0.1 |
5 |
пл. |
|
0.02 |
0.01 |
0.1 |
11.015625 |
0.1 |
50 |
пл. |
В табл. 3 Eps и c постоянные, меняется Delta. d = 11, a = 10, b = 15.
Таблица 3 – Изменение параметра Delta
|
Значение c |
Значение Eps |
Значение Delta |
Значение X |
Значение _max |
Значение |
Значение _max |
|
25 |
0.01 |
1 |
11.015625 |
0.01 |
0.04 |
пл. |
|
25 |
0.01 |
0.1 |
10.996094 |
0.1 |
0.04 |
хор. |
|
25 |
0.01 |
0.01 |
10.996094 |
1 |
0.04 |
хор. |
|
25 |
0.01 |
0.001 |
10.996094 |
10 |
0.04 |
хор. |
В табл. 4 изменяется только значение Eps. d = 11, a = 10, b = 15.
Таблица 4 – Изменение параметра Eps
|
Значение c |
Значение Eps |
Значение Delta |
Значение X |
Значение _max |
Значение |
Значение _max |
|
2 |
1 |
0.01 |
11.25 |
100 |
0.5 |
хор. |
|
2 |
0.1 |
0.01 |
11.09375 |
10 |
0.5 |
хор. |
|
2 |
0.01 |
0.01 |
10.996094 |
1 |
0.5 |
хор. |
|
2 |
0.001 |
0.01 |
10.999756 |
0.1 |
0.5 |
пл. |
В табл. 5 изменяются c и Delta, Eps постоянная. d = 11, a = 10, b = 15.
Таблица 5 – Изменение параметров c и Delta
|
Значение c |
Значение Eps |
Значение Delta |
Значение X |
Значение _max |
Значение |
Значение _max |
|
200 |
0.01 |
1 |
10.996094 |
0.01 |
0.005 |
хор. |
|
20 |
0.01 |
0.1 |
10.996094 |
0.1 |
0.05 |
хор. |
|
0.2 |
0.01 |
0.01 |
11.015625 |
1 |
5 |
пл. |
|
0.02 |
0.01 |
0.001 |
11.015625 |
10 |
50 |
пл. |
Если значение _max - хор., то будем считать, что задача хорошо обусловлена, иначе пл. – плохо.