Economiko Matematicheskoe Modelirovanie Simonov 1
Economiko Matematicheskoe Modelirovanie Simonov 1
Economiko Matematicheskoe Modelirovanie Simonov 1
РОССИЙСКОЙ ФЕДЕРАЦИИ
Федеральное государственное бюджетное
образовательное учреждение высшего образования
«ПЕРМСКИЙ ГОСУДАРСТВЕННЫЙ
НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ»
П. М. Симонов
ЭКОНОМИКО-МАТЕМАТИЧЕСКОЕ
МОДЕЛИРОВАНИЕ
Часть первая
Пермь 2019
УДК 519.862.4:330.542
ББК 65.012.1
С375
Симонов П. М.
С375 Экономико-математическое моделирование [Электронный ресурс]:
учеб. пособие: в 2 ч. / П. М. Симонов; Перм. гос. нац. исслед. ун-т.
– Электрон. дан. – Пермь, 2019. – Ч. 1. – 3,45 Мб; 230 с. – Режим
доступа: http://www.psu.ru/files/docs/science/books/uchebnie-
posobiya/economiko-matematicheskoe-modelirovanie-simonov-1.pdf. –
Загл. с экрана.
ISBN 978-5-7944-3378-4
Теоретические основы экономико-математического моделирования в пособии
излагаются с широким использованием разнообразного математического аппарата.
Принципиальное отличие этих моделей заключается в том, что они носят динами-
ческий характер. При помощи динамических моделей (и статических, если в этом
есть необходимость) обоснованы выводы экономических закономерностей. Глав-
ное в этом учебном пособии – это модели предприятий, которым посвящены с пя-
того по девятый параграфы. Решение ряда экономических задач доведено до чис-
ленных ответов. Для приобретения практических навыков применения моделей в
каждую главу включены различные формы учебных заданий (вопросы, задачи и т. п.).
Пособие предназначено для бакалавров, изучающих применение математиче-
ских моделей и методов в экономических исследованиях.
УДК 519.862.4:330.542
ББК 65.012.1
3
экономическом факультете Пермского государственного национального
исследовательского университета и на факультете прикладной математики и
механики Пермского национального исследовательского политехнического
университета.
Оно содержит изложение теории, относящейся к производству. Материал
расположен по принципу усложнения применяемого математического
аппарата.
Пособие предназначено для студентов, обучающихся по направлениям
подготовки бакалавров «Прикладная математика и информатика» и «Бизнес-
информатика» с целью углубленного изучения учебной дисциплины
«Экономико-математическое моделирование».
Выражаю особую признательность Э.Д.Каданэру и С.Э.Батищевой – ими
частично были написаны введение и параграфы 7, 8 и 10 (по материалам нашей
совместной книги «Математические модели микроэкономики»).
Кроме того, выражаю глубокую благодарность моим бывшим студентам за
возможность использовать превосходно сделанные ими записи лекций автора
по курсам «Экономико-математическое моделирование», «Динамические
модели экономики», «Математический анализ динамических моделей».
Искренне благодарен коллегам за задачи к параграфам и за ценные
замечания – профессору В.П.Максимову, доцентам: А.Б.Бячкову,
Ю.Н.Еленскому, К.Б.Кузнецову, А.Ю.Куликову, В.А.Соколову и
В.А.Шишкину, старшему преподавателю А.Л.Чадову.
4
ВВЕДЕНИЕ
1
Батищева С.Э., Каданэр Э.Д., Симонов П.М. Математические модели микроэкономики:
учеб. пособие. 2-е изд., перераб. и доп. / Перм. гос. ун-т. Пермь, 2006. Гл. 1.
5
научных дисциплин экономическая наука не всегда дает определенные ответы
на поставленные вопросы. Однако ее практическое применение имеет большое
значение.
Основная цель изучения экономической теории состоит в разрешении
проблемы эффективности использования ограниченных ресурсов для
максимального удовлетворения потребностей людей в экономических благах,
обеспечения наивысшей полезности производимых благ.
При изучении рыночной экономики хозяйственная деятельность
рассматривается в рамках действующего законодательства, которое
предоставляет населению страны определенные свободы.
Исходя из уровня изучения экономическую теорию разделяют на две
достаточно самостоятельные отрасли знаний: макроэкономику и
микроэкономику.
• Макроэкономика – наука о народном хозяйстве страны в целом.
Субъектами макроэкономики являются производственный сектор – все
предприятия (фирмы) и иные производители товаров и услуг, домохозяйства –
все население страны, государство (правительство) и заграница. Народное
хозяйство страны выполняет функции производства, обмена, распределения,
перераспределения и потребления экономических благ. В рыночных условиях
могут выполняться функции производства, обмена и потребления.
Распределение и перераспределение экономических благ осуществляет
правительство (исполнительный орган государства), применяя
административные методы.
Экономические блага – это такие товары и услуги, которые производит
народное хозяйство страны для удовлетворения различных потребностей
населения и организаций. Для производства экономических благ затрачивается
труд людей.
• Микроэкономика – наука о хозяйственной деятельности отдельных
субъектов экономики. Субъектами микроэкономики являются отдельные
предприятия (фирмы) и иные производители товаров и услуг, рынки разной
6
продукции, отдельные люди и их объединения.
В настоящем пособии хозяйство рассматривается в разрезе
самостоятельных субъектов – коммерческих организаций (предприятий, фирм и
т.д.), отдельных групп людей и рынков разных видов продукции в условиях
конкуренции. Структуризация хозяйствующих субъектов осуществляется
только в той мере, которая позволяет анализировать условия равновесия,
оптимизации и рационального выбора решений. Изучению подлежат
экономические законы, действующие на уровне таких субъектов.
Товары, услуги и выполняемые работы удовлетворяют разнообразные
потребности человека и общества, они приносят пользу, т.е. их применение
доставляет человеку благо. По количеству и качеству потребляемых человеком
благ определяется его благосостояние.
Создаваемые в экономике блага называют продукцией.
Продукция (продукт) – это вещественные или информационные
элементы, созданные в процессе производства, а также предоставляемые
услуги и выполняемые работы.
Товары – это продукция, предназначенная для купли-продажи.
Услуги или работы – это труд по осуществлению обмена (купли-
продажи), созданию, изготовлению, обработке, перемещению продукции и т.п.
Например, строительные, ремонтные, сельскохозяйственные работы, торговые
или транспортные услуги.
Произведенная экономикой продукция подлежит обмену на деньги или на
другой продукт либо административному распределению для потребления
(использования, применения) в дальнейшем. Основными показателями
результатов производства являются чистый и валовой продукты.
Валовой продукт (продукция) – статистический показатель,
характеризующий объем продукции народного хозяйства, сфер и отраслей
экономики, объединений и предприятий в денежном выражении. Если же
рассматривать этот показатель с точки зрения затрат на производство
продукции, то валовый продукт образует сумму текущих материальных
7
затрат, амортизации и чистой продукции.
Чистый продукт – часть продукта, которая остается при вычитании из
него материальных затрат. Для отрасли чистый продукт – разность между
валовой продукцией и материальными затратами, включая амортизационные
отчисления. Таким образом, чистый продукт отрасли состоит из оплаты
труда и чистого дохода (прибыли) отрасли.
Ведение хозяйства всегда связано с необходимостью использования
ресурсов. В широком понимании ресурс – это источник или запас чего-нибудь.
Для осуществления хозяйственной деятельности требуются разнообразные
экономические ресурсы. К ним относят природные и трудовые ресурсы,
капитал в денежной форме и в натуральном виде и т.д.
Природные ресурсы – это все предметы окружающей природной среды,
которые человек использует для производства разнообразной продукции. Они
включают площади целинных земель, лес, водные богатства рек, озер и морей,
запасы полезных ископаемых. Некоторые из этих ресурсов расходуются
безвозвратно, например нефть. Другие самовоспроизводятся, например,
рыбные запасы. Для восстановления третьих необходимо выполнение
восстановительных работ, например, возобновление лесных богатств.
Экономическая особенность этих богатств состоит в том, что природные
ресурсы не имеют стоимости, так как на их «производство» не затрачивается
труд.
Источником живого и овеществленного труда в производстве выступают
соответствующие ресурсы: трудовые ресурсы и производственный капитал.
Трудовые ресурсы – это физическая сила, знания и профессиональные
навыки работающих людей. Ресурсы труда в разных странах различаются по
уровню профессиональной подготовки. Соотношение групп разных
специальностей определяется специализацией экономик. Структура трудовых
ресурсов в нефтедобывающих странах в значительной степени отличается от
структуры таковых в промышленно развитых странах.
Производственный капитал – это ресурсы, посредством которых и из
8
которых производятся товары, выполняются услуги и осуществляются
работы. Словом «капитал» в широком смысле обозначают такие денежные
средства (их вещественное или информационное воплощение) в экономике,
которые в результате использования труда приносят доход их владельцам.
Каждый субъект рыночной экономики стремится эффективно
использовать ресурсы. При этом он желает получить наибольшее количество
качественной продукции и/или затратить на производство определенного
количества продукции ресурсы, имеющие минимальную стоимость.
Для оценки эффективности хозяйственной деятельности в каждой отрасли
вводят специальные показатели, соответствующие конкретным экономическим
интересам.
При каждом общественно-политическом устройстве страны действуют
соответствующие способы ведения народного хозяйства. В настоящем пособии
рассматривается экономика стран, в которых соблюдается право частной
собственности на денежный и натуральный капитал (станки, оборудование и
т.д.), на предпринимательскую деятельность, на произведенную продукцию и
свободную торговлю, – стран с рыночной экономикой и капиталистическим
общественно-политическим устройством. Владение капиталом и землей дает
возможность получения платы за их использование, получения дохода в форме
процента или ренты. В капиталистическом государстве частная собственность
установлена законодательно.
9
т.д. В рыночной экономике каждое предприятие действует как
самостоятельный хозяйствующий субъект. Вместе с тем все субъекты
экономики находятся в рыночной среде и подчиняются принятым законам
страны.
10
Цена на производимую продукцию устанавливается на товарных рынках.
В первую очередь цена зависит от конкурентной борьбы продавцов в отрасли,
производящей один вид товаров или услуг, в меньшей степени – от издержек
каждого отдельного предприятия.
11
4. Конкурентные отношения на свободном рынке требуют специальных
исследований. Полученные на их основе выводы дают возможность принимать
научно обоснованные решения и аргументировать рыночную стратегию
предприятия в конкурентных условиях. Однако реальные рыночные отношения
могут существенно отличаться от отношений в идеальной модели. Поэтому на
фоне ряда моделей несовершенной конкуренции особо выделяется модель
совершенной конкуренции.
12
Другой подход состоит в изучении обратных зависимостей –
используемых ресурсов производства от производимой продукции. Ресурсы и
продукция могут быть отражены в натуральных или денежных единицах.
Денежная оценка продукции и ресурсов дает возможность численно
соизмерять затраты с результатами производственной деятельности.
Зависимость затрат ресурсов от объемов выпуска продукции отражается в
математической форме функций производственных затрат.
13
Принцип прибыльности
14
Принцип «затраты – выпуск»
2
Американский экономист российского происхождения, создатель теории межотраслевого
анализа, лауреат Нобелевской премии по экономике за 1973 год «за развитие метода
„затраты – выпуск“ и за его применение к важным экономическим проблемам». (https: //
ru.wikipedia.org / wiki / %D0 %9B %D0 %B5 %D0 %BE %D0 %BD %D1 %82 %D1 %8C %D0
%B5%D0%B2,_%D0%92%D0%B0%D1%81%D0%B8%D0%BB%D0%B8%D0%B9_%D0%92
%D0%B0%D1%81%D0%B8%D0%BB%D1%8C%D0%B5%D0%B2%D0%B8%D1%87).
15
рассматриваемых соотношений является технологический закон и его
математические модели в форме производственных функций.
16
Трудовые ресурсы предприятия имеют следующую структуру:
административно-управленческий и вспомогательный персонал предприятия в
целом, административно-управленческий и вспомогательный персонал цехов,
отделов и т.п., производственные и вспомогательные рабочие цехов и других
подразделений.
18
и затраты на использование производственного капитала. Такая
производственная функция применена для исследования рыночного
равновесия предприятия.
Для динамической модели производства в виде передаточной функции
важной характеристикой является воспроизводство производственного
капитала. В процессе производства составные части ОПФ изнашиваются,
оборотный капитал расходуется, их необходимо замещать. Для возмещения
износа предприятия отчисляют некоторую часть полученного дохода в
амортизационный фонд, который предназначен для инвестиций.
3
Советский математик и экономист, один из создателей линейного программирования.
Лауреат Нобелевской премии по экономике 1975 года «за вклад в теорию оптимального
распределения ресурсов». (https://ru.wikipedia.org/wiki /%D0 %9A %D0 %B0 %D0% BD%
D1%82%D0%BE%D1%80%D0%BE%D0%B2%D0%B8%D1%87,_%D0%9B%D0%B5%D0%B
E%D0%BD%D0%B8%D0%B4_%D0%92%D0%B8%D1%82%D0%B0%D0%BB%D1%8C%D0
%B5%D0%B2%D0%B8%D1%87).
4
Американский экономист, лауреат Нобелевской премии по экономике (1970) «за научную
работу, развившую статическую и динамическую экономическую теорию и внесшую вклад в
повышение общего уровня анализа в области экономической науки». Считается
инициатором «неоклассического синтеза» (объединения в одну концепцию неоклассической
микроэкономики и кейнсианской макроэкономики) и одним из основателей
неокейнсианства. (https://ru.wikipedia.org/wiki/%D0% A1% D0% B0% D0% BC% D1% 83%
D1%8D%D0%BB%D1%8C%D1%81%D0%BE%D0%BD,_%D0%9F%D0%BE%D0%BB_%D0
%AD%D0%BD%D1%82%D0%BE%D0%BD%D0%B8).
19
практической деятельности человека. Человек применяет арифметику и
геометрию для обслуживания хозяйственной практики на протяжении
тысячелетий. Слово модель обозначает такой материальный или знаковый
объект, который в процессе познания реальности заменяет объект-оригинал.
В экономических науках применяют знаковые модели. Широко используют
модели в виде схем, диаграмм, таблиц, графиков, математических формул и
т.д. Воспользуемся определением модели из монографии Л.И. Лопатникова5.
Экономико-математическая модель – математическое описание
экономического процесса или объекта, произведенное в целях их исследования и
управления ими: математическая запись решаемой экономической задачи
(поэтому часто термины «модель» и «задача» употребляются как синонимы).
Существуют еще несколько вариантов определения этого термина.
В самой общей форме модель – условный образ объекта исследования,
сконструированный для упрощения этого исследования. При построении
модели предполагается, что ее непосредственное изучение дает новые знания
о моделируемом объекте. Все это полностью относится и к экономико-
математической модели.
Необходимость применения моделей обусловлена тем, что некоторые
объекты (или проблемы, относящиеся к этим объектам) непосредственно
исследовать невозможно, например, экономическую эффективность
строящихся предприятий. В других случаях натурные эксперименты требуют
много времени и средств, а применение моделей существенно сокращает эти
издержки.
Вид модели экономического объекта или явления зависит от ее
назначения. По назначению модели можно разделить на иллюстративные и
аналитические. Иллюстративные модели предназначены для наглядного
объяснения экономических объектов и явлений. К иллюстративным моделям
относятся преимущественно схемы, диаграммы и графики. Аналитические
5
Лопатников Л.И. Экономико-математический словарь: Словарь современной
экономической науки. 5-е изд., перераб. и доп. М.: Дело, 2003. С. 402–405.
20
модели служат для анализа определенных сторон, свойств, характеристик
объектов моделирования. На основе анализа моделей осуществляют синтез
реальных объектов экономики с заданными свойствами. К аналитическим
моделям можно отнести таблицы, графики и математические формулы.
Аналитические модели экономики принято называть экономико-
математическими моделями (ЭММ). Наряду с методами математического
анализа моделей в пособии широко применяется метод графических
доказательств.
С помощью математических соотношений математические модели
воспроизводят такие свойства экономических объектов и явлений, для изучения
которых эти модели предназначены. Для получения новых знаний или
нахождения нужных экономических решений модели позволяют заменять
натурные эксперименты опытами над моделями. Вместе с тем модель не
является точной копией оригинала, она должна адекватно отражать только
такие характеристики объекта моделирования, которые подлежат изучению.
Под адекватностью моделей понимают достаточно правильное отражение
реальных объектов и явлений.
Слово «моделирование» имеет два значения:
- процесс создания таких моделей, которые соответствуют объектам
моделирования и поставленным задачам;
- процесс анализа или решения модели для использования полученных
ответов на поставленные вопросы.
Экономико-математическое моделирование6 – описание экономических
процессов и явлений в виде экономико-математических моделей. (Иногда тем
же термином обозначают также реализацию экономико-математической модели
на ЭВМ, т.е. «искусственный эксперимент» или машинную имитацию,
машинное решение экономико-математической задачи, – однако это может
вводить в заблуждение.)
6
Лопатников Л.И. Указ. соч. С. 411–412.
21
Как и всякое моделирование, экономико-математическое моделирование
основывается на принципе аналогии, т.е. возможности изучения объекта
(почему-либо труднодоступного для исследования) не непосредственно, а через
рассмотрение другого, подобного ему и более доступного объекта, его модели.
В данном случае таким более доступным объектом является экономико-
математическая модель. При построении моделей те или иные теории или
гипотезы благодаря формализации и квантификации становятся обозримыми,
уточняются, и это способствует лучшему пониманию изучаемых проблем.
Моделирование оказывает и обратное влияние на исследователей, требуя
четкости формулировки исследовательской задачи, строгой логичности в
построении гипотез и концепций.
Практическими задачами моделирования являются, во-первых, анализ
экономических объектов; во-вторых, экономическое прогнозирование,
предвидение развития хозяйственных процессов; в-третьих, выработка
управленческих решений на всех уровнях хозяйственной иерархии.
Последнее, впрочем, требует пояснения. Далеко не во всех случаях
данные, полученные из экономико-математического моделирования, могут
использоваться непосредственно как готовые управленческие решения. Гораздо
чаще они служат в качестве «консультирующих» средств, принятие же самих
управленческих решений остается за человеком. Это объясняется чрезвычайной
сложностью экономических и – шире – социально-экономических процессов.
Экономико-математическое моделирование, таким образом, является лишь
компонентом, хотя и очень важным, в человекомашинных системах
планирования и управления народным хозяйством и экономическими
единицами разного уровня.
Процесс создания экономико-математической модели включает ряд этапов
исследования реальных объектов и явлений, выбора методов их формализации
и проверки адекватности модели. В этом процессе используют все научные
методы познания – абстрагирование, логику и выдвижение гипотез, анализ и
синтез, индукцию и дедукцию и т.д.
22
Анализ или решение моделей зависит от их сложности. Как правило,
объекты экономики (предприятие, отраслевой рынок и т.д.) и происходящие в
них явления (производство и реализация продукции, спрос и предложение
товаров и услуг и т.д.) представляют собой сложные системы. На предприятии
сконцентрирован производственный (основной и оборотный) капитал,
трудовые ресурсы, прибыль и т.д. Отношения работников по поводу
использования ресурсов в процессе производства продукции и получения
дохода и прибыли представляют собой сложную систему.
В широком смысле система состоит из совокупности частей,
находящихся в отношениях и взаимосвязях друг с другом и образующих
некоторую целостность, единство. Части системы называют элементами,
подсистемами и т.п. Эти части могут быть различной природы:
экономической, технической, социальной или иной. Они могут
взаимодействовать между собой и с окружающей средой. В экономических
системах отношения реализуются посредством связей между частями.
Существуют связи по передаче материалов, энергии, информации и др.
Действуют отношения распоряжения и подчинения, последовательности
выполнения отдельных работ и т.д.
Сложность системы определяется количеством и многообразием
входящих в нее частей и связей между этими частями и с внешней средой. При
более глубоком изучении систем в понятие сложности следует включать
способы функционирования частей системы и связей между ними,
многообразие их свойств, вид структуры и др.
Вместе с тем один и тот же объект в зависимости от целей его
исследования может быть представлен экономико-математическими моделями
разной сложности.
Под структурой системы понимают её строение, конфигурацию связей
между отдельными частями. Структура модели отражает внутреннюю
организацию объекта исследования, его составные части и их взаимосвязи, а
также связи с окружающей средой. Экономические системы отличаются
23
большим разнообразием структур. Типичными структурами систем являются:
последовательная, параллельная, встречно-параллельная (имеющая обратные
связи) и смешанная.
Сложные системы нередко характеризуются тем, что обладают такими
свойствами, которых нет у отдельных частей. Поэтому механическое
объединение моделей частей не равноценно модели сложной системы в целом.
Важной характеристикой сложной системы является эмерджентность –
наличие новых свойств, которые не присущи ни одному из элементов, входящих
в систему.
Синергетическая связь определяется как связь, которая при
кооперированных (совместных) действиях независимых элементов системы
обеспечивает увеличение общего эффекта до величины, большей, чем сумма
эффекта этих же элементов, действующих независимо. Следовательно, это
усиливающая связь элементов системы.
Принципиально новые свойства, как правило, возникают в экономических
системах с обратными связями. Например, качество управления в системах
экономической кибернетики невозможно измерить отдельно ни у объекта
управления, ни у управляемого объекта. Другой пример: продавцы стремятся
повысить цены на товары, покупатели при этом уменьшают спрос, и продавцы
вынуждены сдерживать рост цен. Поэтому при изучении сложной системы
недостаточно, а иногда и невозможно, пользоваться методом их расчленения на
части с последующим изучением этих частей в отдельности. Именно сложные
объекты экономики представляют наибольший интерес для моделирования.
Такое моделирование может дать результаты, которые нельзя получить иными
способами.
24
применяются в настоящем пособии.
Экономические объекты и явления обладают большой сложностью.
Сложные модели далее будем трактовать в том смысле, что они состоят из
простых моделей. Однако деление на простые и сложные модели не имеет
достаточно четкой границы. Сложность модели обычно зависит от постановки
задачи изучения объекта – оригинала. Поэтому, даже очень сложные объекты
экономики нередко удается отражать при помощи достаточно простых
математических моделей. Если объектом моделирования является зависимость
определенного набора показателей y от вектора внешних переменных x,
причем характер этой зависимости не меняется, то для создания
математической модели применяют принцип черного ящика. «Черный ящик» –
это объект изучения, внутренняя структура которого совершенно не известна.
Для наблюдения доступны только входные переменные и реакция объекта на
эти воздействия.
Такую математическую модель можно интерпретировать как оператор Φ
преобразования «входа», вектора внешнего воздействия x на объект, в
показатель «выхода», векторную реакцию y объекта исследования. Построить
такую модель означает отыскать оператор Φ, связывающий переменные x и y:
y = Φ(x).
Модель отражает поведение объекта так, что, задавая значения «входа» x,
можно получать значения «выхода» y без участия информации о самом
объекте. Компоненты вектора y = (уi), i = 1, 2, …, n, называют эндогенными
величинами, т.е. определяемыми с помощью модели. Внешние по отношению
к моделируемому объекту показатели, компоненты вектора х = (хj), j = 1, 2, …,
m, называют экзогенными величинами, т.е. определяемыми вне модели.
Простые модели отражают реальные объекты в форме экономических
процессов и соотношений разных показателей. При этом особенно большое
значение имеют функциональные зависимости, отражающие причинно-
следственные связи между моделируемыми показателями.
В экономических системах все явления развертываются во времени,
25
значения экономических показателей изменяются, одни события следуют за
другими. По признаку явной зависимости значений переменных от времени
различают динамические и статические модели. Динамическими моделями
являются экономические процессы и динамические звенья. Рассмотрим
указанные виды простых моделей.
1. Экономические процессы отражают математические модели
зависимости значения показателя y от времени t. Она имеет форму функции
времени
y = W (t, A),
где коэффициенты модели обозначены через вектор A. Например,
экспоненциальный рост производства продукции может иметь вид
y = a et,
где a > 0 и > 1 – коэффициенты модели или параметры моделируемого
объекта. Изменение объема выпуска продукции y происходит в зависимости
от времени t. График функции представлен на рис. 2.
2. Экономические функции и соотношения – однозначные зависимости
экономических показателей от воздействующих на них факторов и
однозначные соответствия между этими величинами. Они имеют
математическую форму равенства, неравенства, предпочтения и т.д. Их можно
записать в общем виде
y F(x, A).
Переменная величина y соответствует значениям вектора x при
фиксированном значении вектора коэффициентов A.
y y
y = ax
dy(t)/dt = a x(t),
y = a et y(t) = a x(t)
y x(t) y(t)
t t x
Рис. 4.
Рис. 2. Рис. 3.
26
а) Факторная модель – математическая модель в виде функциональной
зависимости экономического показателя от воздействующих на него
факторов. Факторные модели являются одной из форм соотношения y
F(x,A):
y = F1(x, A).
В этих моделях переменная y является функцией F1 от вектора
экзогенных переменных x. Коэффициенты модели или параметры
моделируемого объекта обозначены через вектор A. Упрощенно зависимость y
= F1(x, A) переменной y от факторов x обозначают как y(x).
Факторными моделями являются, например, производственные функции,
отражающие зависимость выпуска продукции от вектора используемых в
производстве ресурсов; функции производственных затрат, отражающие
зависимость денежных затрат производства от объема выпуска продукции;
функции спроса и предложения в зависимости от цены и т.д.
Пусть однофакторная производственная функция задана формулой
y = ax,
где a > 0 и 0 < < 1 – коэффициенты (параметры) модели. Тогда ее график
имеет вид, представленный на рис. 3.
б) Ограничения – математическая форма соотношений y F(x,A)
экономических показателей и/или функций в виде больше (>), меньше (<), не
больше (≤) или не меньше (≥). Ограничения можно обозначить следующим
образом:
y F2(x, A),
где переменная величина y имеет указанную в формуле взаимосвязь >, <, ≤, ≥ с
вектором переменных величин x; A – вектор постоянных коэффициентов.
Например, в математическом программировании ограничения такого типа
накладывают на допустимые расходы ресурсов.
При необходимости в пособии вводятся другие формы соотношений: y
F(x, A).
3. Динамическая модель отражает зависимость реакции экономического
объекта в виде (вектор-)функции y(t) от экзогенного (вектор-)процесса x(t),
27
причем это зависимая (вектор-)переменная y(t), или её производные, входят в
уравнения. При этом значение y(t) может относиться к периоду времени,
отличному от времени t, к которому отнесены другие значения y(t). Общую
форму динамической модели можно обозначить следующим образом:
W(x, y, A) = 0,
где (вектор-)функционал W задан на множестве процессов (x(t),y(t)), A –
вектор постоянных коэффициентов.
Так заданную динамическую модель будем называть словом неявная
(модель).
Если удалось разрешить уравнение W(x, y, A) = 0 и выражение для y(t)
такое: y(t) = V(x, A)(t), где V – функционал на множестве процессов x(t), то
такую модель назовем явная модель. Таким образом, нелинейный оператор V
математически преобразует функцию x(t) в y(t), подобно преобразованию
входного процесса в реакцию в объекте-оригинале. В любой момент времени
численные значения x(t) и y(t) определяют состояние входа и выхода объекта
моделирования. В случае если явная модель имеет вид y V ( x(t ), A) , где V (, )
– функция, то тогда эту модель назовем кинематической7.
Один из конкретных видов динамических объектов описывается
следующим дифференциальным уравнением:
dy(t)/dt = a x(t), при x(t) 0,
или в другой форме:
y(t) = a x(t), при x(t) 0,
где a 0 – коэффициент модели (параметр моделируемого объекта). На
внешнее воздействие процесса x(t) модель откликается процессом выхода y(t).
Реакция y(t) зависит от входного процесса x(t). Графически модель
представлена на рис. 4.
По способу определения ответов на поставленные вопросы модели делят
на аналитические и численные. В результате аналитического решения модели
7
Кобринский Н.Е., Кузьмин В.И. Точность экономико-математических моделей. М.: Финансы
и статистика, 1981. С. 76.
28
получают ответ в общем виде. Применение численного метода дает конкретный
количественный результат.
По применяемому математическому аппарату можно выделить следующие
методы и соответствующие им модели: алгебраические, дифференцирование,
вычислительные и т.д. В частности, динамические модели делятся на
разностные, дифференциально-разностные, в форме дифференциальных
уравнений и имитационные модели. В пособии объекты экономики в основном
формализованы в виде линейных обыкновенных дифференциальных уравнений
(ЛОДУ), а также в виде линейных дифференциально-разностных уравнений
(ЛДРУ). Для их решения применено операционное исчисление. Кроме этого,
рассматриваются линейные разностные уравнения (ЛРУ).
По степени определенности модели делят на детерминированные и
вероятностные. К вероятностным моделям в пособии относится только
паутинообразная имитационная модель товарного рынка. Остальные модели
детерминированно отображают реальные объекты.
Важным признаком различия моделей является их назначение. По этому
признаку можно выделить прикладные и теоретические модели.
Благодаря использованию прикладных моделей существует возможность
принимать конкретные хозяйственные решения. Эти модели предназначены
для количественного отображения объектов изучения, они воспроизводят
переменные величины в числовой форме. От количественной модели требуется
точность значений отображаемых экономических показателей. К ним
относятся, например, модели оптимального планирования методом линейного
программирования и двойственные задачи. Для прикладных моделей
применяют в основном диаграммы, таблицы, графики и математические
формулы.
Теоретические модели дают качественную характеристику объекта
исследования. К ним относятся, например, модели кибернетического типа,
модели рыночного равновесия фирмы, модели оптимизации прибыли. Модели
качественного отражения должны правильно воспроизводить вид и форму
29
изменения экономических показателей. Например, выпуск продукции может
иметь растущий вид и монотонную форму, а товарный рынок к точке
равновесия может приближаться либо монотонно, либо колебательно. Для
качественной характеристики, как правило, используют схемы, графики и
математические формулы.
Схемы получили широкое применение для отражения структур сложных
объектов, графики – для анализа равновесия, устойчивости и характеристики
экономического роста.
Равновесие товарного рынка обычно изучают при помощи статических
моделей. Системы, которые стремятся к равновесному состоянию, называют
устойчивыми. Критерий устойчивости статических моделей ввел Дж.Р. Хикс
(J.R. Hicks)8.
Вместе с тем известен классический взгляд на устойчивость. Устойчивость
или отсутствие устойчивости выявляют при помощи динамических моделей.
Асимптотическая устойчивость – это способность системы в ответ на внешнее
воздействие давать реакцию, постепенно стабилизирующуюся на постоянном
уровне. Динамические асимптотически устойчивые системы рыночной
экономики принято называть равновесными системами.
30
реальности, принятия хозяйственных решений и т.п., то эта величина
становится экономическим показателем. Запас материалов можно
рассматривать как простую экономическую величину. Доход предприятия –
экономическая величина, которая одновременно является показателем
экономической эффективности работы предприятия. Неизмеримые
характеристики, например, деловая активность предпринимателя или уровень
квалификации рабочего, нельзя выразить экономическими величинами.
Каждый экономический показатель является характеристикой одной из
сторон деятельности субъекта экономики: экономической эффективности
его деятельности, эффективности затрат или использования ресурсов,
достижения цели принятия тех или других решений, характера
экономического процесса и т.д.
Экономический показатель количественно отражает определенную
характеристику либо отдельной экономической величины, либо зависимости
экономических величин. Например, минимальные затраты средств на оплату
ресурсов или зависимость выручки предприятия от затрат ресурсов,
зависимость цены на товар от объемов спроса или предложения на рынке и т.
д. Экономические показатели можно разделить на две большие группы:
измеряемые и вычисляемые. Данные о значении измеряемых показателей
можно получить непосредственно на реальном объекте, например,
численность рабочих на предприятии. Вычисляемые показатели находят по
формулам, например, фондоотдача как частное от деления получаемого
дохода на стоимость ОПФ.
При создании экономических моделей следует различать интервальные и
моментные экономические величины. К интервальным показателям
относят такие, значение которых можно определить только за некоторый
интервал времени (неделю, месяц, квартал, год и т.п.), например, выпуск
продукции на предприятии. Моментными показателями являются такие,
которые количественно могут отражать экономические величины в любой
момент времени; к примеру, запас продукции или сырья на складе (в
31
натуральном или денежном выражении).
Введем ряд показателей, которые применяют в качестве характеристик
моделей, указанных в п. 5.
Динамизм экономических процессов отражают следующие показатели:
– темп роста (скорость или интенсивность роста)
= dy/dt или 1= y/ t;
– индекс роста (темп или скорости прироста)
Ind y(t) = (dy(t)/dt)/y(t) или Ind1 y(t)= ( y(t)/ t)/y(t).
Если, например, задана функция y(t) = V(x(t), A) в виде y = a et, то темп
роста величины y описывает формула = a αet, а индекс роста – Ind y (t ) = α.
Любая статическая факторная модель y = F1(x, A) отражает зависимость
переменной y от факторов x. Характер зависимости переменных y от каждого
из факторов x обычно исследуют при помощи следующих показателей:
– средний показатель λ, который вычисляют по формуле
λ = y/x;
– предельный (маржинальный) показатель μ имеет формулы
μ = dy/dx или μ1 = y/x;
– эластичность Exy (процентный прирост величины y в ответ на
однопроцентный прирост значения x) вычисляют по формуле
Exy =(dy/dx)/(y/x) или Exy (y / x) /( y / x) .
Смысл этого показателя и другие способы вычисления рассматриваются в
следующем разделе.
Показатели λ, μ и Exy могут быть вычислены применительно к любой
точке зависимости y(x).
Исходные данные y и x в формулах λ = y/x, μ = dy/dx (μ1 = y/x) и
Exy =(dy/dx)/(y/x) ( Exy (y / x) /( y / x) ) всегда являются интервальными
величинами.
33
отражающей зависимость выпуска продукции y от вектора ресурсов
производства x.
Показателями отдачи ресурсов являются: средняя производительность
ресурса λ, предельная производительность ресурса μ и эластичность выпуска
по ресурсу Exy . Эти показатели отражают эффект от использования ресурсов
за некоторый период времени, например за год.
Следует подчеркнуть, что в экономической теории понятие отдачи
ресурсов специфично. Это понятие включает, с одной стороны, выпуск
продукции или его денежное выражение (доход, выручка), с другой –
используемые за тот же период времени ресурсы или расходы на их оплату. В
качестве фактора производства вводят не накопленное количество ресурсов, а
их использование в течение некоторого времени, интервальные величины.
Пусть, например, зависимость объема выпуска продукции y от
используемых ресурсов x выражается формулой y = 10x0,5. Тогда при
использовании x = 4 ед. ресурсов в единицу времени выпуск будет произведен
в объеме y = 20 ед. продукции в единицу времени, а при x = 9 ед. – выпуск y =
30 ед. продукции. Графически данная ситуация иллюстрируется рис. 5.
Среднюю отдачу ресурсов λ вычисляют по формуле λ = y/x.
В нашем примере получаем значение показателя средней отдачи ресурсов
λ = y/x = 10 x0,5/x = 10/x0,5.
y
30 y = 10 x 0,5
20
0 4 9 x
Рис. 5.
34
угла наклона луча, выходящего из начала координат и пересекающего график
зависимости y от x на рис. 5 в точке определения показателя.
Предельную отдачу ресурсов μ вычисляют по формуле μ = dy/dx (μ1 =
y / x ). В приведенном примере имеем μ = dy/dx = d(10x 0,5
)/dx = 5/x 0,5
.В
случае x = 4 получаем значение μ4 = 2,5, а в случае x = 9 – значение μ9 1,67.
Предельная отдача ресурсов μ9 в точке x = 9, y = 30 равна тангенсу угла
наклона касательной к графику производственной функции в точке
определения показателя.
Показатель эластичности выпуска по ресурсам выражает процентный
прирост выпуска при однопроцентном приросте использования ресурсов (за
некоторый период времени).
Различают показатели интервальной, дуговой (средней) и точечной
эластичности.
Показатель интервальной эластичности вычисляют в случае дискретно
заданных точек производственной функции как отношение
y относительный прирост y в процентах
Ex .
относительный прирост x в процентах
35
функции имеем значение точечной эластичности
Exy = μ/λ = (dy/dx)/(y/x) = (5/x 0,5)/(10/x 0,5) = 0,5.
Это означает, что в любой точке графика производственной функции на
рис. 5 эластичность y по x равна 0,5.
Точечную эластичность можно также представить в логарифмической
форме: Exy = d (ln y)/d (ln x).
Заметим, что для любой степенной функции вида y = axα справедливо
Exy = (dy/dx)/(y/x) = [d(axα)/dx)]/(axα/x) = (α axα)/(axα-1) = α.
На российских предприятиях широкое применение получили
вычисляемые показатели, в которых участвуют интервальные и моментные
величины, например, производительность труда и фондоотдача.
Производительность труда вычисляют делением размера дохода
(интервальная величина) за определенный период времени на среднюю
численность работников предприятия (моментная величина).
Фондоотдачу определяют делением размера дохода (интервальная
величина) за определенный период времени на среднюю величину стоимости
производственных фондов (моментная величина).
8. Экономико-математические методы
9
Лопатников Л.И. Указ. соч. С. 409–411.
36
экономических дисциплин»10.
Однако хотя тенденция подмечена верно, она, по-видимому, реализуется
еще не скоро. Экономико-математические методы в действительности имеют
общий объект исследования с другими экономическими дисциплинами –
экономику (или шире: социально-экономическую систему), но разный предмет
науки: т.е. они изучают разные стороны этого объекта, подходят к нему с
разных позиций. И главное, при этом используются особые методы
исследования, развитые настолько, что сами они становятся отдельными
научными дисциплинами особого методологического характера. В отличие от
дисциплин, в которых преобладают онтологические аспекты, а методы
исследования выступают лишь в большей или меньшей степени как
вспомогательные средства, в «методологических» дисциплинах, составляющих
значительную часть комплекса экономико-математические методов, методы
сами оказываются объектом исследования. Кроме того, действительный синтез
экономики и математики еще впереди, потребуется немало времени, пока он
осуществится в полной мере.
10
Шаталин С.С. Функционирование экономики развитого социализма. М.: Изд-во Моск. ун-
та, 1982. С. 25.
37
Список литературы к введению
38
§ 1. Решение линейных обыкновенных дифференциальных
уравнений (ЛОДУ), систем ЛОДУ, линейных разностных
уравнений (ЛРУ)
1. Дифференциальные уравнения
39
1.2. Автономное линейное неоднородное обыкновенное дифференциальное
уравнение (ЛНОДУ) (ЛОДУ) первого порядка
x ax f , a 0 ,
имеет общее решение
x(t ) x* ( x0 x* )eat ,
и, значит, c = x0 x* , следовательно,
x(t ) ( x0 x* )eat x* ,
или
f f f
x(t ) x0 eat = x0 eat (eat 1) .
a a a
Полученное решение выражает x(t) как сумму стационарного значения и
отклонения от стационарного значения. Если первоначально x0 > x*
(отклонение от стационарного значения – положительное), то и в дальнейшем
40
всегда x(t) > x* . Аналогично, если первоначально отклонение от стационарного
значения отрицательное, то и в дальнейшем оно отрицательное. Модуль
отклонения от стационарного значения неограниченно возрастает, если a > 0, и
убывает, сходясь к нулю, если a < 0.
При a 0 получаем уравнение x f , откуда простым интегрированием на
отрезке [0, t ] получим решение x(t ) f t x(0) = f t x0 .
1.3. Формула Коши для ЛНОДУ (ЛОДУ) первого порядка
x(t ) a(t ) x(t ) f (t ) , t 0 ,
имеет вид
t t
a ( s ) ds t
a ( ) d
x(t ) e 0
x(0) e s f ( s)ds ,
0
или
t t s t s
a ( s ) ds a ( ) d t
a ( ) d a ( s ) ds t
a ( ) d
x(t ) e 0 x(0) e 0 e
0
0
f (s )ds = e 0 ( x(0) e
0
0
f ( s)ds) .
a ( s ) ds
Функцию X (t ) e 0 назовем фундаментальным решением ЛОДУ
t
a ( ) d
первого порядка. Функцию C (t , s) e s назовем функцией Коши ЛОДУ
первого порядка. Нетрудно увидеть, что
t t s
a ( ) d a ( ) d a ( ) d
C (t , s) e s = e0 e 0 = X (t ) / X (s) .
Таким образом, справедлива формула Коши для ЛНОДУ первого порядка
в общем виде:
t
x(t ) X (t ) x(0) C (t , s) f ( s)ds .
0
41
где a0 (t ) 0 при всех t 0 (так называемый регулярный случай), имеет вид
t
a ( )
t
a1 ( s ) a10 ( ) d
ds t
e s
x(0)
a0 ( s )
x(t ) e 0
f ( s)ds ,
0
a0 ( s)
или вид
t
x(t ) X (t ) x(0) C (t , s) f ( s)ds ,
0
t
a (s)
a10 ( s ) ds
где X (t ) e 0
– фундаментальное решение ЛОДУ первого порядка;
t
a ( )
a10 ( ) d
e s
42
Замечание 3. Для автономного ЛОДУ первого порядка вида (не
разрешённое относительно производной)
( Lx)(t ) a0 x(t ) a1 x(t ) f (t ) , t 0 ,
a1
a1 (t s )
a0
t e
где a0 , a1 const , a0 0 , получаем X (t ) e a0
, C (t , s) C (t s) .
a0
Здесь D a0 a1 – характеристический многочлен, символ
1
дифференциального оператора; W – передаточная функция.
a0 a1
1.4. Автономное линейное неоднородное обыкновенное дифференциальное
уравнение (ЛНОДУ) (ЛОДУ) n-го порядка.
ЛОДУ – линейные обыкновенные дифференциальные уравнения порядка
n, т.е. уравнения вида
43
1 1
Выражение W p – обратный символ, символ правого
L p D p
t
Интегральный оператор Вольтерра Cf C t s f s ds – оператор
0
1 0 ... 0
0 1 ... 0
= = 1.
... ... ... ...
0 0 ... 1
1.5. Методы решения автономных ЛОДУ.
Lx t 0,
x 0 x0 ,
x 0 x1 , – задача Коши при любых x0 , , xn 1 ,
x ( n1) 0 xn1
44
причем задача Коши всегда имеет единственное решение.
x 0 x 0 1 .
x 0 c1 x1 0 c2 x2 0 c1e c2 e c1 c2 1,
1 0 2 0
x 0 c1 x1 0 c2 x2 0 c11 e c2 2 e c11 c2 2 1,
1 0 2 0
c1 1 c2 , 1 c2 1 c2 2 1, 1 c2 1 c2 2 1, c2 2 1 1 1 ,
1 5 1 5
1
1 1 2 5 1 2 5 5 1
c2 2 , c1 1 c2
2 1 1 5 1 5 5 2 5 2 5
2
5 1
= .
2 5
Решение этого уравнения имеет вид
1 5 1 5
5 1 5 1
x t
t t
e 2
e 2
.
2 5 2 5
45
2) – вещественный корень кратности k (резонанс), тогда
x1 t et , x2 t tet , …, xk t t k 1et .
1-й случай
– не является корнем характеристического уравнения (т.е. корень
кратности 0), поэтому x * t t 0 e t Ak t = e t Ak t , где Ak t – многочлен
степени k с неопределенными коэффициентами.
Пример. Найти общее решение уравнения x t 4 x t 3x t e2t , t 0 .
Решение. Сначала решаем однородное ЛОДУ:
x t 4 x t 3x t 0 .
46
Решаем неоднородное уравнение: здесь f t e2t ; 2 2 , – не
c 1 c2 ,
x 0 c1 3c2 2 1, 1 2c2 2 , c2 1 , c1 0 .
1 c2 3c2 2 1,
И, наконец: x t e3t e2t .
2-й случай
– корень характеристического уравнения кратности j, тогда
x * t t j e t Ak t .
корнем один раз (k=0, степень одночлена j = 1 равна кратности корня), поэтому
x * t t j e t A0 t te3t A .
47
6 9t e3t A 4 1 3t e3t A 3te3t A = 2 Ae3t = e3t ,
1 1 1
находим A , откуда x * (t ) te3t . Ответ: x t c1et c2 e3t te3t .
2 2 2
Пример. Решим задачу Коши:
x t 4 x t 3x t e3t , t 0,
x 0 1,
x 0 0.
Подставляем начальное значение: x 0 c1 c2 0 1 ,
c1 1 c2 ,
1 3 3 3
x 0 c1 3c2 0 , 1 2c2 , c2 , c1 1 .
2 1 c2 3c2 2 0, 2 4 4
7 3 3t 1 3t
И, наконец: x t e te .
4 4 2
2) f t e t Pk t cos t Ql t sin t , составим = i .
1-й случай
– не корень характеристического уравнения (т.е. корень кратности 0),
тогда x * t t 0 e t Am t cos t Bm t sin t , где m max k , l .
2-й случай
– корень кратности j, тогда x * t t j e t Am t cos t Bm t sin t ,
m max k , l .
e t e t e t e t
3) f t e Pk t ch t Ql t sh t , где sh t
t
; ch t ,
2 2
где , – вещественные числа.
~ ~ ~ ~
f t et
Pm t e Qm t e
t t
e e P m t e e Q m t
t t t t
P t Q t
k l Pk t Ql t
2 2
48
~ ~ ~ ~
e( )t P m t e( )t Q m t e1t P m t e2t Q m t .
e1t e 2t
Пусть j1 j2 j , тогда x * t t j e t Am t ch t B m t sh t .
~ ~
1.6. Формула Коши для автономного ЛНОДУ (ЛОДУ) n-го порядка.
Для ЛОДУ порядка n
Lx t 0, x 0 1, x 0 0, … , x( n1) 0 0 ; x2 t : Lx t 0, x 0 0,
x 0 1, … , x( n1) 0 0 ; … ; xn t : Lx t 0 , x 0 0 , x 0 0 , … ,
x( n1) 0 1 .
1
Можно показать, что C t xn (t ) .
a0
Пример. Записать формулу Коши для уравнения
Lx t = x t x t f (t ) , t 0 .
1
c1 c2 1, c1 c2 0 . Находим c1 c2 . Первое фундаментальное решение
2
49
1
имеет вид x1 t (et et ) ch t . Аналогично можно показать, что второе
2
фундаментальное решение имеет вид x1 t sh t . Отсюда находим функцию
t
= sh t s f s ds + ch t x 0 + sh t x 0 .
0
a11 a1n
матрица, A , то рассматривается автономная, стационарная
a ann
n1
система ЛОДУ, x col x1, , xn – вектор неизвестных.
d d
Дифференциальный оператор A L , A x f , т.е. выражение
dt dt
Lx f .
Матрица E A L ( pE A L p ) – матричный символ
50
t
x t C t s f s ds X t x 0 .
0
общее решение
частное решение однородного ур-я
неоднородного ур-я
Составим матрицу X t x1 t , , xn t , где x1 t , , xn t – линейно
столбец
независимые частные решения однородного уравнения Lx 0 . Выполняется
аналог вронскиана – det X t 0 для любого t 0 . Вычислим X 0 , X 1 0 и
f t e t Pk t cos t Ql t sin t ,
51
характеристического уравнения; j 0 , если число i совпадает с корнем
кратности j.
1.9. Автономные линейные системы однородных обыкновенных
дифференциальных уравнений (ЛСООДУ) (ЛОДУ) второго порядка:
x1 a11 x1 a12 x2 ,
x2 a12 x1 a22 x2 .
Общее решение указанной системы второго порядка имеет вид
x1 (t ) c1 A1e1t c2 A2 e2t ,
x2 (t ) c1 B1e1t c2 B2 e2t ,
a a12
где 1 , 2 – собственные числа матрицы A 11 , т.е. корни
a21 a22
A1e1t A2 e2t
матрицу X (t ) t 2 t
. Вычислим X 1 0 и X 1 0 X t . Тогда будет
B1e B2 e
1
C t s X t s .
Пример. Решить неоднородную систему ЛОДУ, записать формулу Коши:
x1 2 x1 4 x2 4e2t ,
x2 2 x1 2 x2 .
Решение. Сначала решим однородную систему ЛОДУ
x1 2 x1 4 x2 ,
x2 2 x1 2 x2 .
Характеристическое уравнение системы
52
2 4
= 2 4 =0
2 2
имеет корни 1,2 2i .
2 2i 4 x1
2 y = 0,
2 2i 1
т.е.
(2 2i) x1 4 y1 0,
2 x1 (2 2i) y1 0,
1 i
что равносильно (1 i) x1 2 y1 или y1 x1 .
2
Положив x1 2 , имеем
x 2
X 1 (0) 1 .
1
y 1 i
Получим частное решение
2 2it
X 1 (t ) e .
1 i
Выделив в полученном выражении вещественную и мнимую части,
получим два линейно независимых решения:
2 2cos 2t 2i sin 2t
X 1 (t ) (cos 2t i sin 2t ) = =
1 i cos 2t i sin 2t i (cos 2t i sin 2t )
2cos 2t 2i sin 2t 2cos 2t 2sin 2t
= = + i .
cos 2t sin 2t i (sin 2t cos 2t ) cos 2t sin 2t sin 2t cos 2t
Тогда
2cos 2t 2sin 2t
X 1 (t ) = c1 Re X1 (t ) c2 Im X1 (t ) = c1 + c2 .
cos 2t sin 2t sin 2t cos 2t
2cos 2t 2sin 2t
Составим матрицу-функцию X (t ) = .
cos 2t sin 2t sin 2t cos 2t
53
1
0
2 0 1 1 0 2
Вычислим X (0) = , X 1 (0) = = . Следовательно,
1 1 2 1 2 1
1
2
фундаментальная матрица имеет вид
1
2 0
2cos 2t 2sin 2t cos 2t sin 2t
X (t ) = cos 2t sin 2t sin 2t cos 2t = sin 2t cos 2t .
1 1
2
Отсюда находим матрицу-функцию Коши:
cos 2(t s) sin 2(t s)
C t s = X (t s) = .
sin 2(t s ) cos 2(t s )
Найдем теперь частное решение исходной системы. Имеем
f1 (t ) 4e2t
F (t ) .
f 2 (t ) 0
Таким образом, 2 , m 0 – многочлен нулевой степени. Частное
Ae 2t 2 Ae2t
решение X 2 (t ) системы ищем в виде X 2 (t ) = 2t . Тогда X 2 (t ) = 2 t
,
Be 2 Be
и, подставив в исходную систему найдённое решение, получим
2 Ae2t 2 Ae2t 4 Be2t 4e2t ,
2 t 2 t 2 t
2 Be 2 Ae 2 Be .
0
Таким образом, 4B 4 A 4 , откуда A 0 , B 1. Тогда X 2 (t ) = 2t ,
e
следовательно,
x1 (t ) 2cos 2t 2sin 2t 0
x (t ) = X 1 (t ) + X 2 (t ) = c1 cos 2t sin 2t + c2 sin 2t cos 2t + e 2t .
2
Ответ: x1 (t ) = 2c1 cos2t 2c2 sin 2t , x2 (t ) = c1 (cos2t sin 2t ) +
cos 2t sin 2t
c2 (sin 2t cos 2t ) + e2t . Фундаментальная матрица: , матрица
sin 2t cos 2t
54
cos 2(t s) sin 2(t s)
Коши: C t s = X (t s) = .
sin 2(t s) cos 2(t s)
1.10. Автономные линейные системы однородных обыкновенных
дифференциальных уравнений (ЛСООДУ) (ЛОДУ) n-го порядка
Для решения системы
x1 a11 x1 ... x1n xn ,
.................................. (1)
x a x ... x x ,
n n1 1 nn n
x1 a11...a1n
x ... , A ............ ,
x a ...a
n n1 nn
надо найти корни характеристического уравнения
a11 a12 ... a1n
a21 a22 ... a2 n
0. (2)
..........................................
an1 an 2 ... ann
3 4 2 5 2 0, 1 2, 2 3 1.
Для простого корня 1 2 находим собственный вектор ( , , ) , решая
систему
0,
2 2 0, (6)
2 0
(коэффициенты этой системы равны элементам детерминанта (5) при
2 ). Из (6) находим 2 . Значит, вектор (1, –2, 2) – собственный и
x e2t , y 2e2t , z 2e2t (7)
– частное решение системы (4).
Для кратного корня 1 сначала определим число линейно
независимых собственных векторов. При 1 из (5) получаем матрицу
56
1 1 1
2 1 1 .
2 1 1
Ее порядок n 3 , ранг r 2 . Число линейно независимых собственных
векторов равно m n r 1. Корень 1 имеет кратность k 2 . Так как
k m, то решение надо искать в виде произведения многочлена степени
x C2et C3e2t , y (C1 C2t )et 2C3e2t , z (C2 C1 C2t )et 2C3e2t .
1.11. Частное решение линейной неоднородной системы с постоянными
коэффициентами
xi ai1xi ... ain xn fi (t ), i 1,..., n, (11)
57
можно искать методом неопределенных коэффициентов в том случае, когда
функции fi (t ) состоят из сумм и произведений функций b0 b1t ... bmt m ,
e t , cos t , sin t. Это делается по тем же правилам, что для одного линейного
уравнения с постоянными коэффициентами (см. п. 1.5 и 1.8 со следующим из-
менением). Если fi (t ) Pmi (t )e t , где Pmi (t ) – многочлен степени mi , то
В системе (13) для функций te3t , e3t sin t , e3t cos t числа i
соответственно равны 3, 3 i, 3 i. Поэтому надо отдельно найти частные
решения систем:
x 4 x y te3t , y x 2 y, (14)
58
x 4 x y e3t sin t , y x 2 y te3t cos t. (15)
Для системы ( 1 4 ) i 3 1 2 , s 2, m 1. Согласно (12)
частное решение можно искать в виде
x1 (at 3 bt 2 ct d )e3t , y1 ( ft 3 gt 2 ht j )e3t .
Для системы (15) i 3 i 1,2, s 0, m 1. Частное решение
имеет вид
x2 (kt l )e3t sin t (mt n)e3t cos t ,
2. Разностные уравнения
59
стационарному состоянию из любого начального состояния, т.е. стационарное
состояние – асимптотически устойчиво.
Как видим, между разностными уравнениями и дифференциальными
уравнениями (они рассматривались в первой части §1) – большое сходство, оно
имеет не внешний, а органический характер. Возникает вопрос о
существовании единой математической теории, охватывающей разностные и
дифференциальные уравнения. О такого рода теории можно прочитать в статье
П.Б. Гусятникова в электронном издании «Соросовский образовательный
журнал» (1999. № 10. С. 115–121). В этой теории рассматривается абстрактное
алгебраическое кольцо, вводится формально понятие оператора
дифференцирования кольца, составляется линейное «дифференциальное»
уравнение и затем достаточно легко выписывается решение этого уравнения.
Следствием оказываются формулы общих решений разнообразных
дифференциальных и разностных уравнений.
2.2. Автономное линейное неоднородное разностное уравнение (ЛНРУ)
(ЛРУ) первого порядка
xt axt 1 f , a 1, t 1,2, ,
где f – усиливающая переменная (forcing variable), имеет общее решение
xt x (t ) x* = ca t x* ,
f
где x* – стационарное состояние.
1 a
Частное решение с начальным состоянием x0 имеет вид
xt ( x0 x* )at x* ,
или
f t
xt x0
a
1 a
f
1 a
= xt x0at
f
1 a
1 at .
60
x1 x0 f , при t 1 получаем x2 x1 f = x0 f f = x0 2 f , при t 2
получаем x3 x2 f = x0 2 f f = x0 3 f . И, наконец, получаем xt =
f t x0 .
Комментарий. Как и в случае дифференциального уравнения, общее
решение неоднородного разностного уравнения представляет собой сумму
общего решения однородного уравнения и частного решения неоднородного
уравнения. В качестве частного решения неоднородного уравнения используем
стационарное значение, т.е. такую точку x* , что уравнение справедливо при
xt x* :
x* ax* f , a 0 .
i 0
61
собой общее решение однородного уравнения xt = axt-1, x (t ) x(0)a t . Проверим,
что другое слагаемое Сft есть стационарное состояние неоднородного
уравнения при условии x(0) 0 . Действительно, (1 – aL) Сft = f t при x(0) 0 .
2.4. Автономное ЛНРУ (ЛРУ) первого порядка, формула Коши
xt axt 1 ft , a 0 , a 1 , t 1,2, ,
где { ft } – экзогенная последовательность значений усиливающей переменной,
имеет общее решение, это формула Коши для автономного ЛНРУ (ЛРУ)
первого порядка:
t 1 t
xt = a x(0) Cft = a x(0) a ft i = a x(0) a t i fi .
t t i t
i 0 i 1
Ct s at s .
Комментарий. Как и раньше, общее решение однородного уравнения
имеет вид x (t ) = cat. Чтобы получить частное решение неоднородного
уравнения, положим x0 = 0.
Тогда x1 f1 , x2 = ax1 + f2 = a f1 + f2, x3 = ax2 + f3 = a2f1 + af2 + f3, ...,
t 1 t
xt = a f
i 0
i
t i = a
i 1
t i
fi .
Lx t a0 x t n a1 x t n 1 an 1 x t 1 an x t 0 , t 0,1, 2 ,…,
x t x t 1 .
62
Если взять обыкновенное дифференциальное уравнение, при этом взяв
t 1, получим x(t ) x(t ) и сопоставим ему разностное: x t ax t 0 на
x(t 1) a1 x(t ) 0 , a1 1 a .
x(t+1)
x(t) Δx(t)
x(t)
x(t–1)
t–1 t t+1
x t x t 1 2x t x t 1 ,
n
при этом мы воспользовались биномом Ньютона a b Cnk a k b n k .
n
k 0
63
преобразование Лапласа ( z -преобразование): F p
e
t
pt
x t z x t F z
t
t
1 – дискретное преобразование Лапласа,
производящая функция.
2.6. Методы решения автономных ЛРУ
Lx t 0,
x 0 x0 ,
x 1 x1 , – задача Коши при любых x0 , , xn 1 ,
x n 1 xn 1
причем задача Коши всегда имеет единственное решение.
x 0 x 1 1.
x 0 c1 x1 0 c2 x2 0 c11 c2 2 c1 c2 1,
0 0
x 1 c1 x1 1 c2 x2 1 c11 c2 2 c11 c2 2 1,
1 1
64
c1 1 c2 , 1 c2 1 c2 2 1, 1 c2 1 c2 2 1, c2 2 1 1 1 ,
1 5 1 5
1
1 1 2 5 1 2 5 5 1
c2 2 , c1 1 c2
2 1 1 5 1 5 5 2 5 2 5
2
5 1
= .
2 5
Числа Фибоначчи:
t t
5 11 5 5 11 5
x t .
2 5 2 2 5 2
x1 t cos t arg ,
t
– тоже простой комплексный корень. Тогда
x2 t sin t arg , a bi , a bi .
t
65
b = a + bi
| |
arg
a
0
-b = a – bi
Проверка: t 2 2 t 1 t t 2 2t 2 t 0 .
66
Решение ЛРУ со специальной правой частью
1) f t t Pk t , k – степень многочлена.
1-й случай
– не является корнем характеристического уравнения, поэтому
x * t t Ak t , где Ak t – многочлен степени k с неопределенными
коэффициентами.
2-й случай
– корень характеристического уравнения кратности j, тогда
x * t t j t Ak t .
находим решение x1 1t 1, x2 3t .
Находим А , подставим x t 1 t 1 3t 1 A , x t 2 t 2 3t 2 A в
уравнение:
t 2 3t 2 A 4 t 1 3t 1 A 3t 3t A = A3t t 2 9 12 t 1 3t = A3t 6 3t ,
1 1 1
находим A , откуда x * (t ) t 3t . Ответ: x t c1 c2 3t t 3t .
6 6 6
Пример. Решим задачу Коши:
67
x t 2 4 x t 1 3x t 3t , t 0,1,2, ,
x 0 1,
x 1 0.
Решение. Подставляем начальное значение: x 0 c1 c2 0 1 ,
c1 1 c2 ,
1 3 3 3
x 1 c1 c2 3 3 0 , 1 2c2 , c2 , c1 1 .
6 1 c2 3c2 2 0, 2 4 4
7 3 t 1 t
И, наконец: x t 3 t3 .
4 4 6
2) f t t Pk t cos t Ql t sin t , составим r cos + i sin ei .
1-й случай
r – не корень характеристического уравнения (т.е. корень кратности 0),
тогда x * t t 0 t Am t cos t Bm t sin t , где m max k , l .
2-й случай
r – корень кратности j, тогда x * t t j t Am t cos t Bm t sin t ,
m max k , l .
e t e t e t e t
3) f t t Pk t ch t Ql t sh t , где sh t ; ch t .
2 2
Подставим m , где m max k , l ,
~ ~ ~ ~
f t t P m t e t Q m t e t t e t P m t t e t Q m t
Pk t Ql t Pk t Ql t
2 2
1 2
2.7. Формула Коши для автономного ЛНРУ (ЛРУ) n-го порядка
Для ЛРУ порядка n:
Lx t a0 x t n a1x t n 1 an1x t 1 an x t f t , t 0,1,2, ,
t 0,1,2, .
Функция Коши C (t ) удовлетворяет условию
Lx t 0, x 0 0, x 1 0, , x n 1 1 .
Lx t 0, x 0 1, x 1 0, … , x n 1 0 ; x2 t : Lx t 0, x 0 0,
x 1 1, … , x n 1 0 ; … ; xn t : Lx t 0, x 0 0, x 1 0, , x n 1 1 .
1
Можно показать, что C t xn (t ) .
a0
Пример. Записать формулу Коши для уравнения второго порядка:
Lx t x t 2 x t f t , t 0,1,2, .
Решение. Функция Коши уравнения
Lx t x t 2 x t 0 , t 0,1,2, ,
t
x2 (t ) = sin . Первым фундаментальным решением является функция x1 (t ) =
2
t
cos . Следовательно, формула Коши для уравнения имеет вид
2
t
t s 1 t t
x t 2 sin f ( s) cos x 0 sin x 1 , t 0,1,2, .
s 0 2 2 2
2.8. Решение автономной линейной однородной системы разностных
уравнений (ЛОСРУ) (ЛРУ)
69
xt 1 Axt , t 0,1,2, ,
где А – невырожденная матрица размера n n , x – n-мерный вектор.
Если матрица А приводится к диагональному виду A = S S 1 (для этого
достаточно, чтобы собственные векторы-столбцы были линейно
независимыми), то
xt At x0 S t S 1 x0 c11t u1 c2 2t u2 ... cn nt un ,
где i – собственные числа матрицы А; ui – соответствующие правые
собственные векторы; ci – постоянные, которые удовлетворяют начальному
условию
c1u1 c2u2 ... cnun x0 .
Таким образом, решение представляет собой комбинацию «гармонических
составляющих» – частных решений, проходящих вдоль собственных векторов
матрицы А.
Матрица называется подобной матрице A, если существует
невырожденная матрица S такая, что = S 1 AS . Преобразование A S 1 AS
называется преобразованием подобия (или просто подобием). Отношение «
подобна A » иногда записывается сокращенно: A . Подобие является
отношением эквивалентности, а именно: оно рефлексивно, симметрично и
транзитивно.
Матрица A называется диагонализуемой, если она подобна диагональной
матрице. Имеется простое условие, обеспечивающее диагонализуемость, –
когда все собственные значения различны.
2.9. Формула Коши для систем неоднородных автономных ЛРУ
Формула Коши для системы неоднородных автономных ЛРУ
xt 1 Axt ft 1 , t 0,1,2, , имеет вид
t
xt C
s 1
f
t s s X t x0 , t 1,2, .
общее решение
частное решение однородного ур-я
неоднородного ур-я
70
Здесь Ct ,s Ct s – матрица-функция Коши, Ct = X t – фундаментальная
71
Список литературы к § 1
73
§ 2. Устойчивость решений ЛОДУ, систем ЛОДУ, ЛРУ
, an = const, a0 0 , Lx f , x 0 x0 , , x 0 xn1
n 1
т.е. a0 , – задача Коши.
y t x t , для любого t 0 .
ε
δ
x0 x(t)
δ ε
74
2) любая траектория (решение) y t обладает свойством: y t x t при
t .
Отрицание асимптотической устойчивости означает, что 1) решение
вообще неустойчиво, 2) решение устойчиво, но не асимптотически. В случае 2)
будем говорить, что решение x t просто устойчиво.
Решение x x(t ) , не обладающее свойством устойчивости, называется
неустойчивым.
Утверждение 2.1.1. Для ЛОДУ все решения устойчивы или неустойчивы
одновременно, в том числе все решения и асимптотически устойчивы
одновременно, следовательно, достаточно для ЛОДУ изучать устойчивость
нулевого решения, следовательно, достаточно изучать устойчивость
однородного уравнения.
ЛОДУ в случае устойчивости решений называется устойчивым ЛОДУ, а в
случае неустойчивости решений называется неустойчивым ЛОДУ.
то ek t 0 .
75
2.2.2. Критерий устойчивости стационарных ЛОДУ
неограниченная функция.
Случай Случай
Re λ < 0 Re λ 0
(ас. уст.) (уст.)
• λ4
• λ7 • λ3
λ2 λ5
• • • • • λ1
λ1 λ4
• λ6 • λ2
• λ3
Рис. 1 Рис. 2
76
Случай Случай Re λ 0,
Re λ > 0 но кратные
(неуст.) корни с Re λ = 0
(неуст.)
• λ4
• λ1 • λ7 • λ1 = λ3
λ3 λ6
• • •
λ5
• λ2 • λ8 • λ2 = λ4
• λ5
Рис. 3 Рис. 4
77
Очевидно, что ЛОДУ x( n ) (t ) x(t ) 0 при любом 0 – неустойчивое
уравнение. Действительно, характеристическое уравнение этого уравнения
имеет вид n 0 , откуда 1 1 n 0 . Заметим, что уравнение
a11 a1n
матрица – A .
a ann
n1
78
Определение 2.3.1. Решение x x(t ) системы уравнений устойчиво по
n , так как в n
все нормы эквивалентны.
79
2.4.2. Критерий устойчивости стационарных систем ЛОДУ
12
Гантмахер Ф.Р. Теория матриц. 4-е изд., доп. М.: Наука, Гл. ред. физ.-мат. лит., 1988. Гл.
VI, § 6; Гл. VII, § 7.
13
Мишина А.П., Проскуряков И.В. Справочная математическая библиотека. Высшая алгебра
(линейная алгебра, многочлены, общая алгебра) / под общ. ред. Л.А. Люстерника, А.Р.
Янпольского. М.: Физматгиз, 1962. Гл. II, § 1, 11.
80
Проблема состоит в том, что одному характеристическому числу 1 может
соответствовать несколько клеток Жордана. Сумма размерностей этих клеток
1 0 0
равна n1 . Например, при n1 3 могут быть случаи: 0 1 0 – три клетки,
0 0 1
1 1 0 1 0 0
0 1 0 – две клетки, 0 1 1 – две клетки. В последних двух случаях
0 1 0 1
0 0
x 0 0 0
Пример 1. Рассмотрим систему
, здесь A , 1 2 0 ,
y 0 0 0
клетки Жордана простые, т.е. имеем простые элементарные делители. Решение
системы: x(t ) x(0) , y(t ) y(0) . Очевидно, что решение ограничено, т.е.
система устойчива по Ляпунову.
x y, 0 1
Пример 2. Рассмотрим систему , здесь A , 1 2 0 .
y 0 0 0
Клетка Жордана непростая, поэтому имеем непростой элементарный делитель.
Решение системы имеет вид x(t ) x(0) ty(0) , y(t ) y(0) , т.е. хотя бы одно
решение не ограничено на полуоси [0, ) , значит, система неустойчива по
Ляпунову.
81
2.4.3. Критерий неустойчивости стационарных систем ЛОДУ
2 1
= 2 4 3 8 6 2 3 0 – система неустойчива,
3 4
так как 1 1, 2 5 .
14
Постников М.М. Устойчивые многочлены. М.: Наука, Гл. ред. физ.-мат. лит., 1981. Гл. 1.
15
Мишина А.П., Проскуряков И.В. Указ. соч. Гл. III, § 4, 4, 5.
82
2.5.1. Схема критерия Рауса–Гурвица
83
Многочлен Pn ( ) степени n 1 называют стандартным полиномом16, если
a0 0 и an 0 . В этом случае полином не имеет нулевых корней.
Многочлен Pn ( ) называют устойчивым многочленом (многочленом
(полиномом) Гурвица), если Re 0 для любого корня этого многочлена.
Через 1 ,..., n обозначим последовательные главные (угловые) миноры
матрицы Гурвица H .
a a0
Матрица Гурвица имеет вид H 1 . Далее находим 1 a1 0 ,
0 a2
16
Демидович Б.П. Лекции по математической теории устойчивости: учеб. пособие / предисл.
В.М. Миллионщикова. 2-е изд. М.: Изд-во Моск. ун-та, 1998. Гл. II, § 9.
84
a2 0 . Асимптотическая устойчивость этого уравнения будет при a0 0 ,
a1 0 , a2 0 .
Пример 3. Пусть n 3 , уравнение имеет вид a0 x a1 x a2 x a3 x 0 .
a1 a0 0
a0 0 . Матрица Гурвица имеет вид H a3 a2 a1 . Далее находим
0 a3
0
a1 a0
1 a1 0 , 2 a1a2 a0 a3 0 , 3 det H 2 a3 (a1a2 a0 a3 )a3 0 .
a3 a2
17
Демидович Б.П. Указ. соч. Гл. II, § 9.
85
получаем два корня: k k ik , k 1 k ik . Поэтому
a1 a0
a3 0 , 2 2 a1a2 a0 a3 0 .
a3 a2
a1 a0 0
Здесь матрица Гурвица имеет вид H a3 a2 a1 .
0 a3
0
86
3. Рассмотрим ЛОДУ 4-го порядка: асимптотическая устойчивость
уравнения имеет место тогда и только тогда, когда
a0 0 , a1 0 , a2 0 , a3 0 , a4 0 ,
a1 a0 0
3 a3 a2 a1 a1a2a3 0 0 0 a12a4 a0a32
0 a4 a3
a1 a0 0 0
a a2 a1 a0
Здесь матрица Гурвица имеет вид H 3 .
0 a4 a3 a2
0 0 0 a4
Распишем
a11 a12 a1n
a21 a22 a2 n
det( E A) n a1 n2 an1 an
87
Далее обозначим
aii aij aik
aii aij
A2 , A3 a ji a jj a jk , …, An det A .
i j a ji a jj i j k
aki akj akk
a a12
коэффициентов A 11 и характеристический многочлен
a21 a22
88
a11 a12 a13
имеет матрицу коэффициентов A a21 a22 a23 и характеристический
a a33
31 a32
a1 a0 A1 1
Выпишем 2 = A1 A2 A3 . Тогда условия Льенара–
a3 a2 A3 A2
действительными коэффициентами a0 0 , a1 , … , an 0 .
Определение 2.6.1. Пусть x(t ) – решение разностного уравнения порядка
n , определяемое начальными условиями x(0) x0 ,…, x(n 1) xn 1 , а y (t ) –
18
Араманович И.Г., Лунц Г.Л., Эльсгольц Л.Э. Функции комплексного переменного.
Операционное исчисление. Теория устойчивости. М.: Наука, Гл. ред. физ.-мат. лит., 1968. Ч.
3, § 5.
89
решение того же уравнения при измененных начальных условиях
y(0) y0 ,…, y(n 1) yn 1 . Решение x(t ) называется устойчивым, если для
любого 0 существует такое ( ) 0 , что из совокупности неравенств
| y0 x0 | ,…, | yn 1 xn 1 | следует неравенство | y(t ) x(t ) | при любом
t n.
Если, кроме того, lim | y(t ) x(t ) | 0 , то решение x(t ) называется
t
асимптотически устойчивым.
Если при сколь угодно малом ( ) 0 неравенство | y(t ) x(t ) | не
выполняется для какого-либо решения y (t ) , то решение x(t ) называется
неустойчивым.
Утверждение 2.6.1. Для ЛРУ все решения устойчивы или неустойчивы
одновременно, следовательно, достаточно для ЛРУ изучить устойчивость
нулевого решения, а значит, достаточно изучить устойчивость однородного
уравнения.
характеристического уравнения 1.
a0 a1 ... an 0 ,
w 1 w 1
лежащие в левой полуплоскости. Корни последнего уравнения являются
корнями многочлена
a0 (w 1)n a1 (w 1)n1 (w 1) ... an ( w 1) n b0 wn b1wn1 ... bn ,
w 1 1
где , w .
w 1 1
91
1
0 1
0
λ w
!
где C – биноминальные коэффициенты, причем считается, что
!( )!
C0 1 и C 0 , если или если 0 .
В частности,
b0 a0 a1 ... an , bn a0 a1 ... (1)n an .
После того как многочлен будет построен, останется применить к нему
теорему Рауса–Гурвица либо критерий Льенара–Шипара.
Пример. Для ЛРУ
a0 x t 2 a1 x t 1 a2 x t 0 , t 0,1, 2,... , a0 0 , a2 0 ,
т.е. для характеристического уравнения
a0 2 a1 a2 0 ,
коэффициенты bk равны: b0 a0 a1 a2 , b1 2(a0 a2 ) , b2 a0 a1 a2 – и
условия Гурвица приводят к неравенствам a0 a1 a2 0 , a0 a2 0 ,
a0 a1 a2 0 .
Действительно, в этом случае получаем уравнение
92
a0 a1 a2 w2 2 a0 a2 w a0 a1 a2 0 .
b0 b1 b2
3a0 a2 0,
2
Или в виде: a0 a0 a2 a1a3 a32 0,
a0 a2 | a1 a3 | 0.
Частный случай для n 4: уравнение четвертого порядка
a0 x(t 4) a1x(t 3) a2 x(t 2) a3 x(t 1) a4 x(t ) 0 ( a0 0, a4 0 )
асимптотически устойчиво тогда и только тогда, когда
a0 a4 0,
3a0 3a4 a2 0,
a0 a2 a4 | a1 a3 | 0,
(a0 a4 ) 2 (a0 a4 a2 ) (a1 a3 )(a1a4 a0 a3 ).
n
Теорема Конa19. Если | a || a
i 1
i 0 | , то нулевое решение уравнения
a x(t n i) 0 ,
i 0
i t 0,1,2, , где a0 0, an 0 , является асимптотически
устойчивым.
19
Cohn A. Über die Anzahl der Wurzeln einer algebraischen Gleichung in einem Kreise //
Mathematische Zeitschrift. 1922. Vol. 14, № 1. P. 110–148.
93
Список литературы к § 2
94
Киселев, Г.И. Макаренко, Е.В. Шикин, В.И. Заляпин. М.: Едиториал УРСС,
2005. Т. 3. 240 с.
11. Лихтарников Л.М. Элементарное введение в функциональные
уравнения. СПб.: Лань, 1997. Гл. 4.
12. Мишина А.П. Справочная математическая библиотека. Высшая алгебра
(линейная алгебра, многочлены, общая алгебра) / под общ. ред. Л.А.
Люстерника, А.Р. Янпольского; А.П. Мишина, И.В. Проскуряков. М.:
Физматгиз, 1962. 300 с.
13. Моисеев Н.Д. Очерки развития теории устойчивости. М.; Л.: Гос. изд-во
техн.-теор. лит., 1949. 663 с.
14. Панюкова Т.А. Основы теории дифференциальных уравнений для
экономистов: учеб. пособие. М.: Кн. дом «ЛИБРОКОМ», 2011. 256 с.
15. Постников М.М. Устойчивые многочлены. М.: Наука, Гл. ред. физ.-
мат. лит., 1981. 176 с.
16. Романко В.К. Курс разностных уравнений: учеб. пособие. М.:
ФИЗМАТЛИТ, 2012. 200 с.
17. Сюдсетер К. Справочник по математике для экономистов / пер. с
норвеж. под ред. Е.Ю.Смирновой; К. Сюдсетер, А. Стрём, П. Берк. СПб.:
Экономич. школа, 2000. X + 229 с.
18. Филиппов А.Ф. Сборник задач по дифференциальным уравнениям.
Ижевск.: НИЦ «Регулярная и хаотическая динамика», 2000. 176 с.
95
§ 3. Характеристики экономического развития20
20
Гранберг А.Г. Динамические модели народного хозяйства: учеб. пособие для студ. вузов,
обуч. по спец. «Экон. кибернетика». М.: Экономика, 1985. Гл. 1, § 2.
96
3.1. Темпы роста и прироста, индексы роста дискретные (базисные
и цепные)
Базисные Цепные
(на постоянной базе) (по переменной базе)
t / 0 x(t ) x(0) t / t 1 x(t ) x(t 1)
– абсолютный (цепной) прирост
– абсолютный (базисный) прирост;
(аналог первой производной), часто
x(t )
t / 0 используют обозначения
x(0)
– (базисный) темп роста, (базисный) t / t 1 x(t 1) x(t ) x(t ) x(t 1) ;
97
3.2.Абсолютное (дискретное) ускорение. Первые и вторые разности
.
t-1 t t+1
98
В последнем случае шаблон имеет следующий вид:
.
t-2 t–1 t
t
Далее определим t – относительное ускорение, индекс роста
t / t 1
абсолютного ускорения, темп прироста абсолютного ускорения. Заметим, что
x(t )
t .
x(t )
99
5) показатель Ляпунова
ln | x(t ) |
[ x] lim .
t t
Очевидно, что ˆ (t ) Ind x(t ) для функции x(t ) x(0)e t . В случае
ln | x(t ) | ln | x(t ) |
существования lim получаем, что [ x] lim . Например,
t t t t
[ x] для функции x(t ) x(0)e t .
Справедливы равенства
T T
0 1 2 3… T
1 T
прогрессии:1 2 ... T T .
2
100
Иногда возникает обратная задача: по заданным T , T ,b найти x(T ) и
T
T
T
x(T ) x(0) TT , x(t ) 2 (2 x(0) (T 1)
t 1
T ,b ).
TT,b1 (1 xT )T ,b xT ,
101
T
x(t )
где xT = T / 0 = t 1
.
x(0)
Здесь была использована формула суммы для геометрической прогрессии.
xt (1 g )t x0 . (2)
В случае когда величина xt имеет переменный темп прироста, равный
g x (t ) , и задано начальное значение x0 , то
t
xt (1 g x ( s))x0 .
s 0
102
Если темп прироста g близок к нулю, можно использовать приближенное
равенство:
ln(1 g x (t )) g x (t ) , 1 g e g .
Пользуясь этим равенством и (1), получаем, что
g x (t ) ln xt 1 ln xt .
В непрерывном времени темп прироста определяется иначе. Пусть x(t) –
некоторая величина в модели с непрерывным временем. Ее темпом прироста
(growth rate) называется величина
x(t )
g x (t ) , (4)
x(t )
где x(t ) – производная по времени. Нетрудно видеть, что темп прироста
представляет собой производную по времени от логарифма ln x(t)
(логарифмическую производную).
Видим, что величина x(t) имеет постоянный темп прироста g x в том и
только в том случае, если она изменяется экспоненциально:
x(t ) e g xt x(0) . (5)
При этом если x(0) > 0, g x > 0, то величина x(t) растет, а если x(0) > 0, g x <
0, то x(t) уменьшается.
Если величина x(t) имеет переменный темп прироста g x (t ) и задано
начальное значение x(0) , то
t
g x ( s ) ds
x(t ) e 0 x(0) .
Заметим, что в непрерывном времени переход от величины x(t ) к её
логарифму вполне естествен. Если величина x(t ) имеет темп прироста g x (t ) , то
угловой коэффициент кривой ln x(t ) равен g x (t ) . В частности, если темп
прироста g постоянен, то график ln x(t) представляет собой прямую с угловым
коэффициентом, равным g x (t ) .
103
Примером может служить зависимость логарифма реального ВВП в США
от времени. Этот график напоминал бы на большом промежутке времени
прямую, если бы не «искажения», обусловленные Великой депрессией и
Второй мировой войной.
Можно заметить, что определения темпа прироста в дискретном времени
(3) и в непрерывном времени (4) аналогичны: (3) может быть записано как
xt
g x (t ) t ,
xt
xt
где xt xt 1 xt , t (t 1) t = 1, а (4) – как g x (t ) lim t .
t 0 x
t
x(t )
темп прироста частного величин x(t), y(t) равен разности их темпов
y (t )
прироста:
g x (t ) g x (t ) g y (t ) .
y
прироста произведения xt yt , то
1 g xy (t ) (1 g x (t ))(1 g y (t )) ;
xt
темп роста частного величин xt , yt равен частному их темпов роста:
yt
1 g x (t )
1 g x (t ) .
y 1 g y (t )
105
Упражнения к § 3
x(t )
темп прироста частного величин x(t), y(t) равен разности их темпов
y (t )
прироста:
g x (t ) g x (t ) g y (t ) .
y
прироста произведения xt yt , то
1 g xy (t ) (1 g x (t ))(1 g y (t )) ;
106
xt
темп роста частного величин xt , yt равен частному их темпов роста:
yt
1 g x (t )
1 g x (t ) .
y 1 g y (t )
107
Список литературы к § 3
108
§ 4. Математические методы исследования динамических
моделей21
Структурные схемы экономико-математических моделей: элементарные
экономические звенья, основные схемы соединения звеньев, законы
преобразования передаточных функций. Модели с обратной связью,
экономические мультипликаторы
4.1. Преобразование Лапласа
f t 0 при t 0 .
21
Батищева С.Э., Каданэр Э.Д., Симонов П.М. Математические модели микроэкономики:
учеб. пособие. 2-е изд., перераб. и доп. Пермь: Перм. гос. ун-т, 2006. Гл. 8, 9.
22
Фихтенгольц Г.М. Курс дифференциального и интегрального исчисления. 3-е изд.,
стереотип. Л.: Физматгиз, 1963. Т. III. С. 438.
109
f(t)
t
0
1, t 0,
t – эта функция позволяет записать любую функцию в
0, t 0
зануленном виде.
f(t)
0 t
f(t)∙η(t)
0 t
110
Если функция удовлетворяет аксиомам 1–3, то все они образуют линейное
пространство, кроме того, функции с одинаковым образуют нормированное
пространство с нормой || f || inf | f (t ) | M e t почти при всех t [0, ) .
M 0
f(t)
0 t
1
Следовательно, lim e pt 0 при a Re p 0 . Итак, 1 , Re p 0 .
t p
Значит, изображение гарантировано (задано) в полуплоскости Re p 0 .
111
bi
0 a
L x X p , x t X p .
x t 0 x t pe pt dt 0 x 0 pX p pX p x 0 .
pt
e
0
Итак, L x pX p x 0 , x pX p x 0 .
Пример. Решим задачу Коши
x t x t 1,
x 0 1.
1
Таким образом, x x pX p 1 X p – это алгебраическое
p
1 1
уравнение первого порядка, X p p 1 1 , откуда получаем X p .
p p
Далее перейдем к оригиналам:
1
X p x t , X p x t (t ) 1.
p
1 1 1,
Проверка:
1 1.
112
4.2. Свойства преобразования Лапласа
Af t Bg t AF p BG p .
Это следует из линейности интеграла.
Вычислим для при Re p :
e t L e t e t e pt dt e p t 1
e
p t
dt
0 0
p 0
1 1
= 0 1 .
p p
1
Итак, получили e t при Re p .
p
2. Запаздывание, смещение (аргумента оригинала)
Пусть f – оригинал, 0 – число, тогда справедлива формула
f t e p F p ,
где f t f t t .
f(t) f(t–τ)
0 τ t
113
1 η(t– τ)
0 τ t
полуоси , ) .
Доказательство. Непосредственно вычисляем:
f t e py
f t (t )dt e pt f t dt
0
( t s , t s , dt ds , t s 0 , t s )
f s ds e p F p .
p s
e
0
f t pF p f 0 .
f g t f s g t s ds f t s g s ds ,
0 0
причем f g t F p G p .
Доказательство. Пусть f , g – оригиналы, т.е.
f t Me t , g t Ne t (*).
t
Рассмотрим e pt
f g t dt e f s g t s dsdt
pt
0 0 0
114
В силу оценок (*) двойной интеграл существует, так как существует
двойной интеграл от модуля, поэтому можно поменять порядок
интегрирования. Изобразим область интегрирования, поменяем пределы.
s
0 t s t
pt
pt
e f ( s ) g (t s ) dt ds f s e g t s dt ds
0 s 0 s
( t s , dt d , t s 0 , t )
p s p
f s e g d ds f s e ds e g d
ps
0 0 0 0
Так как переменные разделились, то повторный интеграл равен
произведению одинарных интегралов:
e ps
f s ds e p g d F p G p .
0 0
F p G p
Свойства свертки
1) Коммутативность: f g g f .
2) Ассоциативность: f gh f g h .
3) Дистрибутивность: f g h f h g h.
Все оригиналы с операциями сложения, умножения на число, свертки
образуют алгебру.
115
4.3. Импульсные (дельта-) функции
(t )
0 t
, если t 0,
Характеристические свойства -функции: (t )
0, если t 0;
(t )dt 1 .
116
(t )
1
0 t
0 1
При t 0 видим, что (t ) 0 . В случае t 0 и t 0 имеем .
t
При t 0 получаем значение 0. Итак, хотя производной в точке t 0 не
существует, но левосторонняя производная равна , поэтому можно условно
считать, что (0) .
Таким образом, -функция – бесконечная скорость, возникает в точке
разрыва функции.
V0
117
Преобразование Лапласа -функции
det
L( ) e (t )dt e (t )dt e p 0 1 , т.е. по определению L( ) 1.
pt pt
0
1. Дельта-функция (t )
0, t 0,
(t ) dt 1
1
118
Таблица 2
№ Название Оригинал Изображение
п/п
1. Линейность A f(t) + B g(t) A F ( p) B G( p)
2. Подобие f(λt), λ > 0 1/λ F(p/λ)
3. Затухание eat f(t) F(p – a)
4. Запаздывание f(t – ) η(t – τ), τ > 0 e - p F(p)
5. Дифференцирование оригинала d f(t)/ dt = f '(t) p F (p) – f (0)
6. Интегрирование оригинала t
f ( s )ds F (p) / p
0
( f g ) (t) =
9. Свертка оригиналов t F (p) G (p)
= f (t s ) g ( s )ds
0
119
задачи. Линейные операторы блоков Н выражают передаточными функциями
Hˆ ( p) . Такие блоки часто называют линейными звеньями. Тем самым звено
моделируют с помощью передаточной функции. Если H – линейный оператор,
то в дальнейшем Hˆ ( p) , Ф(p) или W(p) обозначения передаточных функций
оператора.
H оператор
x(t) y(t)
Ĥ(p),Ф(p),
W(p)
X(p) Y(p)
t
1 1
при y(0) 0 получаем y (t ) x( s)ds , Y ( p) X ( p) , ( p) передаточная
0
p p
функция.
x(t) y(t)
= 1/p
X(p) Y(p)
y(0)
Ty(0)
121
С помощью операционного исчисления находим
y Y ( p) , x X ( p) , y pY ( p) y(0) ,
T pY ( p) y(0) Y ( p) X ( p) ,
(Tp 1)Y ( p) X ( p) Ty(0) ,
откуда
1 T
Y ( p) X ( p) y (0) .
Tp 1 Tp 1
1 T
Здесь X ( p ) – вынужденное движение, y (0) – свободное движение.
Tp 1 Tp 1
1
Отсюда видим, что ( p) – передаточная функция инерционного
Tp 1
звена.
Найдем оригинал для функции y(t). По таблице изображений вычислим
1 Tt t
1 (t Ts )
t t
1 T
e , поэтому e , отсюда y (t ) e
T
x( s)ds y (0)e ,T
1
T(p ) T Tp 1 0
T
T
где первое слагаемое это вынужденное движение, частное решение
неоднородного уравнения, а второе слагаемое это свободное движение, общее
решение однородного уравнения. Полученную формулу часто называют
1 T1 (t s )
формулой Коши. Здесь C (t , s) e функция Коши (весовая функция), а
T
X (t ) TC (t ,0) et / T фундаментальное решение.
x 1
Пример. Пусть x const , тогда x , в частности, x 1 . Отсюда
p p
при y(0) 0 получаем
1 1 A B Ap B(Tp 1)
Y ( p) .
Tp 1 p Tp 1 p (Tp 1) p
Приравниваем числители: 1 ( A BT ) p B , откуда находим, что B 1,
A BT T . Окончательно:
122
T 1 T 1
t
t
Y ( p) e 1 1 e y(t ) .
T T
Tp 1 p Tp 1 p
Если входной процесс представляет собой поток материалов, энергии,
денежных средств и т.д., то внутри звена накапливается определенное
количество этих материально-вещественных элементов, равное разности
входного и выходного процессов. Инерционными звеньями моделируют также
реакцию покупателей в ответ на поступление товара в продажу, ввод основных
производственных фондов в ответ на капиталовложение.
Какова реакция системы на единичное воздействие?
При x 1 получаем
t
y(t ) 1 e T
Выход у на x 1 происходит
1 с инерцией, с запаздыванием.
0 T 2Т 3T t
t
x 1 y(t ) 1 e T
1
0 t
123
Иногда используется другая величина: = 3Т – лаг запаздывания. В этом
3T
1
случае y (3T ) 1 e T
1 e3 1 0,95, z (3T ) = x(3T ) y(3T ) =
e3
1
= 0,05 , т.е. 5%. Через лаг времени капитальные вложения будут почти
e3
полностью реализованы.
x 1
инерционное звено
1 пор. и т.д.
Т 2Т
x y
x 1 x(t T ) y(t ) T
0 T t 0 T t
125
системы равна произведению передаточных функций (ПФ) отдельных звеньев.
Ф1 Ф2 … Фn
1
Отсюда находим: Y ( p) = X ( p) =
(T1 p 1)(T2 p 1)
126
1
= X ( p) модель, выраженная через ПФ. Тогда
TT
1 2 p2
(T1 T2 ) p 1
TT
1 2 p (T1 T2 ) p 1 Y ( p) X ( p) , значит,
2
1 2 y (T1 T2 ) y y x
TT
– модель в виде ЛОДУ второго порядка.
Аналогично возникают инерционные звенья третьего и более высокого
порядков.
распределительный узел
x1 x2 xn
y1 y2 yn
127
x x1 ... xn – инфор-
мационные потоки
X1 x1 X2 x2 Xn xn
…
Ф1(p) Ф2(p) Фn(p)
Y1 y1 Y2 y2 Yn yn
+ материальные потоки
Y=Y1+…+Yn y=y1+…+yn
x(t) y(t)
Ф1(р)
X(p) Y(p)
z
Ф2(р)
Z(p)
Выведем зависимость Y ( p) от X ( p) :
Y ( p) 1 ( p) X ( p) Z ( p) ,
Z ( p) 2 ( p)Y ( p) ,
128
Y ( p) 1 ( p) X ( p) 2 ( p)Y ( p) ,
Y ( p) 1 ( p) 2 ( p) Y ( p) 1 ( p) X ( p) ,
1 ( p) X ( p)
Y ( p) .
1 1 ( p)2 ( p)
Находим, что
1 ( p)
( p) .
1 1 ( p) 2 ( p)
Ф2 = а
y(t)
x(t) Ф1 = 1
+
X(p) Y(p)
u(t) Ф2 = k2
U(p)
x(t) Ф1 = k1 y(t)
X(p) Y(p)
129
4.6.2. Пример (декомпозиции) модели накопленного выбытия (износа)
основных производственных фондов (ОПФ)
S(t)
Ф2(p) =1
S(p)
F(t) 1 S(t)
– Ф1 ( p)
F(p) – Tp S(p)
130
4.6.2. Модель ОПФ без учета выбытия
131
4.6.3. Модель ОПФ с учетом выбытия
s(t)
n Здесь:
S(p)
s(t ) n F (t ) ,
S ( p) n F ( p ) ,
1
I(t) I(t)-s(t) F(t) F ( p) I ( p) S ( p) F (0) .
_ 1/p p
F(p)
I(p) I(p)-S(p)
1
Итак, – передаточная функция всей системы, откуда
pn
1 1
F ( p) I ( p) F (0) .
pn pn
вынужд. движ. своб. движ.
Получим модель в виде изображений. Переход к оригиналам:
1
e nt
pn t
I ( p) I (t ) F (t ) e n (t ) I ( )d e nt F (0) – модель в виде интегрального
F ( p) F (t ) 0
y (t )
Здесь F (t ) – ОПФ с учетом выбытия.
133
Упражнения к § 4
u(t)
k
U(p)
v(0)
–
p
v(t) v1(t)
p
V(p) V1(p)
x(t) y(t)
1
инф. поток
Y(p)
X(p) v(t)
2p v2(t)
V(p) V2(p)
v(0)
–
2p
Рис. 1.
Рис. 2.
134
x(0)
–
p
инф. поток
X(p) Y(p)
x(t) z(t)
3
X(p) Z(p)
Рис. 3.
135
Список литературы к § 4
136
12. Колемаев В.А. Математическая экономика: учеб. для вузов. 3-е
стереотип. изд. М.: ЮНИТИ-ДАНА, 2005. 400 с.
13. Колемаев В.А. Экономико-математическое моделирование.
Моделирование макроэкономических процессов и систем: учеб. для студ.
вузов. М.: ЮНИТИ-ДАНА, 2005. 296 с.
14. Конторов Д.С. Основы физической экономики (Физические аналогии и
модели в экономике) / Д.С. Конторов, Н.В. Михайлов, Ю.С. Саврасов. М.:
Радио и связь, 1999. 184 с.
15. Кугаенко А.А. Основы теории и практики динамического
моделирования социально-экономических объектов и прогнозирования их
развития. 2-е изд. М.: Вузовская книга, 2005. 392 с.
16. Кугаенко А.А. Тринадцать тренажеров по управлению социально-
экономическими процессами. М.: Финансы и статистика, 2001. 236 с.
17. Лопатников Л.И. Экономико-математический словарь: Словарь
современной экономической науки. 5-е изд., перераб. и доп. М.: Дело, 2003. 520
с.
18. Математические методы и модели исследования операций: учебник
для студ. вузов / под. ред. В.А. Колемаева. М.: ЮНИТИ-ДАНА, 2008. 592 с.
19. Мирошник И.В. Теория автоматического управления. Линейные
системы. СПб.: Питер, 2005. 333 с.
20. Пантелеев А.В. Теория управления в примерах и задачах: учеб. пособие
/ А.В. Пантелеев, А.С. Бортаковский. М.: Высш. шк., 2003. 583 с.
21. Форрестер Дж. Основы кибернетики предприятия (Индустриальная
динамика) / пер. с англ.; общ. ред. и предисл. Д.М. Гришиани. М.: Прогресс,
1971. 340 с.
22. Царьков В.А. Динамические модели экономики: теория и практика
экономической динамики. М.: Экономика, 2007. 216 с.
23. Шишмарев В.Ю. Основы автоматического управления: учеб. пособие
для студ. высш. учеб. заведений. М.: Изд. центр «Академия», 2008. 352 с.
137
24. Экономико-математический энциклопедический словарь / гл. ред. В.И.
Данилов-Данильян. М.: Большая Российская энциклопедия: Изд. дом «Инфра-
М», 2003. 668 с.
25. Эртли-Каякоб П. Экономическая кибернетика на практике / сокр. пер.
с нем.; под. ред. К.А. Багриновского. М.: Экономика, 1983. 160 с.
138
§ 5. Динамический производственный оператор валовой
продукции (модель динамики основных производственных
фондов с учетом выбытия)
139
износ. Размер производственного капитала при вводе в эксплуатацию принято
учитывать по первоначальной стоимости, а в процессе эксплуатации, кроме
того, принято учитывать по остаточной стоимости с исключением
амортизации.
Из сказанного следует, что, в принципе, могут быть разработаны четыре
типа динамических моделей производства: модели производства валовой
продукции в зависимости от величины остаточной или первоначальной
стоимости капитала и модели производства чистой продукции в зависимости от
величины остаточной или первоначальной стоимости капитала. Для
практического моделирования основной интерес представляет модель
производства валовой продукции в зависимости от остаточной стоимости
капитала и модель производства чистой продукции в зависимости от
первоначальной стоимости основного капитала. В случае моделирования
производства валовой продукции изменение остаточной стоимости капитала
происходит в соответствии с валовыми инвестициями и износом капитала.
Валовые инвестиции, если они превышают износ, увеличивают капитал, но
если меньше износа – не обеспечивают поддержание капитала на прежнем
уровне. Поэтому модель производства валовой продукции может включаться в
модели систем с расширенным, сужающимся и простым воспроизводством
капитала. Чистые инвестиции создают только прирост капитала.
Следовательно, модель производства чистой продукции может включаться в
модели систем только расширенного или простого воспроизводства капитала.
Если износ капитала в единицу времени обозначить s, то в любой момент
остаточная стоимость капитала F1 , находящегося в эксплуатации, определяется
дифференциальным уравнением
dF1 (t ) dt v(t ) s(t ) ,
где v(t ) – прирост валовых капиталовложений в каждую единицу времени.
Капитал F1 (t ) накапливался с момента начала его создания в соответствии с
текущими инвестициями и износом.
140
Для возмещения износа основного капитала в законодательном порядке
устанавливается возмещение амортизации. Любое предприятие обязано
отчислять из получаемого дохода денежные средства по норме амортизации
для полного возмещения износа основного капитала до первоначальной
стоимости. При норме амортизации n процесс начисления амортизационных
средств на возмещение износа равен nF1 (t ) .
В задачах экономики по росту и расширению производства, по
воспроизводству капитала и т.п. валовые инвестиции, как правило, разделяют
на чистые инвестиции и возмещение износа:
v(t ) i(t ) nF1 (t ) ,
где i (t ) – чистые инвестиции в единицу времени. Будем предполагать, что
износ капитала полностью возмещается за счет средств амортизации
( s(t ) nF1 (t ) ). Тогда в формулу для динамики остаточной стоимости капитала
вместо износа можно подставить процесс амортизации. В результате получим
производную остаточной стоимости капитала в зависимости от валовых
инвестиций:
dF1 (t ) dt nF1 (t ) v(t ) .
Отсюда нетрудно получить модель выпуска валовой продукции в
следующем виде:
dy(t ) dt ny(t ) v(t ) .
141
F1 (0)
1
F1 (0)
2
0 t* t
t*
86%
F1 63% 95%
0 T 2T 3T t
142
v(t ) i(t ) Am(t ) i(t ) nF1 (t ) ,
где i (t ) – интенсивность чистых инвестиций, дающих прирост ОПФ; Am(t ) –
интенсивность амортизационных отчислений на поддержание ОПФ,
Am(t ) nF1 (t ) . Далее из уравнения F1(t ) nF1 (t ) v(t ) с учетом равенства
v(t ) i(t ) nF1 (t ) получим уравнение относительно F1 (t ) : F1(t ) i(t ) . Выразим
F1 (t ) через i (t ) :
t
F1 (t ) F1 (0) i( s)ds .
0
143
отчисления растут по линейному закону. В случае i(t ) 1 и произвольного
начального положительного значения ОПФ F1 (0) получаем:
v(t ) 1 nF1 (0) nt , F1 (t ) F1 (0) t .
144
Упражнения к § 5
Упражнение 1. Модель динамики основных производственных фондов
(ОПФ) с учетом остаточной стоимости F имеет вид v i nF , где v , i –
прирост валовых капиталовложений и чистые инвестиции в единицу времени.
Найти функции F (t ) и i (t ) и построить их графики:
1
1) v(t ) 2 , n , F (0) 15 ;
10
1
2) v(t ) et , n , F (0) 10 ;
5
1
3) v(t ) t , n , F (0) 100 ;
2
1
4) v(t ) et , n , F (0) 36 ;
4
1
5) v(t ) e2t , n , F (0) 18 ;
2
1
6) v(t ) 3 , n , F (0) 7 ;
5
1
t 1
7) v(t ) 1 e 5 , n , F (0) 5 ;
5
1
8) v(t ) 0 , n , F (0) 15 ;
10
1
9) v(t ) et , n , F (0) 10 ;
5
1
10) v(t ) t 2 , n , F (0) 100 ;
2
1
11) v(t ) 1 t , n , F (0) 10 ;
2
1
12) v(t ) t 3 , n , F (0) 10 .
2
145
Список литературы к § 5
146
§ 6. Две динамические непрерывные модели предприятия
F(t)
y (t ) 1
– постоянная фондоотдача ;
F (t ) вр.
1
Багриновский К.А. Модели и методы экономической кибернетики. М.: Экономика, 1973. Гл.
1, § 1.
147
y (t ) ден.
u (t ) aF (t ) – интенсивность функции инвестиций , где 0 a 1 –
y (t ) вр.
F (t )
Отсюда получаем y(t ) F (t ) , т.е. u (t ) a F (t ) a F (t ) , или, по-
F (t )
другому,
a
u (t ) y (t ) .
Структурная схема модели имеет следующий вид:
y (0)
p
u(t)
ap
U(p)
148
a
так как u (t ) y (t ) ;
a a 1
Y ( p) I ( p) pY ( p) y (0) y (0) ;
p p
a 1
Y ( p)(1 a) I ( p) y (0) y(0) ;
p p p
1 1 a 1
Y ( p) I ( p) y (0) ;
1 a p p 1 a
1 y (0)
Y ( p) I ( p) ;
1 a p p
1
Y ( p) .
(1 a) p p 1 a
Получили модель для выпуска продукции:
1 1
Y ( p)
I ( p) y (0) ;
1 a p p
1 1 1 1
модель для ОПФ: F ( p) Y ( p) , F ( p) I ( p) F (0) ;
1 a p p
1 1 1
F ( p) .
(1 a) p p 1 a
1
Модель в виде ЛОДУ: y (t ) I (t ) , F (t ) I (t ) .
1 a 1 a
1-й случай. Пусть I (t ) (t ) I ( p) 1, y(0) 0 , тогда
1 1
Y ( p)
1 y (t ) , y (t ) F (t ) .
1 a p 1 a p 1 a 1 a
1
Так как 0 a 1, то 1 , причем y(t ) , когда a 0 .
1 a
149
y(t)
1 a
0 t
т.е.
y (t ) y (0) .
1 a
y(t)
y (0)
1 a
1 a
0 t
t
Если на входе интенсивность капвложений постоянна,
y (t ) то фонды растут по линейному закону.
1 a
y(0) = 0
0 t
y(t)
y (t ) y (0)
1 a
y(0) 0
0 t
151
Замечание. Как показано в работе2, инвестиций должно быть столько,
чтобы функция y (t ) росла не медленнее, чем экспонента. Иначе возникает
некорректность модели.
F(0)
2
Чадов А.Д., Кузнецов К.Б. Об области применимости первой модели Багриновского //
Информационные системы и математические методы в экономике: сб. науч. тр. / Перм. гос.
ун-т. Пермь, 2008. С. 128–139.
3
Багриновский К.А. Модели и методы экономической кибернетики. М.: Экономика, 1973. Гл.
1, § 1.
152
1 1
Y ( p) I ( p) U ( p)) F (0) ; Y ( p) F ( p) ;
pn pn
Y ( p) I ( p) U ( p) F (0) ; U(p) aY ( p) ;
pn
Y ( p) I ( p) aY ( p) F (0) ;
pn
a
Y ( p) 1 I ( p) F (0) ;
pn pn
1 1
Y ( p)
a p n
I ( p) F (0)
p n a p n
I ( p) F (0)
1
pn pn
y (0)
I ( p) .
p n a p n a
1 F (0)
Отсюда Y ( p) F ( p) , F ( p) I ( p) .
p n a p n a
Запишем модель в виде ЛОДУ: y(t ) (n a ) y(t ) I (t ) – устойчивость
этого уравнения зависит от коэффициента n a .
Пример 1. Пусть I (t ) (t ) I ( p) 1, причем F (0) 0 ( y(0) 0 ), тогда
Y ( p) , откуда y(t ) e( a n )t .
p n a p ( a n)
Итак, получили, что решение y (t ) начальной задачи
153
y(t) При a > n в дальнейшем система
развивается сама без внешних
капиталовложений.
e( a n ) t
0 t
y(t)
0 t
y(t)
0 t
1-й сл.
y(t)
2-й сл.
0 t
154
Пример 2. В случае I(t) = 1 рассуждения аналогичны, но действия
сложнее. Имеем I (t ) 1 I ( p) 1/ p , откуда
1 y (0) ( a n ) t
Y ( p)
y (t ) y (0)
n a
e .
p n a p p n a n a
u(t) a y(t)
p
U(p) Y(p)
F(0)
1
y (t ) y (0) 1 e n1t .
n1 n1
1 F(t) y(t)
I(t)
pn T1 p 1
I(p) F(p) Y(p)
a
p
y (0)
p
156
Видим, что a0 = T1 > 0, a1 = T1n + 1 a > 0, a2 = n > 0, т.е. ЛОДУ второго
порядка устойчиво асимптотически.
u(t) y(t)
a
U(p) Y(p)
I(t) 1 F(t)
y(t)
pn T1 p 1
I(p) F(p) Y(p)
инерционное звено
F(0)
выпуска продукции
157
Пусть I (t ) (t ) , y(0) 0 , y(0) 0 , тогда решение уравнения y (t ) будет
пропорционально второму фундаментальному решению y2 , т.е.
y(t ) C (t ) y2 (t ) , где C (t s) – функция Коши этого уравнения.
T1
Из теории устойчивости известно, что ЛОДУ второго порядка
асимптотически устойчиво тогда и только тогда, когда коэффициенты этого
порядка уравнения положительны. В нашем случае a0 = T1 > 0, a1 = = T1n + 1 >
0, причем коэффициент a2 = n – a меняет знак. В зависимости от изменения
этого знака соответственно ведет себя система.
Уравнение асимптотически устойчиво, если и только если a2 = n – a > 0,
т.е. норма выбытия > нормы отчисления на развитие от ОПФ. Графики
фундаментальных решений в этом случае ведут себя таким образом:
y(0) = 0,
y1(t) y(0) = 1, y'(0) = 1
y'(0) = 0
1
y2(t)
0 t
0 t
y2(t)
0 t
158
В случае n = a норма выбытия равна норме отчислений на развитие от
ОПФ.
Уравнение модели примет вид T1 y (T1n 1) y I , отсюда видим, что
( n )t
1
1
фундаментальные решения имеют вид: y1 (t ) 1, y2 (t ) = 1 e
1 =
T1
n
T1
1
n n1
1 e ( nn1 )t , n1 1/ T1 .
y1=1
1
n n1
y2 (t )
1
n n1
1 e ( nn1 )t
0 t
159
y(0)
= y(0) y(0)
1
n n1
1 e ( nn1t ) y (0)
n n1
при t . Отсюда видим,
что Ind y 0 .
y(0)
y(0)
y (0)
n n1
y(0)
0 t
I(p) Y(p)
F(p)
160
Дифференциальным уравнением модели является линейное
дифференциально-разностное уравнение вида
y(t ) ny (t ) a y (t T ) I (t T ), t 0,
y ( ) ( ), если 0.
Здесь – некоторая начальная функция для y .
b1= – a
b1= a1
/(2T)
1/T
-1/T a1= n
0
b1= – a1
161
Упражнения к § 6
7 7
3 5
7) I (t ) e , a , , F (0) 10 ;
t
7 7
1 3
8) I (t ) (t ) , a , , F (0) 2 ;
2 4
1 1
9) I (t ) e , a , , F (0) 10 ;
t
2 2
1 1
10) I (t ) e , a , , F (0) 10 ;
t
2 2
1 1
11) I (t ) 1 , a , , F (0) 2 ;
2 2
1 1
12) I (t ) (t ) , a , , F (0) 2 .
2 2
162
1 1 2
3) I (t ) e , a , n , , y(0) 1 ;
t
10 5 3
1 1 2
4) I (t ) e , a , n , , y(0) 1 ;
t
5 5 3
1 1 2
5) I (t ) e , a , n , , y(0) 1 ;
t
5 5 3
3 1 1
6) I (t ) 10 (t ) , a , n , , F (0) 21;
8 8 3
1 1 1
7) I (t ) 3 , a , n , , F (0) 2 ;
3 8 3
2 1 2
8) I (t ) t , a , n , , F (0) 100 ;
5 8 5
1 1 3
9) I (t ) 3 , a , n , , F (0) 2 ;
3 8 4
2 1 7
10) I (t ) t , a , n , , F (0) 100 ;
5 8 10
1 1 3
11) I (t ) 3 (t ) , a , n , , F (0) 2 ;
3 8 4
1 1 1
12) I (t ) (t ) , a , n , , F (0) 2 .
2 8 2
163
Список литературы к § 6
164
12. Кугаенко А.А. Тринадцать тренажеров по управлению социально-
экономическими процессами. М.: Финансы и статистика, 2001. 236 с.
13. Потапов В.Д. Разработка и инструментальная реализация модели
формирования оптимальных планов текущей деятельности для предприятий с
серийным производством: автореф. дис. … канд. экон. наук: 08.00.13. Пермь,
2006. 27 с.
14. Сиразетдинов Т.К. Динамическое моделирование экономических
объектов. Казань: Изд-во «Фэн», 1996. 224 с.
15. Соколовский Л.Е. Модели оптимального функционирования
предприятия. М.: Наука, 1980. 174 с.
16. Форрестер Дж. Основы кибернетики предприятия (Индустриальная
динамика) / пер. с англ., общ. ред. и предисл. Д.М. Гришиани. М.: Прогресс,
1971. 340 с.
17. Хачатрян С.Р. Методы и модели решений экономических задач: учеб.
пособие / С.Р. Хачатрян, М.В. Пинегина, В.П. Буянов. М.: Экзамен, 2005. 384 с.
18. Чадов А.Д. Об области применимости первой модели Багриновского /
А.Д. Чадов, К.Б. Кузнецов // Информационные системы и математические
методы в экономике: сб. науч. тр. / Перм. гос. ун-т. Пермь, 2008. С. 128–139.
19. Царьков В.А. Динамические модели экономики: теория и практика
экономической динамики. М.: Экономика, 2007. 216 с.
20. Ширяев В.П. Экономико-математическое моделирование управления
фирмой / В.П. Ширяев, И.А. Баев, Е.В. Ширяев. 2-е изд., испр. и доп. М:
КомКнига, 2006. 224 с.
21. Эртли-Каякоб П. Экономическая кибернетика на практике / сокр. пер.
с нем.; под. ред. К.А. Багриновского. М.: Экономика, 1983. 160 с.
165
§ 7. Две динамические дискретные модели развития предприятия
t
Непрерывная модель имеет вид F (t ) I ( s)ds F (0) , т.е. F (t ) I (t ) .
0
F(t) = F(0)
F(0)
0 1 2 3 t
Пусть I (t ) 1 , тогда
F (t 1) F (t ) 1,
F (1) F (0) 1,
– арифметическая прогрессия с разностью один.
F (2) F (1) 1,
166
7.1.2. Модель ОПФ с учетом выбытия
t
Непрерывная модель имеет вид F (t ) e n (t s ) I ( s)ds e nt F (0) , или
0
F (t ) nF (t ) nI (t ),
– модель в виде ЛОДУ.
TF (t ) F (t ) I (t )
F(0)
0
1 2 3 t
167
Найдем F*(t) = С частное решение неоднородного уравнения. В этом
уравнении t = (1 n)t, 1 – n = 1, F (t ) C1 (1 n)t общее решение
однородного ЛРУ.
F (t ) F (t ) F *(t ) C1 (1 n)t C .
Надо найти константы C1, C : C берется из уравнения. Для начала запишем
F (t 1) C1 (1 n)t 1 C (1 n) F (t ) n ,
откуда
F(t + 1) = C1(1 – n)t+1 + C = (1 – n) (C1(1 – n)t + C) + n,
C = (1 – n)C + n,
в итоге C = 1.
Отсюда F(t) = C1(1 – n)t + 1, константу C1 находим из равенства F(0) = С1 +
1, C1 F (0) 1. Получим F (t ) ( F (0) 1)(1 n)t 1 .
F(0)
Если F(0) > 1,
то траектория убывает.
Если F(0) = 1, то F(t) 1.
Если 0 < F(0) < 1,
то траектория возрастает.
1
F (t ) 1
t
0 t
Запишем по-другому формулу решения ЛРУ:
F (t ) ( F (0) 1)(1 n)t 1 F (0)(1 n)t 1 (1 n)t .
Другая запись этой модели:
168
F (t ) F (t ) F (t ) F (t 1) левая (восходящая) разность,
F (t ) F (t 1) nF (t ) nI (t 1) , (n 1) F (t ) F (t 1) nI (t 1) ,
1 n
F (t ) F (t 1) I (t 1) , t =1, 2, …
n 1 n 1
1
Неравенство 1 всегда верно при 0 n 1 , значит, это уравнение
n 1
асимптотически устойчиво.
t
1
Находим F (t ) F (t ) F *(t ) , F (t ) C1 *
, F (t ) C .
n 1
t
1
Пусть I (t ) 1 , тогда F (t ) C1 C . Как и в предыдущем случае, из
n 1
1 n
уравнения для F* C находим C C , откуда C 1. Далее
n 1 n 1
F (0) C1 1 1 C1 1, C1 1 F (0) , откуда
t
1
F (t ) (1 F (0)) 1.
n 1
Ответ аналогичный, только коэффициенты разные.
169
Дискретная модель получается следующим образом. Составим
аппроксимацию y(t ) y(t ) y(t 1) y(t ) правая (нисходящая) разность,
получим
y (t 1) y (t ) I (t ) (*)
1 a
модель в виде ЛРУ, в случае I (t ) const – арифметическая прогрессия,
I (t ) разность.
1 a1
Решим уравнение (*). Рассмотрим однородное уравнение y(t 1) y(t ) 0 .
I(t)=1
0 1 2 3 t
170
Видим 1t , откуда 1 1 , поэтому y*(t) = t1C частное
1 a 1 a
решение неоднородного уравнения.
Общее решение неоднородного уравнения имеет вид
y(t ) y (t ) y *(t ) C1 t Ct C1 Ct ,
поэтому y(t 1) C1 C (t 1) . Из уравнения видим y(t 1) y (t ) C 1,
1 a
откуда C .
1 a
y(0)
0 1 2 3 t
171
7.1.4. Первая модель предприятия. Учет выбытия
n
Ранее получили, что y(t ) y(t ) I (t ) непрерывная модель в
1 a 1 a
виде ЛОДУ первого порядка.
n 1
Введем обозначения: n1 n , 1 . Заметим, что I (t )
1 a 1 a 1 a
мгновенный прирост ОПФ без учета выбытия.
Положим y(t ) y(t ) y(t 1) y(t ) , получим ЛРУ первого порядка
y(t 1) (1 n1 ) y(t ) 1I (t ) , t 0, 1, 2, ... .
Пусть 1 1 n1 0 , в случае 1 0 требуется дополнительное
исследование. Пусть далее I (t ) 1 , тогда y (t ) y (0) 1 1t 1 , t 0, 1, ...
n1 n1
Тогда справедливы следующие утверждения.
Асимптотическая устойчивость уравнения имеет место тогда и только
тогда, когда |1| < 1, что эквивалентно неравенству |1 n1 | 1, которое
n
эквивалентно неравенствам 0 n1 2 , т.е. эквивалентно неравенству
1 a
n 2(1 a) .
1
В этом случае при I (t ) 1 получаем, что y (t ) при t .
n1
172
I (t ) 1 I (t ) 0
y(t)
y(0)
1/n1
y(0) y(t)
0 1 2 3 t 0 1 2 3 t
y(t)
y (t )
0 1 2 3 t
173
Очевидно, что y (t ) монотонно растет при t , т.е. происходит рост
выпуска продукции даже при отсутствии внешних инвестиций.
Просто устойчивость имеет место тогда и только тогда, когда | 1 | 1 , что
эквивалентно равенству n 2(1 a) . Тогда при I (t ) 1 получаем
y(t ) y(0) 1t , при I (t ) 0 получаем y(t ) y(0) , т.е. в отсутствие внешних
инвестиций выпуск постоянный.
174
Неравенство a > n означает = 1 + a n > 1 (|| > 1), это эквивалентно
неустойчивости и это эквивалентно тому, что даже в отсутствие внешних
инвестиций выпуск монотонно растет к бесконечности.
Неравенство a = n означает = 1, это эквивалентно простой
устойчивости решений однородного уравнения, и это означает, что в
отсутствие внешних инвестиций выпуск постоянный: y(t) = y(0).
Общее решение однородного уравнения имеет вид
y (t ) C11t C1 (1 a n)t .
1
y (t ) 1t y (0) 1t y (0)
1 1 1 1 n a n a
1t y (0) (1 1t ) ,
n a
где 1 1 a n .
175
- положительная величина,
y(0) - отрицательная величина,
n a n a
n > a n < a
при при
n a
y(0)
y(0)
0 1 2 3 4 t 0 1 2 3 4 t
Так как 1t et ln 1 , то 1t растет как экспонента при 1 1 , т.е. при a n .
Пусть a = n, тогда уравнение принимает вид y(t + 1) – y(t) = I(t), откуда
при I (t ) 1 выводим y(t) = y(0) + t.
y(0)
0 1 2 t
176
Положим y(t ) y(t ) y(t 1) y(t ) . Для второй производной возможны
три аппроксимации. Например, y "(t ) y(t ) y(t ) y(t 1)
t1 t t+1
t t+1 t+2
t–2 t–1 t
178
Упражнения к § 7
179
Упражнение 7. Записать вторую дискретную модель Багриновского с
учетом выбытия и с учетом дискретного технологического лага производства в
виде уравнения y(t 1) (n 1) y(t ) a y(t 2) I (t 2) , t 2,3,4, , 0 n 1,
0 a 1, 0 1. Исследовать на устойчивость.
180
Список литературы к § 7
181
§ 8. Динамические имитационные модели простых
экономических объектов
182
активам Kп является показателем их маржинальной доходности Eд ,
измеряемой в относительных единицах размерностью [1/год]:
Eд yдс / Kп [1/год] (1)
либо в процентах годовых:
Eд yдс 100%/ Kп [%/год]. (2)
Отношение потока платежей yпр за привлекаемые (заемные) ресурсы K пр
183
капитала. Эти расходы относят к так называемым добавленным затратам.
Поэтому величина Em при увеличении капитала будет уменьшаться.
Другой вид расходов – это затраты на комплектующие и расходные
материалы и другие затраты, содержащие налог на добавленную стоимость.
Эти расходы растут обратно пропорционально времени оборачиваемости
капитала об и пропорционально капиталу и образуют поток перенесенной
стоимости yпc :
yпc Kn / об . (7)
Введем обозначение Eоб 1/ об . Отношение потока перенесенной
стоимости к капиталу также может быть измерено в единицах размерности
[1/год] по формуле
Eоб yпc / Kn [1/год]. (8)
Как видим, величина расходности Eоб служит еще и мерой скорости
оборота капитала.
184
Рис. 1. Блок-схема преобразования капитала предприятия в поток перенесенной стоимости
У пс
ya ( p) Kфн ( p) / сл . (14)
187
Допустим, как и ранее, равенство Kфт (t ) Kфн1(t ) , тогда Kфн ( p) Kфн / p , в
188
преобразуется оператором интегрирования в функцию капитала, поступающего
на вход системы. В результате получим блок-схему на рис. 4.
189
Сопоставив эти формулы, можно утверждать, что в модели кругооборота
капитала постоянная времени инерционного звена равна времени
оборачиваемости капитала: об .
Из этого утверждения следует, что трехзвенную блок-схему на рис. 4
можно заменить эквивалентной однозвенной схемой с коэффициентом
передачи W ( p) 1/ об . Однозвенная блок-схема представлена на рис. 5.
190
Модель налогообложения на основе рентабельности затрат. В блок-
схеме обратной связи прибыль является функцией от независимых векторов yпс
и yдз , а также передаточных коэффициентов трех звеньев. Из блок-схемы
несложно убедиться в справедливости уравнений:
yп pyсп ;
yсп p( yпс yдз ндс yдз ндс yп ) .
191
Рис. 7. Модель налогообложения с маржинальной рентабельностью: yпс – перенесенная
стоимость; yдз – добавленные затраты; yндс – налог на добавленную стоимость; ндс –
выплаты по налогу на добавленную стоимость; yсп – себестоимость продукции; p –
рентабельность по себестоимости продукции; yв yсп (1 p) – выручка от реализации
продукции/услуг (валовой операционный доход); yп – прибыль от реализованной продукции
Из блок-схемы на рис. 7 несложно вывести следующие соотношения:
yсп yпс yдз yдз ндс yп ндс ;
yп yв yсп yпс (1 pм ) ( yпс yдз yдз ндс yп ндс ) .
После решения системы из двух уравнений относительно потока прибыли
yп получим
yпс pм
yп yдз . (23)
1 ндс
Уравнения (22) и (23) неприменимы к банковским услугам, так как
большинство банковских услуг реализуются без НДС.
192
Упражнения к § 8
193
6. Модель налогообложения с маржинальной рентабельностью.
194
Список литературы к § 8
195
§ 9. Моделирование экономической динамики предприятия
(опыт применения методов теории автоматического регулирования для
моделирования экономических объектов)
9.1. Введение
1
См. работы в списке литературы к параграфу 9.
196
Эффективность воспроизводства капитала E , равная отношению прибыли
yП (t ) к текущей величине капитала K (t ) , определяется соотношением
E yП (t ) / K (t ) p / . (2)
Эти два уравнения являются фундаментальными, определяющими
динамику роста экономической системы.
Однако с переходом к рыночной экономике все шире начинает
применяться в практике экономических расчетов такая категория, как
добавленная стоимость. Добавленная стоимость служит исходным показателем
как для расчета цены продукции, доходности капитала, так и для формирования
рентабельности. В связи с этим ниже разработана обобщенная модель с
использованием этого показателя. Это вносит определенную новизну, а также
необходимость разработки иной нетрадиционной системы показателей
эффективности использования ресурсов предприятия, адекватных рыночной
экономике. Последнее обстоятельство не означает отрицания показателей
эффективности, основанных на рентабельности предприятия, тем более что все
показатели оказываются взаимосвязанными.
197
KСП – объем собственных средств (активов) предприятия;
пр Kпр / KСП – коэффициент привлечения ресурсов, тогда
K П KСП K П KСП (1 пр ) . (4)
198
Интегральные характеристики предприятия, такие как ежемесячные доход,
расходы, средняя величина собственного основного и оборотного капитала,
привлеченных ресурсов и другие, являются достаточно стабильными
агрегированными экономическими показателями его деятельности,
характерными для детерминированной системы. Очевидно, что такое
предприятие правомерно рассматривать как детерминированную систему,
параметры которой имеют свою динамику, взаимозависимость и свою
траекторию изменения во времени, которую возможно, в принципе, как
прогнозировать, так и планировать.
На рис. 1 представлена обобщенная блок-схема динамической модели
воспроизводства производственного капитала предприятия (основного и
оборотного), построенная на основе методов, применяемых в теории
автоматического регулирования.
200
звено капитализации чистой прибыли с коэффициентом передачи ЧП ,
преобразующее чистую прибыль yЧП в прирост производственного капитала
предприятия K П .
В блок-схеме в наглядном виде зафиксирована следующая система
уравнений:
производственный капитал ( K П (t ) ) = капитал при t 0 ( K П 0 ) +
добавленная стоимость ( y ДС ).
Особенностью модели является свойство ее саморазвития после подачи
собственного начального капитала K П 0 . При выполнении равенства K П 0 0
система остается в режиме покоя, все другие функции равны нулю. После
подачи на вход модели K П 0 0 система переходит из состояния покоя в
состояние динамического развития. Величина каждого значения будет
изменяться со временем. Характер траектории изменения значения будет
зависеть от параметров операторов, входящих в блок-схему модели.
Для любого значения из блок-схемы можно вычислить аналитические
выражения траектории их изменения во времени.
Не будем останавливаться на процедуре математических выкладок.
Прежде чем перейти к уравнениям значения, введем следующие обозначения:
Д ЧП (1 p ) , (7)
201
где Д – это коэффициент капитализации добавленной стоимости,
показывающий долю добавленной стоимости, направляемую на увеличение
капитала;
Е Д p ДС / об , (8)
yЧП (1 Д )( p ДС / об ) K П 0 (t ) ; (12)
yЦП yПС y ДС ( K П 0 (1 p ДС ) / об ) (t ) ; (13)
K П (t ) K П 0 ( (t ) 1) ; (14)
K П (t ) K П 0 (t ) . (15)
Рассмотрев совместно уравнения (11) и (15), получим
Е Д p ДС / об y ДС / K П (t ) . (16)
ЕП (1 Д ) Е Д .
202
Полученные уравнения (10)–(16) наглядно иллюстрируют
экспоненциальный характер траектории роста экономики предприятия.
Введем обозначение темпа роста капитала предприятия:
dK П (t ) /(dtK П ) . (17)
Из совместного рассмотрения (8), (9), (16) и (17) получим
Д Е Д Д p ДС / об . (18)
Таким образом, темп роста капитала равен произведению коэффициента
капитализации добавленной стоимости и коэффициента добавленной
стоимости, деленному на время оборачиваемости капитала предприятия.
Для понимания экономической динамики важно четко представлять, что
для экономического роста требуется одновременное выполнение двух условий:
1) Д 0 ; 2) Е Д 0 .
Расширенное воспроизводство капитала невозможно реализовать, если
эффективность использования капитала Е Д или коэффициент капитализации
203
Выражение в скобках разложим в ряд Маклорена, оставив первые два
члена, в результате будем иметь
(t ) (1 Д p ДС )t / .
об
(19)
В показателе степени теперь имеем число циклов оборачиваемости
капитала, равное t / об , а в скобках – процент наращивания капитала за один
цикл.
Дальнейшее упрощение с применением разложения в ряд Маклорена уже к
выражению (20):
(t ) 1 Д p ДСt / об 1 Д ЕД t . (20)
Равенство (20) позволяет вычислить рост капитала по «формуле сложных
процентов», а равенство (21) – по «формуле простых процентов». Оба
равенства можно рассматривать как частные случаи, вытекающие из общей
модели воспроизводства капитала в экономических процессах производства и
реализации продукции предприятия.
205
где
K пр ( s) – коэффициент передачи прямого канала;
206
Для вычисления коэффициента затрат Д суммируются расходы на оплату
труда, платежи в налоговые органы, в различные фонды, процентные расходы
за привлеченные ресурсы. Полученная сумма делится на величину добавленной
стоимости. Величина коэффициента капитализации чистой продукции ЧП
зависит от ежемесячных платежей дивидендов, спонсорских расходов и других
платежей, уменьшающих долю чистой продукции, направляемую на
увеличение собственного капитала предприятия. При отсутствии этих платежей
ЧП 1 , при этом коэффициент капитализации добавленной стоимости
становится равным Д 1 p .
207
Список литературы к § 9
208
Приложения
Самостоятельная работа № 1
Решение автономных ЛОДУ, систем ЛОДУ, ЛРУ. Устойчивость
автономных ЛОДУ, систем ЛОДУ, ЛРУ
209
x(t ) x(t ) 4et , x(t ) x(t ) 2 x(t ) et ,
13. 31.
x(0) 4, x(0) 3. x(0) 1, x(0) 0.
x(t ) 3x(t ) 2 x(t ) tet , x(t ) 2 x(t ) 2 x(t ) 0,
14. 32.
x(0) 1, x(0) 2. x(0) 1, x(0) 1.
x(t ) 4 x(t ) 2sin(t ), x(t ) 6 x(t ) 9 x(t ) 0,
15. 33.
x(0) 1, x(0) 0. x(0) 1, x(0) 1.
x(t ) 2 x(t ) 5 x(t ) tet , x(t ) x(t ) 6 x(t ) 2,
16. 34.
x(0) 0, x(0) 0. x(0) 1, x(0) 0.
x(t ) 2 x(t ) 2 x(t ) tet , x(t ) 9 x(t ) 2 t ,
17. 35.
x(0) 0, x(0) 0. x(0) 0, x(0) 1.
x(t ) 2 x(t ) x(t ) et , x(t ) 4 x(t ) 4t ,
18. 36.
x(0) 1, x(0) 0. x(0) 1, x(0) 0.
210
x(t ) x(t ) y (t ), x(0) 2,
8.
y(t ) x(t ) y (t ), y (0) 0.
x(t ) 2 x(t ) y(t ), x(0) 1,
9.
y(t ) x(t ) 2 y (t ), y(0) 3.
x(t ) 4 x(t ) 6 y (t ), x(0) 0,
10.
y(t ) 2 x(t ) 3 y (t ) t , y (0) 0.
x(t ) y (t ) 3et , x(0) 1/8,
11.
y(t ) x(t ) 2e , y (0) 5/8.
3t
x(t ) y (t ) et , x(0) 0,
12.
y(t ) x(t ) e , y (0) 0.
t
x(t ) y (t ), x(0) 2,
14.
y(t ) x(t ), y(0) 5.
x(t ) x(t ) 2 y (t ), x(0) 1,
15.
y(t ) x(t ) 4 y (t ) 1, y (0) 1.
x(t ) 2 y(t ) 3t , x(0) 2,
16.
y(t ) 2 x(t ) 4, y (0) 3.
211
x(t ) 2 x(t ) 4 y(t ) cos(t ), x(0) 0,
22.
y(t ) x(t ) 2 y (t ) sin(t ), y(0) 0.
x(t ) 7 x(t ) y (t ), x(0) 1,
23.
y(t ) 2 x(t ) 5 y (t ), y (0) 1.
x(t ) x(t ) y (t ) 3t 2 / 2, x(0) 0,
24.
y(t ) 4 x(t ) 2 y (t ) 4t 1, y(0) 0.
x(t ) y (t ), x(0) 1,
25.
y(t ) x(t ), y (0) 1.
x(t ) 2 y(t ) 3t , x(0) 0,
26.
y(t ) 2 x(t ) 4t , y (0) 0.
214
x(t 2) 2 x(t 1) 2 x(t ) 8,
10.
x(0) 1, x(1) 0.
x(t 2) x(t ) t ,
11.
x(0) 0, x(1) 1.
x(t 2) 2 x(t 1) x(t ) 1,
12.
x(0) 0, x(1) 0.
x(t 2) 2 x(t 1) x(t ) t ,
13.
x(0) 0, x(1) 0.
x(t 2) 2 x(t 1) x(t ) 32 3t ,
14.
x(0) 0, x(1) 0.
x(t 2) x(t ) sin(t ),
15.
x(0) 0, x(1) 1.
x(t 2) 4 x(t 1) 4 x(t ) 3t ,
16.
x(0) 0, x(1) 0.
x(t 2) 5x(t 1) 6 x(t ) 2 4t ,
17.
x(0) 0, x(1) 0.
x(t 2) 3x(t 1) 4 x(t ) (1)t ,
18.
x(0) 0, x(1) 0.
2 x(t 2) 5 x(t 1) 2 x(t ) cos( t / 3),
19.
x(0) 0, x(1) 0.
x(t 2) 3x(t 1) 2 x(t ) 2t ,
20.
x(0) 0, x(1) 0.
x(t 2) x(t ) 4t ,
21.
x(0) 0, x(1) 4.
x(t 2) 2 x(t 1) 5 x(t ) t ,
22.
x(0) 0, x(1) 0.
x(t 2) 3x(t 1) 2 x(t ) t ,
23.
x(0) 0, x(1) 0.
215
x(t 2) 2 x(t 1) 5 x(t ) t 1,
24.
x(0) 0, x(1) 0.
x(t 2) 2 x(t 1) 2 x(t ) t ,
25.
x(0) 0, x(1) 0.
x(t 2) 2 x(t 1) x(t ) (1)t ,
26.
x(0) 0, x(1) 0.
x(t 2) 9 x(t ) sinh(t ),
27.
x(0) 0, x(1) 0.
x(t 2) 9 x(t ) sinh(t ),
28.
x(0) 0, x(1) 0.
x(t 2) 4 x(t 1) 5 x(t ) t 1,
29.
x(0) 0, x(1) 0.
x(t 2) x(t 1) 2 x(t ) t 1,
30.
x(0) 0, x(1) 0.
x(t 2) 2 x(t 1) x(t ) (1)t ,
31.
x(0) 0, x(1) 0.
x(t 2) x(t ) sinh(t ),
32.
x(0) 1, x(1) 0.
x(t 2) x(t 1) 6 x(t ) 3t ,
33.
x(0) 0, x(1) 0.
x(t 2) 3x(t 1) 2 x(t ) 2t ,
34.
x(0) 0, x(1) 0.
x(t 2) 7 x(t 1) 10 x(t ) 5t ,
35.
x(0) 0, x(1) 0.
x(t 2) 4 x(t 1) 3x(t ) 3t ,
36.
x(0) 0, x(1) 0.
216
Самостоятельная работа № 2
Исследование на устойчивость автономных ЛОДУ, систем ЛОДУ, ЛРУ
2. x( IV ) 4 x 7 x 6 x 2 x 0
3. x 5x 9 x 5x 0
4. x( IV ) 2 x x 2 x 2 x 0
5. x( IV ) 7 x 17 x 17 x 6 x 0
6. x 3x 12 x 10 x 0
7. x( IV ) 5x 18x 34 x 20 x 0
217
23. x 13x 43x 51x 40 x 12 x 0
(V ) ( IV )
29. x 4 x x x x 0
(V ) ( IV )
218
2 1 1 1 2 0 2 1 -2
10. 1 0 1 19. 0 1 1 28. -1 0 0
3 1 2 0 0 1 1 1 -1
0 4 1 2 1 1 0 8 0
11. 0 0 1 20. 1 2 1 29. 0 0 2
0 4 0 1 1 2 2 8 2
0 1 1 1 1 5 3 1 1
12. 1 1 0 21. 0 2 1 30. 1 5 1
1 0 1 0 0 3 1 1 3
1 0 1 2 1 0 21 8 19
13. 0 2 1 22. 1 2 0 31. 18 7 15
0 1 1 1 3 1 16 6 15
2 1 1 1 2 2 1 2 2
14. 1 0 1 23. 1 4 2 32. 2 1 2
3 1 2 1 5 3 3 2 3
3 2 2 3 3 1 2 1 1
15. 3 1 1 24. 3 2 2 33. 1 0 1
1 2 0 1 2 0 1 1 0
0 1 1 0 1 1 2 1 2
16. 1 0 1 25. 1 1 0 34. 1 0 2
2 2 1 1 0 1 2 0 3
0 1 1 4 2 2 2 0 1
17. 1 0 1 26. 1 3 1 35. 1 1 0
2 2 3 3 3 1 3 1 1
0 1 0 2 1 0 1 0 3
18. 0 0 0 27. 0 2 1 36. 1 2 0
0 0 2 0 0 2 0 1 1
219
Упражнение 3. Исследовать на устойчивость нулевое решение ЛРУ
сведением характеристического уравнения и исследовать его методом Рауса–
Гурвица.
1. x(t 3) x(t 2) 8x(t 1) 12 x(t ) 0, t 0,1, 2,...
2. 11x(t 4) 8x(t 3) 8x(t 2) 4 x(t 1) x(t ) 0, t 0,1, 2,...
3. x(t 4) x(t 3) x(t ) 0, t 0,1, 2,...
4. 12 x(t 4) 3x(t 3) 2 x(t 2) 2 x(t 1) 2 x(t ) 0, t 0,1,2,...
5. 7 x(t 4) 4 x(t 3) 30 x(t 2) 4 x(t 1) 3x(t ) 0, t 0,1, 2,...
6. x(t 4) 2 x(t 3) 4 x(t 2) 2 x(t 1) 5x(t ), t 0,1, 2,...
7. x(t 4) x(t 2) 2 x(t 1) 2 x(t ) 0, t 0,1, 2,...
8. x(t 4) x(t ) 0, t 0,1, 2,...
9. x(t 4) 2 x(t 3) 3x(t 2) 2 x(t 1) x(t ) 0, t 0,1, 2,...
10. 2 x(t 4) 5x(t 3) 5x(t 2) 2 x(t ) 0, t 0,1, 2,...
11. x(t 4) 4 x(t 2) 3x(t ) 0, t 0,1, 2,...
12. 4 x(t 3) 2 x(t 2) x(t ) 0, t 0,1,2,...
13. 27 x(t 3) 27 x(t 2) 12 x(t 1) 2 x(t ) 0, t 0,1,2,...
14. 4 x(t 3) 10 x(t 2) 4 x(t 1) 3x(t ) 0, t 0,1,2,...
15. 4 x(t 3) 3x(t 2) x(t 1) x(t ) 0, t 0,1,2,...
16. 4 x(t 3) 2 x(t 2) 4 x(t 1) 3x(t ) 0, t 0,1,2,...
17. 27 x(t 3) 9 x(t 2) 2 x(t ) 0, t 0,1,2,...
18. 4 x(t 3) x(t 2) 2 x(t 1) x(t ) 0, t 0,1,2,...
19. 2 x(t 3) x(t 2) x(t ) 0, t 0,1,2,...
20. x(t 3) 5x(t 2) 3x(t 1) x(t ) 0, t 0,1,2,...
21. 2 x(t 3) x(t 2) x(t ) 0, t 0,1,2,...
22. 27 x(t 3) 9 x(t 2) 2 x(t ) 0, t 0,1,2,...
23. 2 x(t 4) x(t 3) x(t 2) x(t 1) x(t ) 0, t 0,1, 2,...
220
24. 2 x(t 4) 2 x(t 3) x(t 2) 2 x(t 1) x(t ) 0, t 0,1, 2,...
25. 2 x(t 4) 3x(t 3) 4 x(t 2) 2 x(t 1) x(t ) 0, t 0,1, 2,...
26. x(t 4) x(t 3) 4 x(t 2) x(t 1) x(t ) 0, t 0,1,2,...
27. x(t 4) 2 x(t 3) 4 x(t 2) 5x(t 1) 5x(t ) 0, t 0,1,2,...
221
Самостоятельная работа № 3
Вопросы
W(p)
X(p) Y(p)
Y0 ( p)
W ( p)
x(t) y(t)
W(p)
X(p) Y(p)
Y0 ( p)
W ( p)
Y p = W p X p + Y0 p свободное движение определяется функцией
X p ?
222
Вопрос 8. Верно ли, что в модели
x(t) y(t)
W(p)
X(p) Y(p)
Y0 ( p)
W ( p)
W1(p)
x(t ) 1 y (t )
p
y(0)
есть интегрирующее звено?
Вопрос 11. Верно ли, что модель
x(t )
1 y (t )
p
y(0)
223
x(t ) 1 y (t )
pn
y(0)
W1(p)
x(t ) 1 y (t )
p bp c
2
x(t ) 1 y (t )
p 2 bp c
( p n) 2
y 224
можно представить в виде двух параллельно соединенных инерционных
звеньев первого порядка?
Вопрос 17. Верно ли, что реакция y t инерционного звена первого
1 x(t)
y(t)
225
Самостоятельная работа № 4
Модель динамики основных производственных фондов
226
ОГЛАВЛЕНИЕ
Предисловие 3
Введение 5
Список литературы к введению 38
227
4.6. Системы кибернетического типа (обратная связь). Экономический
мультипликатор. Пример экономического мультипликатора с постоянным
лагом…………………………………………………………………………………. 128
Упражнения к § 4…………………………………………………………………… 134
Список литературы к § 4…………………………………………………………… 136
228
§ 9. Моделирование экономической динамики предприятия (опыт
применения методов теории автоматического регулирования для
моделирования экономических объектов)………………………….……......... 196
9.1. Введение………………………………………………………………………… 196
9.2. Капитальные ресурсы и финансовые потоки………………………………… 197
9.3. Динамическая пятизвенная модель предприятия…………………….……… 198
9.4. Упрощенный расчет динамики………………………………………………... 203
9.5. Обоснование модели на основе теории автоматического
регулирования………………………………………………...…………………….. 204
9.6. Возможности практического применения модели………………………...… 206
Список литературы к § 9…………………………………………………………… 208
Приложения………………………………………………………...……………… 209
Самостоятельная работа № 1. Решение автономных ЛОДУ, систем ЛОДУ,
ЛРУ. Устойчивость автономных ЛОДУ, систем ЛОДУ,
ЛРУ……………………………………………………………………….………….. 209
Самостоятельная работа № 2. Исследование на устойчивость ЛОДУ, систем
ЛОДУ, ЛРУ…………………………………………………………………….…… 217
Самостоятельная работа № 3. Вопросы…………………………..……………… 222
Самостоятельная работа № 4. Модель динамики основных производственных
фондов…………………………………………………………………………….…. 226
229
Учебное издание
Экономико-математическое моделирование
Часть первая
Учебное пособие
Редактор Л. А. Богданова
Корректор Л. А. Семицветова
Компьютерная верстка: Л. Н. Калинина, П. М. Симонов, М. В. Дудина
Издательский центр
Пермского государственного
национального исследовательского университета
614990, г. Пермь, ул. Букирева, 15