Будь умным!


У вас вопросы?
У нас ответы:) SamZan.net

тематическое моделирование в экономике Тема 1

Работа добавлена на сайт samzan.net:


PAGE  23

математическое моделирование в экономике

Тема 1. Математические модели в экономике

§1. Общие понятия

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

Практическими   задачами    экономико-математического моделирования являются:

•   анализ функционирования и развития экономических объектов и процессов;

•  экономическое прогнозирование, предвидение развития экономических процессов;

•  выработка управленческих решений  на всех уровнях хозяйственной деятельности.

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

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

- эмерджентность;

- массовый характер экономических явлений и процессов;

- динамичность экономических процессов;

- случайность и неопределенность в развитии экономических явлений;

- невозможность изолировать протекающие в экономических системах явления и процессы от окружающей среды;

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

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

§2. Этапы экономико-математического моделирования

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

- постановка экономической проблемы, ее качественный анализ;

- построение математической модели;

- математический анализ модели;

- подготовка исходной информации;

- численное решение;

- анализ численных результатов, их интерпретации и применение.

Рассмотрим перечисленные этапы экономико-математического моделирования подробнее.

1. Постановка экономической проблемы и ее качественный анализ.

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

2. Построение математической модели.

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

3. Математический анализ модели.

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

4. Подготовка исходной информации.

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

5.  Численное решение.

Этот этап включает разработку алгоритмов численного решения задачи, подготовку программ на компьютере и непосредственное проведение расчетов; при этом значительные трудности вызываются большой размерностью экономических задач. Обычно расчеты на основе экономико-математической модели носят многовариантный характер. Многочисленные модельные эксперименты, изучение поведения модели при различных условиях возможно проводить благодаря высокому быстродействию современных вычислительных средств. Численное решение существенно дополняет результаты аналитического исследования, а для многих моделей является единственно возможным.

6. Анализ численных результатов, их интерпретация и применение.

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

Отметим, что верификация модели – это проверка правильности структуры (логики) модели; валидация модели — проверка соответствия данных, полученных на основе модели, реальному процессу.

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

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

Следует иметь в виду, что далеко не во всех случаях данные, полученные в результате экономико-математического моделирования, могут использоваться непосредственно как готовые управленческие решения. Они скорее могут быть рассмотрены как «консультирующие» средства. Принятие управленческих решений остается за человеком. Таким образом, экономико-математическое моделирование является лишь одним из компонентов (пусть очень важным) в человеко-машинных системах планирования и управления экономическими системами.

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

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

Тема 2. Моделирование экономической динамики

§1. Экономические ряды динамики

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

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

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

уровни.

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

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

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

Если во временном ряду проявляется длительная («вековая») тенденция изменения экономического показателя, то говорят, что имеет место тренд. Таким образом, под трендом понимается изменение, определяющее общее направление развития, основную тенденцию временных рядов. В связи с этим экономико-математическая динамическая модель, в которой развитие моделируемой экономической системы отражается через тренд ее основных показателей, называется трендовой моделью.

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

Предположим, имеется временной ряд, состоящий из п уровней:

у1, у2,…, yn.

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

• тренд, составляющие которого будем обозначать Ut, t = l,2,...,n;

• сезонная компонента, обозначаемая через Vt, t = 1,2,...,n;

•  циклическая компонента, обозначаемая через Ct, t= 1,2,..., n;

• случайная компонента, которую будем обозначать εt, t= 1,2,..., n.

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

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

Тренд, сезонная и циклическая компоненты называются регулярными, или систематическими компонентами временного ряда. Составная часть временного ряда, остающаяся после выделения из него регулярных компонент, представляет собой случайную, нерегулярную компоненту. Она является обязательной составной частью любого временного ряда в экономике, так как случайные отклонения неизбежно сопутствуют любому экономическому явлению. Если систематические компоненты временного ряда определены правильно, что как раз и составляет одну из главных целей при разработке трендовых моделей, то остающаяся после выделения из временного ряда этих компонент так называемая остаточная последовательность (ряд остатков) будет случайной компонентой ряда, т.е. будет обладать следующими свойствами:

