Моделирование распределения фондов минеральных удобрений сельскохозяйственной организации
МИНИСТЕРСТВО
СЕЛЬСКОГО ХОЗЯЙСТВА РОССИЙСКОЙ ФЕДЕРАЦИИ
ФГБОУ
ВПО Российский государственный аграрный университет - МСХА имени К.А.
Тимирязева
Инженерно-экономический
факультет
Кафедра
«Инжиниринг бизнес-процессов»
РАСЧЕТНАЯ
РАБОТА
по
дисциплине «Моделирование экономических систем»
«Моделирование
распределения фондов минеральных удобрений сельскохозяйственной организации»
Выполнила: студентка 21
МАГ/Э Орлова А.В.
Приняла: к.т.н., доцент
Качанова Людмила Сергеевна
Москва,
2014
Содержание
Введение
. Разработка экономико-математической модели
.1 Система переменных экономико-математической модели
.2 Система ограничений экономико-математической модели
.2.1 Группа ограничений по балансу выноса элементов питания
продукцией и внесения их с удобрениями
.2.2 Группа ограничений по границе интервалов прибавки урожайности
.2.3 Группа ограничений по формированию годовых норм удобрений в
ассортименте поставки
.2.4 Группа ограничений по распределению годовых норм удобрений по
срокам внесения
.2.5 Группа ограничений по допустимому удельному весу отдельных
форм удобрений в общей дозе
.2.6 Группа ограничений по суммарному приросту урожайности на
участке
.2.7 Группа ограничений по балансу ресурсов и потребления
удобрений
1.2.8 Группа ограничений по производству продукции
1.2.9 Условия неотрицательности переменных
экономико-математической модели
.2.10 Целевая функция экономико-математической модели
. Подготовка исходной информации
. Решение экономико-математической задачи распределения фондов
минеральных удобрений сельскохозяйственной организации по полям севооборотов и
кормовым угодьям
.1 Формирование отчетов по результатам решения
. Анализ результатов решения
Список литературы
Приложение 1
Введение
Одним из важнейших разделов
математического программирования является линейное программирование.
Методы и модели линейного
программирования широко применяются при оптимизации процессов во всех отраслях
народного хозяйства: при разработке производственной программы предприятия,
распределении ее по исполнителям, при размещении заказов между ис-полнителями и
по временным интервалам, при определении наилучшего ассортимента выпускаемой
продукции, в задачах перспективного, текущего и оперативного планирования и
управления; при планировании грузопотоков, определении плана товарооборота и
его распределении; в задачах развития и размещения производительных сил, баз и
складов систем обращения материальных ресурсов и т. д. Особенно широкое
применение методы и модели линейного программирования получили при решении
задач экономии ресурсов (выбор ресурсосберегающих технологий, составление
смесей, раскрой материалов), производственно-транспортных и других задач.
Начало линейному программированию
было положено в 1939 г. советским математиком-экономистом Л. В. Канторовичем в
работе «Математические методы организации и планирования производства».
Появление этой работы открыло новый этап в применении математики в экономике.
Спустя десять лет американский математик Дж. Данциг разработал эффективный
метод решения данного класса задач - симплекс-метод. Общая идея симплексного
метода (метода последовательного улучшения плана) для решения задач линейного
программирования (ЗЛП) состоит в следующем:
) умение находить начальный
опорный план;
2) наличие признака
оптимальности опорного плана;
) умение переходить к
нехудшему опорному плану.
Математическое моделирование
экономических ситуаций на базе современной вычисли-тельной техники позволяет
автоматизировать сбор и обработку первичной информации, выделить основные
параметры, влияющие на деятельность фирмы, рассчитать различные варианты
деятельности фирмы, обеспечивающие необходимую эффективность производства или
предпринимательства, и на основе этих данных принять решение о выборе
оптимальной стратегии по управлению деятельностью фирмы.
Современный экономист должен хорошо
разбираться в экономико-математических методах, уметь их практически применять
для моделирования реальных экономических ситуаций. Это позволит лучше усвоить
теоретические вопросы современной экономики, повысить уровень квалификации и общей
профессиональной культуры специалиста.
Как показывает опыт, оптимальные
решения бывают лучше решений, принятых традици-онными методами, на 5…15 %
величины критерия, по которому производится оптимизация. Вместе с тем, принятие
оптимальных решений связано со следующими трудностями:
. необходимы знания, без которых
принятие оптимальных решений невозможно;
. требуется специальное прикладное
программное обеспечение. В компьютерных классах ФГБОУ ВПО МГАУ им. В.П.
Горячкина для решения оптимизационных задач используется приложение Поиск
решения (Excel Solver) надстройка MS Excel.
Сервисы и функции MS Excel
незаменимы для научных работников и практиков. Эффективность поддержки принятия
решений определяется грамотной постановкой задачи, подбором средств ее решения,
подготовкой исходных данных, формулировкой целевой функции, условий,
ограничений, интерпретацией результатов и оценкой их надежности.
1. Расчетная часть. Разработка
экономико-математической модели
.1 Система переменных
экономико-математической модели
Участок № 11, x2,
x3 - дозы действующего вещества, соответственно N, P2O5,
и K2O, отнесенные на прирост урожайности озимой пшеницы по
первому интервалу прибавки, кг д. в/га;
х4 - прирост урожайности
озимой пшеницы по первому интервалу прибавки, ц;5, x6, x7
- дозы действующего вещества, соответственно N, P2O5, и
K2O, отнесенные на прирост урожайности озимой пшеницы по второму
интервалу прибавки, кг д. в/га;
х8 - прирост урожайности
озимой пшеницы по второму интервалу прибавки, ц;9, x10, x11,
x12 - дозы в физической массе соответственно аммиачной
селитры, карбамида, суперфосфата и калийной соли для основного внесения под
озимую пшеницу, ц/га;13, x14 - дозы в физической массе
соответственно аммиачной селитры и карбамида для внесения в подкормку под
озимую пшеницу, ц/га;
х15 - общий прирост
урожайности озимой пшеницы, ц.
Участок 216, x17,
x18 - дозы действующего вещества, соответственно N, P2O5,
и K2O, отнесенные на прирост урожайности картофеля по первому
интервалу прибавки, кг д. в/га;
х19 - прирост урожайности
картофеля по первому интервалу прибавки, ц;20, x21, x22
- дозы действующего вещества, соответственно N, P2O5, и
K2O, отнесенные на прирост урожайности картофеля по второму
интервалу прибавки, кг д. в/га;
х23 - прирост урожайности
картофеля по второму интервалу прибавки, ц;24, x25, x26,
x27 - дозы в физической массе соответственно аммиачной
селитры, карбамида, суперфосфата и калийной соли для основного внесения под
картофель, ц/га;28, x29 - дозы в физической массе
соответственно аммиачной селитры и карбамида для внесения в подкормку под
картофель, ц/га;
х30 - общий прирост
урожайности картофеля, ц.
.2 Система ограничений
экономико-математической модели
.2.1 Группа ограничений по балансу
выноса элементов питания продукцией и внесения их с удобрениями
Ограничения обеспечивают
соответствие норм удобрений в единицах действующего вещества приросту
урожайности. Привязка осуществляется в рамках границ интервалы в урожайности,
где исследована зависимость урожая от удобрений, с дифференциацией по элементам
питания.
Ограничения по затратам действующего
вещества удобрений в первом интервале прибавки урожайности озимой пшеницы:
по азоту (ограничения по затратам
1):
х1+3,3х4=0;
по фосфору (ограничения по затратам
2):
х2+3,6х4=0;
по калию (ограничения по затратам
3):
х3+2,9х4=0;
Каждый центнер зерна в первом
интервале прибавки урожайности ячменя требует внесения, в перерасчете на
действующее вещество, 3,3 кг азотных, 3,6 кг фосфорных и 2,9 кг калийных
удобрений.
Ограничения по затратам действующего
вещества удобрений во втором интервале прибавки урожайности озимой пшеницы:
по азоту (ограничения по затратам
5):
х1+4,0х4=0;
по фосфору (ограничения по затратам
6):
х2+4,2х4=0;
по калию (ограничения по затратам
7):
х3+3,7х4=0;
Каждый центнер зерна во втором
интервале прибавки урожайности озимой пшеницы требует внесения, в перерасчете
на действующее вещество, 4,0 кг азотных, 4,2 кг фосфорных и 3,7 кг калийных
удобрений.
Ограничения по затратам действующего
вещества удобрений в первом интервале прибавки урожайности картофеля:
по азоту (ограничения по затратам
15):
х1+0,6х4=0;
по фосфору (ограничения по затратам
16):
х2+0,7х4=0;
по калию (ограничения по затратам
17):
х3+0,6х4=0;
Каждый центнер зерна в первом
интервале прибавки урожайности картофеля требует внесения, в перерасчете на
действующее вещество, 0,6 кг азотных, 0,7 кг фосфорных и 0,6 кг калийных
удобрений.
Ограничения по затратам действующего
вещества удобрений во втором интервале прибавки урожайности картофеля:
по азоту (ограничения по затратам
19):
х1+0,7х4=0;
по фосфору (ограничения по затратам
20):
х2+0,8х4=0;
по калию (ограничения по затратам
21):
х3+0,8х4=0;
Каждый центнер зерна во втором
интервале прибавки урожайности картофеля требует внесения, в перерасчете на
действующее вещество, 0,7 кг азотных, 0,8 кг фосфорных и 0,8 кг калийных
удобрений.
Обобщенная математическая запись
ограниченной данной группы имеет вид
-λ1xrlk+arlkxrk=0,
(kϵK, lϵL, rϵR), (1.1)
где r-номер элементарной культуры, rϵR;
R- множество, составленное номерами элементарных культур; l- номер вида удобрения
(элемента питания), lϵL; L-множество, составленное номерами видов
удобрений (элементов питания); k-номер интервала прибавки урожайности, kϵK;
K-множество, составленное номерами интервалов прибавки урожайности; λ1- содержание элементов питания l-го вида в единице действующего
вещества удобрения l-го вида; arlk- затраты действующего вещества
удобрений l-го видана единицу прироста урожайности r-й элементарной культуры в
k-м интервале прибавки, кг д. в/га; xrlk - искомая доза действующего
вещества удобрения l-го вида, отнесенная на прирост урожайности r-й
элементарной культуры по k-му интервалу прибавки, кг д. в/га; xrk -
искомый прирост урожайности r-й элементарной культуры в k-м интервале прибавки,
ц/га.
.2.2 Группа ограничений по границе
интервалов прибавки урожайности
Ограничения реализуют условия по
пределу прироста урожайности элементарной культуры в выделенном интервале
прибавки.
Обобщенная математическая запись
ограничений данной группы:
λrkхrk≤brk,
(kϵK, rϵR), (1.2)
где brk - верхняя граница
k-го интервала прибавки урожайности r-й элементарной культуры, ц/га; λrk - выход основной продукции с единицы прироста урожайности r-й
элементарной культуры по k-му интервалу прибавки; λrk =1, (kϵK, rϵR).
По величине первого интервала
прибавки урожайности озимой пшеницы (огр. 4):
х4 ≤6.
По величине второго интервала
прибавки урожайности озимой пшеницы (огр. 8):
х8 ≤5.
По величине первого интервала
прибавки урожайности картофеля (огр. 18):
х18 ≤30.
По величине второго интервала
прибавки урожайности картофеля (огр. 22):
х22 ≤60.
.2.3 Группа ограничений по
формированию годовых норм удобрений в ассортименте поставки
Ограничения предназначены для
перехода от суммарной годовой нормы удобрений в единицах действующего вещества
к дозам конкретных форм удобрений в единицах физической массы.
По формированию доз азотных
удобрений (аммиачной селитры и карбамида) для внесения под озимую пшеницу (огр.
9):
х1+х5= 34х9+46х10+34х13+46х14;
или в результате преобразования:
х1+х5 - 34х9-46х10-34х13-46х14=0.
Обе части уравнения определяют
годовые нормы азотных удобрений в единицах действующего вещества.
Левая часть уравнения представляет
собой сумму переменных, обозначающих нормы азотных удобрений, обеспечивающих
прирост урожайности по первому и второму интервалам прибавок.
Правая часть уравнения представлена
суммой произведений переменных, обозначающих дозы аммиачной селитры и карбамида
в основное внесение и в подкормку в физической массе, на технико-экономические
коэффициенты. В качестве технико-экономических коэффициентов использованы
нормативы содержания действующего вещества в единице физической массе
удобрений.
Отличительной особенностью
ограничений 10, 11, 24, 25 является более простая форма записи.
Эта особенность определяется главным
образом упрощением условий по ассортименту и срокам использования фосфорных и
калийных удобрений в настоящей постановке задачи:
предоставление в распоряжение
трудового коллектива простых форм удобрений при отсутствии возможности выбора:
в числе ресурсов представлены по одной форме фосфорного и калийного удобрения;
внесение фосфорных и калийных
удобрений в полном объеме в один срок-перед посевом (основное внесение).
Обобщенная математическая запись
(1.3)
где: f-номер формы
удобрения, fϵF; F-множество, образованное номерами форм удобрений; F1-подмножество
множества F, составленное номерами форм удобрений, отнесенных к l-м виду
удобрений; номер срока внесения удобрений, tϵT; T- множество,
составленное номерами сроков внесения удобрений; Vef - содержание
действующего вещества l-го вида в единице физической массы f-й формы удобрения,
кг д. в/га; xrtf - искомая доза f-й формы удобрения для внесения в
t-й срок вод r-ю элементарную культуру, ц.
По формированию доз
фосфорных удобрений для внесения под озимую пшеницу (огр. 10):
х2+х6-
20х11=0.
По формированию доз
калийных удобрений для внесения под озимую пшеницу (огр. 11):
х3+х7-
40х12=0.
По формированию доз
азотных удобрений (аммиачной селитры и карбамида) для внесения под картофель
(огр. 23):
х16+х20 -
34х24-46х25-34х28-46х29=0.
По формированию доз
фосфорных удобрений для внесения под картофель (огр. 24):
х17+х21-
20х26=0.
По формированию доз
калийных удобрений для внесения под картофель (огр. 25):
х18+х22-
40х27=0.
.2.4 Группа ограничений
по распределению годовых норм удобрений по срокам внесения
Ограничения
предназначены для формирования доз удобрений в ассортименте в соответствии с
требованиями по агросрокам применения удобрений.
По распределению годовых
норм аммиачной селитры и карбамида под зерновые (озимая пшеница) для основного
внесения и в подкормку не менее 20% от годовой нормы азотных удобрений в
единицах действующего вещества (огр. 12):
х13+46х20
≥0,2(34х9+46х10+34х13+46х14);
или в результате
преобразований
,2х13-36,8х14+6,8х9+9,2х10≤0;
Для картофеля (огр.26):
34х28+46х29 ≥0,2(34х24+46х25+34х28+46х29);
,2х28-36,8х29+6,8х24+9,2х25≤0;
Обобщенная математическая запись
ограничений данной группы
, (1.4)
где - Wlt -
предельно допустимая доля внесения в t-й срок удобрения l-го вида в годовой
норме удобрения l-го вида.
.2.5 Группа ограничений по
допустимому удельному весу отдельных форм удобрений в общей дозе
По удельному весу карбамида в дозе
азотных удобрений, вносимой под зерновые в подкормку: не более 50% от общей
дозы в пересчете на действующее вещество (огр. 13):
х14≤0,5(34х13+46х14);
х14-17х13≤0
Под картофель (огр.27):
х29≤0,5(34х28+46х29);
х29-17х28≤0.
Обобщенная математическая запись
ограничений данной группы:
, (1.5)
где- Wlt -
предельно допустимая доля дозы внесения в t-й срок f-й формы удобрения l-го
вида в общей дозе удобрения в рассматриваемый срок внесения этого же вида.
1.2.6 Группа ограничений
по суммарному приросту урожайности на участке
Ограничения выполняют
вспомогательную функцию: позволяют сформировать значение специальной
переменной, характеризующей общую прибавку урожайности по элементарной
культуре, суммированием частных приростов урожайности по интервалам прибавки.
По суммарному приросту
урожайности озимой пшеницы (огр. 14):
х4+х8=х15,
или
х4+х8-х15=0.
По суммарному приросту
урожайности картофеля (огр. 28):
х19+х23=х30,
или
х19+х23-х30=0.
Обобщенная
математическая запись
(1.6)
где выход
продукции с единицы прироста урожайности r-й элементарной культуры; хк=1
(r;
xк- вспомогательная переменная, обозначающая суммарный прирост
урожайности по r-й элементарной культуре, ц.
.2.7 Группа ограничений
по балансу ресурсов и потребления удобрений
Ограничения входят в
состав связующего блока.
По фонду аммиачной
селитры, т. физ. массы (огр. 29):
,1(71(х9+х13)+64(х24+х28))≤12,53,
или
,1х9 + 7,1х13
+ 6,4х24 + 6,4х28 ≤ 12,53.
По фонду карбамида, т.
физ. массы (огр. 30):
,1(71(х10+х14)+64(х25+х29))≤1,5,
или
,1х10 + 7,1х14
+ 6,4х25 + 6,4х29 ≤ 1,5.
По фонду суперфосфата,
т. физ. массы (огр. 31):
,15(71х11 +64х26)≤
20,105,
или
,65 х11 + 9,6
х26 ≤ 20,105.
По фонду калийной соли,
т. физ. массы (огр. 32):
,15(71х12 +64х27)≤
24,5,
или
,65 х12 + 9,6
х27 ≤ 24,5.
Обобщенная
математическая запись
(1.7)
где -
технико-экономический коэффициент, определяющий изменение единиц измерения, =0,1
;
ar-площадь участка под r-й элементаной культурой, га; bа-фонды
удобрений f-й формы, т.
.2.8 Группа ограничений
по производству продукции
Ограничения реализуют
условия по выполнению заданий по производству отдельных видов продукции в
натуре (прирост урожая за счет удобрений). По приросту производства озимой
пшеницы, т (огр. 33):
,1*71х15≥49,8
или
,1х15≥49,8.
По приросту производства
картофеля, т (огр. 34):
,1*64х30≥320
или
,4х30≥320.
Обобщенная
математическая запись
(1.8)
где m-номер вида
интегральной продукции, mM;
М- множество, составленное номерами видов продукции; Pm-подмножество
множества R, составленное номерами элементарных культур, отнесенных к m-му виду
интегральной продукции; -выход
продукции m-вида в пересчете на единицу продукции r-й элементарной культуры; Wr-масштабный
множитель перевода центнеров в тонны; , bm-
минимальный допустимый объем производства продукции m-го вида, т.
.2.9 Условия неотрицательности
переменных экономико-математической модели
хrtlf≥0,
где fϵF, lϵL, tϵT, rϵR;rlk≥0, где
kϵK, lϵL, rϵR;rk≥0, где kϵK, rϵR;r≥0,
где rϵR.
.2.10 Целевая функция
экономико-математической модели
В качестве критерия
оптимальности использован показатель дополнительного чистого дохода,
отнесенного на удобрение.
Целевая функция имеет
видz = - 72420х9-113600 х10- 124250х11-347190х12-
83709 х13-121055х14
+66849,34 х15-65280
х24-102400 х25-112000х26- 312960 х27-75456х28-
109120х29-+48720х30.
Обобщенная математическая
запись целевой функции может быть представлена следующим образом:
(2.9)
где cft -
удельные суммарные затраты на применение в t-й срок 1ц f-й формы удобрения.
Суммарные затраты складываются из затрат на приобретение, транспортировку,
приготовление и внесение удобрений и накладных расходов в расчете на 1 ц
удобрений; сmr-цена 1ц продукции r-й элементарной культуры,
отнесенной к m-му виду интегральной продукции,р.; сr- удельные
затраты на уборку, транспортировку и доработку 1ц продукции r-й элементарной
культуры, ц.
2. Подготовка исходной
информации
модель распределение
фонд удобрение
Полный учет
взаимодействующих факторов, определяющих потребность элементарных культур в
питательных веществах, проблематичен и не является самоцелью в моделировании
процесса использования удобрений. Кроме того, в современном земледелии
действующее вещество удобрений используется растениями не в полном объеме.
Степень использования удобрений определяется набором наиболее существенных почвенно-климатических.
Агротехнических и организационно-экономических нормообразующих факторов.
Многообразие методик
расчета доз затрудняет формирование единого информационного обеспечения,
планирования и анализа применения удобрений. Показатели использования
питательных веществ из удобрений в первый, второй и другие годы, коэффициенты
распределения годовой нормы по срока внесения определяются в соответствии с
действующими справочниками по удобрениям, рекомендациями и данными анализов
научно-исследовательских учреждений, проектно-изыскательных станций химизации,
зональных агрохимических лабораторий, агрохимических центров.
Почвенно-агрохимические
параметры в модели не получают прямой оценки. Опосредованное их влияние на
эффективность удобрений учитывается при нормировании частного эффекта от
удобрения по элементарной культуре.
Затраты на приготовление
и внесение удобрений и уборку дополнительного урожая калькулируются на
основании норм выработки и расценок на соответствующие работы.
Расчет показателей последствия
удобрений, внесенных под предшественники, а также плановые затраты удобрений,
отнесенные на исходную урожайность представлены в таблице 1.
Таблица 1
Расчет показателей
последствия удобрений, внесенных под предшественники, а также плановые затраты
удобрений, отнесенные на исходную урожайность
Участок
|
Площадь, га
|
Нормы внесения удобрений под предшественник, кг д.в./га
|
«Стартовые» нормы под плановую культуру, кг д.в./га
|
|
|
Всего
|
Из них переходит (в форме последействия) на плановый период
|
Исходные
|
Скорректированные на величину последействия удобрений
предшественника
|
|
|
N
|
P2O5
|
K2O
|
N
|
P2O5
|
K2O
|
N
|
P2O5
|
K2O
|
N
|
P2O5
|
K2O
|
1
|
2
|
3
|
4
|
5
|
6
|
7
|
8
|
9
|
10
|
11
|
12
|
13
|
14
|
3
|
71
|
60
|
60
|
40
|
6
|
9
|
6
|
60
|
70
|
30
|
54,0
|
61,0
|
24,0
|
4
|
64
|
80
|
59
|
40
|
8
|
8,85
|
6
|
60
|
40
|
30
|
52,0
|
31,15
|
24,0
|
-
|
135
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
Потребность в удобрениях под исходный урожай (скорректированная
с учетом последействия удобрений предшественника) в расчете на всю площадь
|
В единицах действующего вещества, ц д.в.
|
В пересчете на конкретный ассортимент, т. физ. Веса
|
N
|
P2O5
|
K2O
|
Аммиачная селитра
|
Карбамид
|
суперфосфат
|
Калийная соль
|
15
|
16
|
17
|
18
|
19
|
20
|
21
|
38,34
|
43,31
|
17,04
|
11,28
|
-
|
21,655
|
4,26
|
33,28
|
19,94
|
15,36
|
9,79
|
-
|
19,94
|
3,84
|
71,62
|
63,25
|
32,4
|
21,07
|
-
|
41,595
|
8,1
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Расчет фондов удобрений, отнесенных
на прирост урожая, приведен в таб. 2, выход продукции, отнесенной на исходную
урожайность, определён в таб. 3., расчет обеспечения заданных объемов
производства продукции представлен в табл. 4.
Таблица 2
Расчет фондов удобрений, отнесенных
на прирост урожая, т физ. веса
Вариант
|
Фонды удобрений под урожай планового периода в ассортименте
поставки
|
Потребность в удобрениях под исходный урожай, скорректированная
на последствие удобрения предшественника
|
Подлежит распределению под прирост урожая
|
1
|
Аммиачная селитра
|
карбамид
|
суперфосфат
|
Калийная соль
|
Аммиачная селитра
|
карбамид
|
суперфосфат
|
Калийная соль
|
Аммиачная селитра
|
карбамид
|
Суперфосфат
|
1
|
2
|
3
|
4
|
5
|
6
|
7
|
8
|
9
|
10
|
11
|
12
|
13
|
3
|
33,6
|
1,5
|
61,7
|
32,6
|
21,07
|
-
|
41,595
|
8,1
|
12,53
|
1,5
|
20,105
|
24,5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 3
Расчет показателя выхода продукции,
отнесенной на исходную урожайность
Участок
|
Площадь, га
|
Культура
|
«Стартовая» урожайность, ц/га
|
Выход продукции за счет «стартовой» урожайности по всей площади,
т
|
1
|
2
|
3
|
4
|
5
|
4
|
71
|
Озимая пшеница
|
22
|
156,2
|
5
|
64
|
Картофель
|
110
|
704
|
Всего
|
135
|
-
|
-
|
860,2
|
Таблица 4
Обеспечение заданных объемов
производства продукции
Вариант
|
Задание по производству продукции, т
|
|
Всего
|
Из них
|
|
|
Выполняется за счет исходного урожая
|
Подлежит выполнению за счет прироста урожая
|
|
Ячмень
|
Озимая пшеница
|
ячмень
|
Озимая пшеница
|
Озимая пшеница
|
Картофель
|
1
|
2
|
3
|
4
|
5
|
6
|
7
|
1
|
206,0
|
1024,0
|
156,2
|
704
|
49,8
|
320
|
|
|
|
|
|
|
|
|
3. Решение экономико-математической
задачи распределения фондов минеральных удобрений сельскохозяйственной
организации по полям севооборотов и кормовым угодьям
Разработанная экономико-математическая
модель может быть решена симплексным методом, так как является задачей
линейного программирования. Общий принцип, которых таков: выбирается
неоптимальный опорный план и его параметры варьируются с целью
последовательного улучшения плана, т.е. оптимизации целевой функции при
соблюдении всех ограничений, что дает возможность решать оптимизационные
задачи.
Рассмотренная задача решена с
использованием приложения Поиск решения MS Excel.
Рис.1 Фрагмент рабочего
листа MS Excel c исходными данными по участку №1
Рис.1a Фрагмент рабочего
листа MS Excel c исходными данными по участку №1(продолжение)
Рис.2. Фрагмент рабочего
листа MS Excel c исходными данными по участку №2 и ограничениями связующего
блока по фондам удобрений и заданиям по производству продукции
Рис.2а. Фрагмент
рабочего листа MS Excel c исходными данными по участку №2 и ограничениями
связующего блока по фондам удобрений и заданиям по производству продукции
(продолжение)
По результатам решения
получены дозы внесения удобрений по различным культурам, в установленные сроки,
приросты урожайности культур по интервалам внесения и общий прирост урожайности
и величина дополнительного чистого дохода (см. рис.3.).
Рис.3. План
распределения фондов удобрений сельскохозяйственного предприятия по полям
севооборотов и кормовым угодьям
Внешний вид
экономико-математической модели по распределению фондов минеральных удобрений
сельскохозяйственного предприятия по полям севооборотов и кормовым угодьям
представлен в приложении 1.
.1 Формирование отчетов
по результатам решения
Если средство Поиск
решения нашло решение, то Excel предоставляет возможность получения таких
отчетов, как Результаты, Устойчивость, Пределы.
В отчете Результаты
выводятся исходные и полученные в результате поиска решения значения изменяемых
ячеек и целевой функции, а также сведения об ограничениях задачи.
Отчет Устойчивость дает
основную информацию для анализа чувствительности линейных и нелинейных моделей.
Этот анализ показывает, насколько чувствительно оптимальное решение к небольшим
изменениям параметров модели. Этот тип отчета будет недоступен, если в модели
используются ограничения целочисленности.
Отчет Пределы
представляет собой ограниченный вариант отчета Устойчивость. Здесь показаны
наименьшее и наибольшее значения, которые может принимать каждая изменяемая
переменная целевая функция. Этот тип отчета будет недоступен, если в модели
используются ограничения целочисленности.
Для выполнения анализа
полученного решения формируются следующие отчеты:
по результатам (рис.4)
по устойчивости (рис.5)
по пределам (рис.6)
Рис. 4. Отчет по
результатам
В отчете Результаты,
показанном на рис.4 содержатся следующие данные:
адреса целевой ячейки и
изменяемых ячеек и их имена (если они заданы)
И значения в этих
ячейках до начала выполнения Поиск решения и после завершения.
адреса ячеек, на значения
которых налагаются ограничения, имена этих ячеек (если они заданы), значения в
этих же ячейках, формулы ограничений, статус ограничения (связанное или
несвязанное) и значения разностей.
Значения разностей -
абсолютные разности между вычисленными значениями правых и левых частей
неравенств. Если значение разности для некоторого ограничения равно нулю, то
это ограничение называется связанным, или лимитирующим, поскольку оно
лимитирует найденное решение. Если значение разности для ограничения не равно нулю,
то такое ограничение называется не связанным или не лимитирующим, - найденное
решение не зависит от этого ограничения.
Рис.5. Отчет по
устойчивости
Отчет по устойчивости
показан на рис.5. Такой отчет доступен только для задач, которые не имеют
ограничений целочисленности. В этом отчете содержатся следующие данные.
В таблице Изменяемые
ячейки приведена информация о значениях изменяемых ячеек:
адреса и имена
изменяемых ячеек;
значения этих ячеек,
найденные средством Поиск решения;
нормированная стоимость,
показывающая, насколько изменится значение целевой функции, если на единицу
изменится значение в данной изменяющейся ячейке при условии, что это значение
достигло своей верхней или нижней границы;
целевой коэффициент -
коэффициент, стоящий при данной изменяемой переменной в уравнении целевой
функции;
значения в столбцах
Допустимое увеличение и Допустимое уменьшение показывают, в каких пределах
может изменятся целевой коэффициент при условии, что найденное значение целевой
функции останется неизменным.
В таблице Ограничения
приведена информация об ограничениях:
адреса и имена ячеек, на
значения которых наложены ограничения;
значения в этих ячейках,
найденные средством Поиск решения;
теневая цена показывает,
насколько изменится значение целевой функции, если на единицу изменится
значение правой части данного ограничения при условии, что этот изменение лежит
в пределах, указанных в столбцах Допустимое увеличение и Допустимое уменьшение;
значения правых частей
ограничений;
значения в столбцах
Допустимое увеличение и Допустимое уменьшение показывают, в каких пределах
может изменяться значение правой части ограничения при условии, что найденное
значение целевой функции останется неизменным.
Рис. 6. Отчет по
пределам
Отчет по пределам
показан на рис.6. Такой отчет доступен только для задач, которые не имеют
ограничений целочисленности. В этом отчете показано значение в целевой ячейке.
Для ограничений приведены значения в изменяемых ячейках и возможные значения
целевой функции, если значения в изменяемых ячейках и возможные значения
целевой функции, если значения в изменяемых ячейках достигнут допустимых нижних
и верхний границ (определяются заданными ограничениями).
4. Анализ результатов
решения
Основные результативные
показатели приведены в табл. 5.
Размер дополнительного
чистого дохода от применения удобрений под прирост урожая составил
86439,907руб.
Таблица 5
Основные результативные
показатели использования удобрений по оптимальному плану
Культура
|
Озимая пшеница
|
Картофель
|
Итого
|
|
Количество
|
1
|
10,96
|
-
|
|
Площадь, га
|
71
|
64
|
135
|
|
Стоимость, руб/га
|
700
|
9864
|
-
|
|
Затраты на уборку, руб/га
|
41,54
|
671,3
|
712,84
|
|
Затраты на уборку в расчете на всю площадь, руб.
|
2949,34
|
42963,2
|
45912,54
|
|
Аммиачная селитра, основное внесение
|
0,67
|
0,59
|
-
|
|
Аммиачная селитра, подкормка
|
0,105
|
0,73
|
-
|
|
Карбамид, основное внесение
|
0,127
|
0
|
0,127
|
Карбамид, подкормка
|
0,078
|
0,054
|
-
|
|
Суперфосфат
|
1,92
|
0,41
|
-
|
|
Калийная соль
|
0,805
|
0,19
|
-
|
|
Суммарные затраты на приобретение, транспортировку и внесение
удобрений, р./га
|
941,54
|
8343,3
|
9284,84
|
|
Суммарные затраты на приобретение, транспортировку и внесение
удобрений в расчете на всю площадь, р.
|
66849,34
|
91442,57
|
158291,91
|
|
Дополнительный чистый доход, отнесенный на удобрение, р./га
|
658,46
|
9192,7
|
9851,16
|
|
Дополнительный чистый доход, отнесенный на удобрение на всю
площадь, р.
|
32941,407
|
53498,5
|
86439,907
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
С точки зрения критерия
оптимальности (дополнительного чистого дохода), более эффективным является внесение
удобрений под картофель:
прирост урожайности озимой пшеницы
обеспечивает, но не превышает, задание по гарантированному производству зерна;
производство картофеля превышает
задание на … т. Дальнейший прирост урожайности ограничивается фондами аммиачной
селитры и карбамида.
Список используемой литературы
1. Математическое моделирование экономических процессов в
сельском хозяйстве/под ред. А.М. Гатаулина. -М.: Агропроиздат, 1990.
2. Математическое программирование/Под редакцией Н.Ш.Кремера
- М.: Финстатинформ,1995. -346с.
. Бережная Е.В., Бережной В.И. Математические методы
моделирования экономических систем: учеб. пособие. - М: Финансы и статистика,
2005. -432с.
. Качанова, Л.С. Экономико-математические методы в
экономике: методические рекомендации по выполнению лабораторных работ. - 2-у
изд., перераб.-М.: ФГОУ ВПО МГАУ, 2010, 68с.
. Качанова Л.С., Вуколов М.В. Моделирование распределения
фондов минеральных удобрений сельскохозяйственной организации: методические
рекомендации по выполнению курсовой работы. -М.: ФГБОУ ВПО МГАУ, 2013. - 36с.
Приложение 1