Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Выбор шага интегрирования

Для устойчивости самого метода проведем выбор шага интегрирования h исходя из условия:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

где — собственное значение матрицы Якоби.

Для комплексного значения условие имеет вид:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Собственными значениями матрицы Якоби порядка n называют корни , где , ее характеристического уравнения, определяемого по формуле:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

где А — матрица Якоби динамической модели;

Е — единичная матрица.

Произведем расчет матрицы Якоби по формуле (55), подставляя начальные значения фазовых координат:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Тогда характеристическое уравнение имеет вид:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Вычислим корни характеристического уравнения с помощью программы MathCad, тогда собственные значения матрицы Якоби имеют вид:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Корни характеристического уравнения имеют отрицательные и нулевые значения, что говорит об устойчивости системы.

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Для гидравлической системы рекомендуемый шаг интегрирования h=0.5с. Выполним проверку устойчивости численного метода Эйлера при данном шаге.

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Проверка условий выполняется, следовательно, шаг h=0.5 обеспечит устойчивость метода и приемлемую точность вычислений.

Видео:2. Дифференциальные уравнения с разделяющимися переменными. Часть 1.Скачать

2. Дифференциальные уравнения с разделяющимися переменными. Часть 1.

Решение систем дифференциальных уравнений методом Эйлера

Формула численного интегрирования неявного метода Эйлера имеет вид:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Совместное преобразование двух последних выражений приводит к записи:

где — модифицированная матрица Якоби на k+1 шаге, которая формируется по следующему правилу:

Диагональные элементы матрицы Якоби на k-ом шаге пересчитываются по формуле:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Остальные элементы не изменяются. Для матрицы размерности 3х3 получаем:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

— модифицированный вектор входных воздействий на k+1 шаге, определяемый по формуле:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Решение системы уравнений (66) дает значение фазовых координат на k+1 шаге, то есть в момент времени tk+1.

Алгоритм неявного метода Эйлера с постоянным шагом интегрирования h:

  • 1) задание шага интегрирования h;
  • 2) задание начальных значений фазовых переменных при t0=0;
  • 3) вычисление времени tk+1=tk+h, где k=0,1,2… ;
  • 4) вычисление модифицированных матриц и на k+1 шаге;
  • 5) решение системы уравнений (66) с целью определения в момент времени tk+1;
  • 6) переход к этапу (3) до тех пор, пока в случае устойчивой системы фазовые координаты не достигнут состояния конечного значения .

Начальные значения вектора определяются на основании входных воздействий системы. В качестве начальных значений фазовых переменных берем вектор начальных значений .

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийИнтегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Рисунок 7 — Графики фазовых координат M(n)0, M(n)1, M(n)2, M(n)3

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Рисунок 8 — Переходный процесс гидросистемы

Видео:18+ Математика без Ху!ни. Дифференциальные уравнения.Скачать

18+ Математика без Ху!ни. Дифференциальные уравнения.

6. Интегрирование с переменным шагом. Автоматический выбор шага интегрирования.

До сих пор рассматривался такой процесс решения задачи (1), (2), когда формула Рунге-Кутта применялась с одной и той же величиной шага интегрирования во всей области вычисления решения вне всякой зависимости от характера поведения решения.

В этом случае говорят, что решение задачи получено с постоянным шагом интегрирования. Применение переменного шага интегрирования позволяет учитывать характер поведения решения и уменьшить общее число шагов, сохранив при этом требуемую точность приближенного решения. Тем самым могут быть снижены объем работы и машинное время и замедлен рост вычислительной погрешности.

Имея в распоряжении способы (53), (54), (71), (73) оценки погрешности метода, величину шага интегрирования можно выбирать автоматически в процессе счета.

6.1. Алгоритм выбора с помощью удвоения и деления шага пополам. Пусть rn+1 —оценка локальной погрешности метода на шаге h, допущенной при вычислении приближенного значения решения yn+1 в точке xn + h. Если оценка превосходит некоторую наперед заданную границу e:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, (90)

то считается, что значение yn+1 решения не удовлетворяет предписанной точности и шаг h объявляется неприемлемым. Полученная точка xn + h и значение yn+1 исключаются из рассмотрения. Выбирается новое значение шага

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