•  случайностью колебаний уровней остаточной последовательности;

• соответствием   распределения   случайной   компоненты нормальному закону распределения;

•  равенством математического ожидания случайной компоненты нулю;

•  независимостью значений уровней случайной последовательности, то есть отсутствием  существенной  автокорреляции.

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

Данная проверка осуществляется с использованием  ряда статистических критериев.

§2. Предварительный анализ и сглаживание временных рядов экономических показателей

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

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

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

§3. Расчет показателей динамики развития экономических процессов

Расчет проводится на основе статистического анализа одномерных временных рядов экономической динамики. Для статистического анализа одномерных временных рядов экономических показателей вида

у1, у2,…, yn

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

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

•  по территории ввиду изменения границ региона, по которому собираются статистические данные;

•  по кругу охватываемых объектов по подчинению или форме собственности ввиду перехода, например, части предприятий данного объединения в другое объединение;

•  по временным периодам, когда, например, данные за различные годы приведены по состоянию на разные даты;

•  уровней, вычисленных в различном масштабе измерения;

• уровней ряда из-за различий в структуре совокупности, для которой они вычислены.

Возможны и другие причины несопоставимости.

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

                                   (1)

где yi i-й уровень временного ряда (i = 2,3, ..., n); индекс k = 1,2, ..., n-1 определяет начальный уровень и может быть выбран любым в зависимости от целей исследования: при k = 1 получаются цепные показатели, при h = i-1 получаются базисные показатели с начальным уровнем ряда в качестве базисного и т. д.

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

                                        (2)

В частности, средний абсолютный прирост за весь период наблюдения для данного временного ряда равен

                                          (3)

и характеризует среднюю скорость изменения временного ряда.

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

Коэффициент роста для i-го периода вычисляется по формуле

                              (4)

Коэффициент прироста равен

                     (5)

На практике чаще применяют показатели темпа роста и темпа прироста:

                                  (6)

                (7)

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

Абсолютное значение одного процента прироста определяется как отношение абсолютного прироста Δyi к темпу прироста в процентах .

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

       (8)

Соответственно средний темп прироста определяется как

                        (9)

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

Важной характеристикой временного ряда является также средний уровень ряда. В интервальном ряду динамики с равноотстоящими во времени уровнями расчет среднего уровня ряда производится по формуле простой средней арифметической (здесь и далее суммирование ведется по всем периодам наблюдения):

                                (10)

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

                             (11)

Для моментного ряда с равноотстоящими уровнями средняя хронологическая рассчитывается по формуле:

                (12)

где  п — число уровней ряда.

Средняя хронологическая для моментного временного ряда с разноотстоящими во времени уровнями вычисляется по формуле:

    (13)

Здесь п – число уровней ряда, a ti — период времени, отделяющий i-й уровень ряда от (i+1)-го уровня.

§4. Тренд-сезонные экономические процессы и их анализ

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

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

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

Будем рассматривать тренд-сезонный временной ряд у1,…, уТ порождаемый аддитивным случайным процессом:

уt=Ut+Vtt,                                           (14)

где Ut — тренд;

Vt — сезонная компонента;

εt — случайная компонента;

Т — число уровней наблюдения.

Относительно Ut предполагается, что это некоторая гладкая функция. Сезонная компонента Vt имеет период То: (То = 12 для ряда месячных данных; То = 4 — для ряда квартальных данных).

Кроме того, известно, что То нацело делит Т, т.е. Т = m·То, m — целое число. Очевидно, если То — число месяцев или кварталов в году, то mчисло лет, представленных во временном ряду.

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

Задачи, возникающие при исследовании сезонных временных рядов:

1) определение наличия во временном  ряду тренда;

2) выявление наличия во временном  ряду сезонных колебаний;

3) фильтрация компонент ряда;

4) анализ динамики сезонной волны;

5)  исследование факторов, определяющих сезонные колебания;

