ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
Optimization of extraction and distribution of raw materials
between consumers depending on the period 1 2 Asankulova M. , Djusupbaev A.
Оптимизация добычи и распределения сырья
между потребителями в зависимости от периода
Асанкулова М.1, Жусупбаев А.2
1Асанкулова Майрам /Asankulova Mayram - кандидат физико-математических наук, старший научный сотрудник; 2Жусупбаев Амангельди /Djusupbaev Amangeldi - доктор физико-математических наук, заведующая лабораторией экономико-математических методов, Институт теоретической и прикладной математики, Национальная академия наук Кыргызской Республики, г. Бишкек, Кыргызская Республика
Аннотация: в работе сформулирована математическая модель задачи определения оптимального объема добычи сырья компанией и ее распределение между потребителями в котором затраты на добычу единицы объема сырья и ее транспортировку, затраты на переработку сырья в продукцию и цена реализации единицы объема продукции зависит от периода добычи сырья и его переработки. Для решения сформулированной задачи предложен метод решения.
Abstract: the paper formulated a mathematical model of the problem of determining the optimal amount of extraction of raw materials the company and its distribution among the consumers in which the cost of production volume unit of raw material and its transportation, processing costs of raw materials in the production and selling price of a unit of production volume depends on raw material extraction period and its processing. To solve the above problem is provided a method of solution.
Ключевые слова: математическая модель, транспортировка, объем сырья, метод аппроксимации, договор, компания, предприятия.
Keywords: mathematical model, transportation, the amount of raw materials, the approximation method, a contract company.
УДК 519.8
Постановка задачи. Пусть имеется компания, состоящая из m пунктов добычи сырья А; , ieI= {1,2,...,m}. Объем добываемого сырья на каждом пункте добычи
предполагается неизвестным x' , но ограниченным сверху величиной , ieI,
t=1,2,___, Т.
Сырье, добываемое компанией, согласно договора, должно перевозиться перерабатывающим предприятиям Пк, k=1,2,...,p ассоциации и потребителям Gr, r=1,2, ...,R по единой закупочной цене.
Исходя из договора, компания должна доставлять ассоциации за весь планируемый период сырье в объеме Q0 , в том числе каждому предприятию за
период t, t=1,2,..., Т в объеме не более чем d£, k=1,2,...,p, t=1,2,...,T, а частным
потребителям региона Gr, r=l,2,...,R - в объеме равным , r=1,2,...,R, t=1,2,..., Т.
Сырье, доставляемое предприятиям ассоциации, перерабатывается в продукцию и реализуется в том же периоде t, t = 1,2,...,Т.
Предполагаются, что затраты на добычу единицы объема сырья и ее транспортировка, затраты на переработку сырья в продукцию и цена единицы объема реализуемой продукции зависит от периода добычи сырья и его переработки.
Требуется определить объем добычи сырья для пунктов компании
>0, /е/. 1=1,2,..., Т, план распределения сьфья между предприятиями ассоциации и потребителями региона так, чтобы суммарные затраты на добычу, перевозку сырья были бы минимальными, а для перерабатывающих предприятий ассоциации - чистый доход от реализации готовой продукции из доставляемого сырья после переработки.
Для формулировки математической модели задачи введем следующие обозначения: г -индекс пунктов добычи сырья компании, г е!; /-множество индексов пунктов добычи сырья; г-индекс потребителей сырья региона, г=1,2,.,К; к-индекс перерабатывающих предприятий ассоциации, к=1,2,...,р; /-индекс периода добычи сырья компании, где сырье доставляется предприятиям ассоциации и потребителям региона, 1=1,2,..., Т.
Известные параметры и функции:
-максимальный объем добычи сырья 1-го пункта компании в /-ом периоде,
1=1,2,..., Т, г е/; ((£ - максимально возможный объем сырья, перевозимый компанией к-му предприятию ассоциации за /-ый период по договору, к=1,2, ...,р, 1=1,2,..., Т; Qo -объем сырья, поставляемый компанией предприятиям ассоциации за планируемый
период; - объем сырья, поставляемый компанией г-му потребителю региона за /ый период, г= 1,2,.,К, /=1,2,., Т; (р1{Х1 ) - функция, отражающая зависимость стоимости добываемого сырья от объема добычи г-го пункта в /-ом периоде, ге/,
/=1,2,..., Т; С^, - транспортные расходы на перевозку единицы объема сырья из г-го пункта добычи в к-ое предприятие ассоциации и г-му потребителю региона в /-ом
периоде соответственно, ге/, к=1,2,...,р, г=1,2,...,К, /=1,2,., Т; С £ - оптовая_цена реализации готовой продукции к-ым предприятием в /-ом периоде, к=1,2,...,р,
/=1,2,..., Т; - норма расхода сырья на единицу объема готовой продукции в к-ом
предприятии ассоциации в /-ом периоде, к=1,2,...,р, /=1,2,., Т; ££ - затраты на переработку единицы объема сырья в к-ом предприятии на /-ом периоде, к=1,2,.,р, /=1,2,., Т.
Искомые переменные:
г
Х(к - объем сырья, перевозимый компанией из г-го пункта добычи к-му
предприятию ассоциации перерабатывающих предприятий в /-ом периоде, г е/,
г
к=1,2,.,р, /=1,2,., Т; Хг - объем сырья, перевозимый компанией из г-го пункта
г
добычи г -му потребителю региона в /-ом периоде, 1е/, г=1,2,.,К, /=1,2,., Т; Х{ -
объем сырья, добываемый г-ым пунктом компании в /-ом периоде, ге/, /=1,2,., Т;
- объем готовой продукции к-го предприятии ассоциации в /-ом периоде,
к=1,2,.,р, /=1,2,., Т.
В соответствии с принятыми обозначениями математическая модель задачи определения оптимального объема добычи сырья компанией и ее распределение между потребителями по критерию минимума суммарных затрат на добычу сырья и перевозку запишется в виде. Найти минимум
ад{(¿¿л)(х)} (1)
¡е! г=1 к=1 г=1
Е <+Е 4 = X, е *=1,2-- т (2)
т=1 к=1
0 < х] < а], /е/, t=l,2,..., Т, (3)
0 <ЕХ^ < , k=1,2,..., р, t=1,2,..., Т,
(4)
р т
ЕЕЕ хк=оо,
(5)
/е/ к=1 ]=1
Е хТ = Ь
г=1,2,..., Я, (=1,2,., Т,
(6)
геI
х]к > 0, х'г > 0, 1е1, к=1,2,..., р, г=1,2,..., Я, (=1,2,
(7)
где
х =
14,р ' \х'Ащ\к ' 1 1'2'-' Т
Предполагается, что имеет место условие
т К т
Е ЕК+а0 <Е Е а;-
]=1 Т=1
г=1 не
(8)
Согласно принятым обозначениям математическая модель задачи определения оптимального объема готовой продукции и чистого дохода предприятиями ассоциации в каждом периоде могут быть записана в виде. Найти максимум
D(х, у) = Е («), * = 1,2,...,Т
к=1
iеI
при условиях
К Ук =Е хк, к=1,2,..., р, (=1,2,., Т,
ге/
х]к >0, / еI, к = 1,2,...,р, ] = 1,2,...,Т,
ук > 0, к=1,2,...,р, (=1,2,., Т, (12)
(9)
(10) (11)
где
х = х! , г = 1,2,..,Т, у=у! . I Мщ,р' ' ' ' ' у кк11,р
Из постановки задачи легко заметить, что задача (1)-(7) и (9)-(12) может быть решена последовательно.
Метод решения. Рассмотрим способ решения задачи (1)-(7) в случае, когда функция
(((х]) - выпуклая непрерывно возрастающая по х] е[0, а¿], 1е1, (=1,2,., Т.
Выпуклые функции (((х]), 1е1, (=1,2,., Т заменим кусочно-линейными функциями. Используем способ, приведенный в [1], [2], [3].
Разбиваем интервалы [0, а^] на I\ , (=1,2,..., Травных частей с шагом Ц = а] 1е1, (=1,2,., Т. Построим кусочно-линейную аппроксимацию функции (((х]), 1е1,
(=1,2,..., Т. Переменные х^ заменяем через следующим образом:
К
р
lt
X = £-У, i e I, t = 1,2,...,T, (13)
у=1
где
0 < z*v <h', i e I, t = 1,2,...,T. (14) Преобразовав неравенства (14) в равенства, имеем
И = 4 +¿1, i e I, t = 1,2,...,T (15)
где
С > 0, zl> 0, V=1,2,...,l], ieI, t=1,2,..., Т. Функцию (( (X ) представим приближенно в виде li zt ((X) = £{(((yht)-((((У-Щ)}, ieI, =1,2., Т. (16)
v=i hi
Из системы (13) и (15) получим
-У= И-С, y=1,2,...,l], i e I, t = 1,2,...,T,
it
X = Е(И-С), i e I, t = 1,2,...,T. (17)
у=1
Подставляя (17) в систему ограничений (2), получаем
R p li
Е+Еxk+ЕС = X, e t=1,2,., т (18)
r=1 k=1 У=1
где
0 <С <ht, y=1,2,...,lt, i e I, t = 1,2,...,T. Суммируя по i и t равенства (2) и (13), а (6) по r и t, далее используя (5), получим
T R T l't
Е ЕК+Q =Е £5Х> 0. (19)
t=1 r=1 t=1 ieI у=1
Подставляя значение ((() из (16) в целевую функцию задачи, получим
T p R i \mt
l( x)=££ {(Zctx+z<Xr)+Е^т У,
ieI t=1 k=1 r=1 y=1 hi
где
( _ ((lht)-((((y-1)ht) ц - Ц
-, у = 1,2,...,Ц, ieI, t = 1,2,..,T.
h
- угловые коэффициенты соответствующих звеньев кусочно-линейных
функций (х]) , 1е/, /=1,2,., Т.
Таким образом, окончательно имеем следующую задачу. Найти минимум
т Р К 1
Д х)=¿¿{(2 ¿Л¿х)+Е ^ х,} (20)
¡е! г=1 к=1 г=1 ,=1 ^
К р 1
Ех]т +Ех]к +Е&= а', (=1,2,., Т (21)
т=1 к=1 у=1
Ех]к < dtk, k=1,2,., р, t=1,2,., Т, (22)
ге1
Ех]Т = К, Г=1,2,., R, t=1,2,., Т, (23)
ге1
Р Т
ЕЕЕ х]к=00, (24)
ге/ к=1 !=1 Т ¡1 Т К
Е ЕЕ<=Е ЕК+а, (25)
!=1 ге1 у=1 г=1 т=1
<+£,= к], у=1,2,...,1, г е I, ] = 1,2,...,Т (26) х]> 0, ! 0, г е I, к = 1,2,..., р, т = 1,2,..., К, г = 1,2,...,Т, (27)
0, £1 > 0, У=1,2,...Х, /е/, 1=1,2,., Т. (28)
Задачу (20)-(28) при помощи запрещающих тарифов можно свести к закрытой модели транспортной задачи линейного программирования. Введем дополнительные
переменные х]к > 0, к=1,2,., р, (=1,2,., Т, и обращаем систему неравенств (22) в
равенства.
Определим объем фиктивного поставщика
Е £<к=Е Е< -&.
]=1 к=1 ]=1 к=1
Решив задачу (20)-(28) получим оптимальный план перевозок х*к > 0, / е1,
к=1,2,., р, (=1,2,., Т, х1 > 0, 1е1, г=1,2,., Я, (=1,2,., Т, и объемы добычи
сырья каждого пункта компании на каждом периоде =Е,(=1,2,., Т,
у=1
удовлетворяющий условиям (21)-(28) и доставляющий минимальное значение целевой функции (20).
Далее, используем решение х]к > 0, 1е1, к=1,2,. , р, (=1,2,., Т задачи (20)-(28)
для задачи (9)-(12).
Сформулируем экстремальную задачу следующего вида. Найти максимум
D(y) = Е с у]-Е<х]к), г = 1,2,..,Т (29)
к=1 ¡е1
при условиях
КУк = Е хк, к=1,2,., р, (=1,2,., Т, (30) ук > 0, к=1,2,., р, (=1,2,., Т. (31)
Решив задачу (29)-(31), получим оптимальный план у] > 0, к=1,2,., р, /=1,2,.,
Т, выпуска готовой продукции каждого предприятия ассоциации на каждом периоде и чистый доход ассоциации при договорных условиях работы.
Литература
1. Ланге Э. Г., Жусупбаев А. Комбинаторный метод решения задачи размещения. -Фрунзе, Илим, 1990. -153 с.
2. ХедлиДж. Нелинейное и динамическое программирование. - М.: Мир, 1967. -506 с.
3. Асанкулова М. Методы решения транспортно-производственной задачи. - Бишкек, Илим, 2012. - 159 с.
The parameters of our Universe and its structural elements
Grishchenko S.
Параметры нашей Вселенной и ее структурных элементов
Грищенко С. В.
Грищенко Сергей Васильевич / Grishchenko Sergey - главный редактор, Сайт astronomy3d.ru, г. Санкт-Петербург
Аннотация: заполнение космического пространства однотипными иерархическими звездными структурами позволило представить, в первом приближении, вероятное строение Вселенной. Методом последовательных приближений рассчитали усредненные параметры пространственного каркаса нашей Вселенной и ее структурных элементов: Мегагалактик, мегагалактических колец, протогалактических туманностей, спиральных галактик и пространственного цилиндра звезд. Сформулированы вероятные этапы эволюции нашей Вселенной и ее структурных элементов. Проведен анализ параметров видимых космических объектов на земном небосводе и их вероятная идентификация с известными группировками звезд.
Abstract: filling all space consistent hierarchical stellar structures allowed us to understand, in the first approximation, the likely structure of the Universe. Iterative detailed mathematical analysis has allowed to calculate the parameters of the spatial framework of our Universe and its structural elements: megagalaktik, megagalactic rings, protogalactic nebula, spiral galaxies and the spatial cylinder of the stars. Formulated probable stages of evolution of our Universe and its structural elements. The analysis is made of the parameters of visible space objects in earth's sky and their possible identification with known groups of stars.
Ключевые слова: расчетные параметры, пространственный каркас, Вселенная, эшелон Мегагалактик, слой эшелонов Мегагалактик, эллиптическая галактика, пространственный цилинд, мегагалактическое кольцо, спиральная галактика, протогалактическая туманность, строение Вселенной и ее структурных элементов, этапы эволюции, параметры видимых космических объектов.
Keywords: design parameters, space frame, the Universe, echelon of megagalaxies, layer of echelons of megagalaxie, an elliptical galaxy, space cylinder, megagalactic ring, spiral galaxy, protogalactic nebula, the structure of the Universe and its structural elements, the stages of the evolution, the parameters of visible space objects.