2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки


Правила форума


Посмотреть правила форума



Начать новую тему Ответить на тему На страницу 1, 2  След.
 
 Резольвента самосопряженного оператора
Сообщение22.02.2016, 03:06 


21/02/15
27
Москва
Столкнулся с такой задачей: найти резольвенту следующего интегрального оператора
$\mathscr K (x,y) = 
\begin{cases}
x,&\text{если } 0\leqslant x\leqslant y\leqslant1;\\
y,&\text{если } 0\leqslant y\leqslant x\leqslant1.
\end{cases}
$

Нетрудно выбрать для него систему ортонормальных на $L_2[0,1]$ собственных функций:
$\varphi_{m+1} = \sqrt{2} \sin \left[ (\frac{\pi}{2} + \pi m)x \right], \lambda_{m+1} = \left(\frac{\pi}{2} + \pi m \right)^2, m \in \mathbb N_0$
Поскольку ядро оператора непрерывно и эрмитово, резольвента представляется в виде билинейного ряда:
$\mathscr R (x,y; \lambda) = \sum\limits_{n=1}^{\infty} \frac{\varphi_n(x) \bar{\varphi}_n(y)}{\lambda_n - \lambda} = \sum\limits_{m=0}^{\infty} \frac{2 \sin \left[ (\frac{\pi}{2} + \pi m)x \right] \sin \left[ (\frac{\pi}{2} + \pi m)y \right]}{(\frac{\pi}{2} + \pi m)^2 - \lambda}$
Найти функцию в обычном виде у меня не получилось.

Вопрос такой: есть ли способ вычисления такого ряда и можно ли найти эту резольвенту как-нибудь проще?

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение22.02.2016, 10:34 
Заслуженный участник
Аватара пользователя


31/01/14
11057
Hogtown
Очевидно, с.з. найдены неверно: оператор ограниченный, а они у Вас стремятся к бесконечности (поскольку оператор компактный, то они должны стремиться к 0).

Подсказка: Ваш интегральный оператор можно описать так: $f\to u$, где $u$ решение $u''=-f, u(0)=u'(1)=0$. Поэтому с.з. будут обратными к тому, что Вы нашли.

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение22.02.2016, 15:30 


21/02/15
27
Москва
Red_Herring в сообщении #1101229 писал(а):
Очевидно, с.з. найдены неверно

$\lambda_{m+1}$ - это не с.з. оператора, а его характеристические числа. И в представлении резольвенты $\mathscr R (x,y; \lambda)$ имеются ввиду именно характеристические числа.
Red_Herring в сообщении #1101229 писал(а):
с.з. будут обратными к тому, что Вы нашли

Все верно, $\Lambda_{m+1} = \left\{ \sqrt{2} \sin \left[ (\frac{\pi}{2} + \pi m) x \right] \right\}^{-1}$

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение22.02.2016, 15:39 
Заслуженный участник
Аватара пользователя


31/01/14
11057
Hogtown
А по-вашему, что такое резольвента оператора?

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение22.02.2016, 16:28 


21/02/15
27
Москва
Red_Herring в сообщении #1101287 писал(а):
А по-вашему, что такое резольвента оператора?

$\mathscr R (x,y; \lambda) = \sum\limits_{k=0}^{\infty} \lambda^k {\mathscr K}_{k+1} (x,y)$, где $\mathscr K_{k+1} (x,y)$ - повторные ядра ядра $\mathscr K (x,y)$.
В таком случае решение интегрального уравнения $u = \lambda K u + f$ представляется (в достаточно хороших случаях) в виде $u = f + \lambda R f$, где $R$ - интегральный оператор с ядром $\mathscr R (x,y; \lambda)$

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение22.02.2016, 17:16 
Заслуженный участник
Аватара пользователя


31/01/14
11057
Hogtown
Я допускаю, что такое определение ($(I-\lambda K)^{-1})$) где-то даётся, однако стандартное определение это $(K-\lambda)^{-1}$.

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение24.02.2016, 12:20 
Заслуженный участник


