Контрольная работа по «Прикладной математике»

Министерство  экономического развития Российской Федерации

Московский  Государственный  университет коммерции

Ростовский  институт (филиал) 
 
 
 
 
 
 

Контрольная работа 

по  дисциплине «Прикладная математика» 
 
 
 
 
 

                                                 Выполнила: Бражкина Вера А.

                                                                      Группа КФ-1 З/О

                                          Проверил:   Дерезина Н.П.  
 
 
 
 
 

Ростов-на-Дону

2011г. 

 

Задачa № 2.

Построить на плоскости область решений системы линейных неравенств и найти максимальное и минимальное значения линейной функции в этой области:

Решение

     Построим на плоскости Х12  многоугольник решений рис 1. Для этого в неравенствах системы ограничений и условиях неотрицательности переменных знаки неравенств заменим на знаки точных равенств:

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

      Рис 1. 

     Многоугольником решений задачи является пятиугольник АВСDЕ, координаты точек которого удовлетворяют условию неотрицательности и неравенствам системы ограничений задачи.

Для нахождения точек экстремума построим начальную  прямую F(X)=0=14x1+6x2 и вектор N (14,6) Передвигая прямую F(X)-0 параллельно самой себе в направлении вектора N, найдем точку D,в которой начальная прямая принимает положение опорной прямой. Следовательно, в точке D целевая функция принимает максимальное значение, так как точка D получена в результате пересечения прямых 1 и 5 , то ее координаты удовлетворяют уравнениям этих прямых: 

 

Решив систему уравнений, получим: x1=0; x2=7, а F(X)=14*0+6*7=42.

Для нахождения минимального значения целевой функции  задачи перемещаем начальную прямую в направлении, противоположном вектору N. Начальная прямая займет положение опорной прямой в вершине A. Целевая функция принимает минимальное значение в угловой точке A: 

 

где x1=5/6=0,833, x2=0,833,  а F(X)=14*0,833+6*0,833=16,66.

  Найдем координаты угловых точек:  В, С и Е. Для этого решим следующие системы уравнений: 

            

В результате получим координаты точек В(1;1), C(1;4,66) и E(0;5).

Вычислим  значения целевой функции во всех угловых точках многоугольника решений АВСDЕ: 

A(0,833;0,833), F(X)=14*0,833+6*0,833=16,66 (min),

В(1;1), F(X)=14*1+6*1=20,

C(1;4,66), F(X)=14*1+6*4,66=41.96,

D(0;7), F(X)=14*0+6*7=42(max),

E(0;5), F(X)=14*0+6*5=30. 

Задачa № 28

Для реализации трех групп товаров коммерческое предприятие располагает тремя видами ограниченных материально-денежных ресурсов в количестве b1, b2, b3 единиц. При этом для продажи 1 группы товаров на 1 тыс. руб. товарооборота расходуется ресурса первого вида в количестве a11 единиц, ресурса второго вида в количестве a21 единиц, ресурса третьего вида в количестве a31 единиц. Для продажи 2 и 3 групп товаров на 1 тыс. руб. товарооборота расходуется соответственно ресурса первого вида в количестве  а12, а13 единиц, ресурсов второго вида в количестве а22, а23 единиц, ресурсов третьего вида в количестве а32, а33 единиц. Прибыль от продажи трех групп товаров на 1 тыс. руб. товарооборота составляет соответственно c1, c2, с3 (тыс. руб.).

Определить  плановый объем и  структуру товарооборота так, чтобы прибыль торгового предприятия была максимальной, где:

а11=16, a12=18, a13=9, a21=7, а22=7, а23=2, а31=9, a32=2, а33=3, b1=520, b2=140, b3=810, с1=8, c2=6, c3=4. 

Решение 

Для решения  задачи симплексным методом построим таблицу затрат и объёма ресурсов:

Таблица 1

Виды  материально-денежных ресурсов Норма затрат материально-денежных ресурсов на ед. товарооборота, тыс. руб. Объем ресурсов

b1,2,3,

 
 
1 группа 2 группа 3 группа  
 
Ресурс первого вида 16 18 9 520
Ресурс второго вида 7 7 2 140
Ресурс третьего вида 9 2 3 810
Прибыль, тыс. руб. 8 6 4 max
 

1. Запишем  математическую модель задачи.

Определить  X =(х123), который удовлетворяет условиям

И обеспечивают максимальное значение целевой функции:

Для построения первого опорного плана систему  неравенств приведем к системе уравнений. 

 

В матрице  этой системы уравнений  А(аij) имеет: 
 

 

Векторы A4, A5, A6, - линейно независимы, так как определитель, составленный из компонент этих векторов, отличных от нуля: 

 

