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

высшая математика

.pdf
Скачиваний:
73
Добавлен:
18.02.2016
Размер:
5.88 Mб
Скачать

Пример 2. Пояснить теорему Перрона для случая, когда положи-

тельная матрица

 

A =

1

2

 

 

 

 

.

 

 

 

 

 

 

3

 

Решение.

 

 

4

 

 

Имеем:

 

 

 

1− λ

2

 

= 0 λ2

− 4λ − 5 = 0 λ1 = − 1, λ 2 = 5 λ(A) = 5 .

 

 

 

4

3 − λ

 

Соответствующий λ2 = λ(A) = 5 собственный вектор (см. Л. 15)

(

 

 

 

a12

 

 

 

2

 

 

x A =

 

λ

 

a

 

=

 

4

 

=

c

2

 

c

 

 

 

 

11

 

 

 

 

 

 

 

2

при c = 1.

 

 

 

 

4

 

 

 

 

 

Положительная матрица является частным случаем неразложимой неотрицательной матрицы. Фробениус обобщил теорему Перрона и на эти матрицы.

Теорема 2 (Фробениуса).

Неразложимая неотрицательная матрица А всегда имеет положительное характеристическое число λ(A), являющееся простым корнем характеристического уравнения, причем λ(A) ≥ λi , где {λ i } – совокупность остальных характеристических чисел. Числу λ(A) соответствует собственный вектор x(A с положительными координатами.

Из теоремы 1 вытекает следующий важный результат, который называют предельной теоремой.

Теорема 3.

Пусть c – произвольный неотрицательный вектор (c ≠ 0). Тогда предел

=A nc

vlim

[λ(A)]nn

существует и является собственным вектором матрицы А, соответствующим характеристическому числу λ(A). Этот предел единствен с точностью до скалярного множителя, зависящего от выбора начального вектора с.

30. Простая модель межотраслевого производствен-

ного процесса. Рассмотрим простую модель трех взаимодействую-

3 9 1

щих отраслей, которые будем условно называть «автомобилестроительной», «сталелитейной», «станкостроительной».

Пусть каждая из этих отраслей в любой момент времени может быть описана ее сырьевым ресурсом и ее производительной мощностью.

Введем следующие переменные состояния:

x1(n) – число автомобилей, произведенных к моменту времени n; x2(n) – производственная мощность автомобилестроительных за-

водов в момент времени n;

x3(n) – запас стали в момент времени n;

x4(n) – производственная мощность сталелитейных заводов в момент времени n;

x5(n) – количество (ресурс) станков в момент времени n;

x6(n) – производственная мощность станкостроительных заводов в момент времени n.

Если сделать некоторые допущения относительно экономической взаимосвязи трех отраслей и использовать принцип сохранения: количество материала в момент n + 1 равно количеству материала в момент n минус использованное за период времени (n; n + 1) плюс произведенное за указанный период времени, то получим уравнения, связывающие xi(n + 1) и xi(n) вида:

xi (n + 1) = (1aii )xi (n) + aij x j (n) , i = 1, 2, ..., 6,

(3)

j i

 

где

 

aij ≥ 0 ; i, j = 1, ..., 6.

(4)

Рассмотрение модели (3) , (4) иллюстрирует ту важную роль, которую играют матрицы специального вида, рассмотренные в п.20 и называемые иногда «входными – выходными»ы.

40. Матрицы Минковского-Леонтьева. Рассмотрим теперь специальный класс неотрицательных матриц, определяемый условиями

0 ≤ aij

; i,

j =

 

 

,

(5)

1, n

n

 

 

 

 

 

 

 

 

aij

1 ,

j =

1, n .

(6)

i= 1

Эти матрицы возникают в связи с решением алгебраических уравнений вида

n

 

xi = aij x j + yi , i =

 

,

(7)

1, n

j = 1

 

3 9 2

встречающихся при рассмотрении моделей, очень похожих на предыдущую, для межотраслевых производственных процессов.

Теорема 4. Если 0 ≤

aij

; i, j =

1, n

 

и

n

 

 

 

 

 

 

aij

< 1, j = 1, n ,

j =

1

 

 

 

 

 

то система (7) имеет единственное решение, которое положительно, если все yi положительны.

Если aij > 0, то матрица (I – A)1 также положительна.

! Задания для самостоятельной работы

1. С помощью формулы Фробениуса найти обратные матрицы для следующих матриц:

 

3

 

5

 

7 2

 

 

 

1 2 3

0

 

 

 

 

 

 

 

 

1

 

2

 

3 4

 

 

 

 

 

0 1 2

3

 

 

 

 

 

 

 

а)

 

 

 

 

;

 

б)

 

 

