CC BY
51
Вестник СамГУ — Естественнонаучная серия. 2013. № 9/1(110)
МАТЕМАТИКА
УДК 519.6
ЧИСЛЕННЫЙ МЕТОД РЕШЕНИЯ ОБРАТНЫХ СПЕКТРАЛЬНЫХ ЗАДАЧ, ПОРОЖДЕННЫХ ВОЗМУЩЕННЫМИ САМОСОПРЯЖЕННЫМИ ОПЕРАТОРАМИ
© 2013 С.И. Кадченко1
Разработан новый численный метод решения обратных спектральных задач, порожденных возмущенными самосопряженными операторами по их спектральным характеристикам. Метод был апробирован на задачах для возмущенного оператора Штурма — Лиувилля. Результаты многочисленных расчетов показали его вычислительную эффективность. Найдены простые алгебраические формулы для нахождения собственных значений дискретных операторов. При этом вычисление собственных значений возмущенного самосопряженного оператора можно начинать с любого их номера независимо от того, известны ли собственные значения с предыдущими номерами или нет.
Ключевые слова: дискретные операторы, самосопряженные операторы, собственные значения, собственные функции, обратные спектральные задачи, интегральное уравнение Фредгольма первого рода, некорректно поставленные задачи, задача Штурма — Лиувилля.
Введение
Интерес к обратным задачам все время возрастает в связи с широкой областью их приложения, например, к задачам сейсморазведки, идентификация композитных материалов, проблемам неразрушающего контроля, нелинейных эволюционных уравнений математической физики и др. В статье разработан численный метод решения обратных спектральных задач, порожденных возмущенными самосопряженными операторами по их спектральным характеристикам. Получено интегральные уравнение Фредгольма первого рода, позволяющие восстановить значения возмущающего оператора в узловых точках дискретизации. На основе метода Бубнова - Галеркина найдены простые формулы для вычисления приближенных собственных значений возмущенных самосопряженных операторов.
1 Кадченко Сергей Иванович (kadchenko@masu.ru), кафедра уравнений математической физики Южно-Уральского государственного университета (Национальный исследовательский университет), 454080, Российская Федерация, г. Челябинск, пр. Ленина, 76.
1. Теоретические положения метода
Рассмотрим дискретный полуограниченный снизу оператор Т и ограниченный оператор Р, заданные в сепарабельном гильбертовом пространстве Н. Допустим, что известны собственные значения {^к}&=1 и ортонормированные собственные функции {шк оператора Т, которые занумерованы в порядке невозрастания собственных значений Цк по величине с учетом алгебраической кратности. Пусть {вк}к=1 - собственные значения оператора Т+Р, занумерованные в порядке невозрастания их действительных частей с учетом алгебраической кратности. Для нахождения собственных значений спектральной задачи
(T + = вФ (1)
применим метод Бубнова - Галеркина.
Пусть система функций {wkявляется ортонормированным базисом пространства H. Следуя методу Бубнова - Галеркина, будем искать решение спектральной задачи (1) в виде
п
Фп = ak(п)шк. (2)
k=i
Подставляя (2) в уравнение (1), получим
п п
ak(n)(T + P )wk = ¡3(n)^2 ak (n)wk.
k=i k=i
Здесь ¡3(n) - n-е приближения по Бубнову - Галеркину к соответствующим собственным числам {ek}^=1 оператора T + P. Так как Twk = jkwk, то
пп
У^ ak (n)(jk + P)uk = ak (n)uk.
k=i k=i Коэффициенты {ak(n)}n=i определяются из требования, чтобы левая часть последнего уравнения была ортогональна к функциям {wi }П=1. В результате получаем систему линейных уравнений для нахождения коэффициентов {ak(n)}n=i
п
y^ak(n)^\fi(n) - jk]Ski - (Pwk,wi)} = 0, l = 1,n.
k=i
Приравняв определитель этой системы к нулю, приходим к уравнению
det \\f3(n)E - A\\ = 0,
определяющему приближенные значения первых собственных значений {fík(n)}rn=i оператора T + P. Здесь E - единичная матрица размера n х n, A = \\aki\ni=i, aki = jk Ski + (Pwk ,wi), Ski - символ Кронекера.
Известно [5], что для собственных значений {¡3k(n)}n=i матрицы A справедливы равенства
п
Y,Pi(n)= SpAp, p =1m (3)
k=i
где SpAp - след p-й степени матрицы A.
След p-й степени матрицы A вычисляется по формуле
п p
SpAp = J2 Uajsjr. (4)
p=
j1 j2 jp=i S = i
Здесь г = { * + 1' 5 = р
у 1, в = р.
Используя (3) и (4) при р = 1, получим
п п
У1Дк (п) = акк .
к=1 к = 1
Вводя обозначения £к(п) = Дк — Дк(п), найдем
пп
Дк = ^2\акк + £к (п)]. (5)
к=1 к=1
Запишем уравнение (5) для п — 1
п-1 п-1
^2Дк = ^2\акк + £к (п — 1)]. (6)
к=1 к=1
Вычитая (5) из (6), имеем
Дк = ^к + (Р^к)+ Пк(п), Ук е М, (7)
где
п-1
Пк (п) = £к (п) — [Дк (п) — Дк (п — 1)].
к = 1
Если в (7) подставить £к(п) = Дк — Дк (п), то получим формулы для нахождения собственных значений оператора Т + Р, удобные для численных расчетов
п-1
Дк(п) = Цк + (Р^к,Шк) — ^[Дк(п) — Дк(п — 1)], к = тп. (8)
к=1
Теорема 1. Пусть Т - дискретный полуограниченный снизу оператор, Р -ограниченный оператор, действующие в сепарабельном гильбертовом пространстве Н. Если оператор Т + Р положительно определенный в Н и система координатных функций {фп}^=1 является базисом Н, то метод Бубнова - Галеркина в применении к задаче об отыскании собственных значений спектральной задачи (1), построенный на этой системе функций, сходится. Доказательство. Запишем уравнение (1) в виде
(Т + Р — ХЕ )ф = (Д — Х)ф. (9)
Для дискретного оператора Т + Р существует резольвентный оператор Я\(Т + + Р) = (Т + Р — ХЕ)-1, который вполне непрерывный в Н [3]. Действуя слева на обе части уравнения (9) оператором И\(Т + Р), получим
Ф =(в — Х)Ех(Т + Р )ф.
На основании [4] метод Бубнова - Галеркина в применении к задаче об отыскании собственных значений (9), а следовательно и уравнения (1), сходится. □
Если операторы Т и Р удовлетворяют требованием теоремы 1, то метод Бубнова - Галеркина в применении к задаче для отыскания собственных значений спектральной задачи (1), построенный на системе функций {ф
п}п=1: сходится.
Поэтому с увеличением п точность вычисления собственных значений {Дк(п)}П=1 оператора Т + Р по формулам (8) возрастает.
Нахождение собственных значений возмущенного самосопряженного оператора по формулам (8) можно начинать с любого их номера независимо от того, известны ли собственные значения с предыдущими номерами или нет.
По формулам (8) особенно эффективно находить собственные значения возмущенного самосопряженного оператора с большим номером, что в случае применения метода Бубнова - Галеркина становится затруднительным.
2. Решение обратных спектральных задач
Пусть T — дискретный полуограниченный снизу оператор, а P — ограниченный оператор умножения на функцию p(s), заданные в сепарабельном гильбертовом пространстве L^ [a,,b], где [a,b] — отрезок изменения переменной s. Допустим, что известны собственные значения {¡k }&=1 и ортонормированные собственные функции {ukоператора T, которые занумерованы в порядке невозрастания собственных значений ¡k по величине с учетом их кратности. Предположим, что система функций {uk}|=1 образует базис пространства L2[a,b].
Рассмотрим интегральное уравнение Фредгольма первого рода
Ap = Í K(x,s)p(s)ds = f (x), c < x < d. (10)
•J a
Функции K(x, s) и f(x) такие, что
K(xk,s) = Jk(s), c < Xk < d, f (xk) = ßk - ik.
Пусть ядро интегрального уравнения (10) K(x, s) непрерывно и замкнуто в прямоугольнике П= [a,b] х [c,d], функции p(s) G W^a,^] и f(x) G L2[c,d].
Допустим, что в узловых точках xk отрезка [a, b] вместо точных значений функции f (xk) правой части уравнения (10) известны ее приближенные значения
n — 1
f(xk) = ßk(n) - ik + ^[ßk(n) - ßk(n - 1)]. k=1
Задача решения интегрального уравнения Фредгольма первого рода (10) является некорректно поставленной. Ее приближенные решения могут быть найдены с помощью метода регуляризации Н.А. Тихонова [6].
Численное решение интегрального уравнения (10) определяет приближенные значения {p>i(s)}l=1 функции p(s) в узловых точках s¿, i = 1,I, a = s1 < s2 < ... ... < si = b. Число узловых точек I можно выбрать достаточно большим, чтобы получить необходимую точность при интерполяции функции p(s).
3. Численный эксперимент
Проиллюстрируем метод нахождения значений возмущающего оператора p(s) в узлах дискретизации, на спектральной задаче Штурма — Лиувилля
-u'' + p(s) u = ß u, a < s < b;
cosa u'(a) + sina u(a) = 0; (11)
cosj u'(b) + sinY u(b) = 0, a,Y G R.
Рассмотрим оператор Tw = —ш'', причем функция ш удовлетворяет граничным условиям задачи (11). Нетрудно показать, что оператор T самосопряженный, и его собственные значения являются решением трансцендентного уравнения
[sin a sin(y/fi,a) + yffi, cos a cos(y/r^a)]x x[sin y cos^y^b) — yffi, cos y sin(y^b)]+ + cos a sin^y^a) — sin a cos(y/r^a)]x x[sin y sin(y/fib) + yffi, cos y cos^y^b)] = 0, а собственные функции имеют вид:
wk(s) = Ck{[sin asin(yffika) + cos acos(^ka)\ cos(y/^ks)+ + cos a sin(^f^ka) — sin a cos(yf^ka)] sin(v/^ks)}, к = 1, ж.
Здесь {^k}'kLi — собственные значения оператора T. Постоянные Ck находятся из условия нормировки.
Для проверки полученных формул (8) сравним величины собственных значений спектральной задачи Штурма - Лиувилля (11), найденные по формулам (8), которые обозначим fik(n) и методом Бубнова - Галеркина, которые обозначим /% (n).
Таблица 1
k А(31) Â(31) |А(31) - вк(31)|
1 13,66134 - 2, 55439Í 13, 60139 - 2, 57604Í 0 06375
2 42,92974 - 2, 63944Í 42, 95887 - 2, 62907Í 0 02997
3 92,21746 - 2, 65564Í 92,22718 - 2, 65110Í 0 01035
4 161,78283 - 2, 66092Í 161, 28901 - 2, 65806Í 0 00554
5 250,10059 - 2, 66235Í 250,10377 - 2, 66122Í 0 00338
6 358,66101 - 2, 66367Í 358,66318 - 2, 66291Í 0 00230
7 486,96272 - 2, 66447Í 486,96429 - 2, 66392Í 0 00166
8 635,00474 - 2, 66498Í 635,00593 - 2, 66457Í 0 00126
9 802,78662 - 2, 66534Í 802,78755 - 2, 66501Í 0 00099
10 990,30810 - 2, 66559Í 990,30885 - 2, 66533Í 0 00080
11 1197,56905 - 2,66578Í 1197, 56967 - 2, 66556Í 0 00066
20 3951,18622 - 2,66640Í 3951, 18641 - 2, 66633Í 0 00020
21 4355, 83991 - 2,66642Í 4355, 84007 - 2, 66636Í 0 00018
22 4780, 23281 - 2,66644Í 4780, 23296 - 2, 66639Í 0 00016
23 5224, 36493 - 2,66646Í 5224, 36507 - 2, 66641Í 0 00015
24 5688, 23627 - 2,66648Í 5688, 23640 - 2,66644Í 0 00014
25 6171, 84683 - 2,66649Í 6171, 84695 - 2,66645Í 0 00013
26 6675,19660 - 2, 66651Í 6675, 19671 - 2, 66647Í 0 00012
27 7198,28559 - 2, 66652Í 7198,28569 - 2,66648Í 0 00011
28 7741, 11359 - 2, 66653Í 7741, 11389 - 2,66650Í 0 00010
29 8303, 68121 - 2, 66754Í 8303,68130 - 2,66650Í 0 00011
30 8885, 98783 - 2, 66655Í 8885, 98793 - 2, 66651Í 0 00010
31 9488, 03367 - 2, 66656Í 9488, 03607 - 2,66573Í 0 00254
В табл. 1 приведен пример численных расчетов собственных значений задачи (11) при a = 0, b = 1, а = Pi/5, y = Pi/7, p(s) = s2 + 5s + 1 + (s2 — 3)i, к = ГЖ Расчеты проводились при условии, что 3к (n) — 3 (n — 1) =0 для k,n = 1, 31.
Из табл. 1 видно, что при увеличении порядкового номера к приближенных собственных значений соответствующие величины 13к (31) — Рк (31) | уменьшаются. Исключение составляет последняя строка (к = 31). Если найти Рз\(61), то
\,3з!(31) — З3!(61)\ = 0,00009. Следовательно, увеличение значения |,3з!(31) — /Зз!(31)| связано с погрешностями возникающими при применении метода Бубнова — Га-леркина.
Проведенные многочисленные расчеты для различных значений параметров а, Ь, с, ¿,а, /, р(в) спектральной задачи (11) показали высокую точность и вычислительную эффективность полученных формул (8).
Изменим правую часть уравнения (11) и восстановим приближенные значения функции р(в) в узловых точках {вн}П=!, п = 31. В табл. 2 приведен пример расчетов при а = 0, Ь =1, а = Рг/Ъ, 7 = Рг/7, /(хн) = Рн — Ин — 2 — 3г, к = 1, 31.
Таблица 2
к Эк р(эк) С (Эк )
1 0 00 -0, 325 - - 6, 291г 0 000017
2 0 03 -0, 315 - - 6, 283г 0 000082
3 0 07 -0, 295 - - 6, 288г 0 000020
4 0 10 -0, 269 - 6, 265г 0 000057
5 0 13 -0, 223 - 6, 274г 0 000021
6 0 17 -0,189 - 6, 238г 0 000085
7 0 20 -0,077- 6, 203г 0 000062
8 0 23 0, 063 - 6, 160г 0 000072
9 0 27 0, 227 - 6, 111г 0 000055
10 0 30 0, 413 - 6, 055г 0 000047
11 0 33 0, 617 - 5, 996г 0 000094
21 0 67 2, 949 - 5, 344г 0 000086
22 0 70 3, 150 - 5, 290г 0 000022
23 0 73 3, 400 - 5, 240г 0 000064
24 0 77 3,515 - 5, 194г 0 000081
25 0 80 3, 677 - 5, 153г 0 000053
26 0 83 3, 825 - 5, 116г 0 000003
27 0 87 3,960 - 5, 083г 0 000083
28 0 90 4, 082 - 5, 054г 0 000018
29 0 93 4, 192 - 5, 028г 0 000079
30 0 97 4, 292 - 5, 005г 0 000008
31 1 00 4,364 - 4, 991г 0 000090
Здесь р(вк) — приближенные значения функции р(в) в узловых точках в к.
Величины £н = \/(хн) — а К(хн, в)р(8)с18\ определяют поточечную абсолютную погрешность решения. Невязка приближенного решения р(вн), равна \\Ар — /1| = = 0,000094. Параметр регуляризации а при численном решении интегрального уравнения Фредгольма первого рода (11) методом регуляризации Тихонова вычислялся с помощью метода невязки. В нашем случае а = 0,0000075.
Значения поточечной абсолютной погрешности £н и невязки || АрГ— /1| позволяют сделать вывод о хорошей точности нахождения приближенных значений функции р(в) в узловых точках дискретизации в к.
Заключение
В работе разработан численный метод решения обратных спектральных задач для возмущенных самосопряженных операторов. Найдены простые формулы для вычисления собственных значений дискретных операторов. Результаты многочисленных расчетов показали вычислительную эффективность метода.
Литература
[1] Кадченко С.И. Метод регуляризованных следов // Вестник Южно-Урал. гос. ун-та. Сер.: Математическое моделирование и программирование. 2009. № 37(170). Вып. 4. С. 4-23.
[2] Кадченко С.И., Рязанова Л.С. Численный метод нахождения собственных значений дискретных полуограниченных снизу операторов // Вестник Южно-Урал. гос. ун-та. Сер.: Математическое моделирование и программирование. 2011. № 17(234). Вып. 8. С. 46-51.
[3] Садовничий В.А. Теория операторов. М.: Дрофа, 2004.
[4] Михлин С. Г. Вариационные методы в математической физике. М.: Наука, 1970. 510 с.
[5] Демидович Б.П., Марон И.А. Основы вычислительной математики. М.: Наука, 1966. 659 с.
[6] Васильева А.Б., Тихонов Н.А. Интегральные уравнения. М.: Изд-во МГУ, 1989. 156 с.
Поступила в редакцию 15/IX/2013; в окончательном варианте — 15/IX/2013.
NUMERICAL METHOD FOR SOLVING INVERSE SPECTRAL PROBLEMS GENERATED BY PERTURBED SELF-ADJOINT OPERATORS
© 2013 S.I. Kadchenko2
A new numerical method for solving inverse spectral problems generated by perturbed self-adjoint operators from their spectral characteristics is developed. The method was tested on the problems for perturbed Sturm - Liouville operator. The results of numerous calculations have shown its computational efficiency. The simple algebraic formulas for finding the eigenvalues of discrete operators was found. At that the calculation of eigenvalues of perturbed self-adjoint operator can start from any number, no matter known the eigenvalues from previous numbers or not.
Key words: discrete operators, self-adjoint operators, eigenvalues, eigenfunctions, inverse spectral problems, Fredholm integral equation of the first kind, improperly posed problems, Sturm - Liouville problem.
Paper received 15/IX/2013. Paper accepted 15/IX/2013.
2Kadchenko Sergey Ivanovich (kadchenko@masu.ru), the Dept. of Mathematical Physics Equations, South-Ural State University, Chelyabinsk, 454080, Russian Federation.
