Математические методы. 2

      Содержание

  1. Постановка экономико-математической задачи
  2. Состав переменных и ограничений задачи, их запись
  3. Математическая запись задачи
  4. Анализ результатов решения
  5. I постановка
  6. II постановка
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

1.Постановка  задачи и исходная  информация.

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

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

Критерием оптимальности  или показателем эффективности при составлении рациона в моей задаче служит минимум стоимости рациона,

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

Табл1.Норма кормления дойной коровы

Живая масса, кг Суточный  удой Рацион  должен содержать не менее
Кормовых  единиц, кг Переваримого  протеина, г Каротина, мг
400 8 8.0 840 320
 

Табл2 . Использование корма, их питательность  и стоимость.

Корма и кормовые

добавки

Содержится в 1 кг корма Стоимость 1 кг

Корма, копеек

Корм.ед., кг Перевар. прот., г Каротина,

мг

Комбикорм 0,9 160 2 300
Мука  вико-овсяная 0,94 135 1 300
Отходы  ячменя 0,9 90 1 280
Дерть кукурузная 1,29 74 3 230
Сено  злаковое разнотравное 0,46 41 15 68
Сено  осоковое 0,40 39 5 42
Солома  викоовсяная 0,30 28 2 40
Силос подсолнечн. 0,16 15 15 40
Силос клеверный 0,16 19 20 46
Картофель 0,3 16 0 150
Кормовая  свекла 0,12 9 0 54
Силос разнотравный 0,13 15 10 38
Жом 0,11 7 0 40

Табл.3. Содержание отдельных групп кормов в рационе.

 
 
Группы  кормов
1 постановка

(в % от  общего количества

кормовых  единиц)

2 постановка

( в физическом  весе)

не  менее не более не менее не более
Концентраты - - 0,7 2,5
Грубые - - 6 11
Силос - - 8 30
Корнеплоды - - - 15
Жом 8 - 1 -
 

Удельный  вес комбикорма в концентрированных  кормах (по массе) может быть не более 40%, соломы в грубых - не более 40 %, силоса подсолнечникового во всем силосе - не менее 25%,картофеля в корнеклубнеплодах  не менее 20%.В общем балансе кормовых единиц жом может занимать не менее 8%.

    2.Состав  переменных и ограничений  задачи, их запись

    Состав  переменных:

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

    Х 1  -  комбикорм

    Х 2   - Мука вико-овсяная

    Х 3   - отходы ячменя

    Х 4 - Дерть кукурузная

    Х 5   - Сено злаковое разнотравное

    Х 6   - Сено осоковое

    Х 7 - Солома викоовсяная

    Х 8  - Силос подсолнечн.

    Х 9  - Силос клеверный

    Х 10 - Картофель

    Х 11  - Кормовая свекла

    Х 12  - Силос разнотравный

    Х 13  - Жом

    Дополнительные  переменные:

    Х14_-общее количество кормовых единиц

    Х15 – общее количество переваримого протеина

Основные  ограничения модели – по балансам питательных веществ в рационе записываются на основании исходных данных таблицы 1 и имеют следующий вид:

1) баланс  кормовых единиц в рационе,  кг

0,9 Х  1 +0,94 Х 2+0,9 Х 3+1,29 Х 4+…+0,11 Х 13 ≥  8,0

2) баланс  переваримого протеина, г

160 Х 1 +135 Х 2+90 Х 3+74 Х 4+…+7 Х 13   ≥840

3) баланс  каротина, мг

2Х  1 + 1Х 2+1Х 3+ …+10Х 12≥320

Так как  дополнительные ограничения записываются в первой постановке задачи в %, то сначала  записывается вспомогательное ограничение, позволяющее определить общее содержание кормовых единиц в рационе, и вводится вспомогательная  переменная Х14 – «всего кормовых единиц в рационе»:

0,9 Х  1 +0,94 Х 2+0,9 Х 3+1,29 Х 4+…+0,11 Х 13= Х14

или после  переноса всех переменных в правую часть:

4) 0,9 Х 1 +0,94 Х 2+0,9 Х 3+1,29 Х 4+…+0,11 Х 13  - Х14 = 0

Запишем  вспомогательное ограничение, позволяющее определить общее содержание переваримого протеина в рационе. Введем вспомогательную переменную  Х 15 - «всего переваримого протеина в рационе» и перенесем все переменные в правую часть:

