CC BY
112
45. Современные численные методы решения обыкновенных дифференциальных уравнений: Пер. с англ. / под ред. Дж. Холла, Дж. Уатта. - M.: Мир, 1979. - 312 с.
46. Вержбицкий В.М. Численные методы (математический анализ и обыкновенные дифференциальные уравнения). - M.: Высшая школа, 2001. - 382 с.
47. Арушанян О.Б., Залеткин С.Ф. Численное решение обыкновенных дифференциальных уравнений на Фортране. - M.: Изд-во МГУ, 1990. - 336 с.
48. Gear C.W. The automatic integration of ordinary differential equations // Communications of the ACM. - 1971. - V. 14. - № 3. -P. 176-179.
49. Gear C.W. The Algorithm 407: DIFSUB for solution of ordinary differential equations [D2] // Communications ofthe ACM. - 1971. -March. - V. 14. - № 3. - P. 185-190.
50. Nordsieck А. On Numerical Integration of Ordinary Differential Equations // Mathematics of Computation. - 1962. - V. 16. -№ 77. - Jan. - P. 22-49.
51. Miranker W.L. The Computation Theory of Stiff Diffential Equations. Serie III. № 102. - Roma: IAC. Istitute per Le Applicazioni Del calcolo «Mauro Picone», 1975. - 120 p.
52. Miranker WL. Numerical method for stiff equations and singular perturbation problem. - Reidel: Reidel Publishing Company, 1981.- V. 5. - № 1. - 205 p.
53. Dahlquist G. Convergence and stability in the numerical integrator of ordinary diferential equations // Math. Scand. - 1956. - V. 4. -P. 33-53.
54. Dahlquist G. A special stability problem for linear multistep methods // BIT Numerilal Mathematics. - 1963. - V. 3. - № 1. - P. 27-43.
55. Самарский А.А., Гулин А.В. Численные методы. - M.: Наука, 1989. - 320 с.
56. Влах И., Сингхал К. Машинные методы анализа и проектирования электронных схем: Пер. с англ. - M.: Радио и связь, 1988. - 560 с.
57. Shampine L.F., Thompson S. Stiff systems // Scholarpedia. 2007. URL: www.scholarpedia.org/article/Stiff_systems (дата обращения: 14.04.2010).
58. Mott D., Oran E., B. L. A Quasi-Steady-State Solver for the Stiff Ordinary Differential Equations of Reaction Kinetics // Journal of Computational Physics. - 2000. - № 164. - Iss. 2. - P. 407-428.
59. Cash J. Efficient numerical methods for the solution of stiff initial-value problems and differential algebraic equations // Proc. Royal Society Lond. - 2003. - V. 459. - P. 797-815.
Поступила 15.04.2010 г.
УДК 517.9
КОЭФФИЦИЕНТНАЯ ОБРАТНАЯ ЗАДАЧА ДЛЯ ЛИНЕЙНОГО УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ ЧЕТВЕРТОГО ПОРЯДКА
Т.Д. Асылбеков, М.К. Чамашев
Ошский государственный университет, Кыргызстан E-mail: marat2712@rambler.ru
Методом интегральных уравнений и сжимающих отображений, доказана корректность коэффициентной обратной задачи для линейного гиперболического уравнения четвертого порядка с действительными трехкратными и простыми характеристиками.
Ключевые слова:
Коэффициентная обратная задача, функция Римана, система интегральных уравнений, система интегро-дифференциальных уравнений, корректность, векторная функция, метод сжимающих отображений.
Key words:
Coefficient inverse problem, Riemann function, the integral equation system, integro-differential equation system, correctness, vector function, contraction method.
1. Постановка задачи. Исследование коэффициентных обратных задач является одним из основных направлений теории обратных задач для дифференциальных уравнений. Обзор работ, посвященных таким задачам, приведены, например, в работах [1-3]. Коэффициентные обратные задачи для уравнений в частных производных третьего и четвертого порядков рассмотрены в работах [4-7].
В области D={(x,y):0<x<£, 0<у<Ц рассмотрим линейное уравнение в частных производных четвертого порядка
иххху + « (У>хху + А (УК + А2(У)% +
+71 (х, у)ых + у 2 (х)Ыу + 8 (х, у)и = / (х, у). (1)
Краевые задачи, задачи Коши и Дарбу для уравнения (1) изучены в работах [8, 9].
Пусть C+s(D) означает класс функций, имеющих непрерывные производные
Задача 1. Требуется найти коэффициент у2(х) и решение и(х,у)еС^ ^С+'ф) уравнения (1), удовлетворяющее начальным и краевым условиям:
и(х,0) = т(х), 0 < х < I, (2)
и(0,у) = Х1(y), их(0,у) = Хг(уХ ихх(0, у) = Жэ(уX 0 < у < К
и(х,К) = у(х), 0<х<I, (3)
где а, А, А, 7, 8, /, т (/=1,3), у - известные
функции, причем
«, А, Р2, Yi, 72> S’ f e C(D) z,ye C3[0,¿], Xt e C'[0, A],
(4)
т(0) = т'(0) = т"(0) =
У(0) у'(0) =Ь(у'(0) =Жэ(К). (5)
2. Сведение задачи 1 к системе интегральных уравнений. С этой целью воспользуемся представлением
х у
и (х, у) = Ф^ х, у) + | х,£) Е(£,ц) ёц, (6)
0 0
дающей решение уравнения иххху=Е(х,у), удовлетворяющее условиям (2) и (3), где
1 2 1 2
Ф1(ху) = т(х)+^х Хз(у)-2х 2з(0) +
+хХ2( у) - хХ 2(0) + Х1( у) -%х(0),
а &(х,^)=1/2(х-^)2 - функция Римана.
Перенося в уравнении (1) все младшие члены в правую часть равенства и используя формулу (6), получим интегро-дифференциальное уравнение
и (х, у) = Ф0( х, у) -
х у 0 0
+р2°^,ц)и^п + 71^,'П)% +72^)% + 8(Z,^n)u\dЦ, (7)
ху
где Ф0(х,у) = Ф1(х,у)- |]$(х,у)/(х,у)ёц.
0 0
Дифференцируя (7) по х и у, будем иметь их(х у) = Ф0х(х у) -х у а(^,Ц)ие„ + А(<^)ие +
-/(х,^) +А2(^,П)и|„ +71(%,Ц)и^ + dЦ, (8) 0 0 _+72(^)ил+8(^,П)и
<%, (9)
и y(х у) = ф(х у) -
х а(^, У)иЙу +А(^ у)ие +
J 3( х,^) +А2(^ У) % + Y 1(^ , У) u¡; +
0 _+Y2(^)uy +S(^, У)и
ихх (x, У) =Ф0хх (x, У) -
«(^,п)ие„ +А(<^)ий +
+ Р2(^,Л)и^п + Yi(^rH + ^ (10)
+Y2^K +8(%,л)и
х У
■Jdí¡
ихУ (X У) =Ф0хУ (X У) -
х «(^, У)иЙУ +PЛ, У)ие +
-JX( х,& +А2(^ У)и& + Y i(^, У)и^ +
0 +Y2KK У)и
d$, (11)
иххУ (x, У) =Ф0ххУ (x, У) -
У)иЙУ + А(^ У)и% + + P2(%, У)иЦу + Yl(£, У)и^ + + Y2KK + S(^, У)и
d£.
(12)
Далее интегрируя уравнение (1) по у в пределах от 0 до h и воспользуюсь краевыми условиями (2) и (3), найдем
72( х) = 70( х) -
1
р( х) 0 где
«(x,V)Uxxn+Pl( х,Л)ихх +
+fi2(x,'n)uxn+Yi( x,r)u +S( x,r)u
dr, (13)
Y0(x)=
9( x)
Jf (x,r)dr-?f”(x)
(p( x) =t( x)-щ( x).
Система уравнений (7-13) представляет собой замкнутую нелинейную систему интегро-диффе-ренциальных уравнений относительно
и(х, У), их(х, у), иУ(х, у), Uxx (х, у),
ихУ (x, У), ихху (x, У), Y2(х).
3. Решение системы уравнений методом сжимающих отображений [10]. Для этого систему уравнений запишем в виде операторного уравнения g = Ag,
где g= Cg1,g2,g3,g4,g5,g6,g7> - векторная функция, причем
gi(x) = u(x,y), g2(х) = их(XУ), g3(х) = иу(x, У),
g4 (х) = Uxx (X, УX g5(х) = иху (X, УХ g6(х) = ихху (X, У), gy(X, У) = gy(х) = и(X, У),
а оператор A определен на множестве функций geC(D) и имеет вид
х У
A1 g = g01 - Jd^J^(х,^) X
0 0
S(^,r) gi(^,r) +Yi(^,r) g2(^,r) +
Ptér) g4(^,r) + ^2(^,r) g 5(^,Г) + d a^r) gMr) + g№r) gy(^)
x У
A2 g = g 02 - Jd^jX (X,^) X
0 0
S(^r) g№r) +Yi(^r) g 2(^,r) +
+ Pl(£,r) g4(¿,r) + ^2(^,r) g 5 , Г ) + d^
_+«(^,r) g^r) + gs(^,r) gyOD
A3g = g03 - J^(^,r) X
0
s(%, у)gi(^, у) + Yi(^, У)g2(%, У) +
+ Pi(^ У) g 4 (^, У) + А2(^, У) g5 (^, У) + d^
_+a(Z, У)g6(Z, у) + gз(^, y)g7(^)
A4 g = g 04 -
gi(^) + Yl(4,l) g 2(4,1) +
+ ßl(4,1) g 4(4,1) + ß2(4,1) g 5(4,1) + +a(%,i) g6(^,i) + g3(£,l) g 7(^)
A5 g = go5 - Jöx (X,4) X
dl,
y)gl(4, y) + П(4, y)g2(4, y) +
+ ßl(|, y)gtf, y) + ß2(^, y)g5(|, y) +
+a(4, y^4, y) + g3(4, y)g7(4) _
A6 g = go6 -
S(4, y) g№, y) +Yl(Л, y) g 24, y) + +ßi(|, y) g 4(4, y) + ß2(4, y) gs(4, y) + _+a 4, y)g6 4, y) + g34, y) g7 4 )
d4,
d4,
A7g g07
l
ф( x)
S(4,1) gl(4,1) + Tl(4,l) g 2(4,1) +
+ ßl(4,1) g 4 (4 ,1) + ß2(4,1) g 5(4 ,1) + +a(4,1)g6(4,1)
dl,
где
I, I 2Ä,„
A7 g - g07 <---M1,
a
где L=5H+4M.
Из оценок (14) заключаем, что
II Ag - g (0)|| <в(£, h)M,
(14)
(15)
где в(£, h) = 2L
l
Пусть £ и h такие, что
в{£,h) < l. (16)
Тогда из (15) имеем неравенство lAg-g^KM. Это означает, что при выполнении условия (16) оператор A отображает шар B(g0,M) в себя, то есть
Vg е B(g 0,M): Ag е B(g 0,M).
Теперь покажем, что оператор A является сжимающим оператором. С этой целью рассмотрим два произвольных вектора g®, g(2)eB(g0,M), где
g<i)=(g1(i),g2(,'),g3(i),g4(i),g5(0),g6(‘),g7(i)), i=1,2.
Если учесть неравенство
g3l} g®
■g32) g72) < g
.(l)
■gfl g(7l>
g n g7l) ■
-g^l < 2£||gw - g1 то получим следующие оценки
|Alg® - Alg(2)|< thL|g(1) -g2> |A2g« - A2g(2)| < *2 hL ||g(l) - g^
A3g« - A3g<2>| < ^3L ||g« - g2>|
| A4g (1> - A4g(2) | < h ||g(1> - g(2)
IA5g(1) - A5g<2> < tL ||g« - g2>|
8Г =Ф0(x,y), g2 ’ =Ф0х(x,y), £з =Ф0у(x,y),
я40> = Ф0хх (^ У), g50 } = Ф0ху (x, У), gГ =Ф0хху (x, УХ g70> =70(x>,
компоненты вектора
Я(0)=(Я1(0),Я2(0),Яз(0),Я4(0),Я5(0),Я6(0),Я7(0)).
Норму g определим равенством
|Ы| = тах тах_|g¡ (х, у)|.
11 11 1</<7 (х,у)еВ 1 1
В силу свойств заданных функций_(4) ЗМ>0 V(x,y)еD:^0)||<М. В пространстве C(D) рассмотрим шар Б($,>),М)=^:^-£,>)1\<М}. Тогда VgеB(g<0),M): И<И+М<2М.
Пусть
тах{тах |а|,тах |Д|,тах |А2|,тах |у1|,тах |8} = Н, Vx е[0, £\: |ф(х) >а > 0.
Заметим также, что
Vx е[0,£\ [0,£\: |5(х,^|<£2, \3Х(х,%\< £.
Если VgеB(g<0),M), то AgеC(D) и, кроме того, справедливы неравенства
А g - gol| < 2£3кМЬ, |Аg -go21 < 2£2кМЬ,
Изg - go^ < 2£2кМЬ, |A4g -g0^| < 2£кМЬ,
\А5^ - go51 < 2£2КML, \Аg - g06 | < 2£КМ^
.(2)
\A6g(1) - A6g(2) < £L\\g(1) -i
(l)
.(2)
T(2)
Тогда из (17) получим неравенство
l
|| Ag(1) - Ag (2^| < ^ d(£, h)\\g(l) - g
(2)
(17)
(18)
Отсюда вытекает, что при любых £ и h, удовлетворяющих неравенство (16), в силу (18), оператор А осуществляет сжимающее отображение шара В(^0),М) в себя. Тогда по теореме Банаха в шаре В^0),М) существует единственная неподвижная точка отображения, т. е. существует единственное решение уравнения (15). Таким образом, доказано следующая теорема.
Теорема 1. Если выполняются условия (4), (5) и (16), то система уравнений (7-12) определяет в области D единственное непрерывное решение задачи 1, принадлежащее классу С^^).
СПИСОК ЛИТЕРАТУРЫ
1. Романов В.Г. Обратные задачи математической физики. - М.: Наука, 1984. - 264 с.
1. Кабанихин С.И. Обратные и некорректные задачи. - Новосибирск: Сибирское научное издательство, 2009. - 457 с.
2. Атаманов Э.Р, Мамаюсупов М.Ш. Неклассические задачи для псевдопараболических уравнений. - Фрунзе: Илим, 1990. -100 с.
3. Пакиров С. Обратные задачи для псевдопараболических уравнений // Исследования по интегро-дифференциальным уравнениям. - Бишкек: Илим, 1994. - Вып. 25. - С. 49-53.
4. Аблабеков Б.С. Обратные задачи для псевдопараболических уравнений. - Бишкек: Илим, 2001. - 180 с.
5. Матанова К.Б. Об одной обратной задаче для псевдопараболи-ческого уравнения // Научные труды ОшГУ. - Ош: Редакц.-изд. отдел «Билим», 1999. - Вып. 2. - С. 137-145.
6. Матанова К.Б. О существовании и единственности решения одной обратной задачи // Обратные и некорректные задачи математической физики: Докл. Междунар. конф., посвящ. 75-летию акад. М.М. Лаврентьева, 20-25 авг. 2007 г. - Новосибирск, 2007. - С. 137-145.
7. Джураев Т.Д., Сопуев А. К теории дифференциальных уравнений в частных производных четвертого порядка. - Ташкент: Фан, 2000. - 144 с.
8. Асылбеков Т.Д. Начально-краевые задачи для гиперболических уравнений четвертого порядка: Дис. ... канд. физ.-мат. наук. 01.01.02. - Бишкек, 2002. - 121 с.
9. Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. - М.: Наука, 1968. - 496 с.
Поступила 10.02.2010 г.
УДК 519.17
М-АЦИКЛИЧЕСКИЕ И ДРЕВОВИДНЫЕ ГИПЕРГРАФЫ
В.В. Быкова
Институт математики Сибирского федерального университета, г. Красноярск E-mail: bykvalen@mail.ru
Приводится характеризация двух классов гиперграфов: М-ациклических и древовидных. Установлена связь между этими классами: гиперграф М-ацикличен тогда и только тогда, когда двойственный к нему гиперграф является древовидным. Эта связь дает возможность объединить арсенал известных полиномиальных алгоритмов, позволяющих распознавать принадлежность гиперграфа к указанным классам и строить деревья соединений, деревья декомпозиций и деревья реализаций гиперграфа.
Ключевые слова:
Реализация гиперграфа деревом, дерево соединений гиперграфа, древовидная декомпозиция задач дискретной оптимизации, полиномиальные вычисления.
Key words:
Tree realization hypergraph, tree compounds hypergraph, tree decomposition of problem discrete optimization, polynomial computations.
Введение
Теория гиперграфов - один из разделов современной дискретной математики, имеющий большое прикладное значение. Данная теория располагает группой задач, не имеющих аналогов в теории графов и являющихся в большинстве своем КР-трудными. К этой группе задач, в частности, относятся задачи поиска для гиперграфа декомпозиций и реализаций с заданными свойствами [ 1 -4]. Необходимость практического решения этих задач в условиях, когда размерность задач велика, ставит вопрос о выделении классов гиперграфов, для которых они полиномиально разрешимы. Такие классы, например, образуют М-ациклические и древовидные гиперграфы.
Класс М-ациклических гиперграфов обладает рядом интересных свойств. Основное из них - это существование для связного М-ациклического гиперграфа дерева соединений. Понятие дерева соединений было введено в работе [4] и использовано там для построения монотонных планов соединений и программ полной редукции реляционных баз данных. Между тем, подобные деревья (их на-
зывают деревьями декомпозиции, деревьями клик) широко употребляют при решении задач дискретной оптимизации методами локальной декомпозиции и динамического программирования [2, 3]. Локальная оптимизация, осуществляемая в соответствии с деревом декомпозиции, дает возможность решать многие КР-трудные задачи дискретного программирования с разреженной матрицей ограничений за полиномиальное время [3].
Древовидный гиперграф - гиперграф, для которого существует реализация деревом. Древовидная реализация относится к минимальным (по числу ребер) реализациям гиперграфа. Необходимость построения минимальных реализаций появляется, в частности, при проектировании интегральных и вакуумных схем [1]. Доказано, что в общем случае задача поиска минимальной реализации является КР-трудной [5], но в классе древовидных гиперграфов она может быть решена за полиномиальное время [6-8]. Для древовидных гиперграфов полиномиально разрешимыми являются также задачи распознавания бихроматичности гиперграфа и реализуемости гиперграфа на пло-
