Эконометрия. Методическое пособие
.pdf
ЭКОНОМЕТРИЯ
Длина вектора
Длиной n–мерного вектора a=(a1,a2,...,an) называется число
n
a
(a, a) Rai2 .
i 1
Линейная комбинация векторов a1,a2,....,am одинаковой размерности n с числами c1,c2,...,cm есть вектор атой же размерности
m, вычисленный по формуле:
a=c1a1+c2a2+...+cmam
Линейная зависимость и линейная независимость системы векторов
Система векторов a1,a2,....,am называется линейно-зависимой,
если существует такой набор чисел c1,c2,...,cm, одновременно не
n
равный нулю Rci2=/ 0, с которыми линейная комбинация дает нуле-
i 1
вой вектор. В противном случае система векторов будет линейно-
независимой.
Определение векторного пространства
Векторным пространством En размерности n называется множество векторов одинаковой размерности n, обладающее свойствами:
1)любой вектор из этого множества, умноженный на любое число, принадлежит этому множеству;
2)сумма любых двух векторов из этого множества есть вектор, принадлежащий этому множеству.
Базис векторного пространства
Любая совокупность линейно-независимых векторов из данного
множества, число которых равно размерности n, называется базисом векторного пространства En.
Теорема. Любой вектор из векторного пространства можно представить в виде линейной комбинации векторов базиса, причем единственным образом.
11
О.В.БРАИЛОВА
Тема 3. Система линейных неравенств
Совместимость и несовместимость системы линейных неравенств. Геометрическая интерпретация решения системы линейных неравенств.
Совместимость и несовместимость системы линейных неравенств
Неравенство называется линейным, если оно содержит неиз-
вестные только в первой степени.
a1x1+a2x2+...+anxn>b.
Решением системы линейных неравенств называется такая совокупность чисел c1,c2,...,cn, которая при подстановке их вместо неизвестных x1,x2,....,xn, удовлетворяет каждое неравенство.
Система линейных неравенств называется совместной, если она имеет хотя бы одно решение, иначе система будет несовместной.
Геометрическая интерпретация решения системы линейных неравенств
Если система содержит только два неизвестных x1,x2, то ее решение можно интерпретировать геометрически на плоскости.
Пример 6
Решить систему линейных неравенств геометрически.
3x1 x2 |
6 |
|
|
x2 |
|
|
123–6x1 |
x2 |
6. |
|
|
|
|
–2x1 |
x2 |
0 |
|
|
|
|
|
|
|
|
|
|
3x1+x2=6 |
Полученный |
треугольник |
АВС |
C |
-2x1+x2=0 |
||
является решением системы. |
|
|||||
|
|
|
||||
|
|
|
|
|
B |
-6x1+x2=6 |
|
|
|
|
|
A |
x1 |
|
|
|
|
|
|
|
12
ЭКОНОМЕТРИЯ
Тема 4. Система линейных уравнений
Элементарные преобразования системы m линейных уравнений с n неизвестными. Запись системы в матричном виде. Решение системы линейных уравнений методом исключения неизвестных (метод Жордано-Гаусса, метод Гаусса). Теорема о совместимости системы линейных уравнений (КронекераКапелли).
Элементарные преобразования системы m линейных уравнений с n неизвестными
Пусть дана система m линейных уравнений с n неизвестными:
a11x1 a12 x2 ... a1n xn |
b1 |
||
123 |
a22 x2 |
... a2n xn |
b2 |
a21x1 |
|||
... |
|
am1x1 am2 x2 ... amn xn |
bm |
|
.(1) |
Элементарными преобразованиями системы линейных уравнений называются:
1)умножение обеих частей уравнения на любое не равное нулю число;
2)прибавление к обеим частям уравнения линейной комбинации других уравнений;
3)перестановка местами уравнений.
Для решения системы линейных уравнений используется метод исключения неизвестных, который имеет две разновидности.
Метод Жордано-Гаусса (метод полного исключения)
В качестве первого ведущего неизвестного выберем x1, а в качестве первого ведущего уравнения – любое уравнение, содержащее x1. С помощью элементарных преобразований исключим x1 из всех остальных уравнений.
В качестве второго ведущего неизвестного возьмем x2, а вторым ведущим уравнением будем считать любое уравнение, содержащее x2, кроме первого ведущего уравнения, и т.д.
13
О.В.БРАИЛОВА
Врезультате получим решение системы линейных уравнений (1)
ввиде:
если ранг матрицы коэффициентов системы обозначить через r, то А) при r=n:
123 |
p1 |
|
|
|
|
|
x1 |
|
|
|
|
|
|
x2 p2 |
|
|
|
|
|
|
|
... |
|
|
|
|
|
xn |
pn |
. |
|
|
|
|
|
|
|
|
|
||
В) при r<n: |
r |
|
||||
123r |
r ,r 1 |
r 1 |
rn n |
|
||
x1 |
q1,r 1 xr 1 |
... q1n xn |
p1 |
|
||
x |
q2,r 1 xr 1 ... q2n xn |
p2 |
|
|||
|
|
|
|
... |
|
|
x |
q |
|
x |
... q x |
p |
. |
|
|
|
|
|
|
|
В случае r=n система имеет единственное решение. |
||||||
Если же r<n, получаем бесконечное множество решений:
|
x |
= p |
− q |
|
+ ...+ q |
x |
|
|
1 |
1 |
( 1,r +1 |
|
1,n n ) |
||||
|
x |
= p |
− q |
|
+ ...+ q |
x |
|
|
|
2 |
2 |
( |
2,r +1 |
|
2,n n ) |
||
... |
|
|
|
|
|
|
|
|
x |
= p |
− q |
|
+ ...+ q |
x |
|
||
|
r |
r |
( |
r,r +1 |
r,n n ) .(*) |
|||
Здесь x1,x2,...xr называются базисными переменными, а их совокупность – базисом. Остальные переменные xr+1,xr+2,...,xn свободные.
Совокупность (*) называется общим решением системы. Базисным решением называется такое частное решение систе-
мы линейных уравнений, в которых все свободные переменные равны нулю.
x1 p1 |
r 2 |
n |
|
. |
||
123r |
1 |
|
|
|||
x2 |
p2 |
|
|
|
||
... |
|
|
|
|
|
|
xr pr |
|
|
|
|
||
x |
|
x |
... x |
0 |
|
|
14
ЭКОНОМЕТРИЯ
Пример 7
Решить систему линейных уравнений методом Жордано-Гаусса.
123 |
1 |
2 |
|
3 |
2x1 |
– x2 – x3 |
|
4 |
|
x1 4x2 – 2x3 11
3x – 2x 4x 11
.
С помощью элементарных преобразований приведем матрицу коэффициентов исходной системы к единичной матрице:
2 |
–1 –1 |
|
4 |
1 |
5 |
– 1 |
|
7 |
1 |
0 |
4 |
|
7 |
|
|
|
|
|
|||||||||||
3 |
4 –2 |
|
11 ~ 0 |
–11 1 |
|
–10 ~ 0 |
1 |
–1 |
|
0 |
~ |
|||
3 |
–2 4 |
|
11 |
0 |
1 |
–1 |
|
0 |
0 |
0 |
–10 |
|
– 10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
0 |
0 |
3 |
~ 0 |
1 |
0 |
1 |
0 |
0 |
1 |
1 . |
|
|
|
Таким образом, получим решение:
x1 3 |
|
||
123 |
3 |
|
|
x2 |
1 |
|
|
x |
1 |
. |
|
|
|
|
|
Ответ: x=(3,1,1).
Метод Гаусса (метод неполного исключения)
Этот метод разрешает исключать ведущие неизвестные только из тех уравнений, которые не были ведущими на предыдущих шагах. Тогда система (1) преобразуется к виду:
А) при r=n система имеет единственное решение, легко получаемое из системы треугольного вида:
x1 a12 x2 a13 x3 ... a |
1n xn b 1 |
||
123 |
x2 |
a23 x3 ... a |
xn bn |
|
2n xn b2 |
||
...
15
О.В.БРАИЛОВА
В) при r<n система имеет бесконечное множество решений, в
котором x1,x2,...xr – базисные переменные, а xr+1,xr+2,...,xn – свободные.
Пример 8
Решить систему методом Гаусса:
3x1 + 2x2 + x3 = 52x1 + 3x2 + x3 = 12x1 + x2 + 3x3 = 11.
Приведем матрицу коэффициентов системы к треугольному виду с помощью элементарных преобразований:
3 2 1 |
|
5 1 −1 0 |
|
4 |
1 −1 0 |
|
4 |
||||
|
|
|
|||||||||
2 3 1 |
|
1 |
0 |
5 1 |
|
−7 |
0 |
1 −1 |
|
−5 |
|
2 1 3 |
|
11 |
0 |
1 −1 |
|
−5 |
0 |
0 1 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
. |
Получили систему:
x1 − x2 = 4x2 − x3 = −5x3 = 3
или
x1x2x3
=2
=−2
=3
Ответ: x=(2,-2,3).
Теорема Кронекера-Капелли
Система линейных уравнений совместна тогда и только тогда, когда ранг основной матрицы коэффициентов системы равен рангу расширенной матрицы.
16
ЭКОНОМЕТРИЯ
ГЛАВА 2. ЛИНЕЙНОЕ ПРОГРАММИРОВАНИЕ
Тема 5. Примеры задач линейного программирования
Задача о пищевом рационе. Задача о распределении ресурсов. Задача о пищевом рационе
Имеются четыре продукта питания П1,П2,П3,П4 по цене соответственно C1,C2,C3,C4. Из этих продуктов требуется составить пищевой рацион, содержащий: белков, не менее чем в1 единиц, жиров, не менее чем в2 единиц и углеводов, не менее чем в3 единиц. Содержание элемента в единице каждого продукта указано в таблице:
элемент |
П1 |
П2 |
П3 |
П4 |
белки |
а11 |
а12 |
а13 |
а14 |
жиры |
а21 |
а22 |
а23 |
а24 |
углеводы |
а31 |
а32 |
а33 |
а34 |
Составить пищевой рацион из этих продуктов минимальной стоимости.
Обозначим количество продуктов, входящих в рацион соот-
ветственно x1,x2,x3,x4. Тогда общая стоимость рациона равна
L=c1x1+c2x2+c3x2+c4x4.
Общее количество белков, содержащихся в рационе, составляет
a11x1+a12x2+a13x3+a14x4. Их должно быть не менее в1 единиц. Общее количество жиров рациона есть a21x1+a22x2+a23x3+a24x4,
что не превосходит в2 единиц.
Общее количество углеводов в рационе есть
a31x1+a32x2+a33x3+a34x4, что не должно превосходить в3 единиц. Эти условия представляют собой неравенства-ограничения, кото-
рые накладываются на неизвестные. Возникает следующая задача: найти такой неотрицательный вектор x=(x1,x2,x3,x4)>0, который
удовлетворяет системе неравенств-ограничений
a11x1 + a12 x2 + a13 x3 + a14 x4 ≥в1a21x1 + a22 x2 + a23 x3 + a24 x4 ≥в2a31x1 + a32 x2 + a33 x3 + a34 x4 ≥в3
17
О.В.БРАИЛОВА
и обращает линейную функцию в минимум:
L=c1x1+c2x2+c3x3+c4x4—>min.
Задача о распределении ресурсов
Имеются ресурсы (например сырье, рабочая сила, оборудование) R1,R2,R3 в количестве b1,b2,b3 единиц. С помощью этих ресурсов
производятся три вида товаров T1,T2,T3. Для производства товара Tj необходимо aij единиц ресурса Ri(i=1,2,3, j=1,2,3). Единица ресурса Ri стоит di рублей (i=1,2,3). Единица товара Tj может быть реализована по цене cj рублей (j=1,2,3). По каждому виду това-
ра количество производимых единиц ограничивается спросом. Из-
вестно, что рынок не может реализовать более чем Kj единиц товара Tj(j=1,2,3). Составить план производства товаров, приносящий
максимальную прибыль.
Обозначим через x1,x2,x3 количество товаров T1,T2,T3, которые надо запланировать к производству. Условия спроса налагают ограничения:
1230 |
|
x3 |
|
k3 . |
|
0 |
|
x1 |
|
k1 |
|
0 |
|
x2 |
|
|
k2 |
|
|
|
|||
|
|
|
|
|
|
Ресурсов должно хватить, следовательно, возникают ограничения:
a11x1 a12 x2 a13 x3 |
|
|
b1 |
|||
|
||||||
123a31x1 |
a32 x2 |
a33 x3 |
|
|
b3 . |
|
a21x1 |
a23 x2 |
a23 x3 |
|
|
|
b2 |
Составим прибыль L как функцию от переменных x1,x2,x3.
Себестоимость Sj единицы товара Tj равна
Sj=a1jd1+a2jd2+a3jd3(j=1,2,3).
Прибыль от единицы товара Tj есть
qj=cj-Sj(j=1,2,3).
18
ЭКОНОМЕТРИЯ
Общая прибыль производства составляет
L=q1x1+q2x2+q3x3—>max.
Итак, сформулируем задачу линейного программирования. Найти такой неотрицательный вектор x=(x1,x2,x3)>0, который бы удовлетворял системам неравенств-ограничений
123 |
31 |
1 |
32 |
2 |
33 |
3 |
3 . |
|
a11x1 |
a12 x2 |
a13 x3 |
|
b1 |
||||
a21x1 a23 x2 a23 x3 |
|
b2 |
||||||
|
||||||||
a |
x a |
x |
a |
x |
|
b |
||
k1
k2
k3 .
и обращал бы линейную функцию в максимум:
L=q1x1+q2x2+q3x3—>max.
Тема 6. Основная задача линейного программирования (ОЗЛП)
Постановка ОЗЛП. Допустимое и оптимальное решения ОЗЛП. Геометрическая интерпретация ОЗЛП. Опорное решение. Симплекс-метод решения ОЗЛП.
Постановка ОЗЛП
Рассмотрим задачу линейного программирования с огра- ничениями-равенствами.
Найти такой неотрицательный вектор x=(x1,x2,...,xn)>0, который удовлетворяет системе уравнений-ограничений:
a11x1 a12 x2 |
... a1n xn |
b1 |
123 |
|
|
a21 x1 a22 x2 |
... a2n xn |
b2 (1) |
|
... |
|
am1 x1 am2 x2 |
... amn xn bm |
|
и обращаeт линейную функцию в минимум
L=c1x1+c2x2+...+cnxn—>min.
19
О.В.БРАИЛОВА
Допустимое и оптимальное решения ОЗЛП
Допустимым решением задачи называют любой неотрицательный вектор x=(x1,x2,...,xn)>0, удовлетворяющий системе (1).
Оптимальным решением называют такое допустимое решение, при котором линейная функция L=c1x1+c2x2+...+cnxn обращается в минимум.
Если число неизвестных больше числа линейно-независимых уравнений-ограничений на два, то решение задачи можно рассмотреть на плоскости.
И имеем две свободные переменные, скажем, x1,x2. Выразим остальные базисные переменные x3,x4,...,xn через свободные:
x3 b3 a 31x1 a 32 x2 |
x |
|
|
|
|
|
|
|
|
|||
x4 b4 a 41x1 a 42 x2 |
2 |
|
|
|
|
x |
=0 |
|||||
|
|
... |
|
|
|
|
|
|
4 |
|||
xn bn a n1x1 a n2 x2 |
|
|
|
=0 |
|
|
|
|
|
|||
123 |
x1 |
0, x2 |
0 |
|
|
x |
ОДР |
|
|
|||
|
|
|
|
|
|
3 |
|
|
||||
Получили систему |
|
|
|
|
|
|
||||||
линейных неравенств |
|
|
|
|
|
|
|
|
|
|||
относительно |
двух |
|
|
|
|
|
|
|
|
|
||
переменных x1,x2. Ее |
|
|
|
|
x |
|
|
x |
1 |
|||
можно интерпретиро- |
|
|
|
|
n |
|
|
|||||
|
|
|
|
=0 |
|
|
||||||
вать геометрически |
|
|
|
|
|
|
|
|
|
|||
на плоскости x1,x2: |
|
|
|
|
|
|
|
|
|
|||
Это область допустимых решений. |
|
|
|
|
|
|||||||
Найдем |
среди допустимых |
решений оптимальное. Для это- |
||||||||||
го выразим линейную |
функцию |
через |
свободные |
переменные |
||||||||
L=v0+v1x1+v2x2.
Очевидно, что линейная функция L`=v1x1+v2x2 достигает своего минимального значения при тех же x1,x2, что и L=v0+v1x1+v2x2 .
Построим прямую L`=0 и будем перемещать ее параллельно самой себе в направлении уменьшения. Итак, прямая v x +v x =0
имеет угловой коэффициент -v_1. 1 1 2 2
v2
20