160 Х  1 +135 Х 2+90 Х 3+74 Х 4+…+7 Х 13    =Х 15

5) 160 Х 1 +135 Х 2+90 Х 3+74 Х 4+…+7 Х 13   - Х15 = 0

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

    Для II постановки записываются в физическом весе:

  1. Содержание концентратов в рационе не менее 0,7 кг.

Х1+Х2+Х3+Х4-0,7>0

  1. Содержание концентратов в рационе не более 2,5 кг.

Х1+Х2+Х3+Х4-2,5<0

  1. Содержание грубых кормов в рационе не менее 6 кг.

Х5+Х6+Х7 -6>0

  1. Содержание грубых кормов в рационе не более 11 кг.

Х5+Х6+Х7-11<0

  1. Содержание силоса в рационе не менее 8 кг.

Х8+Х9+Х12 -8>0

  1. Содержание силоса в рационе не более 30 кг.

Х8+Х9+Х12 -30<0

  1. Содержание корнеплодов в рационе не более 15 кг.

Х10+Х11 -15<0

  1. Содержание жома в рационе не менее 1 кг

Х13-1<0

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

  1. Содержание карбамида в рационе, г.

2600Х14-0,1Х15<0

  1. удельное содержание комбикорма в концентрированных кормах:

Х 1 ≤ 0.60 (Х 1 +Х 2+Х 3+Х 4)

0.40Х 1 - 0.60Х 2- 0.60Х 3- 0.60Х 4≤0

  1. удельное содержание соломы в грубых кормах:

Х 7 ≤ 0.60 (Х 5 +Х 6+Х 7)

-0.6Х 5 –0.6Х 6 +0.4Х 7

  1. удельное содержание силоса подсолнечникового во всем силосе:

Х 8≥0.75 (Х 8 +Х 12+Х 9)

-0,75Х 12 -0.75Х 9+0.25Х 8≥0

  1. удельное содержание картофеля в корнеклубнеплодах:

Х 10  ≥ 0.80 (Х 11 +Х 10)

0.20Х 10 –0.80Х 11≥0;

Условие неотрицательности переменных:

Хj≥0         

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

Z = 300Х1 +300Х2 +280Х3 +230Х4 +…+38Х12+ 40Х13       ≥  min

При вводе данных ставится коэффициент «-1», так как находится минимум целевой функции. 
 

3.Математическая  запись модели

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

    SVijXj>Bi  (iÎI1)

                                 jÎJ

    Где j – индекс переменной;

        J – множество,  включающее номера переменных  по видам кормов (кормовых добавок)  в рационе;

        Xj – количество корма (кормовой добавки) j-го вида в рационе;

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

        Bi – суточная потребность животного в i-м питательном веществе;

        i – индекс ограничения;

        I1 – множество, включающие номера ограничений по содержанию питательных веществ.

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

        SVijXj=Υ

                                      jÎJ

    Где Υ – вспомогательная переменная

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

        Для I постановки:

        Bia*Y<S Xj< Bib*Y   (iÎI2), (hÎH)

                                 jÎ h

        Для II постановки:

        Bia<S Xj< Bib   (iÎI2), (hÎH)

                                 jÎ h

    Где H – множество, включающее номера переменных по видам кормов (кормовых добавок) одной группы;

         h – индекс группы  кормов;

         Bia, Bib - минимально и максимально допустимые пределы введения в рацион данной группы кормов;

        I2 - множество, включающее номера ограничений по содержанию отдельных групп кормов в рационе.

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

        Xj< (>) Wij*Xj (iÎI3), (hÎH)

                                      jÎh

    Где Wij – коэффициент пропорциональности;

         I3 - множество, включающее номера ограничений по удельному весу отдельных видов кормов внутри групп.

    Условия неотрицательности переменных:

    Xj>0

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

    Z=∑CjXj  à min

                                   jÎJ

    Где Cj – стоимость (цена приобретения) единицы  корма j-го вида. 

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

    Постановка 1

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

Табл.4. Состав оптимального кормового решения

№№

переменных

 
 
Виды  кормов
Количест

во

корма,

кг

Содержится  в рационе Стоимость

корма и 

рациона,

коп.

Кормовых  ед.,

кг

Перева-

римого

протеина, г

Каро-

тина,