11/05/08
32166
Red_Herring в сообщении #1101321 писал(а):
Я допускаю, что такое определение ($(I-\lambda K)^{-1})$) где-то даётся, однако стандартное определение это $(K-\lambda)^{-1}$.

Оно конечно, однако применительно конкретно к интегральным уравнениям лямбды почему-то принято ставить при операторах (но, правда, и называть не собственными, а характеристическими).

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение24.02.2016, 13:50 
Заслуженный участник
Аватара пользователя


31/01/14
11057
Hogtown
ewert в сообщении #1101722 писал(а):
однако применительно конкретно к интегральным уравнениям

Исторические причины.

А так Ваше выражение максимально простое возможное

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение24.02.2016, 19:54 
Заслуженный участник


10/01/16
2315
Red_Herring

Red_Herring в сообщении #1101732 писал(а):
А так Ваше выражение максимально простое возможное


А чё, ну никак не упростить разве? Разве ряд из первого поста e.mazhnik не есть чей то там ряд Фурье?
Или: решая дифур из Вашего поста (с добавком $\lambda u$) методом вариации постоянной, например, мы разве не получим интегральное представление для резольвенты? Там в какой-то момент появляется нехороший член с $f'$, но потом он благополучно исчезает. Вроде бы...

-- 24.02.2016, 20:57 --

e.mazhnik
Мне кажется, что и повторные ядра вполне считабельны, и ряд - мобыть - свернется

 Профиль  
                  
 
 Posted automatically
Сообщение24.02.2016, 20:09 


20/03/14
12041
 i  Тема перемещена из форума «Помогите решить / разобраться (М)» в форум «Карантин»
по следующим причинам:

- неправильно набраны формулы (краткие инструкции: «Краткий FAQ по тегу [math]» и видеоролик Как записывать формулы);

e.mazhnik
Каждая формула должна начинаться со знака доллара и заканчиваться им. Исправьте по всей теме, пожалуйста.

Исправьте все Ваши ошибки и сообщите об этом в теме Сообщение в карантине исправлено.
Настоятельно рекомендуется ознакомиться с темами Что такое карантин и что нужно делать, чтобы там оказаться и Правила научного форума.

 Профиль  
                  
 
 Posted automatically
Сообщение25.02.2016, 11:16 


20/03/14
12041
 i  Тема перемещена из форума «Карантин» в форум «Помогите решить / разобраться (М)»

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение25.02.2016, 12:01 
Заслуженный участник
Аватара пользователя


31/01/14
11057
Hogtown
DeBill в сообщении #1101811 писал(а):
А чё, ну никак не упростить разве? Разве ряд из первого поста e.mazhnik не есть чей то там ряд Фурье?
Или: решая дифур из Вашего поста (с добавком $\lambda u$) методом вариации постоянной, например, мы разве не получим интегральное представление для резольвенты?


Да, конечно. В первом случае надо произведение синусов разложить в разность синусов и работать с функциями одной переменной $x+y$ и $ x-y$.

Но всё-таки то что записано в оригинальном посте это ядро Шварца $K(1-\lambda K)^{-1}$, а не того, что обычно называется резольвентой.

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение26.02.2016, 03:06 


21/02/15
27
Москва
DeBill в сообщении #1101811 писал(а):
решая дифур из Вашего поста (с добавком $\lambda u$)