и вновь по той же формуле Рунге-Кутта вычисляется значение решения в точке yn+1 в новой точке xn + h (1) . Если опять выполняется условие (90), то шаг снова делится пополам, и вычисления повторяются. Так происходит до тех пор, пока при какой-то величине шага hn оценка локальной погрешности не станет меньше e:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

После этого считается, что решение дифференциального уравнения продолжено до точки xn+1 = xn + hn. Дальнейшее интегрирование производится из точки xn+1 с шагом hn+1, который выбирается описанным ниже способом.

Если оценка локальной погрешности на шаге hn = xn+1 —xn удовлетворяет неравенству

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, (92)

где K —некоторая константа, то считается, что достигнута точность, значительно превышающая заданную, и шаг интегрирования удваиваются:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Если выполняется неравенство

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, (93)

то считается, что полученное в точке xn+1 решение удовлетворяет заданной точности и шаг интегрирования остается без изменения:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Таким образом, на тех участках изменения независимой переменной, где достигается высокая точность приближенного решения, шаг интегрирования возрастает, а там, где точность не достигается, шаг интегрирования сокращается до необходимых для ее достижения значений. Тем самым обеспечивается выбор величины шага в зависимости от характера поведения решения дифференциального уравнения. Константа K обычно полагается равной Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, где n —порядок используемой оценки локальной погрешности метода. Обычно n = s + 1, а s —порядок формулы численного интегрирования.

6.2. Выбор максимальной для заданной точности длины шага.

Так как оценка rn+1 локальной погрешности метода равна с точностью до членов более высокого порядка малости главному члену локальной погрешности метода, то в силу (9)

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Если оценка rn+1 погрешности превосходит заданную границу e:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

то считается, что значение yn+1 решения не удовлетворяет предписанной точности и шаг h объявляется неприемлемым. Полученная точка xn+h и значение yn+1 исключаются из рассмотрения. В этом случае новый размер шага выбирается не последовательным делением пополам, а с помощью соотношения:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

где a находится из условия выполнения равенства

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. (96)

Из (94), (96) следует, что

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений . (97)

Здесь a 6 раз) точность по-прежнему не достигается, то, как правило, это свидетельствует о том, что явный метод типа Рунге-Кутта не подходит для решения данной задачи.

Указанный прием позволяет получить в процессе счета полезную информацию о характере решаемой задачи, обнаружить сильное измельчение шага интегрирования, избежать большого увеличения числа шагов и затрат машинного времени.

Для системы уравнений с проверкой на точность могут вычисляться либо все компоненты решения, либо некоторые из них, в частности одна компонента. При этом контроль точности может вестись поеомпонентно или по норме. Если контроль ведется покомпонентно, то для различных компонент могут использоваться как различные характеристики точности (абсолютная, относительная погрешности, мера погрешности), так и разные допустимые значения погрешности.

Контроль точности по норме означает, что контролируется некоторая норма оценки погрешности Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. Часто используются нормы:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Из вышеизложенного следует, что формулы Рунге-Кутта очень хорошо приспособлены для интегрирования с переменным шагом, т.к. они позволяют легко менять шаг интегрирования и при этом не требуют никаких дополнительных вычислений и преобразований.

Видео:Откуда появляются дифференциальные уравнения и как их решатьСкачать

Откуда появляются дифференциальные уравнения и как их решать

Моделирование мехатронных систем (стр. 6 )

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийИз за большого объема этот материал размещен на нескольких страницах:
1 2 3 4 5 6 7 8 9

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Рис. 3.1. Механизмы продвижения модельного времени

Во многих случаях, при малых значениях временного шага, это не играет существенной роли. В других – может привести к ошибкам моделирования.

Метод Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийцелесообразно использовать в следующих случаях:

● если моделируется непрерывная система, процессы в которой представляют собой непрерывную цепь равнозначных событий;

● если в моделируемой системе моменты появления событий связаны выполнением некоторых условий, касающихся значений переменных системы, в результате чего эти моменты невозможно заранее определить.

Для мехатронных систем реализация принципа Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийявляется основным способом продвижения модельного времени, т. к. основу мехатронных систем составляют механические устройства, обладающие непрерывной динамикой.

Принцип Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, называемый также принципом особых состояний, предполагает, что продвижение модельного времени обусловлено событиями, происходящими в моделируемой системе. Как и в первом случае, модельное время меняется дискретно на величину Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, однако эта величина не привязана к динамическим характеристикам объекта, а представляет собой временной интервал между последовательными событиями в системе. Величина Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийможет иметь произвольное значение, в т. ч. быть равной нулю, если интервал между событиями пренебрежимо мал.

