О решении краевых задач для вырожденных систем линейных интегро-дифференциальных уравнений методом наименьших квадратов Текст научной статьи по специальности «Математика»

УДК 517.518

БО!: 10.14529/ттр150207

О РЕШЕНИИ КРАЕВЫХ ЗАДАЧ ДЛЯ ВЫРОЖДЕННЫХ СИСТЕМ ЛИНЕЙНЫХ ИНТЕГРО-ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ МЕТОДОМ НАИМЕНЬШИХ КВАДРАТОВ

Б.Д. Нгуен, В.Ф. Чистяков

В настоящее время, при анализе сложных электрических и электронных схем, часто встречаются системы, включающие в себя взаимосвязанные дифференциальные, интегральные и алгебраические уравнения. Алгебраические уравнения отвечают за отличие в моделях балансовых соотношений, в частности, законов сохранения или уравнений состояния, системы дифференциальных уравнений описывают динамику процесса. Если процесс обладает последействием, то математическая модель может включать и интегральные уравнения (ИУ). Системы взаимосвязанных дифференци-эльных, алгебраических и интегральных уравнений можно записать в виде векторных интегро-дифференциальных уравнений с матрицей неполного ранга в области определения при старшей производной искомой вектор-функции. Численное решение краевых и начальных задач для таких систем сопряжено с большими трудностями. В данной работе обсуждается метод наименьших квадратов и приведены результаты численных расчетов.

Ключевые слова: вырожденная система; общее решение; интегро-дифферен-циалъные уравнения; краевая задача; метод наименьших квадратов.

Введение

В течение последних сорока лет большое внимание уделяется системам дифференциальных уравнений с матрицей неполного ранга или вырожденным оператором в области определения при старшей производной искомой вектор-функции и численным методам их решения [1-5]. Системы интегро-дифференциальных уравнений с матрицей неполного ранга в области определения при старшей производной искомой вектор-функции и особенно численные методы решения краевых задач для них в прошлое тридцатилетие исследовались фрагментарно. Сейчас это быстро растущая область исследования (см. например, монографию [6], и работы [7-10] с приведенной там библиографией). Ранее изучались только постановки начальных задач. Краевые задачи практически не рассматривались. Методы, применяемые в работах [11, 12] при решении краевых задач для дифференциально-алгебраических уравнений (ДАУ), сложно, а для систем с прямоугольными матрицами коэффициентов невозможно адаптировать к нашим задачам. В данной работе изучается метод наименьших квадратов, применяемый в работе [13] для решения начальных задач для дифференциально-алгебраических уравнений.

1. Постановка задачи

Рассмотрим систему интегро-дифференциальных уравнений (ИДУ)

(Л1 + V )х := А(г)х(г) + Б(г)х(г) + [ к (г, в)х(в№ = ¡, г е т = [а,в], (1)

■) а

где A(t), B(t), K(t, s) — (v x ^-матрицы, x(t) = dx(t)/dt, x(t), f = f (¿)-искомая и заданная вектор-функция соответственно, Л1Х := A(t)x + B(t)x, V—оператор Вольтерра с ядром K(t, s)

Cx(a) + Dx(ß) = a, (2)

где C, D — (m x ^-матрицы, a-заданный вектор. Предполагается, что входные данные достаточно гладкие, и выполнено условие

rankA(t) < min{v, n}, t £ T. (3)

Если v = n, то условие (3) эквивалентно равенству det A(t) = 0,t £ T.

Система (1) называется замкнутой, если v = n, переопределенной, если v > n, и недо-v<n

В работе используются нормы g-мерного вектора b = (bi,b2, ■■■,bq)T, £ Rq, и вектор-функции b(t) = (bi(t), b2(t),bq(t))T, t £ T, вычисляемые по правилам

q rß

\\b\\2E = E b2, \\b\h = ^ma^\b31 , \\b(t)\\l2{T) = ^ \\b(s)\\2E ds, ||b(t)||c(T) = max ||b(t)yj ,

T

t

опускаться, если это не вызывает путаницы. Включения V(t) £ Сг(Т), i > 1, где V(t) — матрица или вектор-функция, означают, что все производные всех ее элементов непрерыв-i V(t) £ C(T).

Запись V (t) £ L2T) означает, что все эле менты V (t) являются функциями суммируемыми T

Под решением задачи (1), (2) мы будем понимать вектор-функцию x = x(t) £ C1(T),

T

вию (2).

Задачи (1), (2), удовлетворяющие условию (3), часто возникают при анализе сложных электрических и электронных схем [14]. Схема обычно изображается в виде направленного графа (цепи), состоящего из ветвей и узлов с фиксированным направлением каждой ветви. В общем случшэ уравнения цбпи состоят из двух частей; уравнения бвлвнев токов (расходов) в узлах (первый закон Кирхгофа), уравнений бэЛШИСе! напряжений (перепадов давлений) в контурах (второй закон Кирхгофа), дополняемых уравнениями ветвей (обобщенный закон Ома). Законы Кирхгофа можно представить в матричной форме с использованием контурной и узловой матриц: A, B соответственно. В результате ряда преобразований (см. [7]) мы получим интегро-дифференциальное уравнение вида

t

(BLM BR)x+(t) ¡x(s)ds+<r =0,

to

где x = x(t)—вектор-функция токов на ветвях, RB, CB - диагональные матрицы активных сопротивлений и величин, обратных емкостям ветвей соответственно; LB - матрица индук-тивностей, которая может быть вырожденной, uc(to) + ев — вектор напряжений на зажимах источников (тока и напряжения). При выписывании уравнений Кирхгофа и матрицы индук-тивностей мы можем получить линейно-зависимые уравнения, и конечная система может быть переопределенной. В зависимости от суточного температурного режима сопротивления, индуктивности, емкости ветвей могут зависеть от времени.

2. Теоремы о разрешимости

Уравнения вида (1) могут обладать сложной внутренней структурой. В таком случае мы можем идти только по пути выделения классов уравнений, для которых доказаны утверждения о существовании и единственности решений, с последующей разработкой численных методов решения.

Определение 1. Если существует оператор

Л k = L0(t) + W)d + ... + Lk (t)( d)\

где Lo(t), Li(t),..., Lk(t) — (v x v)-матрицы из C(T), со свойством

Ak ◦ (Ai + V)y = A(t)y + B(t)y + f K(t, s)y(s)ds Vy € Ck+i(T),

■J a

где A(t), B(t), K(t,s) — (v x п)-матрицы, непрерывные в своих областях определения,

rankA(t) = min{v, n} Vt € T,

то будем говорить, что для системы (1) определен левый регуляризирующий оператор (JIPO). Минимально возможное число к называется индексом системы (1).

Определение 2. Если для пучка (v x п)-матриц AA(t) + B(t) существует (v x v)-матрица P (t) со свойствами

det P (t) = 0 Vt € T, P (t)[AA(t) + B (t)] = A^f^ + (B^) ,

rank ^B^t)) = min{v,n} Vt € T,

