Лекция 10 Системы линейных уравнений

 

Показатели функционирования систем и процессов, экономические показатели часто оказываются взаимозависимыми. Структура связей между такими показателями (переменными) может быть описана с помощью системы (структурных) уравнений [3]. В этих уравнениях присутствуют переменных следующих типов:

- эндогенные, зависимые переменные y, определяемые внутри системы;

- экзогенные, независимые переменные x, значения которых задаются извне, они являются управляемыми, планируемыми;

- предопределенные переменные, включающие в себя как экзогенные переменные за текущий период времени, так и лаговые переменные (т.е. экзогенные и эндогенные переменные за предыдущие периоды времени).

Выделяют следующие виды систем.

1) Системы независимых уравнений, в которых каждая зависимая переменная y_i (i =1,…, n) представлена как функция одного и того же набора независимых переменных x_j (j =1,…, m):

y ₁ = a ₁₁ x ₁ + a ₁₂ x ₂ + …+a _{1 m} x_m + e ₁

y ₂ = a ₂₁ x ₁ + a ₂₂ x ₂ + …+a _{2 m} x_m + e ₂(1)

...............................

y_n = a_{n 1} x ₁ + a_{n 2} x ₂ + …+a_nm x_m + e_n

 

Каждое уравнение этой системы можно рассматривать самостоятельно как уравнение регрессии. В него может быть введен свободный член, коэффициенты регрессии могут быть найдены методом наименьших квадратов (м.н.к).

2) Системы рекурсивных уравнений, в которых зависимые переменные y_i (i= 1 ,…,n) представлены как функции независимых переменных x_j (j= 1 ,…,m) и определенных ранее зависимых переменных y ₁, y ₂ ,…, y_{i- 1}:

 

y ₁ = a ₁₁ x ₁ + a ₁₂ x ₂ + …+a _{1 m} x_m + e ₁

y ₂ = b ₂₁ y ₁ + a ₂₁ x ₁ + a ₂₂ x ₂ + …+a_2m x_m + e ₂(2)

...............................

y_n = b_{n 1} y ₁ + b_{n 2} y ₂ +,…,+b_{nn- 1} y_{n- 1} +a_{n 1} x ₁ + a_{n 2} x ₂ + …+a_nm x_m + e_n

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

3) Системы взаимозависимых уравнений, в которых каждая зависимая переменная y_i (i= 2 ,…,n) представлена как функция остальных зависимых переменных y_k (k ¹ i) и независимых (предопределенных) переменных x_j (j= 1 ,…,m):

y ₁ = b ₁₂ y ₂ + b ₁₃ y ₃ + … + b _{1 n} y_n +a ₁₁ x ₁ + a ₁₂ x ₂ + …+a_1m x_m + e ₁

y ₂ = b ₂₁ y ₁ + b ₂₃ y ₃ + … + b _{2 n} y_n +a ₂₁ x ₁ + a ₂₂ x ₂ + …+a _{2 m} x_m + e ₂(3)

...............................

y_n= b_{n 1} y ₁ + b_{n 2} y ₂ + … + b_{nn- 1} y_{n- 1} +a_{n 1} x ₁ + a_{n 2} x ₂ + …+a_nm x_m + e_n

 

Эта система наиболее распространенная, она получила также название системы совместных, одновременных уравнений. Ее так же называют структурной формой модели (СФМ).

Отдельные коэффициенты при переменных СФМ могут быть равны нулю, что означает отсутствие в уравнении этих переменных. Например, модель динамики цены и заработной платы может быть описана СФМ следующего вида:

 

y ₁ = b ₁₂ y ₂ +a ₁₁ x ₁ + e ₁

y ₂ = b ₂₁ y ₁ + a ₂₂ x ₂ +a ₂₃ x ₃ + e ₂(4)

где y ₁– темп изменения заработной платы;

y ₂– темп изменения цен;

x ₁– процент безработных;

x ₂– темп изменения постоянного капитала;

x ₃– темп изменения цен на импорт сырья.

 

Данная система из двух уравнений содержит две зависимые, эндогенные (y ₁, y ₂) и три независимые, экзогенные (x ₁, x ₂, x ₃) переменные. В первом уравнении отсутствуют переменные x ₂ и x ₃. Это значит, что коэффициенты a ₁₂= 0и a ₁₃= 0.

