27.10.2014 Views

Вычислительная математика - ИСЭМ СО РАН

Вычислительная математика - ИСЭМ СО РАН

Вычислительная математика - ИСЭМ СО РАН

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Заметим, что (6)–(13) – смешанная система неклассических уравнений, в которой (6)<br />

– билинейное уравнение Вольтерра I рода, (7) – билинейное уравнение Вольтерра II рода<br />

относительно переменной x(t) и I рода относительно переменной y(t).<br />

Вопросы существования решения систем уравнений Вольтерра I рода с двумя переменными<br />

пределами интегрирования рассмотрены в монографии А.С. Апарцина [1].<br />

Система (6)–(13) малоизучена.<br />

2. О свойствах двухсекторной модели развивающейся экономической системы<br />

Рассмотрим несколько примеров.<br />

1. Пусть в исходной модели (6)–(13) y(t) = 1. Это означает, что весь произведенный<br />

продукт направляется на развитие самой системы, внешний продукт равен нулю. В этом<br />

случае модель описывается лишь уравнением второго рода. Зададим длину предыстории<br />

равной 1 2 . x(t) =<br />

∫ t<br />

x(s)ds, t ∈ [ 1 2 , 5 ], (14)<br />

2<br />

t− 1 2<br />

x 0 (t) = t, t ∈ [0, 1 ). (15)<br />

2<br />

Решение уравнения с переменными пределами (14) при условии (15) происходит поэтапно.<br />

Отрезок [t 0 , T ] разбивается на части последовательностью точек<br />

{<br />

max z, если a(T ) > z k−1 , z 0 = t 0 ,<br />

z k = a(z)=z k−1<br />

z N = T, если a(t) ≤ z N−1 .<br />

В нашем случае<br />

z 0 = 1 2 , z 1 = 1, z 2 = 3 2 , z 3 = 2, z 4 = 5 2 .<br />

Используя аналитический вид решения на предыстории x 0 (t), решаем классическое уравнение<br />

Вольтерра II рода<br />

x(t) =<br />

∫ t 0<br />

a(t)<br />

x 0 (s)ds +<br />

∫ t<br />

t 0<br />

x(s)ds<br />

на отрезке [z 0 , z 1 ). На следующем этапе, для получения решения на отрезке [z 1 , z 2 ), в<br />

качестве x 0 (t) используется полученное x(t) на [z 0 , z 1 ). Решение (14)–(15) описывается<br />

формулой<br />

⎧<br />

t, 0 t < 1,<br />

2<br />

− ⎪⎨<br />

7 8 et−1/2 + t + 1, 1<br />

t < 1,<br />

2 2<br />

− 7<br />

¯x(t) =<br />

8 et−1/2 + t + 1 + e t−1 (− 11 + 7t), 1 t < 3,<br />

8 8 2<br />

− 7 8 et−1/2 + t + 3 + 2 et−1 (− 11 + 7t) + 8 8 et−3/2 (− 7 16 t2 + 29t<br />

− 143),<br />

3<br />

t < 2,<br />

16 64 2<br />

− ⎪⎩<br />

7 8 et−1/2 + t + 2 + e t−1 (− 11 + 7t) + 8 8 et−3/2 (− 7 16 t2 + 29t<br />

− 143)+<br />

16 64<br />

+e t−2 ( 7<br />

48 t3 − 9 8 t2 + 13t − 11), 2 t < 5.<br />

4 3 2<br />

Из рисунка 1 видно, что в точке стыковки предыстории с основным отрезком z 0 решение<br />

претерпевает разрыв, на [z 0 , z 2 ) решение убывает и на следующих отрезках стремится<br />

89

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!