6)  прогнозирование тренд-сезонных процессов.

тЕМА 3. лИНЕЙНОЕ ПРОГРАММИРОВАНИЕ

§1. Понятие линейного программирования

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

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

          Формы записи задачи линейного программирования:

Общей задачей линейного программирования называют задачу

                              (1)

при ограничениях

                              (2)

                                 (3)

                                   (4)

                                                (5)

- произвольные                                                   (6)

где - заданные действительные числа; (1) – целевая функция; (1) – (6) –ограничения;

- план задачи.

§2. Графический способ решения ЗЛП

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

                                    (7)

                                    (8)

                                               (9)

Дадим геометрическую интерпретацию элементов этой задачи. Каждое из ограничений (8), (9) задает на плоскости некоторую полуплоскость. Полуплоскость – выпуклое множество. Но пересечение любого числа выпуклых множеств является выпуклым множеством. Отсюда следует, что область допустимых решений задачи (7) - (9)  есть выпуклое множество.

        Перейдем к геометрической интерпретации целевой функции. Пусть область допустимых решений ЗЛП – непустое множество, например многоугольник .

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

Найдём частные производные целевой функции по  и

          (10)

          (11)

       Частная производная (10)  ((11)) функции показывает скорость ее возрастания вдоль данной оси. Следовательно,  и скорости возрастания  соответственно вдоль осей  и . Вектор  называется градиентом функции. Он показывает направление наискорейшего возрастания целевой функции:

        Вектор (- ) указывает направление наискорейшего убывания целевой функции. Его называют антиградиентом.

Вектор перпендикулярен к прямым  семейства  

       Из геометрической интерпретации элементов ЗЛП вытекает следующий порядок ее графического решения.

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

§3. Симплексный метод

       Общая идея симплексного метода (метода последовательного улучшения плана) для решения ЗЛП состоит:

  1.  умение находить начальный опорный план;
  2.  наличие признака оптимальности опорного плана;
  3.  умение переходить к нехудшему опорному плану.

Симплексный метод может быть реализован и «вручную», однако лучше воспользоваться компьютерными возможностями, например, возможностями Excel.

§4. Двойственность задач

          Понятие двойственности рассмотрим на примере задачи оптимального использования сырья. Пусть на предприятии решили рационально использовать отходы основного производства. В плановом периоде появились отходы сырья m видов в объемах единиц . Из этих отходов, учитывая специализацию предприятия, можно наладить выпуск n видов неосновной продукции. Обозначим через  норму расхода сырья i-го вида на единицу j-й  продукции, - цена реализации единицы j-й продукции (реализация обеспечена). Неизвестные величины задачи: объемы выпуска j-й продукции, обеспечивающие предприятию максимум выручки.

Математическая модель задачи:

         (1)

        (2)

                                   (3)

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

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

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

Эти требования формализуются в виде следующей ЗЛП.

Требование покупающей организации – минимизация покупки:            (4)

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

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

        (5)

       По смыслу задачи оценки не должны быть отрицательными:

                         (6)

Переменные    называют двойственными оценками или объективно обусловленными оценками.

         Задачи (1)-(3) и (4)-(6) называют парой взаимно двойственных ЗПЛ.

         Между прямой и двойственной задачами можно установить следующую взаимосвязь:

  1.  Если прямая задача - на максимум, то двойственная к ней - на минимум, и наоборот.
  2.  Коэффициенты  целевой функции прямой задачи являются свободными членами ограничений двойственной задачи.
  3.  Свободные члены  ограничений прямой задачи являются коэффициентами целевой функции двойственной.
  4.  Матрицы ограничений прямой и двойственной задач являются транспонированными друг к другу.
  5.  Если прямая задача на максимум, то ее система ограничений представляется в виде неравенств типа . Двойственная задача решается на минимум, и ее система ограничений имеет вид неравенств типа .
  6.  Число ограничений прямой задачи равно числу переменных двойственной, а число ограничений двойственной — числу переменных прямой.
  7.  Все переменные в обеих задачах неотрицательны.

     Теорема 1. Для любых допустимых планов  и прямой и двойственной ЗЛП справедливо неравенство , т.е.

 (7) – основное неравенство теории двойственности.

      Теорема 2. (критерий оптимальности Канторовича)

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

       Теорема 3. (малая теорема двойственности)

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