Решим систему уравнений относительно базисных переменных. 

Функцию цели запишем в виде: 

 

2. Полагая, что свободные переменные x1=0, x2=0, х3=0, получим первый опорный план X1=(0,0,0,520,140,810), F(X1)=0, в котором базисные переменные x4=520, х5=140, x6=810, следовательно, товары не продаются и прибыль равна нулю, а ресурсы не используются.

Заносим первый опорный план 1 в симплексную таблицу 2

 

Таблица 2 

Симплексная таблица  2
 
План
Базисные переменные Ресурсы плана Значения  коэффициентов при переменных

при

 
 
 
 
 
 
 
 
x1
 
x2
 
x3
 
x4
 
x5
 
x6
 
 
I план x4

x5

x6

520

140

810

16

7

9

18

7

2

9

2

3

1

0

0

0

1

0

0

0

1

   32,5

    20

    90

Инд. Строка  
F(X1)
 
0
 
-8
 
-6
 
-4
 
0
 
0
 
0
 
II план x4

x1

x6

200

20

630

0

1

0

2

1

-7

4,429

0,286

4,429

1

0

0

-2,286

0,143

-1,286

0

0

1

45,161

70

1470

Инд. Строка  
F(X2)
 
160
 
0
 
2
 
-1,714
 
0
 
    1,143
 
0
 
II план x3

x1

x6

45,161

7,097

610,645

0

1

0

0,452

0,871

-7.194

1

0

0

0,226

-0,065

-0,097

-2,286

0,143

-1,286

0

0

1

 
Инд. Строка  
F(X3)
 
237,419
 
0
 
2,774
 
0
 
0,387
 
    0,258
 
0
 
 

3. Первый  опорный план I не оптимальный,  так как в индексной строке находятся отрицательные коэффициенты  -8, -6, -4.

4. За  ведущий столбец выберем столбец,  соответствующий переменной х1, так как сравнивая по модулю имеем:

|-8| > {|-6|, |-4|} Рассчитываем  значения  по строкам, как частное от деления 

 и выбираем наименьшее:

 

Следовательно, вторая строка является ведущей

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

5. Формируем  следующую симплексную таблицу. Вместо переменной х5 в план II войдет переменная х1. Строка, соответствующая переменной х1 в плане II, получена в результате деления всех элементов строки х5 плана I на разрешающий элемент РЭ=7.  На месте разрешающего элемента в плане II получаем 1. В остальных клетках столбца х1 плана II записываем нули.

Таким образом в новом плане II заполнены строки х1 и столбец х1. Все остальные элементы нового плана II, включая элементы индексной строки определяется по правилу прямоугольника. Для этого выбираем из старого плана 4 числа, которые расположены в вершинах прямоугольника и всегда включают разрешающий элемент РЭ=8. Во второй вершине по диагонали находится старое значение элемента, например, значение целевой функции F(X1)=0=СЭ, которое указывает на место расположение нового НЭ в новом плане II. Третий элемент А=140 и четвертый элемент В=-8 завершают построение прямоугольника в недостающих двух вершинах и расположены по другой диагонали. Значение нового элемента в плане II находится из выражения:

                                          140* (-8)

НЭ=СЭ - (А - В)/РЭ=0-       7       =160.

Элементы  строки определяются аналогично 

      140* 16

520 -    7        =200, 
 

          16* 7

18 -      7        =2, 
 

          16* 2

9 -       7        =4,429, 
 

         1* 16

0 -      7        =-2,286. 
 

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

6. План II не является оптимальным, так как есть коэффициенты в индексной строке < 0.

7. Аналогично получаем план III. Он оптимален, так как индексные коэффициенты > 0.

Его можно записать так:

X=(7,097;0; 45,161;0; 0; 610,645),                                                                                                                                                                                                                  F( X *)=237,419 тыс. руб.

Следовательно, необходимо продавать товаров первой группы  7,097 ед., а третьей группы 45,161. При этом торговое предприятие получает максимальную прибыль в размере 237,419 тыс. руб. Товары второй группы не реализуются.

В оптимальном  плане среди базисных переменных находится дополнительная переменная x6. Это указывает, что ресурсы третьего вида недоиспользованы на 610,645 ед., так как переменная х6 была введена в третье ограничение.

В индексной  строке оптимального плана в столбцах переменных x2, x4, x5, не вошедших в состав базисных, получены ненулевые элементы, поэтому оптимальный план задачи линейного программирования является единственным.

 

Задача  № 51

Используя вариант предыдущего  контрольного задания  № 28, необходимо:

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

установить  сопряженные пары переменных прямой и  двойственной задач,

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

Решение 

Определим Y = (y1, у2, y3), который удовлетворяет условиям - ограничениям: 

 

