Построение экономико-математической модели

Система переменных данной экономико-математической задачи представлена основными и вспомогательной переменными. За основные переменные принимаются:

- площадь посева j -ой сельскохозяйственной культуры;

- поголовье j -го вида сельскохозяйственных животных;

- объем приобретения j -го вида корма или кормовой добавки;

x1	-	Озимая пшеница товарная
x2	-	Озимая пшеница фуражная
x3	-	Ячмень товарный
x4	-	Ячмень фуражный
x5	-	Соя
x6	-	Горчица
x7	-	Сахарная свекла
x8	-	Подсолнечник
x9	-	Кукуруза на силос
x10	-	Кукуруза на зеленый корм 1 срока сева
x11	-	Кукуруза на зеленый корм 2 срока сева
x12	-	Однолетние травы на сено
x13	-	Однолетние травы на сенаж
x14	-	Однолетние травы на зеленый корм 1 срока сева
x15	-	Однолетние травы на зеленый корм 2 срока сева
x16	-	Многолетние травы на сено
x17	-	Многолетние травы на сенаж
x18	-	Многолетние травы на зеленый корм
x19	-	Озимые на зеленый корм
x20	-	Пар
x21	-	Комбикорм
x22	-	Патока
x23	-	Премикс
x24	-	Молоко на корм
x25	-	Обрат
x26	-	ЗЦМ
x27	-	Поголовье крупного рогатого скота

Вспомогательная переменная () используется для определения суммы производственных затрат по предприятию.

x28

-

Производственные затраты

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

· по использованию ресурсов;

· по выполнению агротехнических требований возделывания сельскохозяйственных культур;

· по выполнению договорных обязательств по реализации продукции;

· по обеспеченности потребностей животноводства кормами собственного производства;

· по определению стоимости производственных затрат.

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

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

,

где - площадь посева j -ой сельскохозяйственной культуры;

- поголовье j -го вида скота;

- затраты ресурсов i -го вида в расчете на 1 га посева j -ой сельскохозяйственной культуры;

- затраты ресурсов i -го вида в расчете на 1 структурную голову j -го вида скота;

- наличие ресурсов i -го вида.

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

,

где - затраты пашни в расчете на 1 га посева j -ой сельскохозяйственной культуры;

- наличие пашни в хозяйстве.

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

1X₁+1X₂+1X₃+1X₄+……+1X₁₇+1X₁₈+1X₁₉+1X₂₀<=5400 i=1.

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

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

,

где - бинарные коэффициенты связи (равны либо 0, либо 1) по i -му виду сельскохозяйственных животных;

- фактическое поголовье i -го вида сельскохозяйственных животных в хозяйстве.

То есть сохранение существующего поголовья крупного рогатого скота обеспечивается через реализацию условия

1X₂₇=450 i=2.

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

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

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

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

,

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

Например, ограничение по насыщению севооборота зерновыми и зернобобовыми культурами может быть записано следующим образом:

1X₁+1X₂+1X₃+1X₄+1X₅+1X₁₉ > =2700 i=3,

1X₁+1X₂+1X₃+1X₄+1X₅+1X₁₉ <=3240 i=4.

Ограничение по насыщению севооборота озимыми:

1X₁+1X₂+1X₁₉ > =1080 i=5,

1X₁+1X₂+1X₁₉ <=1620 i=6.

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

Информация о пределах насыщения севооборота отдельными сельскохозяйственными культурами или группами культур берется из таблицы 21 (столбцы 2 и 4).

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

,

где - площадь посева j -ой озимой культуры;

- площадь посева j -ой сельскохозяйственной культуры;

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

Данное ограничение может быть записано следующим образом

1X₁+1X₂-1X₅-1X₁₀-1X₁₂-1X₁₃-1X₁₄-1X₁₅-0,33X₁₆-0,33X₁₇-0,33X₁₈-1X₂₀<=0 i=12.

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

Вторую группу дополнительных ограничений (с 13 по 16) составляют ограничения по выполнению договорных обязательств по реализации продукции.

Эти ограничения гарантируют, что товарной продукции по тем ее видам, на реализацию которой заключены договора, будет произведено не менее заданных объемов

,

где X_j - площадь посева j -ой сельскохозяйственной культуры или поголовье j -го вида скота;

r_ij - выход продукции i -го вида в расчете на 1 га посева j -ой сельскохозяйственной культуры или 1 структурную голову j -го вида скота;

- минимально необходимый объем производства товарной продукции i- го вида.

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

37,68X₁>=23000 i=13,

31,98X₃>=12000 i=14,

482X₇>=195000 i=15.

Ограничение по реализации молока будет иметь несколько иной вид

-1X₂₄+51,10X₂₇>=20000 i=16.

Поскольку потребность в молоке на корм планируется удовлетворять за счет молока собственного производства, а его объем будет определяться в ходе решения задачи, то в данном ограничении из объема валового производства молока (51,10X₂₇) необходимо вычесть потребность в молоке на корм (X₂₄).

Информация о выходе товарной продукции с 1 га посева соответствующих сельскохозяйственных культур берется из столбца 7 таблицы 24.

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

,

где - выход корма i -го вида в расчете на 1 га посева j -ой сельскохозяйственной культуры;

- содержание питательных веществ i -го вида в единице приобретаемого j -ого вида корма или j -ой кормовой добавки;

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

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

1,05X₁+49,28X₂+10,91X₃+47,69X₄+…….+0,34X₂₄+0,13X₂₅+2,40X₂₆-53,02X₂₇>=0 i=17.

По обеспечению рациона кормления переваримый протеином

0,09X₁+4,08X₂+0,46X₃+3,18X₄+…….+0,05X₂₄+0,04X₂₅+42,60X₂₆-5,41X₂₇>=0 i=18.

В ограничениях 17 и 18 информация о выходе корма (с учетом основной и побочной продукции) берется из столбцов 11 и 12 таблицы 25.

В ограничениях 20-22 и 24-25 – из столбца 11 таблицы 25. Например, ограничение по удовлетворению потребности в комбикорме запишется

1,10X₂₁-8,0X₂₇>=0 i=20

Ограничение по удовлетворению потребности в концентрированных кормах собственного производства будет записано как

1,05X₁+49,28X₂+0,83X₃+37,61X₄+0,41X₅ -6,57X₂₇>=0 i=21.

В ограничении по соломе (23) информация о выходе соломы с 1 га посева зерновых и зернобобовых культур берется из столбца 9 таблицы 25

10,08X₃+10,08X₄+3,84X₅-0,09X₂₇>=0 i=23.

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

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

27,90X₁₅+22,60X₁₈-0,28X₂₇>=0 i=29.

Вспомогательным является ограничение по определению стоимости производственных затрат в целом по предприятию:

,

где - материально-денежные затраты в расчете на 1 га посева j -ой сельскохозяйственной культуры;

- цена приобретения j -ого вида корма или j -ой кормовой добавки;

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

Это ограничение может быть записано следующим образом

10,10X₁+10,10X₂+8,50X₃+8,50X₄+…+0,72X₂₁+…+23,40X₂₇ ‑ X₂₈=0 i=37.

Значения переменных должны быть неотрицательными

, , ,

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

,

где - площадь посева j -ой сельскохозяйственной культуры или поголовье j -го вида скота;

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

Целевая функция может быть записана следующим образом:

Z_max=18,09X₁+13,11X₃+18,82X₅+20,54X₆+53,02X₇+22,34X₈+53,33X₂₇ ‑X₂₈.

В матричном виде экономико-математическая модель по оптимизации отраслевой структуры производства представлена в таблице 8 (при распечатке курсового проекта целесообразно распечатать модель из MS Excel и поместить ее в текст в виде приложения.

Реализация разработанной экономико-математической модели производства осуществляется с помощью надстройки Поиск решения пункта меню Сервис Microsoft Excel.

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

Таблица 8. Экономико-математическая модель по оптимизации отраслевой структуры производства (часть 1)

№

Ограничения

Единица измерения

Переменные

Озимая пшеница

Ячмень

Соя

Горчица

Сахарная свекла

Подсолнечник

Кукуруза

товарная

фуражная

товарный

фуражный

на силос

на зеленый корм

1 срока сева

2 срока сева

x1

x2

x3

x4

x5

x6

x7

x8

x9

x10

x11

Значения по решению = >

1477,2

0,0

595,0

12,4

1080,0

0,0

810,0

675,0

105,9

2,0

1,0

Пашня

га

Поголовье КРС

стр.гол.

Площадь зерновых и зернобобовых min

га

Площадь зерновых и зернобобовых max

га

Площадь озимых min

га

Площадь озимых max

га

Площадь горчицы и сои

га

Площадь сахарной свеклы

га

Площадь подсолнечника

га

Площадь пара min

га

Площадь пара max

га

Обеспеченность озимых предшественниками

га

-1

-1

Объем реализации по договорам:

пшеница

ц

37,68

ячмень

ц

31,98

сахарная свекла

ц

482,00

молоко

ц

Требуется корма всего

ц.к.ед.

1,05

49,28

10,91

47,69

4,25

48,15

65,10

65,10

Требуется переваримого протеина

кг

0,09

4,08

0,46

3,18

0,23

3,37

4,34

4,34

Концентраты всего:

ц.к.ед.

1,05

49,28

0,83

37,61

0,41

в т.ч. комбикорм

ц.к.ед.

прочие концентраты

ц.к.ед.

1,05

49,28

0,83

37,61

0,41

Сено

ц.к.ед.

Солома

ц.к.ед.

10,08

10,08

3,84

Сенаж

ц.к.ед.

Силос

ц.к.ед.

48,15

Зеленый корм во 2-ой половине мая

ц.к.ед.

в 1-й половине июня

ц.к.ед.

во 2-й половине июня

ц.к.ед.

в июле

ц.к.ед.

в августе

ц.к.ед.

65,10

в сентябре

ц.к.ед.

65,10

Патока

ц.к.ед.

Премикс

ц.к.ед.

Молоко

ц.к.ед.

Обрат

ц.к.ед.

ЗЦМ

ц.к.ед.

Производственные затраты

тыс.руб.

10,10

10,10

8,50

8,50

8,90

7,90

44,90

7,20

10,80

10,20

10,20

Zmax

18,09

13,11

18,82

20,54

53,02

22,34

Продолжение таблицы 8 (часть 2).

№	Переменные	Сумма произведений	Тип ограничений	Объем ограничений
Однолетние травы	Многолетние травы	Озимые на зеленый корм	Пар	Комбикорм	Патока	Премикс	Молоко на корм	Обрат	ЗЦМ	Поголовье крупного рогатого скота	Производственные затраты
на сено	на сенаж	на зеленый корм	на сено	на сенаж	на зеленый корм
1 срока сева	2 срока сева
x12	x13	x14	x15	x16	x17	x18	x19	x20	x21	x22	x23	x24	x25	x26	x27	x28	bi
	0,0	0,0	2,3	3,0	126,1	235,4	1,9	2,8	270,0	0,0	1823,9	178,0	125,0	365,4	28,0	450,0	86192,0
																		5400,0	=	5400,0
																		450,0	=
																		3167,4	>=
																		3167,4	<=
																		1480,0	>=
																		1480,0	<=
																		1080,0	<=
																		810,0	<=
																		675,0	<=
																		270,0	>=
																		270,0	<=
	-1	-1	-1	-1	-0,33	-0,33	-0,33		-1									0,0	<=
																		55660,6	>=
																		19028,1	>=
																		390420,0	>=
													-1			51,10		22870,0	>=
	13,95	16,32	27,90	27,90	25,50	28,90	56,54	22,40		1,10	0,76	0,50	0,34	0,13	2,40	-53,02		6627,4	>=
	2,08	1,68	3,2	3,72	4,18	3,57	7,97	2,80		11,00	6,00	5,50	0,05	0,04	42,60	-5,41		13165,6	>=
										1,10						1,44		3601,4	>=
										1,10						8,00		3601,4	>=
																-6,57		0,0	>=
	13,95				25,50											-7,14		0,0	>=
																-0,09		10228,7	>=
		16,32				28,9										-15,12		0,0	>=
																-11,34		0,0	>=
								22,40								-0,14		0,0	>=
							33,90									-0,14		0,0	>=
			27,90													-0,14		0,0	>=
				27,90			22,60									-0,28		0,0	>=
																-0,28		0,0	>=
																-0,14		0,0	>=
											0,76					-3,08		0,0	=
												0,50				-0,20		0,0	=
													0,34			-0,09		0,0	=
														0,13		-0,11		0,0	=
															2,40	-0,15		0,0	=
	5,60	5,40	5,20	5,20	6,70	6,50	6,20	4,62	2,70	0,72	0,15	0,55		0,10	0,00	23,40	-1	0,0	=
Zmax																53,33	-1	50673,3	max