§5. Транспортная задача

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

В т пунктах отправления А1, А2, ...,Ат, которые в дальнейшем будем называть поставщиками, сосредоточено определенное количество единиц некоторого однородного продукта, которое обозначим ai (i = 1, 2, ..., т). Данный продукт потребляется в п пунктах В1, В2, ..., Вп, которые будем называть потребителями; объем потребления обозначим bj (j = 1, 2, ..., п). Известны расходы на перевозку единицы продукта из пункта Ai в пункт Bj, которые равны cij и приведены в матрице транспортных расходов С = (cij).

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

Обозначим количество продукта, перевозимого из пункта Ai в пункт Вj, через xij. Совокупность всех переменных xij для краткости обозначим X , тогда целевая функция задачи будет иметь вид

а ограничения выглядят следующим образом:

Первые условия означают полное удовлетворение спроса во всех пунктах потребления; вторые - определяют полный вывоз продукции от всех поставщиков.

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

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

Если же условие баланса не выполняется, то имеет место открытая транспортная задача.

Тема 4. Нелинейное программирование

§1. Задача нелинейного программирования

Задача математического программирования

max (min) z = f(х);                                  (1)                                                                                                                                                                                                                                                                                                     

φi(x ){≤,=,≥} bi,   (i = l,…, m),                 (2)

x ≥ 0,   x = (x1,... ,xn),                              (3)

в которой либо целевая функция (1), либо ограничения (2-3), либо и то и другое нелинейны, называется нелинейной. Нелинейные задачи составляют широкий класс настолько сложных задач, что до сих пор невозможно разработать общие методы, подобные симплекс-методу в линейном программировании, которые позволяли бы решать любые нелинейные задачи. Но, несмотря на отсутствие универсальных методов, разработаны способы решения отдельных специальных классов задач, и прежде всего задач с выпуклыми (вогнутыми) функциями f(х) и φi(x ).

В теории выпуклого программирования в качестве основной обычно рассматривается задача минимизации выпуклой функции п переменных z = f(х) при ограничениях φi(x ) ≤ 0,   (i = l,…, m),  x ≥ 0, где функции φi(x ) предполагаются выпуклыми.

§2. Метод множителей Лагранжа

Метод множителей Лагранжа является классическим методом решения задач математического программирования (в частности выпуклого). К сожалению, при практическом применении метода могут встретиться значительные вычислительные трудности, сужающие область его использования.

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

Рассмотрим классическую задачу оптимизации

  (4)

Эта задача выделяется из задачи (1)- (3) тем, что среди ограничений  нет неравенств, нет условий неотрицательности переменных, их дискретности, т < п и функции f(х) и φi(x ) непрерывны и имеют частные производные по крайней мере второго порядка.

Классический подход к решению задачи (4) дает систему уравнений (необходимые условия), которым должна удовлетворять точка х*, доставляющая функции f(х) локальный экстремум на множестве точек, удовлетворяющих ограничениям задачи (для задачи выпуклого программирования найденная точка х* будет одновременно и точкой глобального экстремума).

Предположим, что в точке х* функция f(х) имеет локальный условный экстремум и ранг матрицы равен т. Тогда необходимые условия запишутся в виде:

                        (5)

где

            (6)

есть функция Лагранжа; λ1,... , λm – множители Лагранжа.

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

Можно указать следующий порядок решения задачи (5) методом множителей Лагранжа:

1) составить функцию Лагранжа (6);

2) найти частные производные функции Лагранжа по всем переменным х1,... , хп, λ1,... , λm и приравнять их нулю. Тем самым будет получена система (5), состоящая из п + т уравнений. Решить полученную систему (если это окажется возможным!) и найти таким образом все стационарные точки функции Лагранжа;

