- •Федеральное агентство по образованию
- •Численные методы в инженерных расчетах
- •Численные методы в инженерных расчетах. Учебное пособие по выполнению лабораторных работ.
- •Содержание
- •Предисловие.
- •Лабораторная работа № 1 Тема: Методы оценки погрешностей Теоретическое введение.
- •Абсолютная и относительная погрешности.
- •Правильная запись и округление чисел.
- •Вычисление ошибок арифметических действий.
- •Оценка погрешностей значений функций.
- •Определение количества верных цифр по относительной погрешности приближенного числа.
- •Способы приближенных вычислений по заданной формуле.
- •6.1 Вычисления по правилам подсчета цифр.
- •6.3 Вычисления по методу границ.
- •Приближенные вычисления по формулам с использованием инструментального пакета Excel.
- •Практическая часть.
- •Варианты заданий.
- •Контрольные вопросы
- •Лабораторная работа № 2 Тема: Решение уравнений с одной переменной. Теоретическое введение.
- •Постановка задачи решения уравнений. Отделение корней алгебраических и трансцендентных уравнений.
- •Уточнение корня методом половинного деления.
- •Метод хорд.
- •Метод Ньютона (метод касательных).
- •Практическая часть.
- •Варианты заданий.
- •Контрольные вопросы
- •Лабораторная работа № 3 Тема: Численные методы решения систем линейных алгебраических уравнений (слау). Теоретическое введение.
- •3.1. Системы линейных алгебраических уравнений.
- •3.2. Метод Гаусса.
- •3.3. Решение систем уравнений с помощью табличного процессора Excel.
- •Практическая часть.
- •Варианты заданий.
- •Лабораторная работа № 4 Тема: Приближение функций Теоретическое введение.
- •4.1 Постановка задачи аппроксимации функций.
- •4.2. Существование и единственность интерполяционного многочлена.
- •Интерполяционный многочлен Лагранжа.
- •4.4. Организация ручных вычислений по формуле Лагранжа.
- •4.5. Интерполяционные формулы Ньютона.
- •4.5.1. Конечные разности
- •4.5.2. Первая интерполяционная формула Ньютона
- •4.5.3. Вторая интерполяционная формула Ньютона
- •4.7. Погрешность многочленной интерполяции.
- •Практическая часть.
- •Пояснения к выполнению лабораторной работы № 4
- •Контрольные вопросы
- •Варианты заданий.
- •Лабораторная работа № 5 Тема: Численное дифференцирование и интегрирование Теоретическое введение.
- •5.1. Постановка задачи численного дифференцирования.
- •5.2. Численное дифференцирование на основе интерполяционной формулы Лагранжа.
- •Пример 5.2 Составить интерполяционный многочлен Лагранжа для функции, заданной своими значениями на равноотстоящих узлах ( ):
- •5.3 Численное дифференцирование на основе интерполяционной формулы Ньютона.
- •5.4. Постановка задачи численного интегрирования.
- •5.5. Квадратурные формулы Ньютона-Котеса.
- •5.6 Формула трапеций
- •5.7. Формула Симпсона.
- •5.8 Полуэмпирические оценки точности вычислений по квадратурным формулам.
- •Практическая часть.
- •Варианты заданий.
- •Контрольные вопросы
- •Лабораторная работа № 6 Тема: Численное решение обыкновенного дифференциального уравнения. Теоретическое введение.
- •6.1. Постановка задачи.
- •6.2 Метод Эйлера.
- •6.3. Метод разложения в степенной ряд.
- •6.4 Метод Рунге-Кутта.
- •Практическая часть.
- •Пояснения к выполнению лабораторной работы № 6
- •Варианты заданий.
- •Контрольные вопросы.
- •Приложение 1 Табличный процессор Excel
- •Приложение 2 Задания к контрольной работе для заочной формы обучения.
- •Список литературы
- •Составитель Юрий Алексеевич Занора численные методы в инженерных расчетах
Определение количества верных цифр по относительной погрешности приближенного числа.
Количество верных значащих
цифр в приближенном числе и величина
относительной погрешности этого числа
взаимосвязаны. Эта связь со всей
очевидностью вытекает уже из того, что
по величине х,
учитывая формулу (1.6), можно вычислить
абсолютную погрешность
,
которая, как следует из определения
верных значащих цифр, явно влияет на их
количество в приближенном числе. На
практике иногда удобнее пользоваться
правилом, устанавливающим взаимосвязь
количества верных цифр непосредственно
с величиной относительной погрешности.
Рассмотрим этот вопрос применительно
к приближенному числу х,
записанному в форме с плавающей запятой
х=
,
в предположении что мантисса М
удовлетворяет условию 0,1
<1.
У числа в таком виде первый разряд
мантиссы после запятой всегда отличен
от нуля (число х
в этом случае называют нормализованным).
Заметим также, что для числа, записанного
в указанной Фоме, справедливо ограничение:
<
. (1.7)
Итак, имеются приближенное число х и его относительная погрешность . Нужно установить количество верных в строгом смысле значащих цифр в числе х.
Для каждого известного значения можно подобрать такое наибольшее натуральное n, чтобы имело место неравенство:
(1.8)
Тогда
<
=
<
,
т.е.
<
. (1.9)
Сопоставляя теперь (1.7) и
(1.9) и используя определение цифры, верной
в строгом смысле, можно сделать вывод,
что в мантиссе приближенного числа
верны в строгом смысле по крайней мере
цифр после запятой (и, стало быть, все
эти цифры – значащие). Таким образом,
для того, чтобы по заданной величине
относительной погрешности
найти количество верных значащих цифр
в числе
,
достаточно подобрать наибольшее
натуральное
такое, чтобы имело место неравенство
,
а потом полученное значение
уменьшить на единицу.
Пример 1.8.
Пусть x=984,6;
=0,008.
Очевидно, что 0,008
.
Это означает, что число х
имеет по крайней мере одну верную в
строгом смысле цифру (это первая слева
цифра 9). Полученный результат легко
подтвердить, используя определение
цифры, верной в строгом смысле. Вычислим:
х=984,6·0,008=7,8768,
откуда следует, что в числе 984,6 цифра 9
действительно верна в строгом смысле.
Полученное правило в отдельных случаях проявляет завышенную “осторожность” - при выполнении условия (1.8) в числе могут оказаться верными все n цифр. Зависит это от величины первых значащих цифр числа x.
Вычисление ошибок арифметических действий.
Рассмотрим формулы для учета распространения ошибок при выполнении основных арифметических действий.
Таблица 1.1.
# |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Используя приведенные в таблице 1.1 формулы, можно вести пооперационный учет ошибок арифметических действий на компьютере.
.
Оценка погрешностей значений функций.
Вычисления по формулам нередко предполагают нахождение значений различных математических функций. При этом становится актуальным вопрос о методах определения погрешностей значений элементарных функций.
Пусть функция
дифференцируема
в некоторой окрестности приближенного
значения аргумента
,
а
- абсолютная ошибка значения аргумента.
Тогда абсолютная ошибка значения функции
.
Поскольку на практике ошибка
обычно мала по сравнению со значением
,
воспользуемся приближенным равенством:
.
Заменим
на
.
Это означает, что можно принять
(1.10)
Равенство (1.10) позволяет
получить целую серию формул для оценки
предельных абсолютных погрешностей
значений элементарных функций. Например:
.
Тогда
,
или
;
и так далее.
Из формулы (1.10) можно сделать
важное наблюдение – если значение
модуля производной функции
в точке
меньше единицы, то
<
,
т.е. абсолютная ошибка значения функции
оказывается меньше абсолютной ошибки
значения аргумента. Если же
>1,
то значение функции будет иметь ошибку,
большую ошибки аргумента.
В таблице 1.2 приведены формулы для вычисления предельных абсолютных погрешностей значений некоторых функций одной переменной.
Таблица 1.2.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
