x₉

-

Кукуруза на силос

x₁₀

-

Кукуруза на зеленый корм 1 срока сева

x₁₁

-

Кукуруза на зеленый корм 2 срока сева

x₁₂

-

Однолетние травы на сено

x₁₃

-

Однолетние травы на сенаж

x₁₄

-

Однолетние травы на зеленый корм 1 срока сева

x₁₅

-

Однолетние травы на зеленый корм 2 срока сева

x₁₆

-

Многолетние травы на сено

x₁₇

-

Многолетние травы на сенаж

x₁₈

-

Многолетние травы на зеленый корм

x₁₉

-

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

x₂₀

-

Пар

x₂₁

-

Комбикорм

x₂₂

-

Патока

x₂₃

-

Премикс

x₂₄

-

Молоко на корм

x₂₅

-

Обрат

x₂₆

-

ЗЦМ

x₂₇

-

Поголовье крупного рогатого скота

x₂₈

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

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

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

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

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

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

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

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

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

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

,

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

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

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

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

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

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

,

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

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

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

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

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

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

,

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

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

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

1X₂₇=425 i=2.

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

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

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

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

,

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

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

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

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

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

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

1X₁+1X₂+1X₁₉ <=1740 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-го вида.

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

34,74X₁>=23000 i=13,

30,02X₃>=12000 i=14,

400X₇>=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-ого вида сельскохозяйственных животных.

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

0.97X₁+45,44X₂+10,30X₃+44,83X₄+…….+0,34X₂₄+0,13X₂₅+2,40X₂₆-59,14X₂₇>=0 i=17.

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

0,08X₁+3,76X₂+0,43X₃+2,98X₄+…….+0,05X₂₄+0,04X₂₅+42,60X₂₆-5.73X₂₇>=0 i=18.

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

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

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

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

0,97X₁+45,44X₂+0,78X₃+35,31X₄+0,39X₅-6.97X₂₇>=0 i=21.

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

9,60X₃+9,60X₄+3,60X₅-0,10X₂₇>=0 i=23.

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

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

25.90X₁₅+20.90X₁₈-0,30X₂₇>=0 i=29.

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

,

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

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

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

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

10,10X₁+10,10X₂+8,50X₃+8,50X₄+…+0,72X₂₁+…+29.40X₂₇ _ 1X₂₈=0 i=37.

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

, , ,

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

,

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

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

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

Z_max=16,68X₁+12,31X₃+17,64X₅+19,26X₆+44.00X₇+16.76X₈+53,30X₂₇ _X₂₈.

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

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

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

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

№	Ограничения	Единица измерения	Переменные
			Озимая пшеница	Ячмень	Соя	Горчица	Сахарная свекла	Подсолнечник	Кукуруза
			товарная	фуражная	товарный	фуражный					на силос	на зеленый корм
												1 срока сева	2 срока сева
			x1	x2	x3	x4	x5	x6	x7	x8	x9	x10	x11
	Значения по решению = >	1285.8	0,0	1323..1	9.5	858.4	301.6	487.5	725.0	114.5	2,1	1,1
1	Пашня	га	1	1	1	1	1	1	1	1	1	1	1
2	Поголовье КРС	стр.гол.
3	Площадь зерновых и зернобобовых min	га	1	1	1	1	1
4	Площадь зерновых и зернобобовых max	га	1	1	1	1	1
5	Площадь озимых min	га	1	1
6	Площадь озимых max	га	1	1
7	Площадь горчицы и сои	га					1	1
8	Площадь сахарной свеклы	га							1
9	Площадь подсолнечника	га								1
10	Площадь пара min	га
11	Площадь пара max	га
12	Обеспеченность озимых предшественниками	га	1	1			-1					-1
	Объем реализации по договорам:
