КРАТКИЕ СООБЩЕНИЯ
УДК 519.642.8
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ СПЕКТРАЛЬНОЙ ЗАДАЧИ ОБ ЭЛЕКТРИЧЕСКИХ КОЛЕБАНИЯХ В ПРОТЯЖЕННОЙ ЛИНИИ МЕТОДОМ РЕГУЛЯРИЗОВАННЫХ СЛЕДОВ
С.Н. Какушкин
Работа посвящена описанию нового численного метода вычисления значений собственных функций возмущенных самосопряженных операторов, основанного на методе регуляризованных следов. Построена математическая модель вычисления значений собственных функций спектральной задачи об электрических колебаниях в протяженной линии. Разработанные алгоритмы позволяют вычислять значения собственной функции возмущенного оператора незавсимо от того, известны предыдущие значения собственных функции или нет. Получены оценки остатков сумм функциональных рядов «взвешенных» поправок теории возмущений возмущенных самосопряженных операторов, и доказана их сходимость. Для вычислительной реализации метода найдены эффективные алгоритмы нахождения «взвешенных» поправок теории возмущений. Проведенные численные эксперименты вычисления значений собственных функций задачи об электрических колебаниях в протяженной линии показывают, что метод хорошо согласуется с другими известными методами А.Н. Крылова и А.М. Данилевского. Метод регуляризованных следов показал свою надежность и высокую эффективность.
Ключевые слова: задача Штурма - Лиувилля, собственные числа, собственные функции, теория возмущений, метод регуляризованных следов.
В последнее время все большее значение приобретают вопросы математического моделирования нахождения собственных чисел и собственных функций возмущенных самосопряженных операторов [1, 2]. Рассмотрим задачу об электрических колебаниях в протяженной линии, собственные колебания которой описываются собственными значениями задачи Штурма - Лиувилля [3]:
-у"(х, у) + д(х)у(х, у) = уу(х, у), (1)
у'(0,у) - (рг + Р2у)у(0,у) = 0,
у'(1,у) + (рз + Р4у)у(1, у) = 0,
С Ьо С 1о „ Т , ,
где рг = —, р2 = ——, рз = , р4 = ——, С и Ь - коэффициенты емкости и самоиндукции,
Со Ь С0 Ь
рассчитанные на единицу длины провода. Исходя из физического смысла задачи, коэффициенты рг, г = 1, 4 должны быть вещественными. Физически граничные условия означают, что левый конец провода заземлен через сосредоточенную самоиндукцию Ьо и емкость Со, а правый - через сосредоточенную самоиндукцию Ьо и емкость Со. Предполагается, что сосредоточенная самоиндукция и емкость соединены последовательно. Для определенности будем считать, что длина провода равна единице.
В работах [4, 5] разработан неитерационный метод регуляризованных следов (РС) вычисления значений собственных функций возмущенных самосопряженных операторов. Следуя обозначениям метода РС, перепишем уравнение (1) в виде:
(Т + Р)п = уп, п(х) € Вт, (2)
<12и ^п(х)
где Бт = {и | и Є С2(0,1) П С1 (0, 1),-^ Є ^2(0,1), их 0 - (р! + р2у)и(х)
-и(х)
ж=0
-X
Х=1
+ (рз + р4у)и(х) =0 к Т = —А - оператор Лапласа, Р = д(х) - потенци-
Х=1 J
ал, х Є (0,1). Собственные числа Ап невозмущенного оператора Т являются корнями уравнения tg /А = /А———, а соответствующие им собственные функции -и(х) имеют вид: А — я1я2
^п(х) = === (^1 8Ш \/АПх + л/ХПес^^АПх), где ^1 = Р1 + Р2А, ^2 = рз + Р4А. Обозначим
у/ Ап + ^1
через По количество всех неравных друг другу Ап, лежащих внутри окружности Тп0 радиуса 1Апо+1 + Апо 1
рпо = -------------------------------------------------------------------^- с центром в начале координат комплексной плоскости.
Теорема 1. Если Т - дискретный полуограниченный снизу оператор, а Р - ограниченный оператор, действующие в сепарабельном гильбертовом пространстве Н, с областью опре-
2| |Р ||
деления в Б, и для всех натуральных п > По выполняются неравенства дп = —------------- —- <
1 Ап+1 Ап1
1, то значение произведения собственной функции ип(х) на ее сопряженную ип(у), при любых значениях аргументов х, у Є Б, можно найти по формулам:
1 1
ип(х)йп(у) = — (Ап^п(х)^п(у) + ^[ак1}(п,х,у) — а^(п — 1,х,у)]) + е(1)(п, х, у), (3)
уп ' - . '
к=1
где для е^1)(п, х,у) справедливы оценки
|е(1)(п,х,у)| < 2||Р|1 С^пР^-----------, V Є N п = 1,по.
Уп 1 — д
/ 1 \ 2 ____________________________________
Здесь £п = вир ^ у------------—г , |^(х)| < Со V* = 1, ж, д = тахдп.
Лі ^ і=1 1 Ап — АгК п^1
і=п
«:Взвешенные> поправки теории возмущений а^^п, х,у), к = 1, ж, входящие в формулы (3), можно найти, используя следующую теорему.
Теорема 2. Если Т - дискретный полуограниченный снизу оператор, а Р - ограниченный оператор, действующие в сепарабельном гильбертовом пространстве Н, и для всех п > по выполняются неравенства дп < 1, то „взвешенные11, поправки теории возмущений акР)(по,х,у) для любых натуральных к, р и по можно найти по формулам:
по го к
акр) (по, х, у) = — Е £ VI(х) ^7к+1 (у)г*Т)(п,І1,-,ік+1) П "К?т,?т + 1 , (4)
п= 1 Л,...,Л;+1 = 1 т=1
" 0, = п, т = 1, к + 1;
где г(р)(п, ;'ь..., ік+1) = <
)Ш
кг і™ ^ 1 =к + 1;
Л——Лп
(і-1)! <ілг-1 (к—+ ), 0 <1 <к;
Л—Лп 4 П (Л-Літ )у
т = 1
Уі,; = (Р^і, ) - скалярное произведение; I- число совпадений = п, т = 1, к + 1.
Был проведен вычислительный эксперимент по нахождению значений собственных функций спектральной задачи (2). Значения собственных функций ип вычислялись по формулам (3). Суммы функциональных рядов Рэлея - Шредингера приближались четырьмя «:взвешенными> поправками теории возмущений по формулам (4). Значения собственных функций ип спектральной задачи (2), вычисленные методом РС, сравнивались со значениями, найденными методом А.Н. Крылова. Результаты вычисления значений тринадцатой собственной функции приведены в таблице. Первые обозначены - йіз(х), а вторые - «із(х).
Таблица
Значения тринадцатой собственной функции задачи Штурма - Лиувилля (1), вычисленные при р1 = 1, р2 = 0, р3 = 1, р4 = 0 и потенциалом д(х) = х2
і Хі «із(хі) «із (хі) |«із(хі) - «із(хі)| |«із(хі) - «із(хі)| %
і«із(хі)| %
1 0,095238 -1, 284995 -1, 285037 0,000041 0,003226
2 0,142857 0, 859796 0, 859596 0,000199 0,023195
3 0,190476 0, 898039 0, 897731 0, 000308 0,034311
4 0, 238095 -1, 264164 -1, 264311 0,000146 0,011609
5 0, 285714 -0,329073 -0,329143 0, 000070 0,021281
6 0, 333333 1,412499 1,412254 0, 000245 0,017351
7 0, 380952 -0,306695 -0,306993 0, 000298 0,097239
8 0,428571 -1, 274733 -1, 274854 0,000121 0,009478
9 0,476190 0, 880529 0, 880433 0, 000095 0,010857
10 0, 523809 0, 878726 0, 878454 0, 000272 0,031011
11 0, 571428 -1, 276228 -1, 276495 0, 000266 0,020867
12 0, 619047 -0,304623 -0,304712 0, 000089 0,029144
13 0, 666666 1,413653 1,413536 0,000117 0,008295
14 0, 714285 -0,331369 -0,331649 0, 000280 0,084491
15 0, 761904 -1, 264955 -1, 265172 0,000216 0,017116
16 0, 809523 0, 900513 0, 900457 0, 000055 0,006215
17 0, 857142 0, 860215 0, 860080 0,000135 0,015679
18 0, 904761 -1, 287601 -1, 287871 0, 000269 0,020894
19 0, 952380 -0,281363 -0,281518 0,000155 0,055095
Из таблицы видно, что результаты вычисления значений собственных функций методом РС хорошо согласуются с результатами, полученными методом А.Н. Крылова. При этом метод РС показал надежность и высокую эффективность.
Литература
1. Свиридюк, Г.А. Разрешимость задачи термоконвекции вязкоупругой несжимаемой жидкости / Г.А. Свиридюк // Известия высших учебных заведений. Математика. - 1990. -№12. - С. 65-70.
2. Новый метод вычисления первых собственных чисел спектральной задачи гидродинамической теории устойчивости течения вязкой жидкости между двумя вращающимися цилиндрами / В.В. Дубровский, С.И. Кадченко, В.Ф. Кравченко, В.А. Садовничий // Доклады Академии наук. - 2001. - Т. 381, № 3. - С. 320-324.
3. Валеев, Н.Ф. О задаче определения параметров граничных условий оператора Штурма-Лиувилля по спектру / Н.Ф. Валеев, С.А. Рабцевич, Э.Р. Нугуманов // Вестник СамГУ. - 2009. - №6 (72). - С. 12-20.
4. Кадченко, С.И. Вычисление значений собственных функций дискретных полуограни-ченных снизу операторов методом регуляризованных следов / С.И. Кадченко, С.Н. Какушкин // Вестник СамГУ. - 2012. - №6 (97). - С. 13-21.
5. Кадченко, С.И. Алгоритм нахождения значений собственных функций возмущенных самосопряженных операторов методом регуляризованных следов / С.И. Кадченко, С.Н. Какушкин // Вестник ЮУрГУ. Серия: Математическое моделирование и программирование. - 2012. - №40 (299), вып. 14. - С. 71-76.
Сергей Николаевич Какушкин, аспирант, кафедра «Прикладная математика и вычислительная техника>, Магнитогорский государственный университет (г. Магнитогорск, Российская Федерация), kakushkin-sergei@mail.ru.
Bulletin of the South Ural State University. Series «Mathematical Modelling, Programming & Computer Software>,
2013, vol. 6, no. 3, pp. 125-129.
MSC 47A75
Mathematical Modelling of Finding the Values of Eigenfunctions for the Electrical Oscillations in the Extended Line Problem Using the Method of Regularized Traces
S.N. Kakushkin, Magnitogorsk State University, Magnitogorsk, Russian Federation, kakushkin-sergei@mail.ru
This paper describes a new numerical method for computing the values of the eigenfunctions of perturbed self-adjoint operators. The new method is based on the method of regularized traces. A mathematical model for calculating the eigenfunction values of the spectral problem concerning electrical oscillations in the extended line is developed. The elaborated algorithms make it possible to calculate the values of the eigenfunction of the perturbed operator whether the previous values are known or not. We’ve obtained the estimates of functional series residual sum«suspended> the corrections of the perturbation theory of perturbed self-adjoint operators, and proved their convergence. Effective algorithms for finding«suspended» perturbation theory corrections are discovered for the numerical implementation of the method. The numerical experiments on the calculation of the values of a problem on its own electrical oscillations in the extended lines show that the method is consistent with the other well-known methods of A.N.Krylov and A.M.Danilevsky. The method of regularized traces proved its reliability and high efficiency.
Keywords: Sturm - Liouville problem, eigenvalues, eigenfunctions, perturbation theory, the method of regularized traces.
References
1. Sviridyuk G.A. Solubility of the Thermal Convection of Viscoelastic Incompressible Fluid.
Soviet Mathematics (Izvestiya VUZ. Matematika), 1990, vol. 34, no. 12, pp. 80-86.
2. Dubrovckiy V.V., Kadchenko S.I., Kravchenko V.F., Sadovnichiy V.A. A New Method of Calculation of the First Eigenvalues of the Spectral Problem of Hydrodynamic Stability Theory Cramped Viscous Fluid Between Two rotating Cylinders. Doklady Akademii nauk [Doklady Mathematics], 2001, vol. 381, no. 3, pp. 320-324. (in Russian)
3. Valeev N.F., Rabtsevich S.A., Nugumanov E.R. The Problem of Determining the Parameters of the Boundary Conditions of the Sturm - Liouville Operator on the Spectrum. Vestnik Samarskogo Gosudarstvennogo Universiteta. Estestvenno-Nauchnaya Seriya, 2009, no. 6 (72), pp. 12-20. (in Russian)
4. Kadchenko S.I., Kakushkin S.N. The Calculation of the Values of Natural Functions of Discrete Semi-Bounded from Below by the Operators of Regularized Traces. Vestnik Samarskogo Gosudarstvennogo Universiteta. Estestvenno-Nauchnaya Seriya, 2012, no. 6 (97), pp. 13-21. (in Russian)
5. Kadchenko S.I., Kakushkin S.N. The Algorithm for Finding the Values of the Eigenfunctions of Self-Adjoint Operators Perturbed by Regularized Traces. Bulletin of the South Ural State University. Series «:Mathematical Modelling, Programming & Computer Softwares, 2012, no. 40 (299), issue 14, pp. 71-76. (in Russian)
Поступила в редакцию 9 июня 2013 г.