Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Эконометрия. Методическое пособие

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
711 Кб
Скачать

ЭКОНОМЕТРИЯ

Длина вектора

Длиной 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

 

1236x1

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. Примеры задач линейного программирования

Задача о пищевом рационе. Задача о распределении ресурсов. Задача о пищевом рационе

Имеются четыре продукта питания П1234 по цене соответственно 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

1230 x1
0 x2
0 x3

ЭКОНОМЕТРИЯ

Общая прибыль производства составляет

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]