13	Пшеница	ц	34.74
14	Ячмень	ц			30,02
15	сахарная свекла	ц							400
16	Молоко	ц
17	Требуется корма всего	ц.к.ед.	0.97	45,44	10,30	44,83	3.99				44.55	60.27	60,27
18	Требуется переваримого протеина	кг	0.08	3,76	0,43	2,98	0,21				3,12	4,02	4,02
19	Концентраты всего:	ц.к.ед.	0,97	45,44	0,78	35,31	0,39
20	в т.ч. комбикорм	ц.к.ед.
21	прочие концентраты	ц.к.ед.	0,97	45,44	0,78	35,31	0,39
22	Сено	ц.к.ед.
23	Солома	ц.к.ед.			9,52	9,52	3,60
24	Сенаж	ц.к.ед.
25	Силос	ц.к.ед.									44.55
26	Зеленый корм во 2-ой половине мая	ц.к.ед.
27	в 1-й половине июня	ц.к.ед.
28	во 2-й половине июня	ц.к.ед.
29	в июле	ц.к.ед.
30	в августе	ц.к.ед.										60.30
31	в сентябре	ц.к.ед.											60,30
32	Патока	ц.к.ед.
33	Премикс	ц.к.ед.
34	Молоко	ц.к.ед.
35	Обрат	ц.к.ед.
36	ЗЦМ	ц.к.ед.
37	Производственные затраты	тыс.руб.	10,10	10,10	8,58	8,50	8,90	7,90	44,90	7,20	10,80	10,20	10,20
Zmax		16,68		12,31		17,64	19,26	44.00	16.76

№	Переменные	Сумма произведений	Тип ограничений	Объем ограничений
	Однолетние травы	Многолетние травы	Озимые на зеленый корм	Пар	Комбикорм	Патока	Премикс	Молоко на корм	Обрат	ЗЦМ	Поголовье крупного рогатого скота	Производственные затраты
	на сено	на сенаж	на зеленый корм	на сено	на сенаж	на зеленый корм
			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.5	3,3	133.9	256.6	2,0	3,1	290,0	0.0	1823.9	178.0	125,0	365,4	28,0	425,0	76472.4
1	1	1	1	1	1	1	1	1	1									5 800,0	=	5 800
2																1		425,0	=	425
3								1										3480,0	>=	2 900
4								1										3480,0	<=	3 380
5								1										1 288	>=	1 160
6								1										1 288	<=	1740
7																		1 160	<=	1 160
8																		487.5	<=	870
9																		725	<=	725
10									1									290	>=	290
11									1									290	<=	580
12	-1	-1	-1	-1	-0,33	-0,33	-0,33		-1									0,0	<=	0
13																		44669.4	>=	23 000
14																		39720.30	>=	12 000
15																		195000	>=	195 000
16																51,10		21592.5	>=	20 000
17	13,05	15,36	26,10	26,10	24,00	26,52	52.36	20.80		1,10	0,76	0,50	0,34	0,13	2,40	-56.14		12135.2	>=	0
18	1,94	1,58	3,48	3,48	3,94	3,28	7,38	2,60		11,00	6,00	5,50	0,05	0,04	42,60	-5.73		13366.9	>=	0
19										1,10						1.52		3601.4	>=	0
20										1,10						8.47		0,0	>=	0
21																6.95		0,0	>=	0
22	13,05				24,00											7.56		0.0	>=	0
23																-0,10		15736.5	>=	0
24		15,36				26,86										-16.01		0,0	>=	0
25																-12.00		0,0	>=	0
26								20.80								-0,15		0,0	>=	0
27							31,40									-0,15		0,0	>=	0
28			26,10													-0,15		0,0	>=	0
29				26,10			20.90									-0,30		0,0	>=	0
30																-0,15		0,0	>=	0
31																-0,15		0,0	>=	0
32											0,76					-3,26		0,0	=	0
33												0,50				-0,21		0,0	=	0
34													0,34			-0,10		0,0	=	0
35														0,13		-0,11		0,0	=	0
36															2,40	-0,16		0,0	=	0
37	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	=	0
Zmax																53,277	-1	38456.0	max

Анализ результатов решения задачи