то будем говорить, что пучок матриц AA(t) + B(t) имеет индекс один на отрезке T.

Лемма 1. [6] Пучок (п x п)-матри ц AA(t) + B(t) имеет индекс один на о трезке T тогда и только тогда, когда выполнен критерий <,щнг-степень>:

rank A(t) = deg det[AA(t) + B(t)] = const = r Vt € T. (4)

Далее нам потребуется такое утверждение.

Лемма 2. [6] Пусть (v x п)-матрица A(t) € Cl(T), i = 0,1, 2, ■ ■ ■, имеет, постоянный ранг r на T. Тогда существуют квадратные матрицы P(t), Q(t) € Сг(T) соответствующей размерности со свойствами

det P(t) = 0, det Q(t) =0 Vt € T, P(t)A(t)Q(t) = jj) , t € T, (5)

где Er—единичная матрица размерности r.

Следствие 1. Для, (v x п)-матриц A(t) полного ранга, когда, rankA(t) = min{v,n} Vt € T,

P(t), Q(t)

A(t)Q(t) = (Ev 0) , v < n, P(t)A(t) = f^E"^ , n<v, t € T.

Свойства систем (1), имеющих индекс один, определяются следующим утверждением.

Теорема 1. Пусть для системы (1) выполнены условия

1) A(t),B(t),f (t) е C1(T),K(t,s) e C1(T x T);

2) пучок матриц XA(t) + B(t) имеет индекс один на отрезке T. Тогда:

1) система, (1) имеет, индекс один;

2) существуют квадратные матрицы P(t), Q(t) е C1(T) подходящей размерности со свойствами det P(t) = 0, det Q(t) = 0 Ш e T, система,

P (t)(Ai + V )Q(t)u — P (t)f,

имеет, вид

(Er 0 0\ Г

[о о o) |u

E

Ev—r

Bi3 0

>ж a

й va W - ™ - o ■

v — n,

v <n,

(6) (7)

(Е 0) (2)+С? СО Ч£) СО=р=(!) ■

где С(Ь)-некоторый (г х г)-блок,

VцX = Кц (г,8)х(8)й8,

■) а

Кц(Ь, в)—блоки подходящей размерности произведения V(Ь)К(t,s)Q(t), иц — компоненты вектор-функции и, соответствующие блочной структуре матриц коэффициентов.

Более того, если п < V, гапк(Л(Ь)\Б(Ь)) = п Ш € Т, то существуют квадратные матрицы V(Ь), Q(t) € С;(Т) подходящей размерности со свойствами det V(Ь) = 0, det Q(t) = 0 Ш € Т, система,

V (Ь)(Л; + V Ш)и = V (¿)/,

имеет вид

fEr 0^

о о о о,