Необходимым условием реализации моделирования по принципу Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийявляется разработка специальной процедуры планирования событий – так называемого календаря событий [10].

Моделирование по особым состояниям целесообразно использовать, если моделируемая мехатронная система является принципиально дискретной, процессы в системе представляют собой цепь событий, которые распределены во времени неравномерно или интервалы между ними велики, между событиями не происходит изменений в системе.

Обычно зависимость между скоростью изменения модельного времени и скоростью изменения физического времени является переменной и определяется ресурсами компьютера. Однако эта связь может быть и постоянной, что часто весьма желательно, например, при анимации.

Для мехатронных систем достаточно характерным является режим, когда обработка модели должна быть связана с работой реального оборудования. В этом случае говорят, что имеет место моделирование в «режиме реального времени». Режим реального времени – режим обработки данных, при котором обеспечивается взаимодействие вычислительной системы с внешними по отношению к ней процессами в темпе, соизмеримом со скоростью протекания этих процессов. Этот режим обработки данных широко используется в информационно-поисковых системах [21]. Кроме того, моделирование в режиме реального времени актуально при полунатурном моделировании и особенно при использовании моделей в контуре управления реальными техническими системами.

Еще одна проблема в управлении модельным временем связана с тем, что многие технические системы имеют в своем составе компоненты, работающие одновременно, или, как обычно говорят, параллельно. Эти компоненты могут взаимодействовать между собой либо работать независимо друг от друга. Учитывая, что в большинстве случаев моделирование ведется на однопроцессорных ЭВМ, возникает задача не только смоделировать параллельные процессы, но и обеспечить их взаимодействие.

Обычно в таких случаях приходится организовывать квазипараллельные модельные процессы. Одновременные события обрабатываются одно за другим при остановленном модельном времени. Время остается фиксированным до тех пор, пока не будут обработаны все события, привязанные к текущему моменту. В результате два одновременных события выполняются на ЭВМ последовательно, но в один и тот же момент модельного времени, т. е. одновременно с точки зрения системы. После этого модельное время опять оживает и начинает двигаться дальше шагами фиксированной длины (принцип Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений) либо прыгая от события к событию (принцип Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений).

3.2. Алгоритмы численного моделирования
нелинейных динамических систем

Реальные мехатронные объекты и мехатронные системы описываются, как правило, системами нелинейных алгебраических и дифференциальных уравнений. Для большинства задач, представляющих практический интерес, решение их аналитическими методами невозможно. Результаты могут быть получены путем построения приближенных решений с помощью численных методов интегрирования, в частности конечно-разностных методов.

Общая идея численного интегрирования ОДУ

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

заключается в том, что производится дискретизация независимой переменной – времени Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийна интервале Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– и замена ее рядом значений Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений(принцип Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений). Расстояние между двумя соседними значениями Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийназывается шагом интегрирования. В частном случае он может быть постоянным на всем заданном интервале изменения переменной Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. В результате системе дифференциальных уравнений тем или иным способом ставится в соответствие система конечно-разностных уравнений

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

где Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– некоторая вектор-функция, определяемая способом построения метода;

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– количество предыдущих точек, которые используются в методе интегрирования.

Процедура интегрирования предполагает решение полученной системы конечно-разностных уравнений для фиксированных моментов времени Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, начиная с момента Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, для которого определено начальное состояние исследуемой системы Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. Соответственно, решение получается в виде совокупности значений Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийдля заданных моментов времени.

В теории численных методов разработано большое количество различных методов интегрирования, каждому из которых соответствует своя система конечно-разностных уравнений.

Общее представление о них можно получить, разделив их, например, на следующие группы:

● методы явные и неявные;

● методы одношаговые и многошаговые;

● методы первого, второго и т. д. порядков;

● методы с постоянным шагом и методы с автоматическим выбором шага.

3.2.1. Свойства методов численного интегрирования

Основными требованиями, предъявляемыми к численным методам, являются достаточная точность и устойчивость. Дополнительными – универсальность, алгоритмическая надежность, умеренные затраты машинного времени и оперативной памяти ЭВМ [29]. При этом следует учитывать, что практически все эти характеристики имеют смысл только применительно к конкретному объекту моделирования. Поэтому выбор подходящего метода интегрирования может иметь очень большое значение с точки зрения эффективности исследования.