3) из стационарных точек, взятых без координат λ1,... , λm, выбрать точки, в которых функция f(х) имеет условные локальные экстремумы при наличии ограничений. Этот выбор осуществляется, например, с применением достаточных условий локального экстремума. Часто исследование упрощается, если использовать конкретные условия задачи.

§3. Градиентные методы

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

Рассмотрим задачу максимизации нелинейной дифференцируемой функции f(х). Суть градиентного поиска точки максимума х* весьма проста: надо взять произвольную точку хо и с помощью градиента , вычисленного в этой точке, определить направление, в котором f(х) возрастает с наибольшей скоростью, а затем, сделав небольшой шаг в найденном направлении, перейти в новую точку x1. Потом снова определить наилучшее направление для перехода в очередную точку х2 и т. д.

Таким образом, надо построить последовательность точек хо, x1,...  так, чтобы она сходилась к точке максимума х*.

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

Движение из точки xk в новую точку хk+1 осуществляется по прямой, проходящей через точку xk и имеющей уравнение

 

где λk — числовой параметр, от которого зависит величина шага.

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

§4. Теорема Куна - Таккера

В теории нелинейного программирования центральное место занимает теорема Куна — Таккера, обобщающая классический метод множителей Лагранжа на случай, когда в нелинейной задаче, помимо ограничений-равенств, содержатся также и неравенства. В частности, для задачи выпуклого программирования: минимизировать z = f(x) при ограничениях φi(x ) ≤ 0,   (i = l,…, m),  x ≥ 0, где все функции выпуклые, теорема Куна — Таккера устанавливает связь между оптимальным решением задачи и седловой точкой функции Лагранжа для этой задачи:

(7)

Точка (х*, λ* ) называется седловой точкой функции (7), если n-мерная точка х* является точкой минимума функции L(x, λ ), а m-мерная точка λ — точкой максимума функции L(x*, λ* ), так что для всех х и λ выполняется неравенство

(8)

В этом смысл теоремы Куна — Таккера. Сама же теорема формулируется следующим образом.

Теорема  (Куна — Таккера). Предположим, что существует вектор х ≥ 0, такой, что φi(x ) ≤ 0,   (i = l,…, m). Тогда необходимым и достаточным условием оптимальности вектора х*, принадлежащего допустимой области, является существование такого вектора λ* , что для всех х ≥ 0 и λ ≥ 0 имеет место неравенство (8).

ТЕМА 5. дИНАМИЧЕСКОЕ ПРОГРАММИРОВАНИЕ

§1. Основные понятия

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

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

Рассмотрим несколько типичных задач, для решения которых естественным является применение метода динамического программирования.

Задача перспективного планирования. Планируется деятельность группы N промышленных предприятий Пi (i = 1,…, N) на период в t (t = 1,…, Т) хозяйственных лет. В начале периода на развитие системы предприятий выделены какие-то средства К, которые должны быть распределены между предприятиями. В процессе деятельности предприятия вложенные в него средства частично амортизируются. Каждое предприятие за год приносит доход, зависящий от вложенных средств, часть которого отчисляется в фонд предприятий. В начале каждого хозяйственного года имеющиеся средства перераспределяются между предприятиями. Возникает задача определения объема средств в начале каждого года, которые нужно выделить каждому предприятию, чтобы суммарный чистый доход за Т лет был максимальным. Это типичная задача динамического программирования. Здесь процесс принятия решения разбивается на Т шагов. Управление им заключается в начальном распределении и последующих перераспределениях средств иt = {иit}, где иit — объем средств, выделенных i-му предприятию в начале to года. Для описания динамики системы вводится вектор состояния хt = {хit}, где хit — состояние i-го предприятия на начало to года. В свою очередь состояние каждого предприятия хit является вектором, компонентами которого служат трудовые ресурсы, основные фонды, финансовое положение и т.д., т.е. хit = { хikt }. Вектор управления — это функция состояния системы на начало соответствующего года: ut = ut(xt-1). Начальное состояние системы х° может быть заданным. Целевой функцией будет суммарная прибыль объединения за Т лет. Если zt — прибыль за t-й год, то получим задачу

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