и обеспечивает минимальное значение целевой функции:  

 

Таким образом оптимальный план двойственной задачи имеет вид: 

Y° = (0,387; 0,258; 0; 0; 2,774; 0)        Z(Y°) =237,419. 

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

 

Задача № 76

     Поставщики товара - оптовые коммерческие предприятия А1, А2,…, Аm, имеют запасы товаров соответственно в количестве а1, а2,...,

aп, ед. и розничные торговые предприятия В12,…,Вn, подали заявки на закупку товаров в объемах соответственно: b1 b2, b3,...,bп. Тарифы перевозок единицы груза с каждого из пунктов поставки в соответствующие пункты потребления заданы в виде матрицы.

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

     a1=222, b1=125,

     a2=188, b2=75,

     a3=210, b3=200,

     a4=380, b4=380,

           b5=220, 

       

Решение

Проверим  необходимое и достаточное условие  разрешимости задачи:

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

 

Таблица 3. Опорный план I.

      B 

A     

               a       b

B1 B2 B3 B4 B5 Потенциалы  α 

α1=0

b1=125 b2=75 b3=200 b4=380 b5=220
A1 a1=222 23 21 11 8 

2

3 

220

 
α2=-6
A2 a2=188 7 17 5 2 

188

4  
α3=-4
A3 a3=210 2 

125

16 8 4 

85

3  
α4=0
A4 a4=380 3 9 

75

21 

200

8 

105

4  
 
Потенциалы

β

 
β1=6
 
β2=9
 
β3=21
 
β4=8
 
β5=3
 
 

1. Используя  метод наименьшей стоимости, построим первый горный план транспортной задачи.

Среди тарифов из всей таблицы первым наилучшим является с24 =2, поэтому в клетку A2B4 направляем максимально возможный груз. Он равен min{188,380}=188. Тогда х24=188 и из базы A2 вывезен весь груз, а потребность пункта В4, удовлетворена не полностью, на 188 ед. и в него необходимо поставить 380-188=192 ед. товара. Строка А2 выходит из рассмотрения. Из оставшихся тарифов наименьший – с31=2. В клетку A3B1 направляем максимально возможный груз, равный min{210,125}= 125 и х31=220. Тогда в строке A3 остается 210-125=85 ед., а потребность столбца В1 удовлетворена полностью и он выходит из рассмотрения. Из оставшихся тарифов наименьший – с15=3. В клетку A1B5 направляем максимально возможный груз, равный min{222,220}= 220 и х15=220. Тогда в строке A1 остается 222-220=2 ед., а потребность столбца В5 удовлетворена полностью и он выходит из рассмотрения. Среди тарифов из всех оставшихся наилучшим является с34 =4, поэтому в клетку A3B4 направляем максимально возможный груз. Он равен min{85,192}=85 (остатки товара и потребностей). Тогда х34=85 и из базы A3 вывезен весь груз, но потребность пункта В4, удовлетворена не полностью, а на 85 ед. и в него необходимо поставить 192-85=107 ед. товара. Строка А3 выходит из рассмотрения. Далее: наилучшим является с14 =8, поэтому в клетку A1B4 направляем максимально возможный груз. Он равен min{2,107}=2 (остатки товара и потребностей). Тогда х14=2 и из базы A1 вывезен весь груз, но потребность пункта В4, удовлетворена не полностью, а на 2 ед. и в него необходимо поставить ещё 107-2=105 ед. товара. Строка А1 выходит из рассмотрения. Из оставшихся тарифов наименьший – с44=8. В клетку A4B4 направляем максимально возможный груз, равный min{380,105}= 105 (остаток потребности В4) и х44=105. Тогда в строке A4 остается 380-105=275 ед., а потребность столбца В4 удовлетворена полностью и он выходит из рассмотрения. Осталось два столбца и из них наилучший тариф – с42=9. В клетку A4B2 направляем максимально возможный груз, равный min{275,75}= 75  и х42=75. Тогда в строке A4 остается 275-75=200 ед., а потребность столбца В2 удовлетворена полностью и он выходит из рассмотрения. Остатки ресурсов и потребностей равны, и записываются в единственную доступную клетку А4В3, где с43=21, и х43=200. В результате получен первый опорный план, который является допустимым, так как все грузы из пунктов вывезены, потребность удовлетворена, а план удовлетворяет системе ограничений транспортной задачи. Построение опорного плана I окончено.

2. Подсчитаем число занятых клеток таблицы, их -8, а должно быть m+n-1=4+5-1=8. Следовательно, опорный план является невырожденным.

3. Определяем значение целевой функции первого опорного плана 

F(X)=2*8+220*3+188*2+125*2+85*4+75*9+200*21+105*8=7357. 

Проверим  оптимальность опорного плана. 

4.Найдем  потенциалы αi βj, по занятым клеткам таблицы, решая систему уравнений, полагая αi+ βj = cij и α1= 0: 
 

 

5. Занесем рассчитанные потенциалы в таблицу 3.- подсчитаем оценки     свободных     клеток,     полагая     что    для     них Δij=cij – (αi+ βj):

Δ11=23-6=17; Δ12=21-9=12; Δ13=11-21=-10;

Δ21=7-0=7; Δ22=17-3=14; Δ23=5-15=-10; Δ25=4+3=7;

Δ32=16-5=11; Δ33=8-17=-9; Δ35=3+1=4;

Δ41=3-6=-3; Δ45=4-3=1;

Первый  опорный план является не оптимальным, так как Δ13<0, Δ23<0, Δ33<0, Δ41<0,  поэтому переходим к его улучшению. Выбираем максимальную по модулю оценку свободной клетки – Δ23=| - 10| = 10.

6. Для клетки А2В3 построим цикл перераспределения груза. Для этого в перспективную клетку A2B3 поставим знак +, а в остальных вершинах многоугольника чередующиеся знаки -,+,-.

А2В3                         А2В4

   +                   - 
 
 

А4В3                А4В4

    -      +

Затем из чисел xij, стоящих в минусовых клетках, выбираем наименьшее, т е min{188,200}= 188, Прибавляем 188 к объемам грузов, стоящих в плюсовых клетках и вычитаем 188 из xij, стоящих в минусовых клетках. В результате получим новый опорный план II. 

Таблица 4. Опорный план II.

      B 

A     

               a       b

B1 B2 B3 B4 B5 Потенциалы  α 

α1=0

b1=125 b2=75 b3=200 b4=380 b5=220
A1 a1=222 23 21 11 8 

2

3 

220

 
α2=-16
A2 a2=188 7 17 5

+

188

2

      -

4  
α3=-4
A3 a3=210 2 

125

16 8 4 

85

3  
α4=0
A4 a4=380 3 9 

75

21

12

  -

8 

293   +

4  
 
Потенциалы

Β

 
β1=6
 
β2=9
 
β3=21
 
β4=8
 
β5=3
 
 
 

7. Определяем значение целевой функции II опорного плана 

F(X2)=F(X) +(-10)*188=7357- 1880=5477. 

Проверим  оптимальность опорного плана. 

8. Число занятых клеток не изменилось и равно 8 – следовательно план  невырожденный.  

9. Найдем потенциалы αi, βj, по занятым клеткам таблицы, решая систему уравнений, полагая αi+ βj = cij и α1= 0:

 

10. Занесем рассчитанные потенциалы в таблицу 4.- подсчитаем оценки     свободных     клеток,     полагая     что    для     них Δij=cij – (αi+ βj):

Δ11=23-6=17; Δ12=21-9=12; Δ13=11-21=-10;

Δ21=7+10=17; Δ22=17+7=24; Δ24=2+8=10; Δ25=4+13=17;

Δ32=16-5=11; Δ33=8-17=-9; Δ35=3+1=4;

Δ41=3-6=-3; Δ45=4-3=1;

Второй опорный план является не оптимальным, так как Δ13<0, Δ33<0, Δ41<0,  поэтому переходим к его улучшению. Выбираем максимальную по модулю оценку свободной клетки – Δ13=| - 10| = 10.

11. Для клетки А1В3 построим цикл перераспределения груза. Для этого в перспективную клетку A1B3 поставим знак +, а в остальных вершинах многоугольника чередующиеся знаки -,+,-.

А1В3                         А2В4

   +                   - 
 
 

А4В3                А4В4

    -      +

Затем из чисел xij, стоящих в минусовых клетках, выбираем наименьшее, т е min{2,12}= 2, Прибавляем 2 к объемам грузов, стоящих в плюсовых клетках и вычитаем 2 из xij, стоящих в минусовых клетках. В результате получим новый опорный план III.

Таблица 5. Опорный план III.

      B 

A     

               a       b

B1 B2 B3 B4 B5 Потенциалы  α 

α1=0

b1=125 b2=75 b3=200 b4=380 b5=220
A1 a1=222 23 21 11

   +

2

8

        -

3 

220

 
α2=-6
A2 a2=188 7 17 5

 

188

2

     

4  
α3=6
A3 a3=210 2 

125

16 8 4 

85

3  
α4=10
A4 a4=380 3 9 

75

21

10

  -

8 

295   +

4  
 
Потенциалы

Β

 
β1=-4
 
β2=-1
 
β3=11
 
β4=-2
 
β5=3