Результаты реализации экономико-математической модели по оптимизации отраслевой структуры производства приведены в таблице 9.

Таблица 9. Результаты решения задачи по оптимизации отраслевой структуры производства

Номер переменной	Наименование переменной	Значение по решению
x1	Озимая пшеница товарная	1285.8
x2	Озимая пшеница фуражная	0,0
x3	Ячмень товарный	1323.1
x4	Ячмень фуражный	9.5
x5	Соя	858.4
x6	Горчица	301.6
x7	Сахарная свекла	487.5
x8	Подсолнечник	725.0
x9	Кукуруза на силос	114.5
x10	Кукуруза на зеленый корм 1-го срока сева	2,1
x11	Кукуруза на зеленый корм 2-го срока сева	1,1
x12	Однолетние травы на сено	0,0
x13	Однолетние травы на сенаж	0,0
x14	Однолетние травы на зеленый корм 1 срока сева	2,5
x15	Однолетние травы на зеленый корм 2 срока сева	3,3
x16	Многолетние травы на сено	133.9
x17	Многолетние травы на сенаж	256.6
x18	Многолетние травы на зеленый корм	2,0
x19	Озимые на зеленый корм	3,1
x20	Пар	290.0
x21	Комбикорм	0.0
x22	Патока	1 823.9
x23	Премикс	178.0
x24	Молоко на корм	125,2
x25	Обрат	365,4
x26	ЗЦМ	28,0
x27	Поголовье крупного рогатого скота	425.0
x28	Сумма производственных затрат	76472.4

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

Все ограничения выполнены. Анализируя результаты, представленные в столбце «Сумма произведений» можно сделать следующие выводы.

Имеющаяся в хозяйстве пашня используется полностью.

Площадь посева зерновых и зернобобовых (ограничения 3 и 4) берется по минимальной границе и составляет 2 500 га. Площадь озимых культур (ограничения 3 и 4) формируется на уровне 1251,4 га (при нижней и верхней границе соответственно 1 000 и 1 500 га). Максимально возможных значений достигает площадь подсолнечника и суммарная площадь посевов сои и гречихи (соответственно 625 и 1 000 га). Площадь свеклы составляет 649,3 га при максимально возможной - 750 га. Площадь пара берется по минимально допустимой границе в 250 га.

Все озимые культуры размещены по хорошим предшественникам. Излишки площадей по таким культурам отсутствуют.

Таблица 10. Выполнение договорных обязательств по реализации продукции, ц

Продукция	Объемы реализации по плану	Объемы реализации по решению	Объемы реализации сверх плана
Пшеница	23 000,0	43 366,0	20 366,0
Ячмень	12 000,0	12 000,0	0,0
Сахарная свекла	95 000,0	292 188,4	97 188,4
Молоко	20 000,0	20 314,8	314,8

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

Анализ обеспеченности отрасли животноводства кормами показывает, что по оптимальному решению излишек кормов составит 15736.5 к.ед. Этот излишек образуется за счет соломы, которая является побочной продукцией при производстве зерновых и зернобобовых культур. Наблюдается избыток в рационе переваримого протеина (около 63 т), что свидетельствует о необходимости корректировки рационов кормления.

Сумма чистого дохода по оптимальному решению составляет 38456.0 тыс. руб., а уровень рентабельности - 50.3%.

Литература

В.А. Колемаев, О.В. Староверов, В.Б. Турундаевский “Теория вероятностей и математическая сатистика”/ М., 1991.

“Теория Статистики” под редакцией Р.А. Шмойловой/ “ФиС”, 1998.

“Многомерный статистический анализ на ЭBM с использованием пакета Microsoft Excel”/ М., 1997.

А.А. Френкель, Е.В. Адамова “Корреляционно регрессионный анализ в экономических приложениях”/ М., 1987.

И.Д.Одинцов “Теория статистики”/ М., 1998.

А.Н. Кленин, К.К. Шевченко “Математическая статистика для экономистов-статистиков”/ М., 1990.

Размещено на Allbest.ru