- •Лабораторные работы по высшей математике
- •432027, Ульяновск, ул. Северный Венец, 32
- •Оглавление
- •4.1. Постановка задачи ……………………………………………… 39
- •5.1. Постановка задачи …………………………………………….. 50
- •5.6. Задание на лабораторную работу …………………………….. 54
- •Инструкция по технике безопасности
- •Введение
- •Пакет программ лабораторных работ
- •1. Решение систем линейных уравнений методом гаусса
- •1.1. Постановка задачи
- •1.2. Задание на лабораторную работу
- •1.3. Порядок выполнения работы в компьютерном классе
- •1.4. Пример выполнения работы
- •1.5. Вопросы для самоконтроля
- •2. Решение нелинейных уравнений
- •2.1. Постановка задачи
- •2.2. Отделение корней уравнения. Графический метод
- •2.3. Метод половинного деления
- •2.4. Метод Ньютона
- •2.5. Метод хорд
- •2.6. Комбинированный метод
- •2.7. Задание на лабораторную работу
- •2.8. Порядок выполнения работы в компьютерном классе
- •2.9. Пример выполнения работы
- •2.9. Вопросы для самоконтроля
- •3. Вычисление определенных интегралов
- •3.1. Постановка задачи
- •3.2. Методы прямоугольников и трапеций
- •3.3. Метод Симпсона
- •3.4. Оценка погрешностей методов
- •3.5. Задание на лабораторную работу
- •3.7. Порядок выполнения работы в компьютерном классе
- •3.6. Пример выполнения работы
- •3.8. Вопросы для самоконтроля
- •4. Численное интегрирование обыкновенных дифференциальных уравнений первого порядка
- •4.1. Постановка задачи
- •4.2. Метод Эйлера
- •4.3. Метод Рунге-Кутта
- •4.4. Выбор шага интегрирования
- •4.5. Задание на лабораторную работу
- •4.6. Порядок выполнения работы в компьютерном классе
- •4.7. Пример выполнения работы
- •Точное решение задачи Коши
- •4.8. Вопросы для самоконтроля
- •5. Аппроксимация функции на основании экспериментальных данных по методу наименьших квадратов
- •5.1. Постановка задачи
- •5.2. Выбор типа кривой
- •5.3. Метод наименьших квадратов
- •5.4. Подбор параметров линейной функции методом наименьших квадратов
- •5.5. Подбор параметров квадратичной функции методом наименьших квадратов
- •5.6. Задание на лабораторную работу
- •5.7. Порядок выполнения работы в компьютерном классе
- •5.8. Пример выполнения работы
- •Определение коэффициентов нормальной системы уравнений
- •5.9. Вопросы для самоконтроля
- •6. Прикладной математический пакет «mathcad»
- •6.1. О программе
- •6.2. Основные понятия и функции
- •6.3. Операторы математического анализа
- •6.4. Функции и операторы матриц
- •6.5. Создание декартовых графиков на плоскости
- •6.6. Программные блоки
- •Библиографический список
3.2. Методы прямоугольников и трапеций
Заменим исходную криволинейную трапецию ступенчатой фигурой, состоящей из n прямоугольников, опирающихся на частичные отрезки, причем высоты этих прямоугольников равны значениям функции в начальных или конечных точках частичных отрезков (рис. 3.2). Значение площади этой фигуры и будет давать приближенное значение интеграла (3.1). Результат будет тем более точен, чем больше число частичных отрезков разбиения.
Рис. 3.2. Геометрическая иллюстрация методов левых и правых прямоугольников
Если обозначить значения функции в точках деления через , то, очевидно, будут иметь место следующие формулы:
, (3.2)
, (3.3)
где в формуле (3.2) взяты значения функции в начальных точках, а в (3.3) в конечных точках частичных отрезков. Эти формулы называются формулами левых и правых прямоугольников. Из рис. 3.2 хорошо видно, что если брать значения функции в концах отрезков , то приближение по площади будет или с избытком, или с недостатком на этих отрезках. Поэтому для повышения точности вычисления интеграла методом прямоугольников берут значения функции в точках , т. е. в серединах частичных отрезков (рис. 3.3).
Рис. 3.3. Геометрическая иллюстрация метода средних прямоугольников
Тогда по формуле средних прямоугольников
(3.4)
где – значения в точках .
Оставим разбиение отрезка прежним, но заменим теперь дугу линии , соответствующей частичному отрезку, хордой, соединяющей конечные точки этой дуги. Таким образом, заменяем данную криволинейную трапецию n прямолинейными (рис. 3.4).
Рис. 3.4. Геометрическая иллюстрация метода трапеций
Как правило, площадь такой фигуры более точно выражает искомую площадь, чем площадь n-ступенчатой фигуры, составленной из прямоугольников. Из рис. 3.4 ясно, что площадь каждой прямолинейной трапеции, построенной на частичном отрезке, равна полусумме площадей, соответствующих этому интервалу левого и правого прямоугольников. Суммируя все эти площади, получим
(3.5)
Эта формула и носит название формулы трапеций.
3.3. Метод Симпсона
Как и раньше, разобьем на n равных частей, но предположим, что n – четное число: n=2m. Заменим дугу линии , соответствующую отрезку , дугой параболы (поэтому метод и называют еще методом парабол), ось которой параллельна оси ординат и которая проходит через следующие три точки: начальную точку дуги , среднюю точку и конечную (рис. 3.5).
Рис. 3.5. Геометрическая иллюстрация метода Симпсона
Аналитически это означает, что в отрезке данная функция заменяется квадратичной функцией
(3.6)
Коэффициенты выбираются так, чтобы значения обеих функций были равны при соответственно:
(3.7)
Решая систему (3.7), находят коэффициенты . Проведя подобные замены во всех интервалах (рис. 3.5), будем считать, что площадь исходной трапеции приближенно равна сумме площадей получившихся параболических трапеций, которые называются элементарными.
Покажем, что площадь S трапеции, ограниченной какой-нибудь параболой с осью параллельной оси ординат будет выражаться формулой
(3.8)
где – ордината начальной, – ордината средней и – ордината конечной точек дуги параболы.
Предположим сначала, что основанием трапеции служит отрезок оси Ox, симметричный относительно начала координат, (рис. 3.6).
Рис. 3.6. Площадь параболической трапеции
Для площади такой параболической трапеции имеем выражение:
Так как здесь
то непосредственной подстановкой этих значений в формулу (3.8) убеждаемся в ее справедливости. Эта формула справедлива для любой параболической трапеции рассматриваемого вида с основанием 2h, т. к. всегда можно выбрать декартову систему координат xOy, как показано на рис. 3.6, чтобы основание стало симметричным относительно начала координат. Тогда, применяя формулу (3.8) для всех элементарных параболических трапеций и суммируя площади этих трапеций, получим формулу Симпсона
(3.9)