Задача об оптимальном управлении поставками. В различных областях народного хозяйства возникает задача определения момента подачи партии поставки и ее объема. С размещением заказов связаны некоторые фиксированные затраты К, не зависящие от величины заказываемой партии, а зависящие только от факта заказывания. С содержанием материальных ресурсов связаны затраты, пропорциональные остатку нереализованной продукции на конец интервала. Пусть Т — промежуток планирования. Обозначим через νt интенсивность потребления ресурса в t-м интервале. Состояние системы будем описывать величиной остатка нереализованной продукции на конец интервала xt. Начальное x0 и конечное xT состояния системы можно считать заданными. Для бесперебойности потребления поставками нужно управлять. Обозначим через u = {ut} вектор управления, координаты которого суть величины поставок в начале соответствующих интервалов. Очевидно, что вектор управления есть функция состояния на начало интервала. Из множества возможных управлений требуется выбрать такое, при котором достигается минимум издержек на заказывание и содержание материальных ресурсов. Если Stиздержки содержания единицы продукции в t-м интервале, то функция цели примет вид

где

Состояние системы опишется соотношением xt = xt-1+ utvt (t =1,…, Т). На состояние системы может быть наложено ограничение, связанное с надежностью снабжения: , где  — величина некоторого страхового запаса, гарантирующего с заданной надежностью от сбоев в системе. Объединение ограничений, налагаемых на состояние системы и вектор управления, обозначим Ω. Получим задачу

§2.Особенности задач динамического программирования

1. Рассматривается система, состояние которой на каждом шаге определяется вектором xt. Дальнейшее изменение ее состояния зависит только от данного состояния xt и не зависит от того, каким путем система пришла в него. Такие процессы называются процессами без последействия.

2. На каждом шаге выбирается одно решение ut, под действием которого система переходит из предыдущего состояния xt-1 в новое xt. Это новое состояние является функцией состояния на начало интервала xt-1 и принятого в начале интервала решения ut, т.е.

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

4.  На векторы состояния и управления могут быть наложены ограничения, объединение которых составляет область допустимых решений u принадлежит Ω.

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

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

Управление — это воздействие, переводящее систему из начального состояния в конечное. Для многих экономических задач не известно начальное либо конечное состояние, а известна область X0 или Хт,, которой эти точки принадлежат. Тогда допустимые управления переводят точки из области Хо в ХТ.

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

§3. Принципы динамического программирования

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

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

Таким образом, одним из условий применимости метода динамического программирования является возможность разбиения процесса оптимизации решения на ряд однотипных шагов (этапов), каждый из которых планируется отдельно, но с учетом состояния системы на начало этапа и последствий принятого решения. Однако из этого правила есть исключение. Среди всех шагов существует один, который может планироваться без учета последствий. Это последний шаг. Он может быть изучен и спланирован сам по себе наилучшим (в смысле выбранного критерия) образом, поскольку за ним нет больше этапов. Отсюда получаем одну из специфических особенностей динамического программирования: всю вычислительную процедуру программирования целесообразно разворачивать от конца к началу. Раньше всех планируется последний N-й шаг, за ним (N-1)-й и т. д. Но как найти оптимальное управление иN на N-м шаге, если оно определяется не только целью управления, но и состоянием системы на начало этого шага? Сделать это можно на основе предположений об ожидаемых исходах предшествующего, но еще не исследованного этапа, т.е. о значениях xN-1.