В СФМ для нахождения параметров модели b_ij и a_ij (называемых также структурными коэффициентами модели), простой м.н.к. неприменим.

Обычно для определения структурных коэффициентов модели СФМ преобразуется в приведенную форму модели (ПФМ):

y ₁ = d ₁₁ x ₁ + d ₁₂ x ₂ + …+d _{1 m} x_m

y ₂ = d ₂₁ x ₁ + d ₂₂ x ₂ + …+d _{2 m} x_m (5)

...............................

y_n = d_{n 1} x ₁ + d_{n 2} x ₂ + …+d_nm x_m

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

Структурные формы модели могут быть:

- идентифицируемые;

- неидентифицируемые;

- сверхидентифицируемые.

Для того, чтобы СФМ была идентифицируема, необходимо чтобы каждое уравнение системы было идентифицируемо. В этом случае число параметров СФМ равно числу параметров приведенной формы.

Если хотя бы одно уравнение СФМ неидентифицируемо, то вся модель считается неидентифицируемой. В этом случае число коэффициентов приведенной формы модели меньше, чем число коэффициентов СФМ.

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

Если обозначить число эндогенных переменных в i -том уравнении СФМ через Н, а число предопределенных переменных, которые содержатся в системе, но не входят в данное уравнение, через D, то условие идентифицируемости модели может быть записано в виде следующего счетного правила:

 

если D+1 < H уравнение неидентифицируемо;

если D+1 = H уравнение идентифицируемо;

если D+1 > H уравнение сверхидентифицируемо;

 

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

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

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

Поясним это на примере следующей структурной модели.

y ₁ = b ₁₂ y ₂ + b ₁₃ y ₃ + a_{1 1} x ₁ + a ₁₂ x ₂

y ₂ = b ₂₁ y ₁ + a ₂₂ x ₂ + a ₂₃ x ₃ + a ₂₄ x ₄ (6)

y ₃ = b ₃₁ y ₁ + b ₃₂ y ₂ +a ₃₁ x ₁ + a ₃₂ x ₂

 

Проверим каждое уравнение системы на выполнение неоходимого и достаточного условия идентификации.

В первом уравнении три эндогенных переменных: y ₁, y ₂и y ₃ (H =3). В нем отсутствуют экзогенные переменные x ₃и x ₄ (D =2). Необходимое условие идентификации D +1= H выполнено.

Для проверки на достаточное условие составим матрицу из коэффициентов при переменных x ₃и x ₄ (см. таблицу 1). В первом столбце таблицы показано, что коэффициенты при экзогенных переменных x ₃и x ₄взяты из уравнений 2 и 3 системы. Во втором уравнении эти переменные присутствуют и коэффициенты при них равны a ₂₃и a ₂₄ соответственно. В третьем уравнении эти переменные отсутствуют, т.е. коэффициенты при них равны нулю. Так как вторая строка матрицы состоит из нулей, определитель матрицы равен нулю. Значит, достаточное условие не выполнено, и первое уравнение нельзя считать идентифицируемым.

Во втором уравнении две эндогенные переменные: y ₁ и y ₂(H =2). В нем отсутствует экзогенная переменная x ₁(D =1). Необходимое условие идентификации D +1= H выполнено

 

 

Таблица 1 - Матрица, составленная из коэффициентов при x ₃и x ₄.

 

Уравнения, из которых взяты коэффициенты при переменных	Переменные
x 3	x 4
	a 23	a 24

 

 

Для проверки на достаточное условие составим матрицу из коэффициентов при переменных y ₃и x ₁, которые отсутствуют во втором уравнении (см. таблицу 2).

 

Таблица 2 - Матрица, составленная из коэффициентов при переменных y ₃и x ₁.

 

Уравнения, из которых взяты коэффициенты при переменных	Переменные
y 3	x 1
	b 13	a 11
	- 1	a 31

 

 

В третьем уравнении при переменной y ₃коэффициент равен –1, так как эта переменная стоит в левой части уравнения. Действительно, третье уравнение можно записать в виде 0 = b ₃₁ y ₁ + b ₃₂ y ₂ - 1 y ₃ +a ₃₁ x ₁ + a ₃₂ x ₂ и тогда равенство b ₃₃ = –1 становится очевидным.