® + ( 0 VJ СО

Доказательство. В качестве ЛРО можно принять оператор

Vil Vl2N

+ I V21 V22 V3I V32 /

fi"

te)—P (t)f—f

(8)

Q = (Pi(t)\ + Q \ 0 J + \(d/dt)P2(t)

((d/dt)P2(t)) ■ ^(i) P(t)'

(9)

где Р (Ь)—матрица из определения 2. Далее, умножим систем у (1) на матрицу Р (Ь) и произведем замену х(Ь) = Q(t)y(t), где Р(Ь), Q(t) матрицы из равенства (5). Для определенности будем считать, что V < п. Получим

Р (Ь)(Л; + V Ш)у =

= (Е 0)(У2) + (Б;; Б;2)(У2) + $£)(£)=р=<-)

0

Блочная структура в равенстве (10) выписана с учетом разбиения на блоки произведения

р (1)А(гш),

со=у, ^ *=айч ^8)г(8)аз>

Кц(Ь, з)—блоки подходящей размерности произведения Р(Ь)К(t,s)Q(t).

В силу условия 2) теоремы гапк ( Е ^ ] = п} Ш € Т. В силу этого гапкБ22 =

\Б21 Б22 /

V — г Ш € Т. Согласно следствию 1 найдется невырожденная матрица Q со свойством

Б2^ = (Еи-Г О) После замены переменных У = О ^^ V система (10) преобразуется к виду

е^ оЩ+(Б2; вз)£)+е^: I1)©, (11)

А (у2\ ^ /цч (Ег —Б:з\

где у; = V;, У2 = Q )■ .Умножая систему (11) на матрицу I и производя

VГ{)з) V О Еу-г)

