2014 dxdy logo

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

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


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


В этом разделе нельзя создавать новые темы.

Если Вы хотите задать новый вопрос, то не дописывайте его в существующую тему, а создайте новую в корневом разделе "Помогите решить/разобраться (М)".

Если Вы зададите новый вопрос в существующей теме, то в случае нарушения оформления или других правил форума Ваше сообщение и все ответы на него могут быть удалены без предупреждения.

Не ищите на этом форуме халяву, правила запрещают участникам публиковать готовые решения стандартных учебных задач. Автор вопроса обязан привести свои попытки решения и указать конкретные затруднения.

Обязательно просмотрите тему Правила данного раздела, иначе Ваша тема может быть удалена или перемещена в Карантин, а Вы так и не узнаете, почему.



Начать новую тему Ответить на тему
 
 Найти логарифм от матрицы (решил, но получается фигня)
Сообщение11.10.2013, 00:21 


25/06/13
27
Нужно найти натуральный логарифм от матрицы A = $ \left[\begin{array}{rr}x & 1 \\ 1 & x \\\end{array}\right]$ Из курса функана — есть формула, похожая на интегральную формулу Коши: $$ f(A) = {1 \over {2 \pi i }}\oint_C f(z)(zI - A)^{-1} dz,$$ где контур охватывает собственные значения матрицы. Определитель $\left[\begin{array}{rr}x - z & 1 \\ 1 & x - z \\\end{array}\right] = (x -  z)^2 - 1 = 0 $ => Собственные значения $ z =  x \pm  1$
$$(zI - A)^{-1} = \left[\begin{array}{rr}z - x & -1 \\ -1 & z - x \\\end{array}\right]^{-1} = \left[\begin{array}{rr}z - x & 1 \\ 1 & z - x \\\end{array}\right]^{T}  = \left[\begin{array}{rr}z - x & 1 \\ 1 & z - x \\\end{array}\right]$$
Как мне объяснили, интеграл берётся поэлементно, т.е. нужно найти $ \oint_C \ln(z) dz $ и $ \oint_C (z-x) \ln(z) dz $. Их проще всего посчитать через вычеты. Если у нас $x \ne \pm 1 $, то контур не охватывает 0 => никаких особых точек внутри него нет и интеграл равен 0.
Потом я нашёл порядок полюса в точке 0 для $  f_1 = \ln(z)  $ и $  f_{2,3} = (z \pm 1) \ln(z)  $. Он равен порядку нуля $  F_1 = {1 \over \ln(z)} $ и $  F_{2,3} = { 1 \over (z \pm 1) \ln(z) }  $. Посчитал производные в вольфраме, там же посчитал пределы - вышло, что уже у первых производных в $z = 0$ значение не 0, а вообще особая точка. Т.е. выходит порядок полюса первый. Потом по этой
$$ {1 \over (n+1)!} \lim {d^{n-1} \over dz^n }[(z-a)^n f(z)]$$формуле для вычетов в полюсах кратности n нашёл их, и они все получились равны 0.
Т.е. у меня получилась нулевая матрица. Мне кажется, я сделал что-то неправильно, пожалуйста, подскажите где.

 Профиль  
                  
 
 Re: Найти логарифм от матрицы (решил, но получается фигня)
Сообщение11.10.2013, 00:30 
Заслуженный участник
Аватара пользователя


06/10/08
6422
Обратная матрица неправильно найдена (забыли поделить на определитель).

 Профиль  
                  
 
 Re: Найти логарифм от матрицы (решил, но получается фигня)
Сообщение11.10.2013, 00:33 


25/06/13
27
Xaositect

Спасибо, воспользовался формулой по памяти и ошибся. Сейчас пересчитаю, посмотрю, что выйдет.

 Профиль  
                  
 
 Re: Найти логарифм от матрицы (решил, но получается фигня)
Сообщение11.10.2013, 18:38 


25/06/13
27
Проверьте решение, пожалуйста, у меня только одна попытка, чтобы сдать, я на грани вылета.
Итак, матрица получилась $$(zI - A)^{-1} = {1 \over (z - x)^2 - 1 }  \left[\begin{array}{rr}z - x & 1 \\ 1 & z - x \\\end{array}\right]$$
Осталось найти $ \oint_C {\ln(z) \over { (z - x)^2 - 1 } } dz $ и $ \oint_C { (z-x) \ln(z)  \over  (z - x)^2 - 1 }dz $. Будем считать вычеты в точке $ x \pm 1$. Обозначим $f_1(z) =  {\ln(z) \over { (z - x)^2 - 1 } } = g_1(z)/h_1(z)$ Тогда $$\operatorname{Res} [f_1(z), x-1] = {g_1(x-1) \over h_1'(x-1)} = - {\ln(x - 1) \over { 2  } }  $$
$$\operatorname{Res} [f_1(z), x+1] = {g_1(x+1) \over h_1'(x+1)} = {\ln(x + 1) \over { 2  } }  $$
Обозначим $f_2(z) =  { (z - x) \ln(z) \over { (z - x)^2 - 1 } } = g_2(z)/h_2(z)$ Тогда $$\operatorname{Res} [f_2(z), x-1] = {g_2(x-1) \over h_2'(x-1)} = {\ln(x - 1) \over { 2  } }  $$
$$\operatorname{Res} [f_2(z), x+1] = {g_2(x+1) \over h_2'(x+1)} = {\ln(x + 1) \over { 2  } }  $$
Т.к. интеграл по контуру — это сумма вычетов во внутренних точках на $2 \pi i$, то в ответе получим
$$ \ln A = {1 \over 2} \left[\begin{array}{rrrrr} \ln (x^2 - 1) & \ln { {x + 1} \over {x - 1}} \\ \ln { {x + 1} \over {x - 1} } & \ln (x^2 - 1) \\\end{array}\right] $$

 Профиль  
                  
 
 Re: Найти логарифм от матрицы (решил, но получается фигня)
Сообщение11.10.2013, 18:56 
Заслуженный участник
Аватара пользователя


06/10/08
6422
Ответ правильный, вычисления не проверял.

 Профиль  
                  
 
 Re: Найти логарифм от матрицы (решил, но получается фигня)
Сообщение11.10.2013, 19:13 


25/06/13
27
Спасибо большое.

 Профиль  
                  
 
 Re: Найти логарифм от матрицы (решил, но получается фигня)
Сообщение17.10.2013, 23:11 


25/06/13
27
Преподаватель посмотрел на ответ и сказал, что я где-то потерял минус. И по его мнению матрица будет действительной для $|x| < 1$ Т.е. выходит, он считает, что вместо $ \ln {(x^2 - 1)}$ должно быть $ \ln {(1 - x^2)}$. Но ведь если $X = \ln A$ и $a_i$ — собственные числа $A$, а $x_i$ — собственные числа $X$, тогда должно выполняться $e^{x_i} = a_i $. Я нашёл собственные числа этой матрицы, они получились $ \ln {(x-1)}$ и $ \ln {(x+1)}$ Либо я ошибся дважды, либо преподаватель не прав.

 Профиль  
                  
 
 Re: Найти логарифм от матрицы (решил, но получается фигня)
Сообщение17.10.2013, 23:25 
Заслуженный участник
Аватара пользователя


06/10/08
6422
Я находил так: привел матрицу к диагональной, взял логарифмы от собственных значений на диагонали и перешел обратно. У меня получился такой же ответ, как у Вас.

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 8 ] 

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



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

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


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

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