Анализ процесса функционирования технического объекта численными методами всегда сопровождается ошибками в определении характеристик и параметров моделируемого процесса. Эти ошибки обусловлены такими причинами, как неадекватность модели, приближенность исходных данных, свойства используемого метода интегрирования. Первые два фактора возникают на этапе получения исходной модели. Последний зависит от выбранного метода численного интегрирования.

Точность метода можно оценить, проанализировав полную ошибку на каждом шаге интегрирования, однако это достаточно сложная задача, т. к. она предполагает наличие точного решения.

Ошибка интегрирования на Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений-ом шаге включает следующие составляющие:

● ошибка дискретизации Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, обусловленная заменой производных конечными разностями;

● ошибка вычислений Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, связанная с округлением чисел при вычислениях;

● ошибка накопления Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, возникающая вследствие ошибок на предыдущих шагах интегрирования.

Ошибка накопления на Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений-ом шаге равна полной ошибке на предыдущем шаге. Ее оценка связана с исследованием устойчивости численного метода. Если метод устойчив, то Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийсущественно не возрастает и общую ошибку интегрирования можно связать в основном с Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Если даже при небольших погрешностях аппроксимации при каждом шаге интегрирования Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийнакопленная погрешность Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийрастет от шага к шагу, то метод неустойчив и результаты исследований неверны.

Основным источником неустойчивости процесса интегрирования является несогласованность выбора метода интегрирования и метода управления шагом интегрирования с характером решаемой задачи, с особенностями исследуемой системы дифференциальных уравнений. Один и тот же метод может быть достаточно эффективен при решении одной задачи и неприемлем для другой.

Анализ устойчивости метода численного интегрирования для конкретного объекта строится на том, что после дискретизации и алгебраизации его модель превращается в систему разностных уравнений. Устойчивость такой системы можно проверить тем же методом, что и устойчивость обычных дискретных систем. Она зависит от расположения корней характеристического уравнения, полученного для системы разностных уравнений, которое, в свою очередь, определяется выбором формулы интегрирования, шагом интегрирования и собственными значениями матрицы Якоби исходной системы дифференциальных уравнений.

Необходимо рассмотреть влияние ошибок округления, возникающих при реализации методов численного интегрирования на ЭВМ, ограничившись следующим интуитивным рассуждением [23]. Ошибка дискретизации любого устойчивого метода зависит от величины шага интегрирования Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийи стремится к нулю при Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. Следовательно, за счет уменьшения шага мы можем сделать ее сколь угодно малой. Однако чем меньше шаг, тем больше потребуется шагов для решения конкретной задачи, тем больше скажутся на полученном решении ошибки округления. На практике, за счет ограниченной разрядности машинных слов, всегда существует величина шага Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, меньше которой вклад ошибок округления начинает доминировать в суммарной ошибке. Эта ситуация схематически изображена на рис. 3.2.

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Рис. 3.2. Зависимость ошибки интегрирования от величины
шага интегрирования

При малых значениях шага интегрирования велико влияние ошибок округления. В средней части диапазона ошибка растет примерно пропорционально шагу интегрирования, что соответствует методу первого порядка (см. раздел 3.2.5). Превышение шагом значения Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийприводит к неустойчивости процедуры. Значения Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийи Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнениймогут быть найдены исходя из заданной точности интегрирования Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Экстремальную величину шага Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийочень трудно установить заранее, но в задачах, где не требуется слишком высокая точность, необходимый шаг обычно будет значительно больше, чем этот минимум, и основной вклад в полную ошибку будет вносить ошибка дискретизации.

3.2.2. Методы явные и неявные

Процесс формирования математической модели для численного интегрирования обязательно включает этап алгебраизации, который состоит в преобразовании ОДУ в алгебраические уравнения. Он основан на использовании одного из методов численного интегрирования.

Если задано дифференциальное уравнение

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений(3.1)

и начальные условия Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, то очередное значение Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийможет быть получено интегрированием (3.1):

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений(3.2)

Определенный интеграл в (3.2) численно равен площади под кривой Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийна интервале Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений(рис. 3.2). Приближенно эта площадь может быть вычислена как площадь прямоугольника, высота которого равна значению функции Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийна левой границе интервала или значению Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийна правой границе интервала. Очевидно, площади обоих прямоугольников, ограниченных сверху отрезками 1 и 2 на рис. 3.3, будут тем ближе к точному значению интеграла, чем меньше шаг интегрирования Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Рис. 3.3. Иллюстрация шага интегрирования