(иЛ ( Ег М (VЛ пп (ал

замену из = vз, I = I I, мы приведем систему (11) к виду (о), где

\и2/ \—Б2: Еи-т) VV2/

1 (Ь) = Б;;— Б;2Б2:, Б;3 = Б;3. Суперпозиции этих действий дадут нам матрицы V(Ь), Я,(Ь).

Существование форм (7), (8) доказывается аналогично. □

Следствие 2. Пусть выполнены условия теоремы 1 и V < п. Тогда системы разрешимы при любой (Ь), и их общее решение имеет вид

х(Ь,с) = Хг (Ь)с + С (Ь), (12)

с(Ь) = I К(t,s)f(з)ёз + Co(t)f(Ь) + I К+ Со(^(Ь),

■) а ■) а

х(Ь,с) = Хг (Ь)с + Ф), (13)

Ш= [ *к(ь,з)/ Шз + сош (Ь),

■)а

где Хг(Ь), Хг(Ь) — (п х г)-матрицы, К(Ь, з), С0(Ь), К(Ь, з), С0(Ь) — (п х п)-матрицы, гладкие в областях определения, 'ш(Ь)—произвольная гладкая вектор-функция, с-вектор произволь-

Т

Доказательство. Рассмотрим случай замкнутой системы. Из второго блочного уравнения (7) получим

и2 = (Еп-г + У22)-: [—У2: и; + /2 ]. (14)

и2

и: + 1 и: + У/и: = Г, (15)

где У — оператор Вольтерра с некоторым ядром КГ — некоторое алгебраическое выражение относительно f: (Ь), f2 (Ь) и соответствующих операторов Вольтерра, действующих на перечисленные функции.

Здесь учтено, что произведение операторов Вольтерра является оператором Вольтерра и (Еп-Г + У22)-1 = Еп-Г + У22, где У22 —некоторый оператор Вольтерра. Далее, разрешив (15) как обыкновенное дифференциальное уравнение вида

иг + Juг = Ф, Ф = Г — уиг,

получим интегральное уравнение Вольтерра

иг = и (г) с + Шиг + Ф,

где и(¿)—матрицант системы иг = —Juг, Ш = IVУ с—вектор произвольных постоянных, Шд = ! и(г)и-1(8)д(в)й8, Ф = Ш = ШГ, иг = (Ег + Ш)-1[и(г)с + Ф]. (15)

■ ' а

Поставляя и\ в (14), найдем выражение для и2 и после умножения на Q получим вид общего решения системы (13). Полагая в (6) компоненту из произвольной функцией и повторяя вышеприведенные выкладки, получим (12). Следствие доказано. □

В случае переопределенной системы для первых двух блочных уравнений из (8) мы можем повторить рассуждения из следствия и получим формулу (13). Из формулы (8) мы видим, что вектор произвольных постоянных должен удовлетворять еще такому уравнению

С(г)с = ф(Ь), (16)

гт = (у V ) ( (Ег + Ш)-ги(г)

с(1) = ^г ^ ^—(Еп-г + V22)-lV2г(Er + Ш)-ги(г)

т = -(V31 Уз2)х , WT,,

(Er + W )-1Ъ) — (En-r + V22)-1V2i(Er + W

где обозначения в формуле (16) соответствуют обозначениям из формул (6) - (15). Согласно [2], система (16) разрешима относительно c тогда и только тогда, когда выполнено соотно-

Ф(г) = L(t)C-в, (17)

где C-полуобратная матрица C: CC-C = C,

ß ß C = J CT(s)C(s)ds, в = J CT(s)^(s)ds. (18)

a a

Все решения описываются формулой

c = C-e + [Er —C-C]c, (19)

где c—произвольный вектор Rr. При выполнении условия (17), а ему удовлетворяют все функции из образа оператора системы (1) (Л1 + V)z, z € C1(T), общее решение переопределенной системы имеет вид

x(t,c) = Xd(t)c + С(t), (20)

где

X (t) = 0 ( (Er + W)-1U(t) \[F

Xd(t) = Q \—(En-r + V22)-1V21 (Er + W)-1U(t)) [En —CC]

ц/- Bulletin of the South Ural State University. Ser. Mathematical Modelling, Programming

ÖU & Computer Software (Bulletin SUSU MMCS), 2015, vol. 8, no. 2, pp. 81-94

С(*) = ха(г)с-е + ¡\(г, в)/ (в)йв + Со(*)/(*),

■) а

где Х^(£) — (п х й)-матрицы, й = гапк[Ег — С-С], К(£, в), С0(£) — (п х т)-матрицы, гладкие в областях определения. Из (19) следует, что размерность многообразия переопределенных систем не зависит в наших условиях от свободного члена /(¿).

Следствие 3. Если выполнены условия теоремы 1 и условие (19), то краевая задача (1), (2) разрешима тогда и только тогда, когда разрешимы относительно векторов с, с системы алгебраических уравнений:

Zic = a - CZ(а) - DZ(в), Zi = CXr(a) + DXr(в), v < n, (21)

Z2C = a - C((a) - D((e), Z2 = CXr(a) + DXr(в), v = n, (22)

Z3C = a - C((a) - DZ(e), Z3 = CXd(a) + DXd(p), v > n. (23)

Если решение системы c, c единственно, то единственно решение задачи (1), (2). Замечание 2. Случай n < v, rank(A(t)|B(t)) > n Ш € T, в работе не рассматривается. Замечание 3. Пусть найдется матрица P(t) из определении 2, но условие о полноте ранга

не выполнено: rank Al(t) < minív, n} Шt € T, Al(t) = | Ai(t) ) , то можно попытаться найти

\B2(t)J

матрицу Pl(t) с аналогичными свойствами и повторить процесс.

Лемма 3. Пусть в системе (1) A(t), B(t) € CA(T), K(t,s) € CA(T x T), где CA(T), CA(T x T)-пространства вещественно-аналитических функций, и индекс системы равен к.

Тогда существуют матрицы Pl(t) € CA(T), i = 1,2, ••• к, det Pl(t) = 0 Ш € T, соответствующей размерности со свойствами: произведение операторов, построенных по формуле (9) является JIPO. Иначе говоря, можно принять

л, = п а. * = ДТИ + ((¿/адм) ■ ДМ = ^

Для замкнутых систем доказательство приведено в [6]. Для незамкнутых систем алгоритм доказательства почти не меняется.

Если для системы (1) определен ЛРО и v > n, то многообразие ее решений конечномерно. В обратном случае мы этого не можем гарантировать.

Пример 1. Рассмотрим систему

(Ai + V)Х = (0 0) Х ^Y 1) Х + ¡' (gds g0°t)) X(S)dS = 0'' € [0' 4-

где 7-вещественный параметр, g(t)-функция из CA(T). Имеем Х2 = -txi, где (xi Х2)т = x. Отсюда (7 - 1)х1 = 0 ^ dim ker (A1 + V) < ж при 7 = 1. Если 7 = 1, то можно принять x = (-u(t) tu(t))T, где u(t) - произвольная функ цпя из C1[0, 1]. В качестве базиса в пространстве решений можно принять функции фj = (y-tj tj+1^ , j = 0,1, • • •. Иначе говоря, dim ker (Ai + V) = ж. Выбором функции g(t) можно задавать любой наперед

y=1

3. Численный метод

Системы индекса 1 обладают важным свойством.

Теорема 2. Пусть выполнены условия теоремы 1, краевая задача (1),(2) разрешима и для некоторой вектор-функции х€(Ь) € С:(Т), выполнены условия

\\ФШа(т) <е, У/ < е, (24)

где

фе(Ь) = (Л: + У)хе(Ь) — f (Ь), К = Схе(а) + Бхе(в) — а. Тогда, найдутся решение задачи х(Ь) и положительные числа е0, к такие, что

\\х€(Ь) — х(Ь)\\С(Т) < ке, е < во-

Более того, если в краевых условиях (2) С = Еп, Б = О (краевая задача является задачей Коши), и для, некоторой вектор-функции хе(Ь) € С:(Т), выполнены условия

\\Ф&)\\12(т) < е, хе(а) — х(а) = = О, (25)

х(Ь) е:, к:

\\х,(Ь) — х(Ь)\\2Ыт) < к:е, е < еЛ.

Доказательство. Выпишем для невязки у(Ь) = х(Ь) — х€(Ь) краевую задачу

(Л: + У )у(Ь) = фС у (а) + Бу(в) = Не. (26)

Пусть V < п. Полагаем в формуле (12) произвольную вектор-функцию w(t) = О и f (Ь) = фе(Ь). Выписываем систему (21) применительно к задаче (26) Z:c = Ь, Ь = —С((а) — (в).

Система совместна. Из условий (24) и формулы (12) в результате простых оценок (с учетом гладкости) вытекает, что справедливо неравенство \\Ь\\/ < ке, где к—некоторая константа. Согласно формуле (19) общее решение с = Z-Ь + [Ег — Z:-Z:]v, где V— произвольный вектор. Полагая V = О, видим, что вектор констант в этом случае имеет оценку \\\\/ < к:е. Подставляя этот вектор в (12), мы находим вектор функцию у(Ь) и нужное нам решение х(Ь) = у(Ь) + х€(Ь). Для случаев V = п, п < V утверждение доказывается аналогично. Не нужно только оговаривать условие w(t) = О.

Во второй части утверждения мы имеем дело с задачей Коши. Пусть V < п. Из формулы (15) мы видим, что при нулевых начальных данных с = О, так как и (а) = Ег. Следовательно, полагая в формуле (12) w(t) = О и подставляя в нее f (Ь) = ф€(Ь), после несложных оценок получим нужное нам неравенство. Следует только учесть: если € С(Т), е(Ь) € Ь2(Т), то те(Ь) € Ь2(Т), \Ше(Ь)\\Ь2(т) < \Ш\\с(т) МШ^т) [15]. □

Итак, для краевых задач (1), (2), где система имеет индекс 1, существует гарантированная оценка отклонения приближающей вектор-функции от множества решений краевой задачи.

Будем искать приближения к решению краевой задачи (1), (2) в виде полинома

Рг(Ь) = со + с:Ь + ... + а?, г = 1,2,...,

где со, с:,..., с^-векторные коэффициенты, подлежащие определению. По методу множителей Лагранжа [16], построим функцию:

Г в

5 (со,с:,...,сг,Х) = \т)\\2Ыт) + (Х, П)= \\Г (Ь)с — f (Ь)\\Е М + (X, П), (27)