Для каждого возможного исхода xN-1 на (N1)-м этапе находим оптимальное управление на N-м этапе. Такой набор оптимальных управлений, зависящих от возможных исходов предыдущего этапа, называется условно-оптимальным решением u*n(xn-i). Завершив анализ конечного этапа, рассматривают аналогичную задачу для предпоследнего этапа, требуя, чтобы функция цели достигала экстремального значения на двух последних этапах вместе. Это дает условно-оптимальное решение на предпоследнем этапе u*n-1(xn-2), т. е. делаются всевозможные предположения о том, чем кончился предыдущий (N-2)-й шаг, и для каждого из предположений находится такое управление на (N - 1)-м шаге, при котором эффект за последние два шага (из них последний уже оптимизирован) будет максимален. Тем самым мы найдем для каждого исхода (N - 2)-го шага условно-оптимальное управление на (N1)-м и условно-оптимальное значение функции цели на последних двух шагах. Проделав такой поиск условно-оптимальных управлений для каждого шага от конца к началу, найдем последовательность условно-оптимальных управлений u1*(xo), u2*(x1),... ,u*N(xN-1).

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

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

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

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

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

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

§4. Функциональные уравнения Беллмана

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

Дадим математическую формулировку принципа оптимальности для задач с аддитивным критерием оптимальности. Для простоты будем считать, что начальное хо и конечное хТ состояния системы заданы. Обозначим через z1(x0, u1) значение функции цели на первом этапе при начальном состоянии системы хо и при управлении u1, через z2(x1, u2) - соответствующее значение функции цели только на втором этапе, ..., через zN(xN-1, uN)— на N-м этапе, ..., через zn(xn~i,un) — на Nэтапе. Очевидно, что

            (1)

Надо найти оптимальное управление u* = (u1*, ..., u*N), такое, что доставляет экстремум целевой функции (1) при ограничениях u € Ω.

Для решения этой задачи погружаем ее в семейство подобных. Введем обозначения. Пусть ΩN, ΩN-1,N,… ≡ Ω - соответственно области определения для подобных задач на последнем этапе, двух последних и т.д.; Ω – область определения исходной задачи. Обозначим через F1(xN-1), F2(xN-2), … , Fk(xN-k), ... , FN(x0) соответственно условно-оптимальные значения функции цели на последнем этапе, двух последних и т.д., на к последних и т.д., на всех N этапах.

Начинаем с последнего этапа. Пусть xN-1 — возможные состояния системы на начало No этапа. Находим:

(2)

Для двух последних этапов получаем

(3)

Аналогично:

(4)

Выражения (2) - (4) называют функциональными уравнениями Беллмана. Отчетливо просматривается их рекуррентный (возвратный) характер, т. е. для нахождения оптимального управления на N шагах нужно знать условно-оптимальное управление на предшествующих N — 1 этапах и т. д. Поэтому функциональные уравнения часто называют рекуррентными (возвратными) соотношениями Беллмана.




1. ROMrdquo; Выполнил- студент ЭОУС311 Васильев А
2. БЕЛГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ НИУ БелГУ МЕДИЦИНСКИЙ И
3. Формирование культуры общения у детей подготовительной группы 2 слайд ~ Культура общения предусматрив
4. тематика на 20132014 учебный год
5. Реферат- Диагностика и восстановительное лечение больных с заболеваниями периферической нервной системы.html
6. Бухгалтерский учет и анализ в области экономического анализа; формирование практических навыков п
7. Цели учебной практики
8. довании страхового случая
9. История развития внутренних войск
10. 1 Понятие и сущность стратегии 5 1
11. Туризм как составляющая непроизводственной сферы деятельности
12.  Внутриорганизационная карьера предполагает прохождение различных ступеней профессионального роста обуч
13. А доц Гвоздева Т
14. а 700 Пилинг 35 гликолевая кислота с хитозаном 750 П
15. Хорошая политика не отличается от здоровой нравственности Г
16. 24w X 1566h cm 8 Count 22
17. Плюсы и минусы рекламы на транспорте
18. Грант 21 Анализ отраслевых особенностей управления персоналом организации на примере ООО
19. Тема- Экскурсия к избе подготовительная к школе группа Разработчик- Якимова Елена Алексан
20. ТЕМА 3 Защита частных прав Иски1