Подставив в (3.2) приближенные значения интеграла, можно получить две формулы:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений; (3.3)

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. (3.4)

Выражение (3.3) представляет собой формулу явного метода Эйлера. Называется метод явным потому, что неизвестное значение Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийможет быть непосредственно вычислено по известному значению Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийв предыдущей точке.

Формула (3.4) соответствует неявному методу Эйлера. Здесь в правой части выражения используется неизвестное значение Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, поэтому вычислить его непосредственно по этой формуле нельзя.

Более точное значение интеграла (3.2) дает метод трапеций, которому соответствует отрезок 3 на рис. 3.3. Тогда

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. (3.5)

Эта формула относится, очевидно, тоже к неявным.

Для явных методов процедура формирования модели для численного интегрирования ограничивается алгебраизацией исходных дифференциальных уравнений. В частности, формула (3.3) не требует дальнейших преобразований и готова для применения.

Для неявных методов дальнейшие действия зависят от того, какой метод решения системы нелинейных уравнений реализован в данном пакете. Одним из вариантов может быть использование итерационного метода Ньютона [23], который, как известно, обладает наибольшей скоростью сходимости среди практически применяемых методов и в котором многократно решается система линеаризованных алгебраических уравнений.

В этом случае реализуется второй этап подготовки математических моделей для неявных методов, который состоит в линеаризации нелинейных алгебраических уравнений, т. е. в разложении нелинейных функций в ряд Тэйлора и сохранении в результате только линейных членов этого ряда.

Пусть задано нелинейное алгебраическое уравнение

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений(3.6)

где Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– вектор переменных.

Разложение (3.6) в ряд Тэйлора с сохранением только линейных членов дает приближенную замену

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений(3.7)

где Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– начальное приближение, в качестве которого берутся значения переменных на предыдущем шаге интегрирования;

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений;

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– неизвестное значение переменной на шаге интегрирования.

Выражение (3.7) может быть записано как линейное алгебраическое уравнение

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, (3.8)

где Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– вычисляется для известных значений переменных на предыдущем шаге интегрирования;

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Таким образом, процесс численного моделирования в общем случае нелинейных систем неявными методами состоит в формировании и решении на каждом шаге интегрирования системы линейных алгебраических уравнений (СЛАУ)

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений(3.9)

которая включает компонентные и топологические уравнения моделируемой схемы. При этом процедурам алгебраизации и линеаризации подвергаются только компонентные уравнения, т. к. топологические уравнения всегда являются линейными алгебраическими.

Рассмотрим пример, связанный с подготовкой модели для численного решения нелинейного дифференциального уравнения второго порядка:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Первым шагом является сведение данного уравнения к задаче Коши, т. е. к системе уравнений первого порядка за счет введения новой переменной Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Явные формулы метода Эйлера имеют вид

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Неявные формулы запишутся следующим образом:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Для перехода к матричной записи выполним ряд преобразований:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Здесь Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Матричная запись имеет вид

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Формулу (3.7), вообще говоря, необходимо применять итерационно. Решение этого уравнения, найденное для заданного начального приближения Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, следует использовать в качестве очередного приближения в (3.7) и повторять формирование и решение линейного алгебраического уравнения до тех пор, пока два последовательных приближения не станут близкими с заданной точностью. При численном моделировании можно ограничиться только одной итерацией, выбирая достаточно малый шаг интегрирования и учитывая, что при этом значения переменных на предыдущем шаге являются достаточно хорошим приближением.

3.2.3. Выбор между явными и неявными методами
в процедурах моделирования
мехатронных систем

Выбор между явными и неявными методами представляет серьезную проблему. Многие специалисты считают неявные методы более универсальным инструментом для решения задач моделирования технических систем [15, 23]. Следует, однако, заметить, что лишь недавно появились достаточно мощные и универсальные системы автоматизированного моделирования, такие, как, например, MATLAB или МВТУ [17], допускающие выбор явного или неявного метода решения задачи. Раньше использовались либо явные, либо неявные методы, т. к. это требовало разных компонентных моделей.