В общем случае СФМ может быть представлена в виде матрицы коэффициентов при переменных. В этом случае третье уравнение может быть задано вектором (b ₃₁, b ₃₂, -1, a ₃₁, a ₃₂, 0, 0), а вся система одновременных уравнений (6) будет представлена матрицей (7):

 

(7)

 

 

В примерах и задачах СФМ представляют в виде такой матрицы коэффициентов при переменных модели.

Определитель представленной в таблице 2 матрицы не равен нулю, а ранг матрицы равен 2. Значит, достаточное условие выполнено, и второе уравнение идентифицируемо.

В третьем уравнении три эндогенные переменные: y ₁, y ₂ и y ₃ (H =3). В нем отсутствует экзогенные переменные x ₃и x ₄ (D =2). Необходимое условие идентификации D +1= H выполнено.

Для проверки на достаточное условие составим матрицу из коэффициентов при переменных х ₃и x ₄, которые отсутствуют в третьем уравнении (см. таблицу 3). Согласно таблице, определитель матрицы равен нулю (первая строка состоит из нулей). Значит, достаточное условие не выполнено, и третье уравнение нельзя считать идентифицируемым.

 

Таблица 3 - Матрица, составленная из коэффициентов при переменных x ₃и x ₄.

 

Уравнения, из которых взяты коэффициенты при переменных	Переменные
x3	x4
	0	0
	a23	a24

 

 

В эконометрических моделях иногда используются балансовые тождества переменных (например, вида y ₃ = y ₁ + y ₂ + x ₁). Коэффициенты при переменных при этом не требуют оценок и уравнение не надо исследовать на идентификацию, но в проверке на идентификацию всей системы эти уравнения участвуют. Присутствующие иногда в моделях свободные и остаточные члены (а ₀₁, а ₀₂, а ₀₃ ,…e ₁, e ₂, e ₃ ,…) не влияют на решение вопроса об идентификации.

При оценивании коэффициентов структурной модели используется ряд методов. С этими методами можно ознакомиться в рекомендованной литературе [3].

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

 

y ₁= b ₁₂ y ₂ + a ₁₁ x ₁ + e ₁ (8)

y ₂= b ₂₁ y ₁ + a ₂₂ x ₂ + e ₂

 

Для построения модели мы располагаем информацией, представленной в таблице 4.

Таблица 4 - Фактические данные для построения модели

 

n	у 1	у 2	х 1	х 2
	33,0	37,1
	45,9	49,3
	42,2	41,6
	51,4	45,9
	49,0	37,4
	49,3	52,3
Сумма	270,8	263,6
Средн. знач.	45,133	43,930	7,500	10,333

 

Структурную модель преобразуем в приведенную форму модели:

y ₁= d ₁₁ x ₁ + d ₁₂ x₂ + u _1, (9)

y ₂= d ₂₁ x ₁ + d ₂₂ x ₂ + u _2,

где u ₁ и u ₂ – случайные ошибки.

Для каждого уравнения приведенной формы при расчете коэффициентов d можно применить м.н.к..

Для упрощения расчетов можно работать с отклонениями от средних уровней y =y - y _cp и x =x - x _cp (y _cp и x _cp – средние значения). Преобразованные таким образом данные таблицы 4 сведены в таблицу 5. Здесь же показаны промежуточные расчеты, необходимые для определения коэффициентов d_ik. Переменные, означающие отклонение от средних значений изображаются далее жирным шрифтом и курсивом.

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

 

Σ y₁ x₁ = d ₁₁ Σ x₁² + d ₁₂ Σ x₁ x₂ (10)

Σ y₁ x₂ = d ₁₁ Σ x₁ x₂ + d ₁₂ Σ x₂²

 

Таблица 5 - Преобразованные данные для построения приведенной формы модели

n	у 1	у 2	х 1	х 2	у 1 *х 1	х 1 2	х 1 *х 2	у 1 *х 2	у 2 *х 1	у 2 *х 2	х 22
	-12,133	-6,784	-4,500	0,667	54,599	20,250	-3,002	-8,093	30,528	-4,525	0,445
	0,767	5,329	-0,500	5,667	-0,383	0,250	-2,834	4,347	-2,664	30,198	32,115
	-2,933	-2,308	-0,500	-1,333	1,467	0,250	0,667	3,910	1,154	3,077	1,777
	6,267	1,969	2,500	-1,333	15,668	6,250	-3,333	-8,354	4,922	-2,625	1,777
	3,867	-6,541	2,500	-9,333	9,667	6,250	-23,333	-36,091	-16,353	61,048	87,105
	4,167	8,337	0,500	5,667	2,084	0,250	2,834	23,614	4,168	47,244	32,115
