- •Конспект лекций
- •Пермь 1995
- •Лекция №1 Введение. Понятие о численных методах. История развития численных методов.
- •Приближенные числа. Абсолютная и относительная погрешности.
- •Погрешность числа (а).
- •Абсолютная погрешность().
- •Предельная абсолютная погрешность.
- •Относительная погрешность ().
- •Предельная относительная погрешность (a).
- •Основные источники погрешности.
- •Значащие и верные цифры.
- •Верные цифры.
- •Округление чисел.
- •Связь относительной погрешности приближённого числа с количеством верных цифр этого числа.
- •Лекция № 2 Погрешность суммы.
- •Погрешность разности.
- •Погрешность произведения.
- •Число верных знаков.
- •Погрешность частного.
- •Общая формула для погрешности.
- •Обратная задача теории погрешности.
- •Лекция №3 и 4 Интерполяция функций.
- •Постановка задачи.
- •Конечные разности различных порядков.
- •Первая интерполяционная формула Ньютона.
- •Вторая интерполяционная формула Ньютона.
- •Общая характеристика интерполяционных формул с постоянным шагом.
- •Интерполяционная формула Лагранжа.
- •Частные случаи.
- •Лекции №5 Фрмула Ньютона для неравностоящих узлов Разделённые разности
- •Интерполяционная формула Ньютона для неравностоящих значений аргумента
- •Погрешность формулы Ньютона для неравностоящих узлов
- •Интерполяция сплайками
- •Многочлены Чебышева
- •Выбор узлов интерполирования
- •Обратное интерполирование для равноотстоящих узлов
- •Обратное интерполирование для неравноотстоящих точек
- •Общие выводы по задаче интерполяции
- •Оценка погрешности интерполяционной формулы Лагранжа
- •Лекция n 6 и 7 приближенное решение алгебраических и трансцендентных уравнений.
- •Определение корней Графический способ.
- •Аналитический способ.
- •Графическое решение уравнений.
- •Метод половинного деления.
- •Алгоритм метода
- •Метод пропорциональных частей (метод хорд)
- •Алгоритм метода
- •Метод ньютона (метод касательных)
- •Сходимость метода итераций
- •Сходимость метода
- •Лекция №8 и 9 Решение системы линейных уравнений Общая характеристика методов решения систем линейных уравнений
- •Метод Гаусса
- •Трудоёмкость метода Гаусса
- •Достоинства метода
- •Метод итераций
- •Сходимость метода итераций для решения системы алгебраических уравнений
- •Достоинства метода итераций
- •Метод Зейделя
- •Лекция №10 Численное решение систем линейных уравнений
- •Метод Ньютона
- •Сходимость метода Ньютона
- •Теорема о существовании корней и сходимости процесса Ньютона
- •Метод скорейшего спуска (градиентный метод)
- •Градиент функции u
- •Сходимость градиентного метода
- •Лекции №11 и 12 Приближённое дифференцирование.
- •Графический способ дифференцирования.
- •Алгоритм построения графика производной.
- •Лекция n 13
- •Особенности метода Эйлера.
- •Метод Руте-Кутта.
- •Метод Адамса.
- •Лекция №14 Краевая задача. Методы её решения.
- •Метод конечных разностей для линейных дифференциальных уравнений.
- •Метод прогонки.
- •Метод конечных разностей для нелинейных уравнений второго порядка.
- •Лекция №15 методы обработки эксперементальных данных. Постановка задачи
- •Узловые точки
- •Класс функций
- •Критерий согласия
- •Среднеквадратический критерий
- •Минимальный критерий или критерий чебышева
- •Вероятностно-зональный критерий
- •Точность
- •Метод наименьших квадратов постановка задачи
Погрешность разности.
Разность двух приближённых чисел u = x1- x2 . Предельная абсолютная погрешность u разности равна сумме предельных абсолютных погрешностей уменьшаемого и вычитаемого, т.е.
ux1+x2.
Предельная относительная погрешность разности:
,
где А - точное значение абсолютной велечины разности чисел x1 иx2.
Потеря точности, если уменьшаемое и вычитаемое близки друг к другу (x1 x2) и имеют малые абсолютные погрешности, то число А мало, отсюда u велико.
Пример: x1 = 47,132;x2 = 47,111 - пять верных цифр.
u = 0,0005 + 0,0005 = 0,001
x1 0,00001; x2 0,00001 в 500 раз превышает предельную
относительную погрешность разности по сравнению с исходными данными.
В таких случаях меняют вычислительные схемы.
Погрешность произведения.
Относительная погрешность произведения нескольких приближённых чисел, отличных от нуля, не превышает суммы относительных погрешностей этих чисел.
Пусть u = x1 · x2 · ...· xn.
Предположим, что числа положительны, тогда:
ln u = ln x1 +ln x2 + ...+ln xn
используя приближённую формулу:
ln x d ln x = x/x исходим:
Беря выражение по абсолютной величине, получим:
т.к. точное число близко к точному, то
,
где Аi - точный сомножитель xi
i - относительная погрешность сомножителей xi
- относительная погрешность произведения
12n
Следствие: предельная относительная погрешность произведения равна сумме предельных погрешностей сомножителей, т.е.
ux1x2xn
Зная n произведение u, можно определить его предельную абсолютную погрешность u по формуле:
u = |u|·u
Пример: x1 = 12,3; x2 = 35,48;x1 = 0,05; x2 = 0,005
ux1x2
т.к. произведение u = 436,404, то
u = u·u = 436,404·0,0041 = 1,8
Результат имеет два верных знака и результат записывают:
u = 436 2
Частный случай.
u = k · x, где k - точный множитель, отличный от нуля. Тогда:
u = u и u = |k|·x
т.е. при умножении числа на точный множитель k относительная предельная погрешность не меняется, а абсолютная предельная погрешность увеличивается в |k| раз.
Число верных знаков.
Если все сомножители имеют m верных знаков и число их не более 10, то число верных знаков произведения на одну или две единицы меньше m.
Если количество верных цифр в сомножителях различно, то число верных цифр в произведении на одну или две единицы меньше числа верных знаков в наименее точном сомножителе.
Пример: u = 93,87 · 9,236
относительная погрешность u = 1/2 · (1/9 + 1/9) · 1/1000 =
= 1/9 · 10-3 < 1/2 · 10-3,
отсюда u имеет не менее трёх верных цифр.
Погрешность частного.
Пусть U=
,
тогда ln u=ln x-ln y
,
отсюда
Относительная погрешность частного не превышает суммы относительных погрешностей делимого и делителя. И как следствие:
Пример: Найти частное, число верных знаков в нём x=25,7; y=3,6
u =7,14
=>
В широком смысле (верные цифры те, у которых абсолютная погрешность не превышает 1) две верные цифры, а в узком - одна.
7.14
0.11
Правило определения верных знаков частного.
Если x
и y имеют m
верных цифр, и если
и
- их первые значащие цифры, то за предельную
относительную погрешность частного u
может быть принята величина:
Отсюда правило:
и
,
то u имеет (m-1) верных
знаков;
или
,
то u имеет (m-2) верных
знаков;
Относительная погрешность корня.
,
отсюда
,
отсюда:
Предельная относительная погрешность корня m-ой степени в m раз меньше предельной относительной погрешности подкоренного числа.
Пример:
т.к.
,
отсюда
число а имеет четыре верных знака (в широком смысле), отсюда
а = 3,513.
Погрешность степени.
,
отсюда ln u
= m ·ln
x, отсюда
,
отсюда
.
Вычисление без точного учёта погрешности.
При массовых вычислениях, когда не учитывают погрешности каждого отдельного результата, рекомендуется пользоваться следующими правилами подсчёта цифр:
При сложении и вычитании младший разряд результата должен быть наибольшим десятичным размером, выражаемых верными последними цифрами исходных данных;
При умножении и вычитании в результате надо сохранять столько значащих цифр, сколько их имеет приближённое данное с наименьшим числом верных значащих цифр;
При возведении в квадрат или в куб приближённого числа в результате надо сохранять столько значащих цифр, сколько верных значащих цифр имеет основание степени;
При извлечении квадратного или кубического корней приближённого числа в результате надо сохранять столько значащих цифр, сколько верных цифр имеет подкоренное число;
Во всех промежуточных результатах следует сохранять на одну цифру больше, чем рекомендуют предыдущие правила. В окончательном результате запасная цифра отбрасывается;
Если берётся логарифм выражения, состоящего из приближённых чисел, то необходимо выявить число верных цифр в приближённом числе, имеющем наименьшее их число, и воспользоваться таблицей логарифмов с числом десятичных знаков, на единицу большим. В окончательном результате последние значащие цифры отбрасываются;
Если данные можно брать с произвольной точностью, то для получения результата с k верными цифрами исходные данные надо брать с таким числом, которое согласно предыдущим правилам обеспечивает (k+1) верную цифру в результате.