В современных перспективных системах автоматизированного моделирования, пригодных для моделирования мехатронных систем, применяются, как правило, неявные методы численного интегрирования. Неявные методы лучше приспособлены для решения систем дифференциальных и алгебраических уравнений, к тому же, они более устойчивы. В результате, несмотря на большие затраты машинного времени на каждом шаге интегрирования, связанные с необходимостью решения СЛАУ, общие затраты могут быть значительно меньше за счет увеличения шага интегрирования и уменьшения общего количества шагов.

Рассмотрим эту особенность неявных методов на примере явного и неявного методов Эйлера [21], определяемых формулами (3.3) и (3.4), соответственно.

Применим указанные формулы для численного интегрирования простейшего линейного дифференциального уравнения

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Характеристическое уравнение данной динамической системы имеет вид

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

где Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– постоянная времени системы.

Единственный полюс системы находится в левой полуплоскости, следовательно, исходная система устойчива. Соответственно, любое решение уравнения при Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийстремится к нулю.

Разностное уравнение, соответствующее численному решению явным методом Эйлера, запишется как

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Известно, что условием устойчивости полученного разностного уравнения является

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийили Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Это означает, что выбор Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийможет качественно изменить вид решения, превратив устойчивый процесс в неустойчивый.

Таким образом, на шаг интегрирования наложено очевидное ограничение: он ограничен постоянной времени системы, иначе дискретная аппроксимация станет неустойчивой. Если система имеет несколько постоянных времени, то подобное ограничение связывает шаг интегрирования с самой маленькой постоянной времени системы.

Переход к методам более высокого порядка мало меняет картину. Для метода Рунге–Кутты 4-го порядка требование устойчивости ограничивает шаг величиной

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

или, в более общем виде:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

где Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– максимальное собственное значение матрицы Якоби [29].

Применение неявного метода Эйлера к той же системе дает

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

где ограничение на величину шага выглядит по-другому:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

что позволяет выбрать шаг любой величины, ориентируясь только на требуемый уровень погрешности.

3.2.4. Многошаговые методы интегрирования

До сих пор мы имели дело с методами, зависящими только от Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийи не использующими никаких иных предыдущих значений переменной. Такие методы называются одношаговыми и могут быть представлены в общем виде как

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

с соответствующей функцией Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. Представляется вполне вероятным, что можно добиться большей точности, если использовать информацию о нескольких предыдущих точках Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийИменно так поступают в многошаговых методах.

Вернемся к задаче Коши

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

при этом рассмотрим лишь один большой и важный класс многошаговых методов, который возникает на основе следующего подхода. Если подставить в приведенное дифференциальное уравнение точное решение Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийи проинтегрировать на отрезке Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, то получим

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

где в последнем выражении предполагается, что Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– полином, аппроксимирующий Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. Чтобы записать этот полином, предположим, как обычно, что Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– приближения к решению в точках Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийи узлы расположены равномерно с шагом Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. Таким образом, Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– полином степени Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, удовлетворяющий условиям Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийЭтот полином можно явно проинтегрировать, что ведет к методу

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений.

Для Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийполином Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийесть константа, равная Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, и мы получаем обычный метод Эйлера. Если Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, то Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– линейная функция, проходящая через точки Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийи Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, т. е.

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений

Интегрируя ее от Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийдо Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, получаем следующее выражение:

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, (3.10)

которое соответствует двухшаговому методу интегрирования, поскольку использует информацию в двух предыдущих точках. Аналогично, если Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, то Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийявляется квадратичным полиномом, а соответствующий трехшаговый метод имеет вид

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. (3.11)

Методы, соответствующие формулам (3.10) и (3.11), называются методами Адамса–Бишфорта.

Процесс можно продолжить, используя интерполяционный полином все более высокого порядка. При этом получаются все более громоздкие формулы, но принцип остается тот же.

Многошаговые методы порождают проблему, которая не возникала при использовании одношаговых методов. Интегрирование начинается с начального значения Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, но при Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийдля счета, например, по формуле (3.11) необходима информация о значении функции в точках Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, которая принципиально отсутствует. Обычный выход из положения состоит в использовании какого-либо одношагового метода того же порядка точности, пока не будет набрана необходимая информация.

Заметим также, что многошаговыми могут быть и неявные методы. В этом случае в формулы входят значения Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, которые могут быть определены только неявно и найдены в результате решения системы алгебраических уравнений. Методы этой группы обычно называются методами Адамса–Моултона.

На практике часто используют совместно явную и неявную формулы, что приводит к методам, известным как методы прогноза и коррекции [23].

