Исходя из маркетинговых исследований рынка, можно считать, что в рассматриваемом регионе или региональном окружении перерабатывающее предприятие сформировало некоторый портфель заказов на закупку сырья (при этом в одних случаях поставщики могут обеспечить предприятие сырьем на 100 %, в других случаях может иметь место дефицит сырья (в том числе и из-за фактора сезонности)). На часть готовой продукции имеются договорные обязательства с бизнес-партнерами, в том числе и госзаказ (больницы, школы, детские дома, войсковые части и т.д.), зависящий от сезонного фактора (осень, зима, лето, весна), и определена номенклатура молочных изделий, востребованных потребителями.
Введем обозначения для всего технологического комплекса молочного комбината через T = {T1, T2, ..., Tn}, обеспечивающего выпуск продукции k-й номенклатуры (кефир, сыр, йогурт и т.д.) на основании соответствующей технологии Tk, (k = 1, 2, ..., n).
Максимальная размерность вектора T составляет порядка 120-150 компонентов. Каждая технология Tk разбивается на технологические процессы или подтехнологии (этапы, составляющие), их количество порядка 20-30. Каждый такой процесс j характеризуется регламентированным технологией временем протекания tjk, являющимся постоянной величиной (k-технология, j-процесс данной технологии); максимально возможной производительностью (загрузкой) оборудования Rjk (например, отношением времени работы оборудования к времени рабочей смены); финансово-экономическими издержками - затратами на j-м процессе технологии Tk-Rik; экологическими издержками j-го процесса технологии Tk - величиной ejk и т.д. Причем, эти величины могут быть детерминированными либо стохастическими величинами. Пусть V - максимально возможный объем молока, который может поставляться производителями на комбинат. В условиях ограниченности валового объема переработки V1, обусловленной имеющимися мощностями комбината, и дефицита исходного сырья V, выражаемого ограничениями V < V1 (а часто и при V << V1), возникает оптимизационная задача оптимального выбора номенклатуры продукции и ее объемов выпуска, которые обеспечат максимальную прибыль комбинату, при выполнении гарантированных поставок продукции по госзаказу и по договорным обязательствам предприятия перед постоянными партнерами-потребителями. Разработка метода решения этой задачи и является целью нашего исследования.
Сформулируем математическую модель задачи. Пусть X = (x1, x2, ..., xn) - вектор объемов производимой продукции (xk - объем выпуска продукции k-й номенклатуры) - искомая величина. Задан вектор V = (v1, v2, ..., vn) с компонентами vk - удельными расходами молока, требующегося на изготовление единицы продукции k-й номенклатуры. Очевидно, что vk - детерминированная величина.
Задан вектор удельного дохода D = (d1, d2, ..., dn) от реализации готовой продукции (dk - соответствующая характеристика единицы продукции k-й номенклатуры).
Задан вектор Z = (z1, z2, ..., zn) - вектор издержек при изготовлении i-й продукции.
Величины D и Z можно считать стохастическими, зависящими от спроса, сезонности продаж, поставок сырья, оптовых скидок, затрат различного вида и т.д., т.е., в этой трактовке D и Z - векторы математического ожидания дохода и издержек.
Можно задать вектор экологических затрат на единицу готовой продукции k, компоненты которого связаны с суммарными экологическими издержками j-го процесса ejk. Обозначим вектор удельных экологических затрат E = (e1, e2, ..., en) (где ek - соответствующая характеристика в расчете на единицу выпускаемой продукции k-й номенклатуры). Эти затраты можно в ряде случаев считать как детерминированными, так и стохастическими, т.е., E - математическое ожидание вектора экологических издержек.
Экологические издержки связаны с неполной переработкой сырья, выбросами отходов в канализацию и в воздух, что влечет за собой существенные штрафные санкции за загрязнение окружающей среды, а также установку (часто дорогостоящих) систем и оборудования для очистки.
Следует иметь в виду, что величины компонент искомого вектора X имеют ограничения снизу и сверху, т.е., rk ≤ xk ≤ Rk, где rk - суммарные обязательные поставки готовой продукции по договорным обязательствам с постоянными партнерами и поставки по госзаказу;
,
где 
- производительность оборудования для изготовления изделия k, Пk - максимально возможный объем продукции k, который в принципе можно реализовать (спрос на k-ю продукцию).
С учетом сказанного имеем:
(1)
(2)
(3)
(4)
(5)
(6)
Rk ≤ xk ≤ Rk, k = 1, 2, ..., n, (7)
где vk, gk, mk - удельные затраты (на единицу выпускаемой продукции) молока, электроэнергии и транспортные издержки соответственно, а G, M - соответствующие квоты.
Следует заметить, что ограничения (1), (2), (4) в определенных случаях могут иметь нелинейный вид (например, в (4) нелинейность может появиться в зависимости от каких-либо условий в процессах технологии Tk, например, температуры, давления и других физических и химических процессов):
(8)
где I - множество индексов процессов с линейной зависимостью удельных расходов сырья, а J - с нелинейным характером зависимостей, причем I ∪ J ≡ {1, 2, .., n}. Функция f(xk) определяется эмпирическим путем, исходя из статистических данных и их последующей обработки. Например, в некоторых случаях она аппроксимируется полиномиальной регрессионной зависимостью
где коэффициенты a, b, c определяются на основе статистических оценок.
Итак, в достаточно общем случае имеем задачу (1)-(7) математического программирования с тремя критериями в виде (1)-(3), и ограничениями (4)-(7), учитывающими стохастические и нелинейные процессы. В последних случаях реализации следует заменить детерминированные величины dk, zk и ek в (1)-(3) на значения их математических ожиданий, а условие (4) - на ограничение (8).
В первом приближении модели поставленная трехкритериальная задача линейного программирования (1)-(7) может быть сведена к однокритериальной в силу свойства аддитивности отражаемых экономических характеристик, т.е., в случае явной связи между критерием, описывающим максимизацию дохода, критерием, минимизирующим издержки комбината, и выделенным критерием, описывающим минимизацию экологических издержек (этот критерий следует учитывать, если суммарные экологические издержки сопоставимы с суммарными издержками на выпуск продукции); в противном случае при незначительной величине экологических издержек критерий (3) можно не рассматривать и задача превращается в двухкритериальную). В большинстве случаев рассматриваемые проблемы сводятся к постановкам задач линейного программирования, так как доход и издержки, как правило, линейно зависят от объема выпускаемой продукции. Это имеет место для комбинатов, средних по мощности, выпускающих стандартную продукцию. Для решения задачи (1)-(7) нами были использованы методы сведения к одному критерию [1] либо алгоритм решения двухкритериальной задачи, разработанный Е.Н. Остроухом и Л.Ч. Гамкрелидзе [2].
Следует заметить, что в случае функционирования больших предприятий (холдингов, объединений), выпускающих большой объем продукции, в том числе и нестандартной, как функция, описывающая доход D, так и функция, описывающая издержки Z, может иметь нелинейный вид, причем, доход описывается полиномиальной функцией второго порядка, а издержки - полиномиальной функцией третьего порядка, вследствие чего критерии (1) и (2) будут иметь следующий вид:
(1*)
(2*)
Для нахождения коэффициентов 
в соотношениях (1*) и (2*) используется теория временных рядов. Создаются диаграммы зависимости данных от времени для отображения тренда, сезонных изменений и выбросов. Если данные изменяются со временем, то они преобразовываются для сохранения постоянства дисперсии, при этом для анализа экономических данных нами были использованы логарифмические преобразования с проверкой постоянства дисперсии для исходных и логарифмированных данных с течением времени. Отфильтруем регулярные составляющие временных рядов (убрав случайный шум), отметим, что они являются либо трендом, либо сезонной составляющей. Тренд представляет собой общую систематическую линейную или нелинейную компоненту, которая может изменяться во времени. Сезонная составляющая - периодически повторяющаяся компонента. В нашем случае имеем модели временного ряда, в которых амплитуда сезонных изменений увеличивается вместе с трендом, т.е., модели с мультипликативной сезонностью. Периодическую зависимость определяем как корреляционную зависимость порядка k между каждым i-м элементом ряда и (i - k)-м элементом. Ее можно измерить с помощью автокорреляции, т.е. корреляции между членами ряда. Сезонные составляющие временного ряда могут быть найдены с помощью коррелограммы.
Авторами был проведен ряд расчетов для решения задачи в постановке (1), (4)-(7) и (1), (4)-(6), (8) соответственно симплекс-методом и методами градиентного спуска и Бокса [3, 4]. Результаты расчетов дают возможность руководителям предприятия оперативно реагировать на изменение внешних условий, например, сезонности, погодных условий, спроса и других факторов, т.е., изменять и номенклатуру и объем выпускаемой продукции в целях оптимизации производства.
Для решения задачи в постановке (1*), (2*), (4)-(7) в случае перехода к одному критерию (весовые коэффициенты каждого из двух критериев выбираются на основании экспертных оценок) предложен генетический алгоритм, представляющий собой ядро, на основе которого можно строить алгоритмы для задач различного класса и различной размерности. Система работает со значениями типа double, позволяя задавать необходимую точность вычислений. Разработан собственный декодер, наследуемый от стандартного класса декодера. Алгоритм одинаково хорошо решает как задачи на максимизацию, так и на минимизацию. Благодаря указанию точности, можно настроить систему на работу с целыми числами. Система реализует довольно гибкую систему ограничений, которые могут быть как линейными, так и нелинейными. Система является объектной, следовательно, все параметры и критерии, требуемые для алгоритма, реализуются с помощью классов объектов. Структуры, требующие наибольшей гибкости, реализованы в виде интерфейсов. Это позволяет создавать собственные классы, наследуемые от этих структур. Используя данную возможность, можно гибко настраивать алгоритм для каждой конкретной задачи.
Использование предложенной модели в целом позволит решать ряд задач оптимизации работы холдингов и объединений, выпускающих как традиционную, так и инновационную продукцию, что поможет руководству оперативно и квалифицированно принимать оптимальные решения в плане стратегии предприятий и объединений в различных ситуациях.
Следует отметить, что в зависимости от особенностей задачи (нелинейные функции, стохастика) можно использовать, например, алгоритмы квадратичного, стохастического программирования [5], метод синтеза критериев [6, 7]. Решения задач указанных классов реализованы на основе разработанного пакета программных средств многокритериальной оптимизации.
Работа поддержана ЦКП «Научные и научно-педагогические кадры инновационной России».
Рецензенты:
Потетюнко Э.Н., д.ф.-м.н., профессор кафедры «Теория упругости» Южного федерального университета, г. Ростов-на-Дону;
Чистяков А.Д., д.т.н., профессор, зав. кафедрой «Экономика и менеджмент» Донского государственного технического университета, г. Ростов-на-Дону.
Работа поступила в редакцию 08.07.2011.