;

 

 

 

 

 

 

− 2 − 3

 

3 2

 

 

 

3 0 1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

3

 

5 4

 

 

 

 

 

2 3 0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 − 2

 

3

4

 

 

 

1 − 3

5 − 7

 

 

 

 

 

2 1

 

− 4 3

 

 

 

 

 

3

− 5

7 −

1

 

 

 

 

в)

 

 

 

;

г)

 

 

 

 

 

 

 

3 − 4

 

− 1 −

 

2

 

 

 

5

− 7

1 −

3

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 3

 

 

2 − 1

 

 

 

 

 

7

− 1

3 −

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Проверить теорему Фробениуса для следующих неразложи-

мых неотрицательных матриц:

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

3

 

 

1

 

1

3

 

 

5

4

 

 

 

1

3

а)

 

4 5 6

 

; б)

 

1 5 1

 

; в)

;

г)

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

0

1

5

 

 

 

3

 

1

 

 

 

 

8

9

 

 

 

 

4

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

3. Доказать, что из условия Ax ≥0 при всех x ≥0 следует, что

A≥ 0.

4.Показать, что из A B ≥ 0 следует, что λ(A) ≥ λ(B) (A B означает, что A B ≥0).

3 9 3

Лекция 46*

Однородные функции. Теоремы о неявных функциях

Дано определение однородной функции, приведена теорема Эйлера об однородных функциях, рассмотрены теоремы о неявных функциях.

10. Однородные функции. В лекции 48 при рассмотрении однородных уравнений первого порядка дано понятие однородной функции двух переменных. Дадим определение однородной функции нескольких переменных: функция нескольких переменных называется однородной функцией этих переменных степени m, если при умножении этих переменных на произвольную величину t функция умножается на tm, т.е. имеет место тождество:

f (tx, ty) = tmf (x, y) или f (tx, ty, tz) = tmf (x, y, z)

(1)

(или f (tx() = f (tx1, ..., tx n ) = t m f (x1, ..., xn ) = t m f (x() )

при любых допустимых значениях переменных x, y, z, t (x1, ..., xn, t). Число m может быть любым фиксированным вещественным числом.

Если, например, m =

1

, то t m = t и t должны быть положительными.

2

 

 

Примеры однородных функций:

а) однородной функцией третьей степени является любой однородный многочлен от х и у третьей степени, т.е. функция вида f (x, y) = ax3 + bx2y + cxy2 + dy3;

б) дроби

x 3

+

y 3

,

xy

 

,

x +

y

суть однородные функции

x 2

+

y 2

x 2 +

y 2

x 2 +

y 2

 

 

 

 

степеней соответственно 1, 0 и –1.

Дифференцируя тождество (1) по t и применяя правило дифференцирования сложной функции, получим, полагая u = tx и v = ty:

xf

(u, v) + yf

(u, v) = mt m 1 f (x, y) .

 

 

u

v

 

Полагая t = 1,

находим

 

 

xf x(x, y ) +

yf y(x, y ) = mf (x, y ) ,

(2)

что выражает следующую теорему Эйлера об однородных функциях.

3 9 4

Сумма произведений частных производных однородной функции на соответствующие переменные равна произведению самой этой функции на степень ее однородности.

Естественно, считаем, что функция f (x, y) имеет непрерывные частные производные по соответствующим переменным, которые

использовались при доказательстве.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для функции y =

