Коэффициенты прямых материальных затрат. Объемы конечной продукции. Изменение удельной условно-чистой продукции, страница 2

X2=0,12 x1 +0,05 x2 +0,04 x3 +0,10 x4 +200;

X3=0,15 x1 +0,10 x2 +0,07 x3 +0,05 x4 +300;

X4=0,12 x1 +0,10 x2 +0,05 x4 +100.

Поскольку экономико-математическая модель МОБ сводится к системе уравнений из n-линейных уравнений (в данном случае из четырех), то к ней применимы обычные методы решения этих систем, которые можно разделить на точные (прямые) и итерационные (приближенные).

Прямые методы позволяют найти точное решение за конечное число шагов.

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

В любом итерационном или приближенном методе решение начинается с изобретения начального или нулевого приближения к решению. В методе межотраслевого баланса начальное приближение – это объемы конечной продукции, т. е. нулевым приближением будет являться:

Xi(0)=Yi

Суть итерационного метода заключается в последовательном приближении к точному решению. По известным x(0) определяем первое приближение, т. е. улучшенное решение:

1.  Для метода простой итерации:

X1(1)=a11x1(0)+a12x2(0)+…+a1nxn(0)+Y1;

X2(1)=a21x1(0)+a22x2(0)+…+a2nxn(0)+Y2;

………………………………………

Xn(1)=an1x1(0)+an2x2(0)+…+annxn(0)+Yn.

По первому приближению – второе и т. д.

Алгоритм матричный: X(k)=AX(k-1)+Y.

Решение будет продолжаться до тех пор, пока не будет выполнено следующее условие: |Xi(k)-Xi(k-1)| ≤ ε, i=1,n.

ε – некое малое наперед заданное число.

Чем меньше величина  ε, тем большее число итераций требуется для выполнения последнего условия, тем точнее расчет.

Решение идет с точность до ε.

2.  Для метода Зейделя:

X1(1)=a11x1(0)+a12x2(0)+a13x3(0)+…+a1nxn(0)+Y1;

X2(1)=a21x1(1)+a22x2(0)+a23x3(0)+…+a2nxn(0)+Y2;

X3(1)=a31x1(1)+a32x2(1)+a33x3(0)+…+a3nxn(0)+Y3;

………………………………………………….

Xn(1)=an1x1(1)+an2x2(1)+an3x3(1)+…+annxn(0)+Yn.

Аналогично определяется 2, 3, ……, k – приближение. Решение продолжается до тех пор, пока не будет выполнено условие: |Xi(k)-Xi(k-1)| ≤ ε, i=1,n.

Систему уравнений (1) решим при помощи компьютера итерационными методами..

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

ε

0,01

0,001

0,0001

X1

460,3031

460,3043

460,3044

X2

307,8256

307,8267

307,8268

X3

440,4488

440,4503

440,4503

X4

195,8086

195,8097

195,8097

Число итераций

10

12

13

Из таблицы видно, что необходимые значения найдены за 25 итераций.

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

ε

0,01

0,001

0,0001

X1

460,303

460,3044

460,3044

X2

307,8262

307,8268

307,8268

X3

440,4498

440,4504

440,4504

X4

195,8094

195,8097

195,8097

Число итераций

7

9

10

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

Изобразим графически зависимость количества итераций от точности расчета.

Таким образом, можно сделать вывод, что необходимые значения объемов конечной продукции методом Зейделя были найдены за меньшее число итераций (9), чем методом простых итераций (12).

В нашей задаче это условие выполняется при ε = 0,001. Использование метода Зейделя дает следующие значения:

X1 =460,3044

X2 =307,8268

X3 =440,4504

X4 =195,8097

Задание 2

Межотраслевой баланс можно представить в виде таблицы, где каждая отрасль представлена дважды.

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

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