Студопедия.Орг Главная | Случайная страница | Контакты | Мы поможем в написании вашей работы!  
 

Типы дифференциальных уравнений



Линейное

Пример. Решить уравнение y¢¢¢ = cos x.

Решение. Интегрируя, находим

y¢¢ = ò cos x dx = sin x + C1,

y¢ = ò (sin x + C1)dx = - cos x + C1x + С2,

y = ò (- cos x + C1x +C2)dx = - sin x + C1x2/2 +C2x+C3.

Итак, общее решение

y = - sin x + C1x2/2 +C2x+C3.

В математической экономике большое применение находят линейные дифференциальные уравнения, и поэтому мы рассмотрим решение таких уравнений. Дифференциальное уравнение (14.1) называется линейным, если имеет вид:

рo(x)y(n)(x) + р1(x)y(n- 1)(x) +... + рn - 1(x)y¢(x) + рn(x)y(x) = f(x), (14.2)

где рo(x), р1(x),..., рn(x), f(x) - данные функции. Если f(x) º 0, то уравнение (9.2) называется однородным, в противном случае - неоднородным. Общее решение уравнения (9.2) есть сумма какого-либо его частного решения y(x) и общего решения соответствующего однородного уравнения:

рo(x)y(n)(x) + р1(x)y(n- 1)(x) +... + рn - 1(x)y¢(x) + рn(x)y(x) = 0. (14.3)

Однородное

Если коэффициенты рo(x), р1(x),..., рn(x) постоянные, то уравнение (14.2) принимает вид:

рoy(n)(x) + р1y(n- 1)(x) +... + рn - 1y¢(x) + рny(x) = f(x) (14.4)

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

Соответствующее уравнению (9.4) однородное уравнение выглядит так:

рoy(n)(x) + р1y(n- 1)(x) +... + рn - 1y¢(x) + рny(x) = 0. (14.5)

Без ограничения общности можно положить рo = 1 и записать уравнение (9.5) в виде

y(n)(x) + р1y(n- 1)(x) +... + рn - 1y¢(x) + рny(x) = 0. (14.6)

Решение уравнения (9.6) будем искать в виде y = e kx, где k - постоянная. Имеем: y¢ = k e kx, y¢¢ = k2 e kx,..., y(n) = kn e kx. Подставляя полученные выражения в (9.6), будем иметь:

e kx (kn + р1kn-1 +... + рn-1k + рn) = 0.

Т.к. e kx ¹ 0, то

kn + р1kn-1 +... + рn-1k + рn = 0. (14.7)

Характеристическое

Равенство (14.7) есть алгебраическое уравнение с неизвестным k. Оно называется характеристическим уравнением для дифференциального уравнения (9.6). Характеристическое уравнение есть уравнение n -й степени, следовательно, оно имеет n корней, среди которых могут быть кратные и комплексные. Если k1, k2,..., kn - действительные и различные корни уравнения (9.7), то - частные решения уравнения (14.7), а общее имеет вид

y = .

Дифференциальное уравнении 2 порядка

Рассмотрим подробно линейное однородное дифференциальное уравнение второго порядка с постоянными коэффициентами:

y¢¢ + рy¢ +qy = 0. (14.8)

Его характеристическое уравнение имеет вид

k2 + рk + q=0 (14.9)

и в зависимости от значения дискриминанта D = р2 - 4q возможны три случая.

1. Если D>0, то корни k1 и k2 уравнения (9.9) действительны и различны, тогда общее решение имеет вид:

y = c1 exр(k1x) + c2 exр(k2x).

2. Если D = 0, т.е. корни k1 и k2 действительные и равные, то общее решение находится по формуле:

y = (c1 + c2x) exр (k1x).

3. Если D<0, то корни комплексные, k1 = a + bi, k2 = a - bi, где i - мнимая единица. Тогда общее решение таково:

y = (c1 cos bx+c2 sin bx) exр (ax).

Пример. Решить уравнение y¢¢ - y = 0.