3.2.5. Порядок метода интегрирования

Главный вопрос при использовании любого численного метода состоит в оценке точности приближенных вычислений Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. В разделе 3.2.1. уже отмечалось, что существуют два источника погрешности при выполнении шага интегрирования:

● ошибка дискретизации, возникающая в результате замены дифференциального уравнения (3.1) разностной аппроксимацией (3.2);

● ошибка округления, накопившаяся при выполнении арифметических операций.

При этом доминирующей является, как правило, ошибка дискретизации.

Будем считать, что все вычисления проводятся точно. Интуитивно ясно, что при Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийошибка дискретизации также должна стремиться к нулю, и это действительно имеет место для любого метода. Однако зависимость скорости уменьшения ошибки от скорости уменьшения Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийдля разных методов интегрирования различна.

Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений,

называемую глобальной ошибкой дискретизации, где Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийи Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений– расчетное и точное значение переменной Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийв момент Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, соответственно. Отметим, что Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийзависит от величины шага, поскольку предполагается, что значения Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийвычисляются при определенном Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. Воспользуемся также стандартным обозначением Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийдля величины, стремящейся к нулю при Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийс той же скоростью, что и Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. В общем случае будем говорить, что функция Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийравна Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, если при Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийвеличина Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийограничена.

Проведя анализ, можно показать [23], что для метода Эйлера Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений. Это обычно выражают утверждением, что метод Эйлера имеет первый порядок. Практическим следствием данного факта является то, что при уменьшении Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийприближенное решение будет все более точным и при Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравненийбудет сходиться к точному решению с линейной скоростью по Интегрирование с переменным шагом в чем смысл такого подхода к решению дифференциальных уравнений, т. е. при уменьшении шага вдвое ошибка дискретизации уменьшится примерно в два раза. Столь медленная сходимость служит препятствием для использования простых методов первого порядка.

📽️ Видео

Дифференциальные уравнения, 1 урок, Дифференциальные уравнения. Основные понятияСкачать

Дифференциальные уравнения, 1 урок, Дифференциальные уравнения. Основные понятия

Дифференциальные уравнения. 11 класс.Скачать

Дифференциальные уравнения. 11 класс.

Интеграл: Азы интегрирования. Высшая математикаСкачать

Интеграл: Азы интегрирования. Высшая математика

Дифференциальные уравнения, 2 урок, Дифференциальные уравнения с разделяющимися переменнымиСкачать

Дифференциальные уравнения, 2 урок, Дифференциальные уравнения с разделяющимися переменными

Дифференциальные уравнения с разделенными переменными. 11 класс.Скачать

Дифференциальные уравнения с разделенными переменными. 11 класс.

Решение дифференциальных уравнений. Практическая часть. 11 класс.Скачать

Решение дифференциальных уравнений. Практическая часть. 11 класс.

Математика без ху!ни. Интегралы, часть 1. Первообразная. Дифференцирование и интегрирование.Скачать

Математика без ху!ни. Интегралы, часть 1. Первообразная. Дифференцирование и интегрирование.

13. Как решить дифференциальное уравнение первого порядка?Скачать

13. Как решить дифференциальное уравнение первого порядка?

Определенный интеграл. 11 класс.Скачать

Определенный интеграл. 11 класс.

Решение дифференциальных уравнений ДИФФУРЫСкачать

Решение дифференциальных уравнений ДИФФУРЫ

ИНТЕГРАЛ С НУЛЯ | определенный интеграл | ТАБЛИЦА ИНТЕГРАЛОВ | сумма РиманаСкачать

ИНТЕГРАЛ С НУЛЯ | определенный интеграл | ТАБЛИЦА ИНТЕГРАЛОВ | сумма Римана

Смысл интеграла и производной. В помощь студентуСкачать

Смысл интеграла и производной. В помощь студенту

Частное решение дифференциального уравнения. 11 класс.Скачать

Частное решение дифференциального уравнения. 11 класс.

Дифференциальные уравнения с разделяющими переменными. 11 класс.Скачать

Дифференциальные уравнения с разделяющими переменными. 11 класс.

Найдем интеграл из дифференциального уравнения!Скачать

Найдем интеграл из дифференциального уравнения!

7. Линейные дифференциальные уравнения первого порядка. Метод Бернулли.Скачать

7. Линейные дифференциальные уравнения первого порядка. Метод Бернулли.
Поделиться или сохранить к себе: