Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Информатика УМК - 1_excel--2101.pdf
Скачиваний:
73
Добавлен:
16.02.2016
Размер:
1.3 Mб
Скачать

12

4. Отчет по работе

Результаты выполнения заданий 13.

Литература: [1], с. 113-116, [6], с. 208-215.

Работа 3. ПЛАНИРОВАНИЕ ВЫПУСКА ПРОДУКЦИИ

1. Цель работы

Научиться использовать возможности Excel в задачах управления и планирования.

2. Основные теоретические положения

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

изменилась.

 

 

 

 

 

Таблица 7

 

 

 

 

 

 

 

 

Объемы

Производственное потребление

Прогнозируе-

Отрасли

производства

отраслей за предыдущий период

мый конечный

 

 

отраслей

1

2

3

спрос

 

1

600

250

100

160

2000

 

2

1000

150

500

0

2000

 

3

800

0

300

400

3000

 

2.1.Математическая постановка задачи

Для решения поставленной задачи можно использовать балансовую модель Леонтьева. Она представляет собой систему уравнений, каждое из которых выражает требование равенства (баланса) между количеством продукции, производимой отдельным экономическим объектом, и совокупной потребностью в этой продукции. В рассматриваемой задаче экономическая система состоит из трех отраслей.

Пусть Хi – величина, равная суммарному выпуску продукции отрасли i;

xij – количество продукции отрасли i, необходимое для того, чтобы отрасль j произвела Xj единиц своей продукции;

X 1 = а11 Х1 + а12 Х2 + а13 Х3 +Y1; X 2 = а21 Х1 + а22 Х2 + а23 Х3 +Y 2 ; X 3 = а31 Х1 + а32 Х2 + а33 Х3 +Y 3.

13

Yi – количество продукции отрасли i, оставшееся для внешнего потребления (конечная продукция).

Тогда взаимосвязь отраслей в процессе производства и потребления отдельного

продуктаХi (i=1, 2, 3) можетбытьописанаввидеследующихуравнений:

 

X 1 = x11 + x12 + x13 + Y1;

(1)

X 2 = x21 + x22 + x23 + Y 2 ;

X 3 = x31 + x32 + x33 + Y 3 .

Используем понятие коэффициентов прямых затрат (технологических коэффициентов) aij:

aij = хij – количество продукции отрасли i, необходимое для того, чтобы

X j

отрасль j произвела одну единицу своей продукции.

Тогда xi j = aijXj и система уравнений (7) будет иметь следующий вид:

(2)

Или в матричной форме

 

a11

a12

a13

 

Х=АХ+Y,

(3)

 

 

 

где

 

a22

a23

 

– матрица прямых затрат,

 

A = a21

 

 

 

a31

a32

a33

 

 

 

Х – вектор-столбец выпуска продукции в предыдущем периоде.

Yвектор-столбец конечного спроса в предыдущем периоде.

2.2.Решение задачи

2.2.1. Определение вектора конечной продукции за предыдущий период

По условию задачи известны объемы производства каждой из отраслей за предыдущий период (суммарный выпуск продукции отрасли i): X1 = 600, X2 = 1000,

X3 = 800 и значения xij (i, j = 1, 2, 3):

х11 = 250;

х12

=100;

х13

=160;

х21 =150;

х22

= 500;

х23 = 0;

х31 = 0;

х32

=300;

х33

= 400.

Отсюда, используя (1), можно определить значения Yi , i = 1, 2, 3, конечной продукции каждой из отраслей за предыдущий период.

Y1 = 600 250 100 160

=90;

 

Y2

=1000 150 500 0 =350;

(4)

Y3

=800 0 300 400 =

100.

 

Таким образом, вектор конечной продукции за предыдущий период найден

Y = (90, 350, 100).

14

Для определения вектора выпуска продукции Х при заданном конечном прогнозируемом векторе спроса Y = (2000, 2000, 3000) надо решить систему

уравнений (3), из которого следует, что

 

Х = (Е-А)-1Y,

(5)

где Е – единичная матрица.

 

Матрица S=(E-A)-1 – называется матрицей полных затрат.

 

2.2.2. Определение коэффициентов прямых затрат

 

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

определим

коэффициенты прямых затрат aij:

а11 =

250

=0,417;

а12 =

 

100

 

=0,1;

а13 =

160

=0,2;

 

 

 

 

 

 

 

 

 

 

600

1000

800

 

 

 

 

 

 

 

 

 

 

а21 =

150

 

=0,25;

а22 =

500

 

=0,5;

а23 =

0

 

=0;

 

 

 

 

 

 

 

 

 

 

 

600

1000

 

 

 

 

 

 

 

 

 

800

 

а31

 

=

2

 

=0;

а32 =

 

300

 

=0,3;

а33 =

 

400

=0,5.

 

 

 

1000

800

 

 

600

 

 

 

 

 

 

Таким образом, матрица коэффициентов прямых затрат будет иметь вид

 

0,417

0,1

0,2

 

 

 

0,25

0,5

0

 

(6)

А=

.

 

0

0,3

0,5

 

 

 

 

 

2.2.3. Проверка продуктивности матрицы

Все элементы матрицы А неотрицательные, А 0.

Для того чтобы система уравнений (5) имела единственное неотрицательное решение при любом неотрицательном векторе Y, необходимо, чтобы матрица А была продуктивной. Экономический смысл продуктивности состоит в том, что существует такой план выпуска продукции, при котором каждая отрасль сможет произвести необходимое количество конечной продукции. Можно показать, что для продуктивности матрицы А необходимо и достаточно, чтобы все главные миноры матрицы (Е-А) были положительными числами строго меньше единицы. Кроме того, известно: если сумма элементов каждого из столбцов неотрицательной квадратной матрицы А положительна и строго меньше единицы, то все главные миноры матрицы (Е-А) положительны и строго меньше единицы.

Суммы элементов каждого столбца матрицы А (6) соответственно равны:

0,417+0,25+0 =0,667;

0,1+0,5 +0,3 =0,9;

0,2 +0 +0,5 =0,7.

Следовательно, в силу вышесказанного, матрица А продуктивна, выражение (5) имеет смысл и вектор Y неотрицателен. Следовательно, для нахождения плана выпуска продукции Х можно воспользоваться формулой (5).

15

2.2.4. Вычисление матрицы Е - А

Вычислим матрицу (Е - А):

 

1

0

0

 

 

0,417

0,1

0,2

 

1

0417

0 0,1

0 0,2

 

 

 

0

1

0

 

 

0,25

0,5

0

 

 

 

10,5

0 0

 

=

(Е А) =

 

 

= 0 0,25

 

 

0

0

1

 

 

 

0

0,3

0,5

 

 

0 0

0 0,3

10,5

 

 

 

 

 

 

 

 

 

(7)

 

 

 

 

 

 

 

0,583

0,1

 

0,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,25

0,5

 

0

 

 

 

 

 

 

 

 

 

 

=

 

.

 

 

 

 

 

 

 

 

 

 

 

0

0,3

 

0,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.2.5. Вычисление обратной матрицы (Е-А)-1

Известно, что матрица В-1 называется обратной по отношению к квадратной матрицеВ, еслипроизведениеВ* В-1 = Е (Е– единичнаяматрица).

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

В1 =

1

(Bij )Т .

(8)

det B

Здесь (Bij) – матрица, полученная из элементов Bij, а Bij

– алгебраические

дополнения элементов матрицы.

 

Bij=(-1)i+j Mij,

(9)

где Mij – минор элемента aij (минор – это такой определитель, который получается из матрицы вычеркиванием строки и столбца, на пересечении которых стоит данный элемент).

Вычислим значения алгебраических дополнений элементов матрицы (Е - А). Обозначим для простоты вычислений Е - А=В

В11

 

=(1)1+1

0,5

0

=0,25;

 

 

В12 = (1)1+2

 

0,25

0

 

= 0,125 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,3

0,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0,5

 

 

 

 

 

 

 

=

 

1+3

 

0,25

0,5

 

=0,075;

В

 

 

2+1

 

0,1 0,2

 

 

=0,11;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13

(1)

 

 

 

 

 

 

 

 

 

=(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

 

 

 

 

 

 

 

 

0

0,3

 

 

 

 

 

 

 

 

 

 

 

 

 

0,3 0,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В22

=(1)2+2

 

0,583

0,2

 

=0,291 ;

 

 

 

2+3

 

0,583

0,1

 

=0,175;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0,5

 

В23=(1)

 

 

 

 

 

0

 

0,3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В31

=

3+1

 

0,1

0,2

 

=

0,175;

В32

= −

 

3+2

 

0,583

 

0,2

 

=0,1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)

 

 

0

 

0,3

 

 

 

( 1)

 

 

 

 

 

 

0,5

 

 

 

 

 

 

0

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,583

0,1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В33 =

(1)

3+3

 

 

=0,276.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,25

0,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, (ЕА)=(

 

 

 

)=

 

 

0,25 0,125 0,075

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,11

0,291

 

0,175 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,1

0,05

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,276

 

 

 

 

 

 

 

 

 

 

 

16

2.2.6. Вычисление транспонированной матрицы

Поменяв в матрице [Е-А] строки и столбцы местами, получаем

 

 

 

= (Вij )T =

 

0,25

0,11

0,1

 

 

 

 

(Е А)Т

 

0,125

0,291

0,05

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

0,075

0,175

0,276

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.2.7. Вычисление определителя матрицы [Е-А]

 

 

 

 

 

 

 

Вычислим определитель, применив разложение по первой строке

 

det (E A) =0,583

 

0,5

0

 

(0,1)

 

 

0,25

0

 

+(0,2)

 

0,25

0,5

 

=

 

 

 

 

 

 

 

 

0,3

0,5

 

 

 

 

0

0,5

 

 

 

 

 

0

0,3

 

 

=0,583(0,5 0,5 +0 0,3) +0,1(0,25 0,5 0 0) 0,2(0,25 0,3 +0 0,5) =0,118.

2.2.8. Вычисление матрицы прямых затрат S

По формуле S=(E-A)-1=B-1= de1t B (Bij )T =

 

1

 

0,25

0,11

0,1

 

 

 

 

0,125

0,291

0,05

 

 

=

 

 

 

=

0,118

 

 

0,075

0,175

0,276

 

 

 

 

 

 

 

0,25

0,11

 

0,1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,118

 

0,118

 

0,118

 

2,118

0,93

0,847

 

 

 

 

 

 

0,125

 

0,291

0,05

 

 

 

 

 

 

=

 

 

 

 

 

 

=

1,059

2,466

0,424

.

0,118

 

0,118

 

0,118

 

 

 

 

 

0,635

1,479

2,253

 

 

0,075

 

0,175

 

0,276

 

 

 

 

0,118

 

0,118

 

0,118

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.2.9. Определение вектора выпуска продукции Х

Зная S и Y, вычислим X по формуле:

 

 

2,118

0,93

0,847

 

 

2000

 

 

 

1,059

2,466

0,424

 

 

2000

 

Х=S Y =

 

 

.

 

 

0,635

1,479

2,253

 

 

3000

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

Х1 = 2,118*2000+0,93*2000+0,847*3000 =8637;

Х2

=1,059* 2000+2,466*2000+0,424*3000 =8322;

Х3

= 0,635*2000+1,479*2000+2,253*3000 =10985.

Таким образом, вектор выпуска продукции в следующем периоде при заданном векторе конечной продукции Y = (2000, 2000, 3000) равен

X = (8637, 8322, 10985).

Очевидно, что с использованием матричных операций в Excel процедура вычислений в балансовой модели существенно упрощается.

17

3. Порядок выполнения работы

Задание. Реализовать балансовую модель в электронной таблице (ЭТ) Excel.

3.1. Выполнение задания

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

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

3.1.1.Создать блок исходных данных. В ячейки А2:D5 ввести исходные данные из таблицы задания.

3.1.2.В ячейках B6:D8 разместить формулы для вычисления технологических коэффициентов:

в ячейку В6 ввести формулу для вычисления первого коэффициента =B3/$A$3 и скопировать ее в ячейки В7:B8;

Таблица 8

 

A

B

C

 

D

1

 

БАЛАНСОВАЯ МОДЕЛЬ

 

Объём про-

 

 

 

 

 

2

изводства

 

Потребление отраслей

3

600

250

100

 

 

160

4

1000

150

500

 

 

0

5

800

0

300

 

 

400

6

Вычисление

=В3/А$3

=С3/А$4

 

=D3/А$5

 

технологиче-

 

 

 

 

 

7

= В4/А$3

=С4/А$4

 

=D4/А$5

ских коэф-

 

 

 

 

 

 

8

фициентов

= В5/А$3

=С5/А$4

 

=D5/А$5

9

 

Проверка продуктивности матрицы А

10

 

=СУММ(B6:B8)

=СУММ(C6:C8)

 

=СУММ(D6:D8)

 

 

 

=ЕСЛИ(A11=ИСТИНА;"Решения нет";"Матрица

11

=ИЛИ(B10>=1;C10>=1;D10>=1)

 

 

продуктивна")

12

Единичная

1

0

 

0

13

 

 

 

матрица

0

1

 

0

 

0

0

 

14

 

 

1

15

Вычисление

=B12-B6

=C12-C6

 

 

=D12-D6

16

=B13-B7

=C13-C7

 

 

=D13-D7

 

Е-А

 

 

 

 

 

17

=B14-B8

=C14-C8

 

 

=D14-D8

 

 

 

18

Вычисление

=МОБР(B15:D17)

=МОБР(B15:D17)

 

 

=МОБР(B15:D17)

19

обратной

=МОБР(B15:D17)

=МОБР(B15:D17)

 

 

=МОБР(B15:D17)

20

матрицы

=МОБР(B15:D17)

=МОБР(B15:D17)

 

 

=МОБР(B15:D17)

21

Спрос на

2000

План выпуска

 

=МУМНОЖ(B18:D20;B21:B23)

22

будущий

2000

 

=МУМНОЖ(B18:D20;B21:B23)

 

период

 

продукции

 

 

 

23

3000

 

=МУМНОЖ(B18:D20;B21:B23)

 

 

 

 

 

18

аналогично в ячейку С6 ввести формулу =C3/$A$4 и скопировать ее в ячейки С7:С8;

вячейкуD6 ввестиформулу=D3/$A$5 ископироватьеевячейкиD7:D8.

3.1.3.В ячейках В10:D10 разместить формулы для подсчета суммы значений элементов по столбцам:

в ячейку В10 ввести формулу =СУММ(В7:В9);

скопировать формулу в ячейки С10:D10.

3.1.4.В строке 11 размещаем формулы для проверки продуктивности матрицы технологических коэффициентов:

в ячейку А11 ввести формулу =ИЛИ(В10>=1;C10>=1;D10>=1).

Эта формула проверяет содержимое ячеек В10:D10. Если хотя бы в одной из этих ячеек значение больше единицы (т.е. сумма значений элементов хотя бы в одном столбце превышает единицу), то в ячейке А11 будет записано значение ИСТИНА». В противном случае – значение «ЛОЖЬ»;

в ячейку C13 ввести формулу

=ЕСЛИ(А11=”ИСТИНА”;”Нет решения”;”Матрица продуктивна”).

Эта формула проверяет содержимое ячейки А11 и если сумма элементов хотя бы одного столбца превысила единицу, выводит сообщение “Нет решения”, а в противном случае – “Матрица продуктивна”.

3.1.5.В строках 12 –14 разместить единичную матрицу Е.

3.1.6.В строках 15 – 17 произвести вычисление матрицы Е-А:

в ячейку В15 поместить формулу =В12-В6;

скопировать формулу в ячейки В16:D17.

3.1.7.В строках 18 – 20 разместим формулы для вычисления матрицы, обратной матрице Е-А:

активизировать ячейку В18;

зажав левую клавишу мыши, выделить диапазон ячеек В18:D20, где будет размещена обратная матрица;

щелкнуть по пиктограмме Мастер функций fx;

в первом окне Мастера функций в поле Категория выбрать

Математические;

в поле Функция среди расположенных по алфавиту функций найти функцию МОБР;

щелкнуть по кнопке Ок и перейти во второе окно Мастера функций;

в поле Массив ввести адрес матрицы Е-А: диапазон ячеек В15:D17;

чтобы формула была введена во все ячейки выделенного диапазона, следует одновременно нажать клавиши Ctrl+Shift+Enter; после этого в ячейки В18:D20 будет введена формула =МОБР(В15:D17).

3.1.8.В строках 21 – 23 поместим формулы для вычисления плана выпуска продукции:

в ячейках В21:В23 разместить значения спроса на будущий период согласно заданию (табл. 8);

выполнить команды Вставка Функция;

19

выделить ячейки D21:D23, в которых будет размещена формула перемножения элементов матрицы, обратной Е-А, и вектора-столбца спроса;

в категории Математические Мастера функций выбрать функцию МУМНОЖ;

во втором окне Мастера в поле Массив1 ввести адрес обратной матрицы: диапазон ячеек В18:D20;

в поле Массив2 ввести адрес вектора-столбца спроса: диапазон ячеек В21:D23;

чтобы формула была введена во все ячейки выделенного диапазона, следует одновременно нажать клавиши Ctrl+Shift+Enter. После этого в ячейки D21:D23 будет введена формула =МУМНОЖ(В18:D20;B21:B23).

 

 

 

 

 

 

 

Таблица 9

 

 

 

 

 

 

 

 

 

 

A

 

B

 

 

C

D

 

1

 

 

БАЛАНСОВАЯ МОДЕЛЬ

 

 

2

Объём

 

 

Потребление отраслей

 

 

производства

 

 

 

 

 

 

 

 

 

 

 

3

600

 

250

 

100

160

 

4

1000

 

150

 

500

0

 

5

800

 

0

 

300

400

 

6

Вычисление

 

0,417

 

0,1

0,2

 

7

технологических

 

0,25

 

0,5

0

 

8

коэффициентов

 

0

 

 

0,3

0,5

 

 

 

 

 

9

Проверка продуктивности

матрицы А

 

 

10

 

 

0,667

 

0,900

0,700

 

11

ЛОЖЬ

 

 

 

Матрица продуктивна

 

 

 

 

 

 

 

 

12

Единичная

 

1

 

0

0

 

13

 

0

 

1

0

 

матрица

 

 

 

14

 

0

 

0

1

 

 

 

 

 

 

 

 

 

 

-0,1

 

 

15

Вычисление Е-А

 

0,583

 

 

-0,2

 

16

 

-0,25

 

0,5

0

 

17

 

 

0

 

 

-0,3

0,5

 

 

 

 

 

 

 

 

18

Вычисление

 

2,113

 

0,930

0,845

 

19

обратной

 

1,056

 

2,465

0,423

 

20

матрицы

 

0,634

 

1,479

2,254

 

 

 

 

 

 

 

 

 

 

 

 

 

21

Спрос на

 

2000

 

 

План выпуска

8619,72

 

22

 

2000

 

 

8309,86

 

будущий период

 

 

 

продукции

 

23

 

3000

 

 

10985,92

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

3.2. Самостоятельная работа

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

Таблица 10

№ варианта

Объем производства

Потребление отраслей

Спрос

1

1000

200

250

120

1000

600

150

50

40

1000

 

800

0

250

700

1000

2

800

150

250

140

2000

1000

200

500

150

2000

 

600

100

0

250

2000

 

 

 

 

 

 

3

600

100

100

0

1500

800

150

200

250

1500

 

1000

200

200

500

1500

 

 

 

 

 

 

4

800

0

250

120

1500

1200

200

300

40

1500

 

1000

150

350

700

1500

 

 

 

 

 

 

5

1000

150

250

140

2000

1500

200

350

150

2000

 

800

100

100

250

2000

 

 

 

 

 

 

6

1500

250

140

150

2000

1000

200

250

150

2000

 

800

100

120

200

2000

 

 

 

 

 

 

7

1000

150

250

100

1500

900

200

150

150

1500

 

800

100

0

250

1500

 

 

 

 

 

 

8

800

50

250

120

1000

1100

200

400

40

1000

 

900

150

350

700

1000

 

 

 

 

 

 

9

1000

200

250

120

1000

1500

150

400

40

1000

 

1000

50

350

700

1000

 

 

 

 

 

 

10

1000

200

250

120

2000

800

250

150

40

2000

 

800

100

350

300

2000

4. Отчет по работе

Распечатки таблицы вычислений.

Литература: [2], c. 509-515.