Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Kursach_16.doc
Скачиваний:
9
Добавлен:
29.08.2019
Размер:
583.17 Кб
Скачать

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-я модификация.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]