Решение. Характеристическое уравнение имеет вид k2 - 1 = 0, корни которого k1 = 1, k2 = -1 действительны и различны. Общее решение:

y = c1 e x + c2 e -x.

Разностные уравнения

На практике простейшие разностные уравнения возникают при исследовании например величины банковского вклада. Эта величина является переменной Yx, представляющей сумму, которая накапливается по установленному закону при целочисленных значениях аргумента x. Пусть сумма Yo положена в банк при условии начисления 100 r сложных процентов в год. Пусть начисление процентов производится один раз в год и x обозначает число лет с момента помещения вклада (x = 0, 1, 2,...). Обозначим величину вклада по истечении x лет через Yx. Мы получаем

Yx = (1+r)Yx-1.

Если начальная сумма составляет Yo, мы приходим к задаче отыскания решения полученного разностного уравнения, подчиненного начальному условию Yx = Yo при x = 0. Полученное разностное уравнение содержит Yx и значение этой переменной на один год раньше, т.е. Yx-1; в данном случае аргумент x явно не входит в разностное уравнение.

Вообще говоря, обыкновенное разностное уравнение устанавливает связь между значениями функции Y = Y(x), рассматриваемой для ряда равноотстоящих значений аргумента x, но можно без ограничения общности считать, что искомая функция определена для равноотстоящих значений аргумента с шагом, равным единице. Таким образом, если начальное значение аргумента есть x, то ряд его равноотстоящих значений будет x, x+1, x+2,... и в обратном направлении: x, x-1, x-2,.... Соответствующие значения функции будем обозначать Yx,Yx+1,Yx+2,... или Yx, Yx-1, Yx-2,.... Определим так называемые разности различных порядков функции Yx с помощью следующих формул:

Разности первого порядка

D Yx = Yx+1 - Yx,

D Yx+1 =Yx+2 - Yx+1,

DYx+2 = Yx+3 - Yx+2,

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

Разности второго порядка

D2Yx= DYx+1 - D Yx,

D2Yx+1= D Yx+2 - DYx+1,

D2Yx+2= D Yx+3 - DYx+2,

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

Разности третьего порядка

D3Yx= D2Yx+1 - D2Yx,

D3Yx+1= D2Yx+2 - D2Yx+1,

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

Обыкновенным разностным уравнением называется уравнение, связывающее значения одного независимого аргумента x, его функцииYx и разностей различных порядков этой функции DYx, D2Yx, D3Yx,.... Такое уравнение можно записать в общем виде следующим образом:

j(x, Yx, DYx, D2YxD3Yx, DnYx) = 0, (14.15)

которое по форме аналогично дифференциальному уравнению.

Порядком разностного уравнения называется порядок наивысшей разности, входящей в это уравнение. Разностное уравнение (10.1) часто удобнее записать, пользуясь не разностями неизвестной функции, а ее значениями при последовательных значениях аргумента, то есть выразить DYx, D2Yx, D3Yx,... через Yx, Yx+1, Yx+2,.... Уравнение (10.1) можно привести к одной из двух форм:

y(x, Yx, Yx+1,...,Yx+n) = 0, (14.16)

x(x, Yx, Yx-1,...,Yx-n) = 0. (14.17)

Общее дискретное решение Yx обыкновенного разностного уравнения n -го порядка представляет функцию x (x = 0, 1. 2,...), содержащую ровно n произвольных постоянных:

Yx = Y(x, C1, C2,..., Cn).

Пример 14.7. Паутинообразная модель. Пусть рынок какого-либо отдельного товара характеризуется следующими функциями спроса и предложения:

D = D(P), S = S(P).

Для существования равновесия цена должна быть такой, чтобы товар на рынке был распродан, или D(P) = S(P).

Цена равновесия `P задается этим уравнением (которое может иметь множество решений), а соответствующий объем покупок-продаж, обозначаемый через `X, - следующим уравнением:

