
Контрольная работа
Контрольная работа для студентов заочной формы обучения предназначена для проверки знаний по теме «Численное интегрирование и дифференцирование» и состоит из двух задач. Первая задача состоит в построении разностной схемы наивысшего порядка точности, что относится к численному дифференцированию. Вторая задача состоит в нахождении величины определенного интеграла и определения эффективного порядка точности этой величины, что относится к численному интегрированию.
Построение разностных схем
Дана
таблично заданная функция на равномерной
сетке с шагом
,
требуется построить разностную схему,
приближающую производную некоторого
порядка в одной из точек сетки с наивысшим
порядком точности.
Пример № 1 решения задачи
Функция задана таблицей на равномерной сетке
,
.
Требуется
построить схему для производной первого
порядка в точке
–
?
Решение задачи состоит из двух этапов. Этап первый состоит в разложении в ряд Тейлора в исследуемой точке всех имеющихся значений функции и правильного использования свойства равномерности сетки. Второй этап состоит в выделении из построенной на первом этапе системы линейных уравнений искомой величины и подборе коэффициентов, при использовании которых сумма всех уравнений, умноженных на эти коэффициенты, даст схему наивысшего порядка точности.
Перейдем
к первому этапу. Исследуемой точкой
является точка
,
следовательно, необходимо записать
разложение в ряд Тейлора значений
функции в точках
,
и
.
В точке получаем
.
,
,
,
.
В точке получаем
.
,
,
,
.
В точке получаем
.
,
,
,
.
Переходим
ко второму этапу решения задачи. Сначала
требуется в каждом уравнении выделить
искомое неизвестное. По условиям задачи
искомой является величина
.
Запишим все уравенения так, чтобы эта
величина содержалась в левой части:
,
,
.
Теперь
требуется подобрать такие коэффициенты
,
и
,
что линейная комбинация полученных
уравнений с этими коэффициентами
даст схему максимального порядка
точности. Запишем эту линейную
комбинацию:
.
Представим полученное уравнение в более удобном для нас виде, приведя подобные члены по степеням в правой части
.
Теперь преобразуем выражения, полученные в квадратных скобках, вынося общие множители за знак этих скобок. В итоге получаем следующее равенство
.
Остается записать систему линейных уравнений для коэффициентов , и , решение которой позволит нам получить схему максимального порядка точности, то есть требуется обнулить квадратные скобки у максимального количества степеней величины , начиная с младшей степени.
Имеется всего три неизвестных, то есть система может состоять только из трех уравнений. Первое уравнение есть уравнение нормировки — перед искомой величиной должен стоять коэффициент, равный единице. Еще два уравнения должны «уничтожать» первую и вторую степень величины , то есть в итоге мы должны получить схему третьего пордка точности.
.
Эту систему можно решать разными способами, в нашем случае разумно применить метод Гаусса. На первом шаге совершим следующие действия: к третьей строке прибавим первую, результат запишем в качестве новой третьей строки, и из второй строки вычтем первую, результат запишем в качестве новой второй строки. Получим следующую систему
.
Теперь
несложно вычислить все коэффициенты,
начиная со значения числа
:
.
Подставим
полученные коэффициенты
,
и
в последнее равенство для искомой
величины
,
обозначив коэффициент при третьей
степени величины
символом
:
.
Последнее выражение является окончательным ответом.