Да, так и правда получилось. Правда, я опять, возможно, все усложняю.
Рассмотрим интегральное уравнение $u = \lambda K u + f$. Путем несложных преобразований получим краевую задачу $\{ u'' + \lambda u = f'', u(0) = 0, u'(1) = 0 \}$.
Найдем функцию Грина $G(x,y)$ краевой задачи. Для это выберем два линейно независимых решения, удовлетворяющих разным граничным условиям: $u_1 = \sin \sqrt{\lambda} x$ и u_2 = \cos \sqrt{\lambda} (x-1)$. Определитель Вронского равен $\omega = -\sqrt{\lambda} \cos \sqrt{\lambda}$.
Тогда $
G(x,y) = \frac{-1}{\sqrt{\lambda} \cos \sqrt{\lambda}}
\begin{cases}
\sin \sqrt{\lambda}x \cos \sqrt{\lambda} (y-1),&\text{если } 0 \leqslant x \leqslant y;\\
\sin \sqrt{\lambda}y \cos \sqrt{\lambda} (x-1),&\text{если } y \leqslant x \leqslant 1.
\end{cases}
$
Общее решение: $u(x) = \int\limits_{0}^{1} G(x,y) f'' (y)\,dy$. Из него следуют условия разрешимости: $\{f(0) = u(0) = 0, f'(1) = u'(1)= 0\}$.
Применив интегрирование по частям получим $u(x) = G(x,y) f'(y) \bigg|_0^1 - \int\limits_{0}^{1} G'_y(x,y) f'(y)\,dy = - \int\limits_{0}^{1} G'_y(x,y) f'(y)\,dy$.
Ядро последнего интеграла $G'_y (x,y)$ разрывно в $x$, поэтому разобьем его на два непрерывных интервала $G'_y (x,y) = {G_1}'_y (x,y) + {G_2}'_y (x,y)$ и проинтегрируем по частям:
$u(x) = -f(y) {G_1}'_y (x,y) \bigg|_0^x - f(y) {G_2}'_y (x,y) \bigg|_x^1 + \int\limits_0^1 G''_{yy} (x,y) f (y)\,dy = {f(x) +\int\limits_{0}^{1} G''_{yy} (x,y) f(y)\,dy}$.
e.mazhnik в сообщении #1101304 писал(а):
решение интегрального уравнения $u = \lambda K u + f$ представляется (в достаточно хороших случаях) в виде $u = f + \lambda R f$, где $R$ - интегральный оператор с ядром $\mathscr R (x,y; \lambda)$

В соответствии с этим определением $\mathscr R (x,y; \lambda) = 
\frac{1}{\sqrt{\lambda} \cos \sqrt{\lambda}}
\begin{cases}
\sin \sqrt{\lambda}x \cos \sqrt{\lambda} (y-1),&\text{если } 0 \leqslant x \leqslant y;\\
\sin \sqrt{\lambda}y \cos \sqrt{\lambda} (x-1),&\text{если } y \leqslant x \leqslant 1.
\end{cases}$

-- 26.02.2016, 03:51 --

Red_Herring в сообщении #1101971 писал(а):
Но всё-таки то что записано в оригинальном посте это ядро Шварца $K(1-\lambda K)^{-1}$, а не того, что обычно называется резольвентой.

Здесь имелось ввиду $(1 - \lambda K)^{-1}$ ? Или я чего-то не понимаю.

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение26.02.2016, 06:12 
Заслуженный участник
Аватара пользователя


31/01/14
11057
Hogtown
e.mazhnik в сообщении #1102172 писал(а):
Здесь имелось ввиду

Во первых "имелось в виду", а во вторых "ввиду того что K интегральный оператор $(1-\lambda K)^{-1} $ интегральным оператором быть не может". А вот если отнять от него $1$, так пожалуйста.

 Профиль  
                  
 
 Re: Резольвента самосопряженного оператора
Сообщение26.02.2016, 08:57 


21/02/15
27
Москва
Red_Herring в сообщении #1102177 писал(а):
Во первых "имелось в виду"

Извиняюсь, нелепая опечатка. Стараюсь такого не допускать. Впрочем, не обособлять вводные слова тоже не есть хорошо.
Red_Herring в сообщении #1102177 писал(а):
во вторых "ввиду того что K интегральный оператор $(1-\lambda K)^{-1} $ интегральным оператором быть не может". А вот если отнять от него $1$, так пожалуйста.

Да, само собой. Просто встречал в литературе такое обозначение резольвенты, что как бы подразумевает, что еще $1$ должна быть. Вы и сами так написали в первый раз:
Red_Herring в сообщении #1101321 писал(а):
Я допускаю, что такое определение ($(I-\lambda K)^{-1}$) где-то даётся

Так я правильно понимаю, что ядро оператора $(1 - \lambda K)^{-1} - 1$ обычно именуется ядром Шварца?

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 19 ]  На страницу 1, 2  След.

Модераторы: Модераторы Математики, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group