Математика изменения

Производная, интеграл и задача с начальными условиями без лишнего формализма.

Производная как скорость изменения

Для моделирования не нужно начинать с большого курса высшей математики. Сначала достаточно понимать, что производная связывает текущее состояние с тем, как оно изменится в ближайший момент. Если скорость ротора выросла на 4 рад/с за 0,2 с, среднее угловое ускорение составило 20 рад/с². Производная dω/dt описывает такой темп в пределе малого интервала.

Это приближение пригодно не только для вычислений, но и для чтения графика. Горизонтальный участок означает почти нулевую производную. Крутой подъём означает большое положительное ускорение. Нулевое значение самой скорости не означает нулевого ускорения: неподвижный ротор может как раз начинать разгон.

Интеграл как накопленное изменение

Если ускорение известно, изменение скорости равно его интегралу по времени. Начальную скорость нужно прибавить отдельно. Для постоянного ускорения α = 10 рад/с² за 0,3 с прирост скорости равен 3 рад/с. При начальной скорости 5 рад/с результат составит 8 рад/с.

Геометрически интеграл — ориентированная площадь под графиком. Отрицательное ускорение уменьшает скорость, поэтому площадь под осью вычитается. Для тока интеграл по времени даёт заряд, а для мощности — энергию. Проверка единиц помогает не перепутать физический смысл: ватт, умноженный на секунду, даёт джоуль.

Как появляется дифференциальное уравнение

Предположим, момент двигателя M постоянен, а сопротивление вращению пропорционально скорости с коэффициентом b. Тогда ускорение зависит от текущей скорости:

Мы не знаем функцию ω(t) заранее, но умеем вычислить её наклон в любой известной точке. Это и есть идея численного интегрирования: находить последовательность состояний по правилу изменения. Чем больше скорость, тем больше вязкое сопротивление и тем меньше ускорение. Когда M = bω, разгон прекращается.

Чтобы решение было однозначным, нужно задать начальное условие. Уравнение описывает и пуск из покоя, и переход с уже вращающегося состояния. Запись ω(0) = 0 выбирает именно пуск. Пропуск начальных условий делает описание задачи неполным.

Вектор состояния

Когда неизвестных несколько, их собирают в вектор. Для двигателя удобно взять y = [i, ω]. Функция правой части получает время и обе величины, а возвращает две производные в том же порядке. Порядок нельзя менять между описанием уравнений, начальными условиями и чтением результата.

def rhs(t, y):
    current, speed = y
    di = (12 - 2 * current - 0.1 * speed) / 0.02
    dw = (0.1 * current - 0.002 * speed) / 0.002
    return [di, dw]

Результат функции — не следующий ток и не следующая скорость. Это их производные: А/с и рад/с². Шаг интегрирования и вычисление нового состояния выполняет численный метод.

Линейность и суперпозиция

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

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

Практическая работа

Для J = 0,002 кг·м², b = 0,002 Н·м·с/рад и M = 0,1 Н·м найдите начальное ускорение при ω₀ = 0 и при ω₀ = 20 рад/с. Затем найдите равновесную скорость. Ответы: 50 и 30 рад/с²; равновесие 50 рад/с. Объясните, почему при разных начальных условиях равновесие одно и то же.

Что запомнить

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