Сумма	0,002	0,001	0,000	0,002	83,102	33,500	-29,001	-20,667	21,755	134,417	155,334

Подставляя рассчитанные в таблице 5 значения сумм в (10), получим:

 

83,102= 33,5 d ₁₁ - 29,001 d ₁₂ (11)

-20,667= -29,001 d ₁₁ + 155,334 d ₁₂

 

Решение этих уравнений дает значения d ₁₁ = 2,822 и d ₁₂ = 0,394. Первое уравнение приведенной формы модели примет вид:

 

y ₁= 2,822 x ₁ + 0,394 x ₂ + u₁ (12)

 

Для нахождения коэффициентов d _2k второго приведенного уравнения можно использовать следующую систему нормальных уравнений:

 

Σ y ₂ x ₁ = d₂₁ Σ x ₁ ² + d₂₂ Σ x ₁ x ₂ (13)

Σ y₂ x₂ = d₂₁ Σ x ₁ x ₂ + d₂₂ Σ x ₂²

Подставляя рассчитанные в таблице 5 значения сумм в (13), получим

21,755 = 33,5 d ₂₁ - 29,001 d ₂₂

134,417= -29,001 d ₂₁ + 155,334 d ₂₂. (14)

 

Решение этих уравнений дает значения d ₂₁ =1,668 и d ₂₂ =1,177. Второе уравнение приведенной формы модели примет вид:

 

y ₂= 1,668 x ₁ + 1,177 x ₂ + u_2. (15)

 

Для перехода от приведенной формы к структурной форме модели найдем x ₂ из второго уравнения (15) приведенной формы модели:

 

x ₂ = (y ₂ - 1,668 x ₁) / 1,177.

 

Подставим это выражение в первое уравнение (12) приведенной модели найдем структурное уравнение

y ₁= 2,822 x ₁ + 0,394 (y ₂ - 1,668 x ₁ ) / 1,177 =

= 2,822 x ₁ + 0,335 y ₂ - 0,558 x ₁ = 0,335 y ₂ + 2,264 x ₁

Таким образом, b ₁₂ = 0,335; a ₁₁ = 2,264.

Найдем x ₁ из первого уравнения (12) приведенной формы модели:

 

x ₁ = (y ₁ - 0,394 x ₂ ) / 2,822.

 

Подставим это выражение во второе уравнение (15), приведенной модели, найдем структурное уравнение:

 

y ₂= 1,177 x ₂ + 1,668 (y ₁ - 0,394 x ₂ ) / 2,822 =

= 1,177 x ₂ + 0,591 y ₁ - 0,233 x ₂ = 0,591 y ₁ + 0,944 x ₂

Таким образом, b ₂₁ = 0,591; a ₂₂ = 0,944.

 

Свободные члены структурной формы находим из уравнений:

 

а ₀₁= y _1,cp - b ₁₂ y _2,cp - a ₁₁ x _1,cp =

=45,133 – 0,335 * 43,93 –2,264* 7,5 = 13,436;

а ₀₂= y _2,cp - b ₂₁ y _1,cp - a ₂₂ x _2,cp =

=43,93 – 0,591* 45,133 - 0,944 * 10,333= 7,502.

 

Окончательный вид структурной модели:

 

y ₁= a ₀₁+ b ₁₂ y₂ + a ₁₁ x ₁ + e ₁= 13,436 + 0,335 y ₂ + 2,264 x ₁ + e _1; (16)

y ₂= a ₀₂+ b ₂₁ y ₁ + a ₂₂ x ₂ + e ₂= 7,502 + 0,591 y ₁ + 0,944 x ₂ + e _2.

 

 

 

 

Контрольные вопросы

1. Что представляет собой система линейных уравнений?

2. Особенности системы независимых уравнений. Методы решения.

3. Особенности системы рекурсивных уравнений. Методы решения.

4. Особенности системы взаимозависимых уравнений. Методы решения.

5. Структурные формы модели системы взаимозависимых уравнений.

6. Косвенный метод наименьших квадратов.