2014 dxdy logo

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

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


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


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



Начать новую тему Ответить на тему На страницу Пред.  1, 2, 3, 4  След.
 
 Re: Формула Иенсена
Сообщение04.01.2018, 18:59 
Заслуженный участник
Аватара пользователя


22/06/12
2129
/dev/zero
Появилась одна идея. Для того, чтоб её реализовать, нужно вычислить
$$
I = \int \limits_0^{2 \pi} \ln (1 + a^2 - 2 a \cos \varphi) \ \mathrm d \varphi, \qquad 0 \leqslant a < 1.
$$
Его можно посчитать по вычетам.
$$
I = \dfrac{1}{2 i \pi} \oint \limits_{|z|=1} \dfrac{\ln\left(1 + a^2 - a \left( \dfrac{1}{z} + z \right)\right)}{z} \ \mathrm dz.
$$
Дальше хочется разложить в ряд Лорана подынтегральную функцию, но с агрегатом вида $z + \dfrac{1}{z}$ я могу лишь сделать $\dfrac{z^2 + 1}{z}$ и дальше из-под логарифма вынести член $\ln z$ с получением
$$
\ln\left(1 + a^2 - a \left( \dfrac{1}{z} + z \right)\right) = \ln \left(z(1 + a^2) - az^2 - a \right) - \ln z.
$$
Тут нужно будет фиксировать ветви функции. Достаточно ли фиксировать, скажем, главную и затем взять разрез по лучу $[0, +\infty)$ для такого преобразования?

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


27/12/17
1411
Антарктика
Попробуйте рассмотреть любую регулярную ветвь функции $\ln(a-\xi)$ в области $\left\lvert\xi\right\rvert>1$, а затем подставить туда $\xi=e^{i\varphi}$. Лобовой способ с выбором более очевидной функции вряд ли подойдет.

UPD. Еще надо будет какой-то разрез сделать, чтобы никакой замкнутый контур не обходил точку $a$

-- 04.01.2018, 21:18 --

StaticZero в сообщении #1281250 писал(а):
но с агрегатом вида $z + \dfrac{1}{z}$

Этот агрегат уже разложен в ряд Лорана в нужной Вам окрестности особой точки 0

 Профиль  
                  
 
 Re: Формула Иенсена
Сообщение04.01.2018, 19:46 
Заслуженный участник
Аватара пользователя


22/06/12
2129
/dev/zero
Попробую так. Если сесть на окружность $|z| = \varepsilon$, то $\xi = z + \dfrac{1}{z}$, $|\xi| = 2\varepsilon |\cos \varphi|$, где $\varphi = \operatorname{Arg} z$.

Теперь $\ln \left(1 + a^2 - a \xi\right) = \ln (1 + a^2) + \ln \left(1 - \dfrac{a \xi}{1 + a^2}\right) = \ln(1 + a^2) + \sum \limits_{k=1}^\infty \dfrac{(-1)^k}{k} \left(\dfrac{a}{1 + a^2}\right)^k \xi^k$. А теперь вычет считать в нуле...

-- 04.01.2018, 20:08 --

У меня выходит $\mathop{\operatorname{res}}\limits_{z = 0} \dfrac{1}{z} \left(z + \dfrac{1}{z}\right)^k = \binom{k}{k/2}$, $k$ чётное, или $0$, если $k$ нечётное. Получаем
$$
\dfrac{1}{2 i \pi} \oint \limits_{|z| = 1} \dfrac{\ln \left(1 + a^2 - a \xi \right)}{z} \ \mathrm dz = \ln (1 + a^2) + \sum \limits_{k=1}^\infty \binom{2k}{k} \dfrac{1}{2k} \left(\dfrac{a}{1 + a^2}\right)^{2k}.
$$
Верно?

 Профиль  
                  
 
 Re: Формула Иенсена
Сообщение04.01.2018, 20:21 
Заслуженный участник
Аватара пользователя


24/02/12
1842
Москва
У Вас же функция неоднозначная. Где разрезы и интегралы по ним?

-- 04.01.2018, 20:24 --

Давайте так. Если бы функция не имела ни нулей, ни полюсов в круге, Вы смогли бы доказать формулу Иенсена для нее?

 Профиль  
                  
 
 Re: Формула Иенсена
Сообщение04.01.2018, 20:30 
Заслуженный участник
Аватара пользователя


22/06/12
2129
/dev/zero
ex-math в сообщении #1281280 писал(а):
Давайте так. Если бы функция не имела ни нулей, ни полюсов в круге, Вы смогли бы доказать формулу Иенсена для нее?

Нет, не смог бы. Но я хочу пока разобраться со случаем, когда в разложении
$$
f(z) = \dfrac{\prod \limits_\text{нули} (z - z_k^+)}{\prod \limits_\text{полюсы} (z - z_k^-)} g(z),
$$
где $g(z)$ не имеет ни нулей, ни полюсов в круге, $g = 1$. Для этого надо вычислить $\ln |z - z_0|$, где $z_0$ — полюс или нуль. Для этого и нужен тот интеграл, который я пытаюсь взять (он вещественный).

 Профиль  
                  
 
 Re: Формула Иенсена
Сообщение04.01.2018, 20:33 
Заслуженный участник
Аватара пользователя


24/02/12
1842
Москва
А попробуйте. Пусть $f$ голоморфна и не имеет нулей.

 Профиль  
                  
 
 Re: Формула Иенсена
Сообщение04.01.2018, 20:49 
Заслуженный участник
Аватара пользователя


22/06/12
2129
/dev/zero
ex-math в сообщении #1281280 писал(а):
У Вас же функция неоднозначная. Где разрезы и интегралы по ним?

