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
Межотраслевой баланс можно представить в виде таблицы, где каждая отрасль представлена дважды.
Представленная строкой отрасль выступает как поставщик продукции всем отраслям, включая себя.
Представленная столбцом отрасль выступает как потребитель, включая себя саму.
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.