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

УДК 517.948

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

Е.В. Табаринцева

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

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

Введение

В статье рассматривается задача с обратным временем для полулинейного дифференциально-операторного уравнения.

Данная задача поставлена некорректно, поэтому основными вопросами при ее исследовании являются вопросы построения устойчивого приближенного решения и оценки погрешности приближенного решения.

Для линейных некорректно поставленных задач общая теория регуляризуемости разрабатывалась в статьях Л.Д. Менихеса и В.А. Винокурова (см., напр., [1, 2]).

В работах В.К. Иванова, В.Н. Страхова, их учеников и последователей (см., напр., [3]) была создана теория, и разработан аппарат для получения оценок погрешности методов приближенного решения линейных некорректно поставленных задач на компактных множествах (классах корректности). В этой теории естественным образом вводятся понятия оптимального и оптимального по порядку метода приближенного решения.

Соответствующие понятия были введены и для нелинейных некорректно поставленных задач (см., напр., [4]).

Пусть Н — гильбертово пространство, М С и , а С[Н] — пространство непрерывных отображений, действующих в Н.

Рассмотрим операторное уравнение

Аои = /; и € Н; / € Н, (1)

Ао € С[Н] — взаимно-однозначный оператор.

Предположим, что при / = /о существует точное решение ио уравнения (1), которое принадлежит множеству М, но точное значение правой части нам не известно, а вместо него дано приближенное значение / € Н такое, что У/о — /|| < 5. Требуется по исходым данным задачи М и / определить приближенное решение уравнения (1) и оценить его уклонение от точного решения.

Определение 1. Семейство операторов (Т : 0 <5 < 5о} будем, называть методом приближенного решения уравнения (1) на множестве М, если для любого 5 € (0; 5о] оператор Т непрерывно отображает пространство Н в Н и Т / ^ ио при 5 ^ 0 равномерно на множестве М при условии ЦЛоио — / || < 5 .

Рассмотрим следующую величину, характеризующую точность метода (Т : 0 <5 < 5о} на множестве М :

Д(Т) = йир(||и — Т/1| : и € М, ЦЛои — /|| < 5}.

Обозначим

ш\(г; М) = 8ир(||и1 — и21 : и\,и2 € М, ЦЛои — Лои2|| < т} модуль непрерывности оператора, обратного к Ло , на множестве ЛоМ.

Определение 2. Метод (Т : 0 <5 < 5о} будем называть оптимальным по порядку на множестве М , если существует число к такое, что для любого 5 € (0; 5о]

Д(Т) < кш1 (5; М).

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

Различные подходы к приближенному решению нелинейных некорректно поставленных задач предложены и исследованы, например, в монографии [5].

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

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

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

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

1. Задача с обратным временем для дифференциальнооператорного уравнения

Пусть Н - гильбертово пространство, Л - линейный неограниченный положительно определенный самосопряженный оператор с областью определения ^(Л), плотной в Н.

Рассмотрим задачу вычисления элемента р € Н такого, что решение задачи Коши

Ни

— = —Ли + /(г,и(г)); г € (*о; Т), (2)

и(го) = р, р € н, 0 < го < т,

удовлетворяет условию и(Т) = х. Здесь / : [0; Т] х Н ^ Н - отображение, удовлетворяющее условию Липшица по переменной и и условию Гельдера по переменной г , т.е. существуют постоянные Ь> 0 , М > 0 , 0 < а < 1, такие, что

||/(иьг) — /(и2,г)||я < £||и1 — и21|н + м|г1 — г2|а

для всех г1, г2 € [0; Т], и1, и2 € Н.

Зафиксируем число г > 0. Рассморим задачу Коши

Ни

— = —Ли + /(г,и(г)); г € (0;Т), (3)

и(0) = ио, ио € Н.

Обозначим р = и(го). Рассмотрим множество

М = (р € Н : ||ио|| < г}.