■¡а

5(t) = A(t)pi(t) + B(t)pi(t) + / K(t, s)pi(s)ds - f (t),

2

L2(T )

P(t)||| dt,

где (A, П)-скалярное произведение в пространстве Rm, c =

(СЛ (Ao\

ci , A = Ai

ci \AmJ

F(t) = A(t)Pi(t) + B(t)Pi(t) + / K(t, s)Pi(s)ds,

Pi(t) = (En tEn ... tiEn), П = Cx(a) + Dx(e) - a = C [c0 + cía + ... + ciai] + D[c0 + cie + ... + ciei] - a,

Будем искать коэффициенты co,ci,... ,ci, доставляющие минимум функции (27). Необходимым условием минимума функции многих переменных является равенство нулю ее частных производных первого порядка по независимым переменным. В функции (27) такими независимыми переменными являются коэффициенты co, ci,... ,ci, A. Выпишем систему линейных алгебраических уравнений порядка n(k + 1) + m относительно неизвестных co, ci,..., ci, A. Имеем

(OS/Oco = 0, dS/dci = 0,

dS/da = 0, [ dS/dA = 0.

С учетом вида функции (27) система выглядит таким образом:

2 гв FT(t) F(t) dt\ c i (dn /dc) \ = 2 rP

(2 /ав FT(t)F(t)dt) c + (dÜ/dc)A = 2 ¡в FT(t)f (t)dt, (d Ü/dc)T c = a.

(28)

В векторном виде система (28) имеет вид:

Ч a) = (¿ $ ® = ^

(29)

ГР ГР

Л = 2 Ет(г)Е(гщ А = дП/де, в = 2 Ет(г)/(г)М.

■) а ■) а

Решение системы (29) может быть осуществлено любым из известных методов (например, методом Гаусса). Подставляя найденные значения с = (со, е\,..., ег), А, получим функции Рг(г) = со + С1* + ... + е{Ьг - приближения к точным решениям задачи (1), (2).

Лемма 4. Пуст,ь V > и, и задача (1), (2) имеет, единственное решение. Тогда, существует единственное решение системы линейных алгебраических уравнений (29). Иначе говоря,

det A = 0.

(30)

t

в

t

Доказательство. Подставим в задачу (1), (2) полином pi(t) = c0 + c\t +... + Cit%, i = 1, 2,..., и вычислим соответствующие ему f (t) = fi(t), a = ai. В силу единственности решения исходной задачи в системах (22), (23) решения единственны. Следовательно, решения краевых задач с входными данными f (t) = fi(t), a = ai единственны. Это возможно тогда и только тогда, когда справедливо неравенство (30). □

Замечание 4. В случае v < n единственность решения задачи (1), (2) невозможна из-за того, что решение зависит от произвольных функций.

Оценим отклонение приближающего полинома от решения задачи Коши. В основу доказательства положим поиск многочлена, отклонение которого от решения мы можем оценить.

Теорема 3. Пусть выполнены условия теоремы 2, краевая задача (1), (2) является задачей Коши и решение x(t) единственно. Тогда, найдется решение системы (29) и соответству-

pi

S (c°, C1, . . . ,Ci,X) < Kn i)2) , Wpi(tt) - X(t)Wl(T) < (i2 . Kn i)2) ,

где ki = const, = const.

Доказательство. Рассмотрим функционал (27) на многочлене l^i(t), у которого i^i(t) есть

x(t)

\\ii(t) - X(t)\\c{T) = O(1/(i • In i)), S(ca, ci,...,5i,X) = O(1/(i2 • (In i)2)),

где co,ci,..., Ci—коэффициенты пол инома li(t). Так как та поли номе pi(t) достигается минимум, то

S (co, ci,...,ci,X) < S(co, ci,...,ci, X). Из второго утверждения теоремы 2 следует вторая оценка утверждения. □

Замечание 5. При определенных требованиях на гладкость решения задачи (1), (2) можно брать другие пробные многочлены, например, интерполяционный многочлен Лежандра, и получать другие оценки скорости сходимости.

4. Численные эксперименты

В качестве иллюстрации приведем несколько примеров. Пример 2. Пусть задана краевая задача

(0 20 X+(01) X+¡1 (г 0) x(s)ds=(1+^— —osr) е м>

(10 x(0)+G 0 x«=(0>

Индекс системы равен 1. Точное решение этой краевой задачи: x(t) = (sint cost)T. Результаты расчетов примера приведены в табл. 1. Здесь Err = \\(sin(t) — x\(tj), cos(t) — x2(tj))T||r

Из табл. 1 мы видим, что при увеличении степени i полинома Pi(t), погрешности расчетов уменьшаются.

Таблица 1

Результаты расчетов для примера 2

i = 3 i = 5

tj Xxitj ) X2 (tj ) Err tj Xi(tj) X2 (tj ) Err

0 0,0624 0,9758 0,0624 0 0,0001 0,9999 0,0001

0,3142 0,3669 0,9533 0,0579 0,3142 0,3095 0,9500 0,0011

0,6283 0,6117 0,8193 0,0239 0,6283 0,5878 0,8088 0,0003

0,9425 0,7928 0,6011 0,0162 0,9425 0,8080 0,5884 0,0010

3,1416 -0,0624 -0,9758 0,0624 3,1416 -0,0001 -0,9999 0,0001

Пример 3. Рассмотрим систему из примера 1 при g(t) = 0, t G [0,1] и ненулевом свободном

(о о)с о)*=г 7e')-tg[о, 1].

Система имеет индекс 2, если 7 = 1 и разрешима при любой f (t) G C2[0,1]. Решение системы единственно, и краевые условия не нужны, x(t) = (e1 e2t^T.

В настоящее время для любой конкретной стандартной разностной схемы можно указать значение параметра 7, при которых схема расходится. Результаты расчетов для примера 2 при 7 = 2 приведены в табл. 2. Здесь Err = ||(et — x1(tj), e2t — x2(tj))T||7.

Таблица 2

Результаты расчетов для примера 3

i = 3 i = 5

tj X1(tj) X2 (tj ) Err j xi(tj) X2 (tj ) Err

0 0,9795 1,0100 0,0205 0 0,9999 1,0000 0,1029 x 10"3

0,1 1,1456 1,2245 0,0405 0,1 1,1057 1,2214 0,5105 x 10"3

0,2 1,2844 1,4832 0,0630 0,2 1,2215 1,4918 0,0540 x 10"3

0,3 1,4068 1,8046 0,0569 0,3 1,3495 1,8222 0, 3362 x 10"3

1,0 2,7388 7,3791 0,0205 1,0 2,7184 7,3890 0,1029 x 10"3

Пример 4.

'1 1\ ( г 1\ г (1 + в Г

2 2| х + 12г + 1 2| х + I г о I х(в)йв = /3(г) г е [о, 1]

чз з) \ зг у ^ \ г о,

с краевым условием

(1 1) х(0) + (3 2) х(1) = 50.

Индекс системы равен 1. Решение краевой задачи единственно. В качестве решения был принят векторный многочлен

X3(t) =(!) + § t + Q t2 + Q t3

и вычислен соответствующий свободный член /з(t). Построена и решена система (29). Вычисленные коэффициенты многочлена рз (t) отличаются от коэффициентов Хз (t) на величину, не превышающую 1О_ 12.

Заключение

В работе получены условия разрешимости недоопределенных, переопределенных и замкнутых систем вида (1), включающих как частный случай ДАУ. Проанализирована возможность применимости метода наименьших квадратов для решения краевых задач (1), (2). Результаты численных экспериментов хорошо соответствуют теоретическим оценкам.

Литература

1. Сидоров, H.A. Исследование непрерывных решений задачи Коши в окрестности точки ветвления / H.A. Сидоров // Известия вузов. Математика. - 1976. - № 9. - С. 99-110.

2. Бояринцев, Ю.Е. Регулярные и сингулярные системы линейных обыкновенных дифференциальных уравнений / Ю.Е. Бояринцев.- Новосибирск: Наука, 1980. - 222 с.

3. Sviridyuk, G.A. Linear Sobolev Type Equations and Degenerate Semigroups of Operators / G.A. Sviridyuk, V.E. Fedorov. - Utrecht; Boston; Köln; Tokyo: VSP, 2003.

4. Свиридюк, P.A. Инвариантные пространства и дихотомии решений одного класса линейных уравнений типа Соболева / P.A. Свиридюк, A.B. Келлер // Известия вузов. Математика. - 1997. № 5. - С. 60-68.

5. Келлер, A.B. Численное решение задачи оптимального управления вырожденной линейной системой обыкновенных дифференциальных уравнений с начальными условиями Шоуолтера-Сидорова/ A.B. Келлер // Вестник ЮУрГУ. Серия: Математическое моделирование и программирование. - 2008. - № 27 (127), вып. 2. - С. 50-56.

6. Чистяков, В.Ф. Алгебро-дифференциальные операторы с конечномерным ядром / В.Ф. Чистяков. - Новосибирск: Наука. Сибирская издательская фирма РАН, 1996. -278 с.

7. Чистякова, Е.В. Методы исследования и решения вырожденных интегро-дифференциальных уравнений и их приложения: дис. ... канд. физ.-мат.наук / Е.В. Чистякова. - Иркутск, 2006.

8. Чистяков В.Ф. О разрешимости линейных интегро-алгебраических уравнений и численных методах их решения / В.Ф. Чистяков // Сибирский математический журнал. -2013. - Т. 54, № 4. - С. 932-946.

9. Булатов, М.В. Об одном семействе вырожденных интегродифференциальных уравнений/ М.В. Булатов, Е.В. Чистякова // Журнал вычислительной математики и математической физики. - 2011. - Т. 51, № 9.- С. 1665-1673.

10. Фалалеев, М.В. Вырожденные интегро-дифференциальные операторы в банаховых пространствах и их приложения / М.В. Фалалеев, С.С. Орлов // Известия вузов. Математика. - 2011. - № 10. - С. 68-79.

11. Clark, Kenneth D. Numerical Solution of Boundary Value Problems in Differential-Algebraic Systems / D. Kenneth Clark, Linda R. Petzold // SIAM Journal on Scientific and Statistical Computing archive. -1989. - V. 10, № 5. - P. 915-930.

12. Marz, R. On Difference and Shooting Methods for Boundary Value Problems in Differential-Algebraic Equations / Marz R. // ZAMM Journal of Applied Mathematics and Mechanics: Zeitschrift fur angewandte Mathematik und Mechanik. - 2006. - V. 64, № 11. - P. 463-473.

13. Чистяков, В.Ф. Применение метода наименьших квадратов для решения линейных дифференциально-алгебраических уравнений / В.Ф. Чистяков, Е.В. Чистякова // Сибирский журнал вычислительной математики. - 2013. - Т. 21, № 1. - С. 81-95.

14. У ШсЖОВ, Е.И. Статическая устойчивость электрических систем / Е.И. У ШсЖОВ. - Новосибирск: Наука, 1988.

15. Маслов, В.П. Операторные методы / В.П. Маслов. - М: Наука, 1973.

16. Бертсекас, Д. Условная оптимизация и методы множителей Лагранжа / Д. Бертсекас. -М.: Радио и связь, 1987.

17. Бахвалов, Н.С. Численные методы / Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков. -М.: Наука, 1987.

Банг Дык Нгуен, аспирант, кафедра «Вычислительная техника:», Иркутский государственный технический университет(г. Иркутск, Российская Федерация), ducbang@mail.ru.

Виктор Филимонович Чистяков, доктор физико-математических наук, Институт динамики систем и теории управления СО РАН (г. Иркутск, Российская Федерация), chist@icc.ru.

Поступила в редакцию 13 марта 2015 г.

MSC 65R20 DOI: 10.14529/mmpl50207

On Boundary Value Problems for Singular Systems of Linear Integro-Differential Equations Method of Least Squares

B.D. Nguyen, Irkutsk State Technical University, Irkutsk, Russian Federation, ducbang@mail.ru,

V.F. Chistyakov, Institute for System Dynamics and Control Theory of Seberian Branch of Russian Academy of Sciences, Irkutsk, Russian Federation, chist@icc.ru

At present, in the analysis of complex electrical and electronic circuits, the system often includes interconnected differential, integral and algebraic equations. Algebraic equations are responsible for the difference of balance relations in the models, in particular, the conservation laws or equations of state, the system of differential equations describing the dynamics of the process. If the process has afteraction, the mathematical model can include integral equation (IE). The systems of interconnected differential, algebraic and integral equations can be written in the form of vector integro-differential equations with a matrix at the highest derivative of a searched vector-function of not full rank in the domain. Numerical solution of boundary and initial problems for such systems conjugates with great difficulty. In this paper we discuss the least squares method and the results of numerical calculations.

Keywords: degenerate system; general solution; integro-differential equations; boundary value problem; least squares method.

References

1. Sidorov N.A. Study of Continuous Solutions of the Cauchy Problem in a Neighborhood of the Branch. Russian Mathematics (Izvestiya VUZ. Matematika), 1976, vol. 20, issue 9, pp. 77-87.

2. Boyarintsev Yu.E. Reguljarnye i singuljarnye sistemy linejnyh obyknovennyh differencial'nyh uravneniy [Regular and Singular Systems of Linear Ordinary Differential Equations]. Novosibirsk, Nauka, 1980. (in Russian)

3. Sviridyuk G.A., Fedorov V.E. Linear Sobolev Type Equations and Degenerate Semigroups of Operators. Utrecht, Boston, Koln, Tokyo, VSP, 2003. DOI: 10.1515/9783110915501

4. Sviridyuk G.A., Keller A.V. Invariant Spaces and Dichotomies of Solutions of a Class of Linear Equations of the Sobolev Type. Russian Mathematics (Izvestiya VUZ. Matematika), 1997, vol. 41, issue 5, pp. 57-65._

5. Keller A.V. Numerical Solution of Optimal Control Problem Degenerate Linear System of Ordinary Differential Equations with Initial Showalter-Sidorov Conditions. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2008, no. 27 (127), issue 2, pp. 50-56. (in Russian)

6. Chistyakov V. F. Algebro-differentsial'nye operatory s konechnomernym yadrom [Algebraic-Differential Operators with Finite-Dimensional Kernel]. Novosibirsk, Nauka, Siberian Publishing Company of the RAS, 1996.(in Russian)

7. Chistyakova E.V. [Research Methods and Solving Degenerate Integro-Differential Equations and Their Applications. Candidate's Dissertation in Mathematics and Physics]. Irkutsk, 2006.(in Russian)

8. Chistyakov V.F. On the Solvability and Numerical Methods for Solution of Linear Integro-Algebraic Equations. Siberian Mathematical Journal, 2013, vol. 54, issue4, pp. 746-758. DOI: 10.1134/S0965542511090065

9. Bulatov M.V., Chistyakova E.V. On a Family of Singular Integro-Differential Equations. Computational Mathematics and Mathematical Physics, 2011, vol. 51, issue 9, pp. 1558-1566. DOI: 10.1134/S0965542511090065

10. Falaleev M.V., Orlov S.S. Degenerate Integro-Differential Operators in Banach Spaces and Their Applications. Russian Mathematics, 2011, vol. 55, issue 10, pp. 59-69.

11. Clark K.D., Petzold L.R. Numerical Solution of Boundary Value Problems in Differential-Algebraic Systems. SIAM Journal on Scientific and Statistical Computing archive, 1989, vol. 10, issue 5, pp. 915-930.

12. Marz R. On Difference and Shooting Methods for Boundary Value Problems in Differential-Algebraic Equations. ZAMM Journal of applied mathematics and mechanics: Zeitschrift fur angewandte Mathematik und Mechanik, 2006, vol. 64, issue 11, pp. 463-473.

13. Chistyakov V. F., Chistyakova E.V. Application of the Least Squares Method for Solving Linear Differential-Algebraic Equations. Numerical Analysis and Applications, 2013, vol. 6, issue 1, pp. 77-90. DOI: 10.1134/S1995423913010102

14. Ushakov E.I. Staticheskaya ustoychivostV elektricheskikh sistem [Static Stability of Electrical Systems]. Novosibirsk, Nauka, 1988.

15. Maslov V.P. Operatornye metodi [Operator Methods]. Moscow, Nauka, 1973.

16. Bertsekas D.P. Conditional Optimization and Lagrange Multiplier Methods. N.Y., Academic Press Inc. 1987.(in Russian)

17. Bahvalov N.S., Zhidkov N.P., Kobelkov G.M. Chislennyye metody [Numerical Methods]. Moscow, Nauka, 1987. (in Russian)

Received March 13, 2015