В этом небольшом вступительном разделе мы кратко обсудим, откуда возникают дифференциальные уравнения и кому необходимо их решение. Этот вопрос не возникает при подготовке специалистов-прикладников, например инженерного профиля, – изучение прикладной области неизбежно сопровождается математическим описанием динамики изучаемых в ней процессов и дифференциальные уравнения воспринимаются как естественный инструмент исследования. Иначе обстоит дело при подготовке специалистов по общей и прикладной математике, так сказать «чистых математиков» – изучение аналитических и численных методов решения систем алгебраических, трансцендентных и дифференциальных уравнений часто воспринимается просто как «упражнения для ума» – преподавание ведется «чистыми» математиками и у студента зачастую остаются весьма смутные представления о прикладном значении изучаемых математических курсов (это особенно заметно в подготовке математиков в педагогических институтах). Если школьные задачи все же включают содержательную часть типа наполняемых жидкостями емкостей или движущихся в разных направлениях автомобилей, то студент-математик чаще всего просто получает задание решить конкретную систему уравнений, не представляя зачем и кому это может понадобиться (кроме необходимости получить оценку). Поэтому мы сочли необходимым хотя бы кратко остановиться на одном из прикладных аспектов изучаемого курса.
Под процессом понимают изменение некоторой величины во времени и пространстве – это может быть изменение температуры металла, нагреваемого в проходной печи перед прокаткой, изменение химического состава вещества в химическом реакторе, изменение прибыли предприятия, изменение координат движущейся ракеты, изменение численности определенного вида животных в некотором регионе и т.п. Процесс начинается и протекает благодаря внешним воздействиям на реализующий его объект. Эти воздействия могут быть целенаправленными (управляющими), приложенными для обеспечения желаемого характера протекания процесса – изменение подачи топлива в горелки нагревательного устройства, изменение подачи реагентов в химический реактор, изменения объема используемых оборотных средств предприятия, изменение тяги ракетного двигателя или положения рулей самолета, изменение факторов, влияющих на прирост численности животных и т.д.
Другой тип воздействий на процесс называют возмущениями (нежелательными воздействиями) – колебания силы ветра, отклоняющего летательный аппарат от расчетного курса, изменения теплотворной способности используемого топлива в нагревательных устройствах, колебания спроса на продукцию предприятия и др.
По-видимому, одной из важнейших и труднейших областей интеллектуальной деятельности является управление процессами, которое состоит в определении и последующей реализации таких управляющих воздействий, которые обеспечили бы желательное или по возможности близкое к нему течение процесса, подвергающегося действию нежелательных и часто непредсказуемых и неконтролируемых возмущений. Но для определения необходимых управляющих воздействий необходимо знать, как процесс реагирует на эти воздействия – другими словами, для вычисления управлений необходимо знать их взаимосвязь с управляемыми (выходными) величинами и эта взаимосвязь должна быть выражена в виде математических соотношений, которые называют математической моделью процесса. Таким образом, математическое моделирование процессов является неотъемлемой составной частью процесса управления.
На рисунке мы представили примерную общую схему системы управления процессом:
устройство управления содержит сведения о математической модели управляемого процесса, в каждый момент времени получает информацию о действующем возмущении и отклонении выходной величины от задания, вычисляет управляющее воздействие для уменьшения рассогласования между заданным и фактическим значением выхода процесса и подает его на вход объекта управления.
С общей задачей управления связаны три основные ее составляющие:
1) Задача идентификации математической модели управляемого процесса возникает, если математические соотношения, связывающие входы и выходы процесса неизвестны. Эти соотношения могут быть определены теоретически с использованием основных законов соответствующей прикладной области или экспериментально (что чаще всего и случается, по крайней мере, в технических системах). Рассмотрим некоторые примеры теоретического решения задачи идентификации.
а) Пусть груз массой mперемещается в вязкой жидкости и центрируется в нейтральном положении (x=0) пружинами с линейной характеристикой. Управляющее воздействие – сила F(t),приложенная к грузу, выходная – перемещение груза x. В соответствии с законами механики входное воздействие уравновешивается сопротивлением деформируемых пружин kпx, возникающей во время движения силой вязкого трения r (r – коэффициент вязкого трения) и силой инерции m (m – масса тела).
Уравнение процесса движения тела имеет вид:
m +r +kпx=F(t).
и представляет собой линейное дифференциальное уравнение с коэффициентами, которые могут быть постоянными или зависящими от времени.
б) Пусть к источнику электрического тока подключена электрическая цепь из включенных последовательно омического сопротивления, конденсатора и катушки индуктивности. Напряжение на конденсаторе будем считать выходом процесса, Е(t) – управлением, внутреннее сопротивление источника пренебрежимо мало. В соответствии с законами Кирхгофа сумма падений напряжений в контуре уравновесит ЭДС источников и уравнение процесса формально не отличается от случая механического движения в предыдущем примере:
L +R +(1/C)x=E(t).
в) Пусть математическая модель строится для определения высоты подъема ракеты при различных углах запуска к горизонту. Обозначим горизонтальную координату полета х, вертикальную у, в момент старта x(0)=f(0)=0. Будем считать, что полет происходит в одной вертикальной плоскости (при отсутствии бокового ветра).
Основой для составления математической модели является в этом случае второй закон Ньютона d(mv)/dt=F,где m(t) – масса ракеты, изменяющаяся во времени из-за расхода топлива, v–вектор скорости с горизонтальной x(1)(t) и вертикальной y(1)(t) составляющими, модулем
v(t)=[(x(1)(t)2+y(1)(t)2]1/2
и углом к горизонту
q(t)=arctg(y(1)(t)/x(1)(t)).
F– сумма действующих на ракету сил: силы тяги T(t), силы гравитации mg и силы сопротивления crsv2, пропорциональной плотности воздуха r, поперечному сечению ракеты s и квадрату скорости. Запишем уравнение движения ракеты через координаты х и у с учетом того, что силы тяги и сопротивления действуют вдоль оси ракеты и в сумме дают T–crsv2:
m(1)x(1)+mx(2)=(T–crsv2)cosq,
m(1)y(1)+my(2)=(T–crsv2)sinq–mg.
Или:
x(2)=(1/m)(T–crsv2/2)cosq–(m(1)/m)x(1),
y(2)=(1/m)(T–crsv2/2)sinq–(m(1)/m)y(1)–g.
Получили систему двух нелинейных дифференциальных уравнений второго порядка с начальными условиями v(0)=0, q(0)=0. В системе только один свободный параметр q0 и его изменение будет определять траекторию.
Если моделируется полет простого снаряда с заданной начальной скоростью, отсутствием тяги и изменения массы, то уравнения значительно упрощаются:
x(2)=(1/2m)(–crsv2/2)cosq,
y(2)=(1/2m)(–crsv2/2)sinq–g при v(0)=v0, q(0)=q0.
г) Пусть моделируется задача о популяции двух видов по типу хищник-жертва. Скорости изменения численности жертв х(1) и хищников у(1) определяются их численностями х и у и вероятностью встречи между собой, то есть произведением ху:
dx/dt=ax+bxy, dy/dt=cy+dxy (a>0, b<0, c<0, d>0)
c некоторыми начальными численностями x(0)=x0, y(0)=y0. Эти уравнения известны под названием уравнений Лотки-Вольтерра.
2) Задача анализа. Если математическая модель процесса идентифицирована, она позволяет осуществить имитационные исследования процесса – задаваясь различными функциями управления, можно решением описывающих процесс дифференциальных уравнений получить функции, описывающие реакцию процесса на эти управления. Эта задача является вспомогательной для решения основной задачи – задачи управления.
Задача анализа может рассматриваться в различных постановках.
Если заданы значения выходной величины процесса и всех ее производных до (n–1)-йвключительно (n – порядок старшей производной в дифференциальном уравнении процесса)в некоторый момент времени, принимаемый начальным, при отсутствии управляющего воздействия – это задача определения свободного движения процесса. Она рассматривается на полубесконечном интервале времени и имеет смысл, если хотя бы одна из n–1производных не равна нулю – в противном случае движение отсутствует и выход процесса постоянен. Эта же задача может решаться и при наличии управляющего воздействия – в этом случае решение дифуравнения представляет собой комбинацию из собственного движения и вынужденного под действием управления. Обе подгруппы задач с начальными условиями носят название задачи Коши для дифференциальных уравнений.
Если процесс анализируется на ограниченном интервале времени и часть условий заданы на левой границе интервала, а часть на правой (общее количество дополнительных условий должно быть равно порядку уравнения), мы имеем дело с так называемой граничной задачей определения такого решения, которое удовлетворяет заданным условиям на обеих границах.
3) Задача управления. В задаче управления ставится задача найти управление, переводящее процесс из произвольного (известного) состояния, характеризуемого значениями выходной функции и ее производных в заданное конечное состояние, тоже заданное значением выходной функции и ее производных. Очевидно, что при избыточном количестве дополнительных условий для получения единственного решения необходимо определить, какое из множества возможных управлений надо считать лучшим, чем остальные. В этой постановке мы приходим к классу задач оптимального управления. Например, при управлении движением ракеты можно в одних случаях считать наилучшей такую управляющую функцию, которая обеспечит вывод ракеты в заданную точку с наименьшим количеством израсходованного топлива (управление, оптимальное по затратам ресурсов); в других случаях наилучшим может считаться управление, выводящее ракету в заданную точку за кратчайшее время (управление, оптимальное по быстродействию).
Управляющее устройство и управляемый объект вместе представляют собой систему, описываемую общей системой дифференциальных уравнений. Параметры объекта (им соответствуют коэффициенты уравнения движения), как правило, не поддаются целенаправленному изменению в процессе управления – они определяются его конструкцией и могут сами дрейфовать в процессе старения объекта во время эксплуатации.
Но параметры другой составляющей системы – управляющего устройства – могут изменяться и это обстоятельство привело к разработке способа управления, основанного на изменении структуры системы в процессе ее движения – это так называемые системы с переменной структурой. Проиллюстрировать это можно следующим примером. Пусть в электрическом контуре зарядки конденсатора с последовательно включенной индуктивностью и омическим сопротивлением у нас есть возможность изменять это сопротивление. Если мы сделаем его маленьким, процесс будет колебательным и его выход на заданный уровень будет длительным. Если сделать его настолько большим, что процесс зарядки станет апериодическим, то длительность процесса тоже будет большой. Но можно применить следующий прием – при большом отклонении напряжения на емкости сделать сопротивление маленьким, что обеспечит движение в автоколебательном режиме с большой скоростью, а при приближении к заданному уровню переключиться на большое сопротивление и небольшое оставшееся рассогласование отработать уже в апериодическом режиме – мы получим процесс зарядки с коротким временем перехода в заданное состояние.
Выводы.Итак, в математике процесс – это функция, описывающая изменение значений управляемых величин процесса в зависимости от управляющих и возмущающих воздействий. Задачи, связанные с решением дифференциальных уравнений, чаще всего возникают при математическом моделировании и исследовании на модели динамики процессов в различных объектах или системах. При этом входное управляющее или возмущающее воздействие на систему задается в виде некоторой, в общем случае произвольной, функции времени в правой части неоднородного дифференциального уравнения. Вызванное этим воздействием изменение процесса на выходе системы называют вынужденным движением.
Изменение выходной величины во времени при снятом входном воздействии зависит от конструктивных особенностей системы и начальных условий и носит название собственного движения, определяемого как решение однородного уравнения без правой части. Для линейных систем результирующее движение представляет собой сумму собственного и вынужденного движений в предположении, что возмущающее воздействие приложено в момент времени t=0 и отсутствует в предшествующие моменты, а для t=0 заданы значения выходной координаты и ее производных.
(r – коэффициент вязкого трения) и силой инерции m 
(m – масса тела).
+R 
+(1/C)x=E(t).