Предположим, что при заданном х € Н существует точное решение р € Н поставленной обратной задачи, принадлежащее множеству М.

Элемент х € Н нам не известен, а вместо него дано приближенное значение х<5 € Н такое, что ||х — Х<5II < 5. Требуется по исходым данным задачи определить приближенное решение р^ задачи с обратным временем и оценить его уклонение от точного решения. Задача Коши (2) равносильна интегральному уравнению

t

и(г) = е-А(^о)р + У е-А^-т) / (т,и(т ))Нт (4)

to

(см., напр., [6]).

Наряду с задачей Коши (2) рассмотрим задачу Коши для соответствующего линейного уравнения

Ни

- = —Ли г € (го; т), (5)

и(го) = р, р € н, 0 < го < т.

Обозначим х = и(Т).

Выполняется следующая лемма [7]

Лемма 1. Пусть р1,р2 € Н, И1(ж,г), П2(ж,г) — соответствующие решения задачи (2),

И1(ж,г), ^2(ж,г) — решения задачи (5). Тогда при каждом г € [го; Т] выполняются неравен-

ства

е~ЬТеЬТ ||^1 — ^21 < 11и 1 — и2 || < еЬТ ||^1 — ^21.

Рассмотрим величины:

ш(М, 5) = 8ир( 11р 1 — р21 : р1, р2 € М, ||х1 — х21 < 5} — модуль непрерывности для нелинейной обратной задачи,

ш(М,5) = 8ир(11р 1 — р21| : р1, р2 € М, ||х1 — х2|| < 5} — модуль непрерывности для линейной обратной задачи.

Справедлива следующая лемма.

Лемма 2. Справедливы неравенства