Так?

(Оффтоп)

Изображение

 Профиль  
                  
 
 Re: Формула Иенсена
Сообщение04.01.2018, 20:53 
Заслуженный участник
Аватара пользователя


24/02/12
1842
Москва
Точка ветвления не в нуле, а в $a$.

-- 04.01.2018, 20:56 --

А дифференцировать Ваш интеграл по $a$ пробовали?

 Профиль  
                  
 
 Re: Формула Иенсена
Сообщение04.01.2018, 20:59 
Заслуженный участник
Аватара пользователя


22/06/12
2129
/dev/zero
ex-math, а тут не понял. Внутри круга $|z| = 1$ у функции $1 + a^2 - a \left( z + \dfrac{1}{z} \right)$ корней нет, они находятся только на единичной окружности. Если обойти эти точки маленькими дугами, загнутыми внутрь круга, то останется только точка $z = 0$.

 Профиль  
                  
 
 Re: Формула Иенсена
Сообщение04.01.2018, 21:05 
Заслуженный участник
Аватара пользователя


24/02/12
1842
Москва
А как же $z=a$?

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


22/06/12
2129
/dev/zero
ex-math в сообщении #1281293 писал(а):
А как же $z=a$?

А, нет, я дурак: тут оба корня вещественные. Понял.

-- 04.01.2018, 21:15 --

Тогда разрез надо брать по лучу $[a, +\infty)$ и контур будет такой же, как на картинке, только окружностью будем обходить точку $a$, а ноль останется внутри, верно?

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


24/02/12
1842
Москва
Дифференцируйте по параметру.
А потом придется вернуться к случаю голоморфной функции без нулей.

-- 04.01.2018, 21:18 --

Неужели не хочется от логарифма избавиться?

 Профиль  
                  
 
 Re: Формула Иенсена
Сообщение04.01.2018, 21:45 
Заслуженный участник
Аватара пользователя


22/06/12
2129
/dev/zero
ex-math, хочется. Но один раз-то нужно продраться, чтоб что-то понять, с многозначными функциями пока очень не очень.

Ладно.
$$
I(a) = \dfrac{1}{2 \pi} \int \limits_0^{2 \pi} \ln (1 + a^2 - 2 a \cos \varphi) \ \mathrm d\varphi.
$$
$$
I'(a) = \dfrac{1}{2 \pi} \int \limits_0^{2 \pi} \dfrac{2 a - 2 \cos \varphi}{1 + a^2 - 2 a \cos \varphi} \ \mathrm d \varphi.
$$
$$\begin{align*}
J_1(a) &= \dfrac{a}{i \pi} \oint \limits_C \dfrac{\mathrm dz}{z\left(1 + a^2 - a \left( z + \dfrac{1}{z} \right) \right)} = \dfrac{a}{i \pi} \oint \limits_C \dfrac{\mathrm dz}{z(-a) (z - a) \left(z - \dfrac{1}{a} \right) } = -2 \left(1 + \dfrac{1}{a \left(a - \dfrac{1}{a} \right)} \right) = \\ &=\dfrac{2 a^2}{1 - a^2}.
\end{align*}
$$
$$\begin{align*}
J_2(a) &= \dfrac{1}{2 i \pi} \oint \limits_C \dfrac{\left(z + \dfrac{1}{z} \right)\mathrm dz}{z(-a) (z - a) \left(z - \dfrac{1}{a} \right) } = \dfrac{1}{2 i \pi} \oint \limits_C \dfrac{\left(z^2 + 1\right)\mathrm dz}{z^2 (z - a) \left(1 - a z\right) } = - \dfrac{1 + a^2}{a^2} + \dfrac{1 + a^2}{a^2 (1 - a^2) } = \dfrac{1 + a^2}{1 - a^2}.
\end{align*}
$$
$$
I'(a) = J_1(a) - J_2(a) = \dfrac{2 a^2 - 1 - a^2}{1 - a^2} = -1.
$$
$$
I(a) = -a + Q, \qquad I(0) = 0 = Q, \qquad I = -a.
$$

-- 04.01.2018, 21:49 --

$C$ — единичная окружность

 Профиль  
                  
 
 Re: Формула Иенсена
Сообщение04.01.2018, 21:54 
Заслуженный участник
Аватара пользователя


24/02/12
1842
Москва
У меня не так. Проверьте вычет в нуле во втором интеграле -- там полюс второго порядка.

-- 04.01.2018, 21:57 --

А, не, в первом $z$ не хватает в числителе.

-- 04.01.2018, 21:59 --

И во втором.

 Профиль  
                  
 
 Re: Формула Иенсена
Сообщение04.01.2018, 22:07 
Заслуженный участник
Аватара пользователя


22/06/12
2129
/dev/zero
$$
\mathop{\operatorname{res}} \limits_{z = 0} \dfrac{z^2 + 1}{z^2 (z-a)(1 - az)} = \left. \dfrac{\mathrm d}{\mathrm dz} \dfrac{z^2 + 1}{(z-a)(1 - az)}\right|_{z = 0} = \left. - \dfrac{(z^2+1)\left \lbrace{ (1 - az) + (-a)(z - a) \right \rbrace }}{(z-a)^2(1 - az)^2}\right|_{z = 0} = - \dfrac{1 + a^2}{a^2}.
$$

-- 04.01.2018, 22:08 --

ex-math в сообщении #1281313 писал(а):
А, не, в первом $z$ не хватает в числителе.

-- 04.01.2018, 21:59 --

И во втором.

Простите?

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

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



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

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


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

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