мг

X1 Комбикорм 2,9831 2,6848 477,2903 5,9661 447,4597
X4 Шрот льняной 0,9944 1,0142 284,3855 -- 120,3169
X6 Сено овсяное 6,5986 3,2333 362,9243 32,9931 211,156
X7 Солома овсян. 2,828 0,8767 39,5917 11,3119 50,9037
X9 Силос подсол. 29,9688 4,795 449,5313 449,531 599,375
X14 Корм.морковь 3,3399 0,4676 23,3794 100,197 100,1976
X15 Жом 5,7128 0,6284 39,9898 -- 114,2566
Итого -- 13,7 1677,0923 600 --
Задаваемая  потреб-

ность в питательных 

веществах

-- 13,7 1500 600 --
Превышение  потреб-

ности (отклонение)

-- 0 177,0923 0 --
 

Табл5. Состав рациона по группам кормов, кг

Группы

кормов

По  условию

По решению
Не  менее Не более Всего

(сумма)

Отклонение

от нижней

границы

Отклонение

от верхней  границы

концентраты 1,37 3,699 3,699 2,329 0
грубые 2,055 4,11 4,11 2,055 0
силос 2,466 4,795 4,795 2,329 0
корнеклуб-

неплоды

0,411 2,055 0,4676 0,0566 1,5874
прочие -- 1,644 0,6284 -- 1,0156
 
 
 
 
 
 
 
 

Табл.6.Состав и структура кормового  рациона по группам  кормов.

Ед.

измерения

 
 
Показатели
Концентраты Грубые  
Силос
Корнеклуб-

неплоды

Прочие  
 
итого
Кг Кормовые  единицы 3,699 4,11 4,795 0,4676 0,6284 13,7
% 27 30 35 3,4 4,6 100
Г Перевар. Протеин 761,67 402,51 449,53 23,379 39,989 1677,0923
% 45,4 24 26,8 1,4 2,4 100
Мг Каротин 5,9661 44,305 449,53 100,19 -- 600
% 1 7,4 74,9 16,7 -- 100
коп. Стоимость 567,77 262,05 599,37 100,19 114,25 1643,6655
% 34,5 15,9 36,5 6,1 7 100
 

Табл.7 .Стоимость единицы  питательного вещества, коп.

Группа  и вид корма 1 кг кормовых

единиц

1 г перевари-

мого  протеина

1 мг каротина

Концентрированные

153,495 0,745 95,167
Комбикорм 166,664 0,938 75,00
Шрот  льняной 118,632 0,423 --

Грубые

63,761 0,651 5,914
Сено  овсяное 65,307 0,582 6,40
Солома  овсяная 58,063 1,286 4,50

Силос (подсолн.)

125 1,333 1,333
К/клубн (морковь) 214,281 4,286 1
Прочие (жом) 181,821 2,857 --
 

Табл.8.Двойственные оценки ограничений.

№

ограничения 

 
 
Название

ограничения

Единица

измерения

 
Объем

Ограничения

(BI)

 
 
Сумма
 
 
Откло-

нение

 
Двойственная

оценка

Баланс кормо-

вых единиц

кг 13,7 13,7 0 -113,3379
Баланс перева-

римого протеина

 
г
1500 1677,0924 177,0924 -0
Баланс каротина мг 600 600 0 0,1515
Общее количество КЕ кг 0 0 0 -68,4783
Общее содержание пер. протеина кг 0 0 0 -0
Концентратов  не менее 0 0 2,329 2,329 -0
Концентратов  не более кг 0 0 0 28,5679
Грубых кормов не менее кг 0 2,055 2,055 -0
Грубых кормов не более кг 0 0 0 119,69
Силоса не менее кг 0 2,329 2,329 -0
Силоса не более кг 0 0 0 71,0227
Корнеплодов не менее  
кг
0 0,056 0,0566 -0
Корнеплодов не более  
кг
0 -1,5874 1,5874 0
Жома не более Кг 0 -1,0156 1,0156 0
Содержание  карбамида  
Кг
0 -251,5639 251,5639 0
Содержание  комбикорма  
кг
0 0 0 -47,087
Содержание  соломы  
кг
0 0 0 2,6654
Содержание  корм. моркови  
кг
0 2,004 2,004 -0
Содержание  картофеля  
кг
0 -0,501 0,501 0