- •8.1.1 Блок-схема………………………………………………………….…..24
- •1 Введение
- •2 Постановка задачи
- •3 Вывод системы дифференциальных уравнений.
- •4 Задание на курсовую работу
- •5.3 Методы численного интегрирования
- •5.3.1 Метод левых прямоугольников
- •5.3.2 Метод правых прямоугольников
- •5.3.3 Метод центральных прямоугольников
- •5.3.4 Метод трапеций
- •5.3.5 Метод Симпсона
- •6 Численная реализация решения системы дифференциальных уравнений
- •6.3.1 Блок-схема
- •6.3.2 Код программы
- •7.1 Реализация в пакете Excel
- •7.2 Реализация в программе MathCad. Метод наименьших квадратов
- •8 Численное интегрирование
- •8.1.1 Блок-схема
- •8.1.2 Код программы
- •8.2 Реализация в программе MathCad
- •9 Выводы
- •10 Список литературы
5.3 Методы численного интегрирования
Вычисление интегралов основано на замене подынтегральной функции интерполяционным многочленом.
5.3.1 Метод левых прямоугольников
Подынтегральная функция заменяется на каждом частичном интервале полиномом нулевой степени, то есть константой, равной значению функции в левой границе частичного интервала, где .
Интеграл будет приближенно равен сумме площадей прямоугольников. При уменьшении шага точность вычисления интеграла повышается.
(7)
5.3.2 Метод правых прямоугольников
Подынтегральная функция заменяется на каждом частичном интервале полиномом нулевой степени, то есть константой, равной значению функции в правой границе частичного интервала, где .
Интеграл будет приближенно равен сумме площадей прямоугольников. При уменьшении шага точность вычисления интеграла повышается.
(8)
5.3.3 Метод центральных прямоугольников
Подынтегральная функция заменяется на каждом частичном интервале полиномом нулевой степени, то есть константой, равной значению функции в центре частичного интервала, где .
Интеграл будет приближенно равен сумме площадей прямоугольников. При уменьшении шага точность вычисления интеграла повышается.
(9)
5.3.4 Метод трапеций
Подынтегральная функция заменяется на каждом частичном интервале интерполяционным полиномом первой степени, то есть графически аппроксимирующая функция является кусочно-линейной. Интеграл будет приближенно равен сумме площадей прямоугольных трапеций. При уменьшении шага точность вычисления интеграла повышается.
(10)
5.3.5 Метод Симпсона
Метод основан на аппроксимации подынтегральной функции отрезками парабол. Разбив отрезок на четное число равных отрезков длиной , рассматриваем отрезки длиной - . На каждом отрезке заменяем подынтегральную функцию полиномом второй степени. В этом случае интеграл будет приближенно равен сумме площадей криволинейных трапеций
, (11)
где .
6 Численная реализация решения системы дифференциальных уравнений
6.1 Реализация в программе MathCad. Метода Рунге-Кутта.
Начальные
условия:
Параметры
цепи:
Вычислим
шаг:
функция,учитывающая
переключение ключа
Зададим
начальные условия:
Итераионные
формулы:
Рисунок 2 - График зависимости силы тока
от времени
Рисунок 3 - График зависимости напряжения
от времени
Таблица 1. Результат работы программы MathCAD
6.2 Реализация в программе MathCad. Метод Эйлера 1-я модификация.
Начальные
условия:
Параметры
цепи:
Вычислим
шаг:
функция,учитывающая
переключение ключа
Зададим
начальные условия:
Итераионные
формулы:
Для
упрощения выражений введем обозначения:
Таблица 2. Результат работы программы MathCAD
Рисунок 2 - График зависимости силы тока от времени
Рисунок 3 - График зависимости напряжения от времени
6.3 Реализация в программе С++. Метода Эйлера 1-я модификация.