
- •Математичні моделі в економіці Вступ
- •Модель міжнародної торгівлі або лінійна модель обміну.
- •Теорія невід’ємних матриць.
- •Нерозкладні матриці. Основні поняття.
- •2. Властивості нерозкладних матриць.
- •Теорема Фробеніуса-Перрона
- •Примітивні матриці. Стійкість
- •Аналіз продуктивності моделі «затрати-випуск»
- •Коефіцієнти трудових витрат.
- •Порівняльна статика моделі Леонтьєва
- •Еластичність
- •Модель динамічного міжгалузевого балансу.
- •Модель фон Неймана
- •Послідовність цін буде стаціонарною траєкторією цін тоді і тільки тоді , коли і початковий вектор справджують нерівність µ . (8)
- •Динамічна рівновага в моделі Неймана.
- •Виробничі функції, загальні положення
- •Виробнича функція Коба-Дугласа
- •Еластичність заміни факторів
- •Виробнича функція ces
2. Властивості нерозкладних матриць.
1. Якщо матриця А нерозкладна, що вона не може мати ні нульових рядків ні нульових стовпчиків.
▼ Справді,
якщо j
стовпчик нульовий, то S={j}
ізольована множина. Якщо нульовий і-ий
рядок, то S=
ізольована
підмножина.
▲
2.
Якщо матриця А
нерозкладна, а вектор
,
то вектор
▼Властивість очевидна, оскільки матриця не має нульових рядків.
▲
3.
Нехай А
- нерозкладна матриця, вектор y
0
, а
і нехай z=Ay,
T
.
Якщо при цьому
і
,
то
▼Доведення. Від супротивного.
Якщо
позначити
і припустити
,
то
,
а
.
Візьмемо
і
тому
і
.
Оскільки
,
то це можливо лише тоді, коли
.
Отже S
-
ізольована, що суперечить, що А
- нерозкладна.
▲
4.
Нехай у
- невід’ємний ненульовий вектор (
),
а вектор z=(I+A)y.
Позначимо через
число координат вектора у,
.
Тоді
при
і
при
.
Крім цього, якщо А
- нерозкладна,
,
то у нерівності
будемо мати
.
▼Доведення.
Перш за все, оскільки
,
то
.
Від
супротивного. Припустимо, що серед
координат вектора x
є
нульові, тобто множина
непорожня.
За властивістю 3) знайдеться
такий що,
, але тоді неможлива нерівність
.
▲
5.
Якщо А
- нерозкладна матриця розмірності n*n,
то
,
тобто всі елементи матриці
є строго додатніми.
▼Доведення.
Дійсно за властивістю 4)
маємо,
що
.
Досить в якості y
взяти
▲
6. Якщо А - нерозкладна, то для будь-якої пари індексів (i,j),
знайдеться
натуральне число m
таке, що
(через
позначено
відповідний елемент матриці
).
▼Справді, запишемо розклад за формулою бінома Ньютона.
(2)
Нехай
спочатку
.
Серед чисел
обов’язково є додатне : в протилежному
випадку, як впливає з формули (2), (i,j)-ий
елемент матриці
дорівнює нулю, а це суперечить властивості
5). Якщо ж i=j,
то слід аналогічно до попереднього
розглянути матрицю
,
яка за властивостями 2) і 6) є додатною
7.
Для того, щоб матриця А
була нерозкладною, необхідно і достатньо,
щоб для будь-яких індексів i,j
існувала послідовність
така
що
,
,
8.
Якщо матриця А
нерозкладна, а
число
з
,
то в матриці
не
можна бути ні нульових рядків, ні нульових
стовпчиків.
Теорема Фробеніуса-Перрона
Теорема
Фробеніуса-Перрона
Нерозкладна
матриця A
має додатнє власне значення
a
таке,
що модулі всіх решта власних чисел не
перевищують
a
:
k|
≤
,
k
= 1,…,m.
Числу a відповідає єдиний (з точністю до скалярного множника) власний вектор XA, всі координати якого ненульові і одинакові за знаком:
Sign(XA)I
= Sign(XA)j
,
i,j = 1
n.
Тобто вектор XA можна вибрати додатнім XA > 0.
Зазначимо,
що число
a
називається
числом Фробеніуса, а вектор XA
– вектором Фробеніуса матриці A
( Axa
=
axa;
a
–
розв’язок рівняння det(A-
I)
= 0 ).
Теорема
2 .
Якщо квадратна n*n
матриця A
невід’ємна, то вона має невід’ємне
власне число
a
≥
0 таке, що для будь-якого іншого її
власного числа
має місце нерівність |
| ≤
a.
При цьому існує невід’ємний власний
вектор XA
≥
0, який відповідає
.
Для
формулювання наступної теореми введемо
позначення:
Теорема 3 (Про оцінки для Фробеніусового числа невід’ємної матриці А).
Якщо квадратна n*n матриця A ≥ 0, то для її Фробеніусового числа мають місце такі нерівності:
≤
a
≤
R
s
≤
a
≤
S
(3)
Якщо матриця A ще й нерозкладна, то всі нерівності в (3) строгі. Вийняток лише у випадку, коли = R, s = S.
Приклад
Для матриці A
=
найти
число і вектор Фробеніуса і оцінити
число Фробеніуса згідно теореми 3.
Характеристичне
рівняння:
= 0
2 – 1,5 + 0,5 = 0
1
= 1
2
= 0,5
A
= 1
(A)
= {0,5; 1}
Вектор XA шукаємо з системи: Ax = Ax або (A- AI)x = 0.
Тоді, наприклад XA = (3,2).
r = 0,9; R = 1,1; s =1; S = 1.
0,9 ≤ A ≤ 1,1
A = 1.