2е^о(1-еЬТV г~° 5Т < ш1(5, М) < 2е^о(1+еЬТ)г5Т.

Доказательство. Пусть р1,р2 € М. По лемме 1 выполняется оценка

||«1(г) — «2(г)| < е^||е"^(и° — ио)| (0 < г < т), (6)

в частности, при г = го получаем неравенство

||р1 — р21 < е“°Не-»(и1 — и§)||. (7)

Из леммы 1 следует также

11е_^(ио — ио)| < еЬТеЬТ|и1(г) — и2(г)| (8)

в частности, при г = Т получаем неравенство

11е-АТ(ио— ио)| <еЬТеЬТ ||х1— хр2|. (9)

Рассмотрим элементы р1, р2, принадлежащие множеству М, т.е. ||и°|| < г, ||и°|| < г. Рассмотрим число 51 = е^Те 5, множество М = (0; г), где В = е-^° и оператор Ло =

е-А(Т-to). Обозначим

ш(5,М) = 8ир{||,г1 — ^21 : *1, *2 € М, ||Ло(^1 — ^2>| < 5}—

модуль непрерывности оператора Л-1 на множестве ЛоМ. Вычислим величину ш(5, М) по схеме, предложенной в [3]:

~ Т—£° £°

ш(5, М) = 2г 5 т . (10)

Рассмотрим элементы *1 = е-^°ио, *2 = е-^°и2, принадлежащие множеству М, так как

||В-1*11 = ||еА°*1|| = 11^11 < г,

11В—1 *21 = ||еА° *2|| = 11 и° | < г.

Из (9) следует ||Ло(*1 — *2)|| = ||е-АТ(ио — ио)| < еЬТеЬТ5 = 51. По определению модуля

непрерывности оператора Л-1, обратного к Ло, на множестве ЛоМ

о

11*1 - *2ІІ = Не-А*°(«о - «0)|| < й(5,М) = 2еьтг1-Т51Т

Следовательно, в силу (7) с учетом определения числа 51

г. -і 4° ^

ІІРі. - р2|| < 2еь*°г1 т £оТ . (11)

Следовательно, ^(5, М2) < 2е^*°(1+еЬТV т ° 54°.

Оценим ^1 (5, М2) снизу. Пусть р>1,р>2 Є М. Неравенство (8) при Ь = го запишем в виде

ІІР1 - Р2ІІ> е-^4° ||е-А*°(«1 - и§)||. (12)

Рассмотрим элементы *1 = е-А*°и°, *2 = е-А*°«0, принадлежащие множеству М = В£(0; г), где В = е-А*°. Рассмотрим оператор Ао = е-А(Т-*°). Вычисляя величину (2)(5, М), находим

~ Т—4° 4°

ш(5, М) = 2г~^5 т . (13)

Пусть *1 = е-^°ио, *2 = е-^°и° - произвольные элементы из М, т.е. ||и°|| < г, ||ио| < г, причем |Ло(*1 — *2)|| < е-ЬТ5 = 51.

Обозначим и1(г) решение задачи Коши

Ни

— = —Ли + /(г, и) (0 < г < т)

и(0) = ио1 ;

и2(г) - решение задачи Коши

Ни

— = —Ли + /(г, и) (0 < г < т)

и(0) = и°;

р1 = и1(го); р2 = и2(го).

По определению элементы р1, р2 € М, причем из неравенства (6), записанного при

г = т,

||х1 — х2| < еЬТ||е-^°(и° — и°)|| < еЬТ51 = 5.

Таким образом, из условий *1, *2 € М и ||Ло(*1 — *2)| < 51 следует р1, р2 € М и

||х1 — х2|| < 5. Следовательно, из (12)

^1(5, М2) = 8ир(11р 1 — р21 : р1,р2 € М2; ||х1 — х21 < 5} >

ТТ Т Т Т_£° ^°

> е-2^е 8ир( : *1, *2 € М, ||Ло(*1 — * — 2)||< 51 = й(5, М) = e-2toe 2г^~51Т ,

т.е. с учетом определения числа 51

ш1(5, М2) > 2е^°(1-еТТ)2гТ—°5Т.

□

2. Метод приближенного решения задачи с обратным временем

Обозначим через (Ед : Л > 0} разложение единицы, порожденное оператором Л. Пусть Ла - линейный ограниченный оператор в Н, определяемый формулой

а

Лаи = J ЛНЕди.

о

Вместо неустойчивой обратной задачи для уравнения (2) рассмотрим задачу вычисления элемента ра = иа(го) € ЕаН, где иа(г) удовлетворяет условиям

Ниа (г)

= —Л«иа(г) + Еа/(г,иа(г)), г € (го,т);иа(т) = Е«х. (14)

Рассмотрим задачу Коши

Ниа (г)

Нг

= —Лаиа(г) + Еа/(г,иа(г)), г € (го,т); иа(го) = ра. (15)

Так как Аа - ограниченный самосопряженный оператор в Н, то задача Коши (15) равносильна интегральному уравнению

Выполняется

Теорема 1. Для любого элемента х € Н существует элемент ра € ЕаН, такой, что решение и(г) задачи Коши (15)удовлетворяет условию иа(Т) = Еах-

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

3. Оценка погрешности метода проекционной регуляризации

Обозначим ра = иа(го), где иа(г) - решение задачи (14); Ра = иа(£о), где Ма(£)- решение задачи (14) с приближенными исходными данными. В качестве приближенного решения задачи с обратным временем рассмотрим элемент ^а^ = Ра(<5)Х<5 при подходящем выборе зависимости а = а(5).

Рассмотрим величину

характеризующую точность построенного метода приближенного решения задачи (2) на классе корректности М.

Воспользуемся неравенством

і

и“(£) = е-А“(і-іо)£ара + е-А“(і-т)Еа/(т,иа(т ))^т.

/

(16)

іо

Амі(а, 5) =вир{У^а - р|| : р Є М, Ух - ха1| < 5},

Дм (а, 5) < Аі (а) + А2(а, 5),

где

Д2(а, 5)

йир ура - рау;

\\x-xs К

Ді(а) = вир ||ра - р||.

^ем

Оценим величину Д2 (а, 5). Рассмотрим функцию ^а(£) = в(*-*°)аиа(г). Функция ^а(£), очевидно, удовлетворяет дифференциальному уравнению

где да(£, V) = Еае(* *°)а/(£, е (* *°)а^) : [£о; Т] х Н ^ Н - отображение, удовлетворяющее условию Липшица по переменной V и условию Гельдера по переменной г, г € [0; Т].

По построению функция Vа (г) удовлетворяет также условиям

Vа (£о) = р«;

(18)

(19)

(20)

іо

Из равенства (20) при Ь = Т с учетом (19) следует равенство

т

еа(Т-*0)£аХ = е-(А“-ае)(*-*0)ЕаРа + ^ е-(А“-аЕ)(Т-т)Еада(т, Vа(r))^Г. (21)

£°

Из (21) следует

т

е-(А“-ае)(*-*0)ЕаРа = еа(Т-*0)ЕаХ - ^ е-(А“-аЁ)(т^Е^аМа(т))^Т. (22)

£°

Подставляя выражение (22) в равенство (20) , имеем

т

V01 (Ь) = е(А“-аЁ)(Т-*)еа(Т-*°)Еах - ^ е(Аа-аЕ)(т-)Еада(т, vа(r))^т. (23)

Аналогично, функция v0(^) удовлетворяет уравнению

т

г?(Ь) = е(Аа-аЕ)(т-)еа(т-*°)ЕаХй - ^ е(Аа-аЕ)(т-)Еада(т, vа(т))^т. (24)

I С' - '' '-ПаУаУ'1 , ^

I

Из (23) и (24) следует неравенство

т

1К - vа(^)У < уе(А“-аЁ)(Т-*)Уеа(Т-*°)Ух - Хй|| + ^ Не^а-^МЩК - vа(т)Уdт. (25)

Из неравенства (25) с учетом леммы Гронуолла и неравенства

||е(А“-аЕ)(Т-*)П < тах е(Л-а)(т-) < 1 0<Л<а

следует

|К - v?(t)|| < еа(т-*°)еь(Т-*)5,

откуда при Ь = Ьо

||ра - Ра(Ь)11 < еа(т-*°)еь(т-*°)5.

Следовательно, для величины Д2(а, 5) имеем оценку

Д2 (а, 5) < еа(т-*°) еЬТеЬТ 5. (26)

Оценим величину Д1(а).

Рассмотрим функцию и(Ь), удовлетворяющую интегральному уравнению (4) и функцию иа(Ь) , удовлетворяющую (16) и функцию и(Ь) = Еаи(Ь). Рассмотрим равенство

т

й(Ь) = еАа(т-4)ЕаХ - / еАа(т-)Еа/(т,й(т))^г. (27)

Из (4) следует равенство

т

Х = е-А(т-*°)р + / е-А(т-т)/(т,й(т))^т. (28)

*°

Подставляя выражение (28) в равенство (16) и учитывая, что еА“(т 4)Еа = еА(т 4)Еа, имеем равенство

т

(Ь) = е А(с *°)ЕаР е а(с тЕа/(т,и(т))^т + У еА(т с)Еа(/(т,и(т)) - /(т,иа(т)))^т.

(29)

с°

С учетом равенства (27) из (29) следует

т

*(Ь) - и(Ь) = Еаи(Ь) - и(Ь) + J еА(т- с)Ео(/(т, и(т)) - /(т, иа(т)))^т. (30)

t

Из (30) с учетом неравенства ||еА(т с)Ео,|| < еа(т t) следует

т

е -а(т - - и(Ь)|| < е -а(т - ^ЦЕаф) - и(Ь)|| + е -а(т- т)|иа(т) - и(т)|Ит. (31)

t

Из (31) в силу леммы Гронуолла

||иа(Ь) - и(Ь)Н < еь(т^ЦЕа^Ь) - и(Ь) Н. (32)

Рассмотрим функцию va(t), удовлетворяющую условиям Н?;а (Ь)

= -Аиа(Ь) + /^а(Ь)), Ь € (0,Т); vа(0) = ЕаЗД. (33)

По лемме 1

1К(Ь)|| < еьт||е-^(Еа - Е)иоН. (34)

Оценим норму разности проекций:

t

Н(Е - ЕаК(Ь) - (Е - Еа)и(Ь)Н < ^ е-а^-т)К(т) - и(т)Нйт. (35)

о

Так как в силу леммы 1 |^а(Ь) - и(Ь)Н < е^т||е-^Еаио||, то из (35) следует

Н(Е - ЕаК(Ь) - (Е - Еа)и(Ь)Н < ^е-а^. (36)

Из неравенств (36) и (34) следует оценка

Н(Е - Еа)и(Ь)Н < Н(Е - Eа)vа(Ь) - (Е - Еа)и(Ь)Н + Н(Е - ЕаК(Ь)|| < (1+ ^)е^ге-а (37)

Из неравенств (37) и (32) следует, что

Д:(а) = 8ир(Нр - Еар|| : р € М} < Д:(а) < ге^°.

Выберем зависимость а = а* (5) из условия

еа(т-^5 = Ге-а^ (38)

Из (38) следует , что оценка погрешности метода проекционной регуляризации на множестве М с выбором параметра регуляризации из условия (38) имеет вид

Дм (а* (5), 5) < 2геьт (1 + ЬТ)е-а*°. (39)

с

Из леммы 2 и неравенства (39) следует

Теорема 2. Метод проекционной регуляризации оптимален по порядку на множестве M.

Другие методы решения нелинейных обратных задач рассмотрены, например, в статьях

[9, 10].

Работа проводилась при финансовой поддержке РФФИ, проект 12-01-00117

Литература

1. Менихес, Л.Д. О регуляризуемости некоторых классов отображений, обратных к интегральным операторам / Л.Д. Менихес // Математические заметки. - 1999. - Т. 65, № 2.

- С. 222-229.

2. Менихес, Л.Д. Об одном достаточном условии регуляризуемости лмнейных обратных задач / Л.Д. Менихес // Математические заметки. - 2007. - Т. 82, № 2. - С. 242-246.

3. Иванов, В.К. Об оценке погрешности при решении линейных некорректных задач / В.К. Иванов, Т.И. Королюк // Журнал вычислительной математики и математической физики. - 1969. - Т. 9, № 1. - C. 30-41.

4. Танана, В.П. Methods for solution of nonlinear operator equations/ В.П. Танана. - Utrecht: VSP, 1997.

5. Васин, В.В. Некорректные задачи с априорной информацией /В.В. Васин, А.Л. Агеев.

- Екатеринбург: Наука, 1993.

6. Хенри, Д. Геометрическая теория полулинейных параболических уравнений. /Д. Хенри.

- М.: Мир, 1985.

7. Табаринцева, Е.В. Об оценке модуля непрерывности одной нелинейной обратной задачи/ Е.В. Табаринцева // Труды ИММ УрО РАН. - 2013. - Т. 19, № 1. - C. 253-257.

8. Танана, В.П. Об одном подходе к приближению разрывного решения некорректно поставленной задачи / В.П. Танана, Е.В. Табаринцева // Сибирский журнал индустриальной математики. - 2005. - Т. 8, № 1(21). - С. 129-142.

9. Табаринцева, Е.В. Об оценке погрешности метода квазиобращения при решении задачи Коши для полулинейного дифференциального уравнения / Е.В. Табаринцева // Сибирский журнал вычислительной математики. - 2005. - Т. 8, № 3. - C. 259-271.

10. Танана, В.П. О методе приближения кусочно-непрерывных решений нелинейных обратных задач/ В.П. Танана, Е.В. Табаринцева // Сибирский журнал вычислительной математики. - 2007. - Т. 10, № 2. - C. 221-228.

Елена Владимировна Табаринцева, кандидат физико-математических наук, доцент, кафедра «Общая математика>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), eltab@rambler.ru.

Bulletin of the South Ural State University. Series «Mathematical Modelling, Programming & Computer Software>,

2013, vol. 6, no. 3, pp. 85-94.

MSC 47J06

On Error Estimate of an Approximate Method to Solve an Inverse Problem for a Semi-Linear Differential Equation

E.V. Tabarintseva, South Ural State University, Chelyabinsk, Russian Federation, eltab@rambler.ru

An inverse problem for a semi-linear differential-operator equation in a Hilbert space is considered in the paper. The projection regularization method is used to get a stable approximate solution to the nonlinear ill-posed problem. The regularization parameter is chosen referring to the Lavrentev scheme. A sharp error estimate of the considered method on a correctness class defined by means of a nonlinear operator is obtained. The value of the continuity module for the corresponding problem on the correctness classes plays an important role in the investigation of the methods for the solution of ill-posed problems in order to state their optimality. The linear operators are used, as a rule, to define the correctness classes. The two-sided estimate of the continuity module for the nonlinear inverse problem on the correctness class defined by a nonlinear operator is obtained in the present work. The obtained estimate of the continuity module is used to prove the order-optimality of the projection regularization method on the analyzed correctness class.

Keywords: inverse problem; a method of approximate solution; continuity module; error estimate; semilinear rquation.

References

1. Menikhes L.D. Regularizability of Some Classes of Mappings that are Inverses of Integral Operators. Mathematical Notes, 1999, vol. 65, no. 1-2, pp. 181-187.

2. Menikhes L.D. On a Sufficient Condition for Regularizability of Linear Inverse Problems. Mathematical Notes, 2007, vol. 82, no. 1-2, pp. 242-246.

3. Ivanov V.K., Korolyuk T.I. Error Estimates for Solutions of Incorrectly Posed Linear Problems. USSR Computational Mathematics and Mathematical Physics, 1969, vol. 9, no. 1, pp. 35-49.

4. Tanana V.P. Methods for Solution of Nonlinear Operator Equations. Utrecht, VSP, 1997.

5. Vasin V.V., Ageev A.L. Inverse and Ill-Posed Problems with a Priori Information. Utrecht, VSP, 1995.

6. Henry D. Geometric Theory of Semi-Linear Parabolic Equations. Berlin, Springer, 1981.

7. Tabarintseva E.V. On an Estimate for the Modulus of Continuity of a Nonlinear Inverse Problem. Trudy Instituta Matematiki i Mekhaniki UrO RAN [Proceedings of the Institute of Mathematics and Mechanics], 2013, vol. 19, no. 1, pp. 253-257. (in Russian)

8. Tanana V.P., Tabarintseva E.V. On an Approximation Method of a Discontinuous Solution of an Ill-Posed Problem. Sibsrskiy zhurnal industrial’noy matematiki [Journal of Applied and Industrial Mathematics], 2005, vol. 8, no. 1, pp. 129-142. (in Russian)

9. Tabarintseva E.V. On Error Estimation for the Quasi-Inversion Method for Solving a Semi-Linear Ill-Posed Problem. Sibirskiy zhurnal vychislitel’noy matematiki [Numerical Analysis and Applications], 2005, vol. 8, no. 3, pp. 259-271. (in Russian)

10. Tanana V.P., Tabarintseva E.V. On a method to approximate discontinuous solutions of nonlinear inverse problems. Sibirskiy zhurnal vychislitel’noy matematiki [Numerical Analysis and Applications], 2007, vol. 10, no. 2, pp. 221-228. (in Russian)

Поступила в редакцию 28 апреля 2013 г.