`X = D (`P) = S(`P).

Динамическая модель получается при наличии запаздывания спроса или предложения. Простейшая модель в дискретном анализе включает неизменное запаздывание или отставание предложения на один интервал:

Dt = D (Pt) и St = S (Pt-1).

Это может случиться, если для производства рассматриваемого товара требуется определенный период времени, выбранный за интервал. Действие модели таково: при заданном Pt-1 предшествующего периода объем предложения на рынке в текущем периоде будет S (Pt-1), и величина Pt должна установиться так, чтобы был куплен весь объем предложенного товара. Иными словами, Pt и объем покупок-продаж Xt характеризуются уравнением:

Xt = D (Pt) = S (Pt-1).

Итак, зная исходную цену Po, с помощью этих уравнений мы можем получить значения P1 и X1. Затем, используя имеющуюся цену P1, из соответствующих уравнений получим значения P2 и X2 и т.д. В общем изменение Pt характеризуется разностным уравнением первого порядка (одноинтер­вальное отставание):

D (Pt) = S (Pt-1).

Решение можно проиллюстрировать диаграммой, представленной на рис.5, где D и S - соответственно кривые спроса и предложения, а положение равновесия (со значениями `P и `X) соответствует точке их пересечения Q. Цена в начальный момент времени равна Po. Соответствующая точка Qo на кривой S дает объем предложения в период 1. Весь этот предложенный объем товара раскупается при цене P1, заданной точкой Q1 на кривой D с той же ординатой (X1), что и Qo. Во второй период времени движение происходит сначала по вертикали от точки Q1 к точке на кривой S, дающей X2, а затем по горизонтали - к точке Q2 на кривой D. Последняя точка характеризует P2. Продолжение этого процесса и дает график паутины, показанный на рис. 5. Цены и объемы (покупок - продаж) в последовательные периоды времени являются соответственно координатами точек Q1, Q2, Q3,... на кривой спроса D. В рассматриваемом случае последовательность точек стремится к Q. При этом точки поочередно располагаются на левой и правой стороне от Q. Следовательно, и значения цены Pt стремятся к `P, располагаясь поочередно по обе стороны от`P. Точно так же обстоит дело и с объемами покупок - продаж (X t).

X S

(D; S)

x2 Q2

Q

`x

x3 Q3

x1 Q0 Q1

D

O P0 P2 `P P3 P1 P

Рис. 5.

Решение можно получить алгебраически для случая линейных функций спроса и предложения: D = a +aP, S = b +bP. Значения равновесия `P и `X будут заданы уравнениями

`X = a +a`P = b +b`P,

то есть

`P = (a - b)/(b - a), `X = (ba - ab)/(b - a). (14.18)

Дискретная динамическая модель задается уравнением

X t = a +aP t = b +bP t-1. (14.19)

Ищем сначала решение, дающее равновесие. Для этого положим P t = `P, X t = `X для всех значений t:

`X = a +a`P = b +b`P. (14.20)

Получаем те же значения `P и `X, что и в (10.4). Следовательно, если в каком-либо периоде существовали цены и объемы, обеспечивающие равновесие, то в динамической модели (10.5) они сохранятся и в последующих периодах.

Вычтем уравнение (10.6) из (10.5) и положим р t = P t -`P, x t = X t -`X. Тогда

x t = aр t = bр t-1. (14.21)

Уравнения (14.21) аналогичны (14.19), за исключением того, что они описывают отклонения от уровней равновесия (теперь уже известно, что таковые существуют). Оба эти уравнения являются разностными уравнениями первого порядка. Положим c = b/a и подставим его в уравнение (14.21), так что разностное уравнение относительно р t будет

р t = c р t-1. (14.22)

При данном значении р o в момент t = 0 из (10.8) получаем решение:

р t = рo c t,

или

P t = `P + (Po - `P) c t.





Дата публикования: 2015-01-10; Прочитано: 265 | Нарушение авторского права страницы | Мы поможем в написании вашей работы!



studopedia.org - Студопедия.Орг - 2014-2024 год. Студопедия не является автором материалов, которые размещены. Но предоставляет возможность бесплатного использования (0.015 с)...