Вычислительная математика - ИСЭМ СО РАН
Вычислительная математика - ИСЭМ СО РАН
Вычислительная математика - ИСЭМ СО РАН
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