f (x()

n переменных x1,

...(,

xn

формулировка

теоремы Эйлера следующая: если функция

 

y =

 

 

f (x )

является в обла-

сти X Rn дифференцируемой однородной степени m, то

x(

X спра-

ведливо равенство:

 

 

f

′(x() (

 

 

 

f (x()

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mf

(

 

 

 

n

 

x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x ) =

 

(

x

=

 

x j

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

j = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

Если m = 0, то, положив в тождестве (1)

 

t

=

 

 

t =

 

 

 

 

, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

y

 

 

z

 

 

 

 

 

 

 

 

 

 

x1

 

f (x, y) = f 1,

 

 

или

f (x, y, z) = f

 

1,

 

,

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2

 

 

 

 

 

x

n

 

 

 

 

 

 

 

 

 

 

 

 

или

f (x1, ..., x n ) =

 

f

1,

 

 

, ...,

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

 

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

Иногда однородную функцию нулевого измерения называют просто однородной.

Примеры однородных функций в экономике. Рассмотрим ос-

новные типы макроэкономических производственных функций.

1) Степенная производственная функция (функция типа КоббаДугласа) имеет вид

Y = AK α Lβ ,

где A, α , β – положительные числа. Здесь Y – национальный доход, К – капитал, L – труд. Отметим, что при α + β = 1, 0 < α < 1 получаем классическую функцию Кобба-Дугласа Y = AK α L1− α .

Очевидно, что степенная производственная функция является однородной степени m = α + β , которая определена для K ≥0 и L ≥0.

2) Функция с постоянными пропорциями имеет вид

Y= A min K , L ,

a b

3 9 5

где A, a, b – положительные числа. Вектор col (a; b) характеризует удельные затраты капитала и труда, необходимые для выпуска продукции в количестве А.

Эта функция определена для K 0, L 0 и является однородной первой степени.

3) Функция с постоянной эластичностью замены (CES) есть

Y = (AK − ρ + BL− ρ )mρ ,

где A, B, ρ, m – положительные числа.

Функция CES однородная, степени m и определена для K > 0, L > 0.

20. Существование неявных функций. Вначале

рассмотрим одно уравнение

 

F (x, y) = 0

(3)

и укажем те условия, при которых оно определяет единственным образом y как функцию от x, непрерывную и имеющую производную.

Теорема 1.

 

Пусть x = x0 и y = y0 – решение уравнения (3),

т.е.

F (x0, y0) = 0;

(4)

пусть F(x, y) и ее частные производные первого порядка по x и y – непрерывные функции при всех x и y, достаточно близких к x0 и y0, и, пусть, наконец, частная производная Fy(x, y) отлична от нуля при x = x0, y = y0. Тогда существует при всех x, достаточно близких к x0, одна определенная функция y(x), удовлетворяющая уравнению (3), непрерывная, имеющая производную и удовлетворяющая условию: y(x0)=y0.

Доказательство. Пусть для определенности Fy(x, y ) >0 при x = x0, y = y0. Т.к. по условию эта производная непрерывна, то она будет положительной и при всех значениях x и y, достаточно близких к x0 и y0, т.е. существует такое положительное число l, что F (x, y) и ее частные производные непрерывны и

 

Fy(x, y ) >0

(5)

при всех x и y, удовлетворяющих условию

 

 

x x0

 

l ,

 

y y0

 

l .

(6)

 

 

 

 

Далее функция F (x0, y) переменной y обращается в нуль при y = y0 в

силу (4), и возрастающая функция от y в промежутке (y0 l; y0 + l), в силу

(5), (6). Таким образом, число F (x0,

y0 l) < 0, а число F (x0, y0 + l) > 0.

Из непрерывности функции F (x, y)

следует, что F (x, y0 l) < 0, а

3 9 6

F (x, y0 + l) > 0 при всех x, достаточно близких к x0, т.е. существует такое положительное число l1, что

 

 

 

 

F (x,

y0 l) < 0 и F (x,

y0 + l) > 0

(7)

при

 

x x0

 

l1 . Пусть

m = min{l, l1}.

Тогда можно утверждать,

что

 

 

выполнены неравенства (5) и (7), если x и y удовлетворяют неравен-

ствам

 

y y0 l .

 

x x0

m ,

(8)

Если взять какое-нибудь определенное x, лежащее в промежутке (x0 m; x0 + m), то F (x, y) как функция от y будет в силу (5) возрастающей функцией в промежутке (y0 l; y0 + l), и в силу (7) – разных знаков на концах этого промежутка. Следовательно, она будет обращаться в нуль при одном определенном значении y из этого промежутка. В частности, если x = x0, то в силу (4) это значение y будет y = y0. Таким образом, доказано существование в промежутке (x0 m; x0 + m) определенной функции y(x), являющейся решением уравнения (3) и удовлетворяющей условию y(x0) = y0. Другими словами, из предыдущих рассуждений следует, что при любом фиксированном x (x0 m; x0 + m) уравнение

(3) имеет единственный корень, лежащий внутри промежутка (y0 l; y0 + l). Покажем теперь, что найденная функция y(x) будет непрерывной при x = x0. Действительно, при любом заданном малом положитель-

ном εчисла F (x0, y0 ε) и F (x0, y0 + ε) будут в силу (7) разных знаков, а,

следовательно, будет существовать такое положительное I, что F (x, y0

l) и F (x, y0 + l) – разных знаков, если только x x0 < I , а поэтому при

x x0 < I корень уравнения (3), т.е. значение найденной функции y(x),

удовлетворяет условию

y y

0

ε

 

 

, что и доказывает непрерывность

y(x) при x = x0.

Покажем теперь существование производной y(x ) при x = x0. Пусть x = x x0, y = y y0 есть соответствующее приращение y. Следовательно, x = x0+ x и y = y0+ y удовлетворяют уравнению (3), т.е. F (x0 + x, y0 + y) = 0, и в силу (4) имеем

F (x 0 + ∆x, y0 + ∆y) F (x0 , y0 ) = 0 .

Принимая во внимание непрерывность частных производных,

можно переписать это равенство так ( Л.42):

 

 

[Fx0 (x0 , y0 ) + ε1]x +

[Fy0 (x0 , y0 ) + ε2 ]y = 0 ,

(9)

где ε и ε → 0, если x и y

0 и где обозначены через

Fx0 (x0 , y0 ) и

1

2

 

 

3 9 7

Fy0 (x0 , y0 ) значения производных при x = x

, y = y . Из доказанной ранее

 

 

 

 

 

 

 

 

 

 

0

 

 

0

непрерывности y(x) при x=x0 следует, что y

0, если x 0.

Из уравнения (9) получаем

 

 

 

 

 

 

 

 

 

 

 

 

y

Fx0 (x0 , y0 ) + ε1

 

 

 

= −

 

 

 

 

 

 

.

 

 

x

F (x

0

, y

0

) + ε

 

 

 

 

y0

 

 

 

 

2

 

 

Переходя к пределу при x

 

0, имеем

 

 

y(x0 ) = −

Fx0 (x0 , y0 )

.

 

 

 

 

 

 

 

 

F

(x

0

, y

0

)

 

 

 

 

 

 

 

 

y0

 

 

 

 

 

 

 

Таким образом, доказаны непрерывность и существование производной функции y(x) только при x = x0. Если взять какое – либо другое значение x из промежутка (x0 m; x0 + m) и соответствующее значение y из промежутка (y0 l; y0 + l), являющееся корнем уравнения (3), то для этой пары значений x, y опять выполнены все условия теоремы 1, и в силу доказанного y(x) будет непрерывной при указанном значении x.

Теорема 1 доказана. Отметим, что совершенно так же, как и выше, формулируется и

доказывается теорема о существовании неявной функции z (x, y), определяемой уравнением

Φ (x, y, z) = 0 .

Рассмотрим теперь систему m функциональных уравнений

Fi (x(, y() = 0 , i =

 

,

(10)

1, m

где x( R n , y( R m . Требуется отыскать решение системы уравнений (10) как совокупность функций y = ϕ (x() i = 1, m , т.е. таких функций

ϕi (x() , которые при подстановке в систему (10) вместо yi обращают ее

втождество. i i

Матрицу частных производных левых частей системы (10) называют матрицей Якоби:

 

 

 

Fj (x(, y()

 

 

 

 

 

 

 

 

 

 

 

 

 

(

(

 

 

 

,i =

1, m

,

yi

J (x

, y) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j = 1, m

 

 

 

 

 

 

 

 

а ее определитель det J (x(, y()

– якобианом.

 

3 9 8

Имеет место

 

 

 

 

 

 

 

Теорема 2.

 

 

 

 

 

 

 

Пусть функции Fi( x(, y() , i =

 

 

 

дифференцируемы по y(

в окре-

1, m

стности точки ( x(0;y(0 ) , причем матрица J( x(, y() непрерывна в точке

( x(0;y(0 ) . Тогда, если F( x(0 , y(0) =

0 , i =

 

 

и det J( x(0 , y(0 ) 0 ,

то для

1, m

i

 

 

 

что для всех x( Oδ( x(0 )

 

любого ε> 0 существует такое δ> 0,

суще-

ствует единственное решение системы (10) yi = ϕi( x() , i = 1, m , удовлетворяющее условию y( Oε( y(0 ) , причем это решение непрерывно и дифференцируемо для всех x( Oδ( x(0 ) .

Здесь Oδ(x(0 ) ( Oε(y(0 ) ) – δ–окрестность точки x(0 (ε– окрестность точки y(0 ).

! Задания для самостоятельной работы

1. Проверить следующие функции на однородность и для однородных функций записать теорему Эйлера:

 

 

 

 

 

 

 

 

y

x

 

 

1)

x + 2y z ; 2) yz 2xy + x 2 + 3z 2 ; 3) e

x

+ e

z

;

 

 

 

y

 

x

 

 

 

 

y

 

 

4)

cos ln

 

 

+ sin

 

; 5)

y + xy + yz ; 6) tg

 

 

z .

 

z

 

 

 

x

 

 

x

 

 

 

 

 

 

 

2. Используя теоремы 1 и 2, проверить существование неявных функций в следующих уравнениях и системах уравнений:

1) x 3 2x 2 y 2 + 5x + y 5 = 0 , y

 

x = 1 = 1;

 

 

 

 

 

 

 

 

2)

x +

xy

y =

 

 

2 + 1, y

 

x = 2 = 1;

 

 

 

 

 

 

 

 

 

 

 

 

3)

ex sin y +

ey

cos x =

1, y

 

x = 0 =

0 ;

 

 

 

 

 

 

 

 

 

 

 

 

x 2 +

xy

1

+ y

2

 

y 2

= 1,

 

 

 

 

 

 

 

 

 

 

 

4)

 

 

 

 

 

 

 

1

 

 

 

 

 

y1

 

x = 0 = 1 , y2

x = 0 = 2 ;

 

 

 

 

 

 

 

 

 

 

=

cos 2,

 

 

sin x ey1 + cos y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1 =

x +

tg y2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

5)

 

 

 

 

 

 

 

π

 

 

 

y1

 

 

= 1 , y2

 

x = 0 =

.

 

 

 

 

 

 

 

 

 

 

 

 

x

2

+

2

=

 

 

 

x = 0

 

4

 

ln

 

y1

 

4

 

y2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 9 9

Лекция 59*

Интегрирование дифференциальных уравнений с помощью рядов

На примере линейного дифференциального уравнения второго порядка иллюстрируется метод нахождения его решения в виде степенного ряда.

10. Разложение решения в степенный ряд. Этот прием особенно удобен в применении к линейным дифференциальным уравнениям. Проиллюстрируем его применение на примере уравнения второго порядка. Рассмотрим однородное дифференциальное уравнение второго порядка

y′′ + p(x )y′ + q(x )y = 0 ,

(1)

где коэффициенты p(x), q(x) представляются в виде рядов, расположенных по целым положительным степеням х, так что уравнение (1) можно переписать в виде

y′′ +

(a

+ a x +

a x 2

+ ...)y′ +

(b

+ b x +

b x 2

+ ...)y = 0 .

(2)

 

0

1

2

 

0

1

2

 

 

Решение этого уравнения также будем искать в виде степенного ряда

 

 

 

 

 

 

 

 

 

 

 

 

 

y =

ck x k .

 

 

 

 

 

(3)

Подставляя выражение для

k = 0

 

 

 

 

получаем:

y и его производных в (2),

 

 

 

 

 

 

k (k 1)ck x k 2

+

ak x k

kc k x k 1

+

bk x k ck x k

=

0 .

(4)

k = 2

 

k = 0

k = 1

 

k = 0

k = 0

 

 

 

Перемножив степенные ряды, собирая подобные члены и приравнивая к нулю коэффициенты при всех степенях х в левой части (4), получаем ряд уравнений:

x 0

2 1c2 +

a0c1 + b0c0 = 0 ,

 

 

x 1

 

3 2c3 +

2a0c2 + a1c1 + b0c1 + b1c0

= 0 ,

(5)

 

 

 

 

 

 

 

x 2 4 3c4 + 3a0c3 + 2a1c2 + a2c1 + b0c2 + b1c1 + b2c0 = 0 .

− − − − − − − − − − − − − − − − − − − − − −

Каждое последующее из уравнений (5) содержит одним искомым коэффициентом больше, чем предыдущее. Коэффициенты c0 и c1 остаются произвольными и играют роль произвольных постоянных. Первое из уравнений (5) дает c2, второе – c3, третье – c4 и т.д. Вообще из (k + 1)-го уравнения можно определить ck+2 , зная c0, c1, ... , ck+1.

4 0 0