Модераторы: Poseidon

Поиск:

Ответ в темуСоздание новой темы Создание опроса
> [Линал] Жордановы формы 
:(
    Опции темы
JAPH
Дата 7.12.2007, 21:33 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 312
Регистрация: 8.1.2007
Где: Ленобласть

Репутация: 12
Всего: 23



Кстати, вполне может оказаться, что решений у системы (A-lE)h=x нет (Например, когда x достигает максимально возможной высоты, то есть (A-lE)^(m-1) x <> 0, где m - кратность собственного числа в минимальном аннуляторе A. Тогда ясно, что для любого вектора h верно (A-lE)^m h = 0, но (A-lE)^m h = (A-lE)^(m-1) (A-lE)h = (A-lE)^(m-1) x <> 0. противоречие), или же решение линейно-зависимо с ранее полученными.

Цитата
Если не получится представить в виде едеениц и нулей. ну допустим
2 0 1
3 4 0. Тут нету единичной м-цы. Ну или любой другой пример. 

smile Честно?
меняем строки местами
3 4 0
2 0 1
делим первую на четыре
0.75 1 0
   2   0 1
Если строго следовать алгоритму Гаусса, сомнений не будет, что ранговая единичная подматрица всегда существует, другое дело, что вычисления не всегда удобные..

Ripper, а обязательно искать сначала Жорданову форму, потом базис, в котором она реализуется? И обязательно через присоединённые вектора? А то я, понимаешь ли, гну свою линию, а, может, задание такое... smile По-моему, схема с башнями проще - сразу находится и Жорданова форма, и соответствующий базис. Если что непонятно в определениях, алгоритме, .. - спрашивайте. Могу рассказать теорию так, как нам давали, с примерами, конечно (хотя давали без примеров smile), может, вам так удобнее будет..

Это сообщение отредактировал(а) JAPH - 7.12.2007, 21:51


--------------------
Что непонятно - спрашиваем smile
PM MAIL ICQ   Вверх
cardinal
Дата 7.12.2007, 22:18 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Инженер
****


Профиль
Группа: Экс. модератор
Сообщений: 6003
Регистрация: 26.3.2002
Где: Германия

Репутация: 15
Всего: 99



Цитата(Ripper @  7.12.2007,  18:36 Найти цитируемый пост)
Но вопрос про свободные переменные (см. выше) я так и не совсем уверен. Т.е. ты тут выбрал x2=1, но тут то ладно. а если система 3х3 хотя бы.. как что выбирать.. )

Тут без разницы какая система. Если других условий кроме x2=-x3 например нет, то тогда просто говоришь, что x2=1, x3=-1. x2=0 делать нельзя!


--------------------
Немецкая оппозиция потребовала упростить натурализацию иммигрантов
В моем блоге: Разные истории из жизни в Германии

"Познание бесконечности требует бесконечного времени, а потому работай не работай - все едино".  А. и Б. Стругацкие
PM   Вверх
Ripper
Дата 7.12.2007, 23:01 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Lonely soul...
**


Профиль
Группа: Участник
Сообщений: 920
Регистрация: 30.6.2004
Где: г. Москва

Репутация: 6
Всего: 23



JAPH, чесно говоря эту схему с башнями я плохо понимаю, из за обилия новых слов )
Да, у нас вроже обязательно сначала форму, потом базис. А потом ещё и ф-ии от матриц, чуствую гемора будет...  smile 
Просто у нас ещё все по разному эти присоединенные ищут, там (a-le)^2h=0 например, незнаю можно так или нет. 

Блин smile
Цитата

1. строим оператор B = A-tE. Находим ФСР СЛАУ Bx = 0, т.е. базис Ker(B). Располагаем эти вектора в первой строке. Если получили k векторов, то goto 2. Иначе строим базис Ker(B^2), т.е. находим ФСР BBx=0, из которого выбираем вектора, линейно независимые с ранее полученными, и располагаем их во второй строке (над первой). И так далее, пока не наберём ровно k векторов.


Беру число t=3 найденное из х.у. кратности 2. решаю  (A-tE)x=0,  нахожу (-1,-1,1). (Я могу ошибатся в числах, уж че то много опечаток постоянно. мне главное суть). Значит дальше я ищу (A-tE)^2 x  =0? На сайте где я ссылку кинул, там (A-tE)h=x, и меня это смущает. Так, ну решил я (A-tE)^2x=0 вроде получил (0,1,0). 
Почему распологаем по строкам. разве не (0,1,0)* брать надо? Не в столбцах?
Цитата

2. Берём последнюю строку, ко всем векторам в ней применяем B, записываем результат в предыдущую строку вместо того, что там было, при необходимости дополняем их до прежней длины строки векторами, которые там были, но так, чтобы набор был л/нез.
После всех таких преобразований получим: в первой строке собственные вектора, так что для любого s из первой строки Bs = 0 <=> As = ts. Во второй строке корневые вектора высоты два, для любого вектора p из второй строки BBp = 0; по построению Bp = s - вектор под p, значит, (A-tE)p = s, т.е. p - как раз присоединённый вектор. То есть алгоритм "схема с башнями" даёт те же вектора, просто с другого конца.

Что значит применить к последней строке В? Мы получили
-1,-1,1
0,1,0.
Последняя строка 0,1,0 (кстати я бы это все перевернул к столбцам т.е. транспонировал изначально решения ...) 
И что значит к 0,1,0 применить В? 
И чем отличается вообще алгоритм башнями от присоединенных? (a-te)^2x=0 и (a-te)h=x? 
Как узнать сколько переменных, которым мы должны присвоить будем константы. или как там, определитель единичный. N-r? тогда это 2 в примере с неоднородной системой а не 1. И как то все совсем подругому :( 
решение неоднородной = решение однородной + частное решение. ? И что брать в качестве вектора, частное?
Как то можно это все решать не приводя к единичной м-цы, а то как то меня это в ступор вводит smile Новые методы

Добавлено через 1 минуту и 29 секунд
Да, я нашел базис, решая ур-ия (a-te)^1 или в ^2 *x=0. 
и получил 
-1 0 -1
-1 -1 0
1 1 1
Но почему то J=C^-1*A*C не выходит. C^-1 точно верно нашел(проверил). 
Строки C находил из решений ур-ий для t=-2,t=3(кратность 2). 



--------------------
"Он знает: надо смеяться над тем, что тебя мучит, иначе не сохранишь равновесия, иначе мир сведет тебя с ума" - Над кукушкиным гнездом
PM MAIL ICQ   Вверх
JAPH
Дата 7.12.2007, 23:41 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 312
Регистрация: 8.1.2007
Где: Ленобласть

Репутация: 12
Всего: 23



C - матрица перехода - составляется по столбцам.
Цитата
И чем отличается вообще алгоритм башнями от присоединенных? (a-te)^2x=0 и (a-te)h=x? 
Цитата
После всех таких преобразований получим: в первой строке собственные вектора, так что для любого s из первой строки Bs = 0 <=> As = ts. Во второй строке корневые вектора высоты два, для любого вектора p из второй строки BBp = 0; по построению Bp = s - вектор под p, значит, (A-tE)p = s, т.е. p - как раз присоединённый вектор. То есть алгоритм "схема с башнями" даёт те же вектора, просто с другого конца.


Не совсем так вы поняли, что есть строки. Вектора (столбцы координат) располагаются в схему.
Например, вектор (1 1 -1)=s,(0 1 0)=p расположатся так:

p
s

, т.е. один под другим, s в первой строке/ступеньке, p во второй. Схема сначала строится снизу вверх, путем нахождения базиса Ker(B), достраивания его до Ker(B^2), и т.д. Затем проходят сверху вниз: берут верхний вектор, p, применяют к нему B, то есть вычисляют Bp, записывают под p.

p
Bp

(если бы нижняя ступенька была больше верхней, то помимо B-образов верхней ступеньки мы бы дописали вектора, оказавшиеся там на проходе снизу-вверх, линейно-независимые с B-образами.)

По сути, p есть вектор, присоединённый к Bp. То есть схема с башнями строит сначала последний присоединённый вектор, потом спускается вниз до собственного вектора. А вам предлагают наоборот.


--------------------
Что непонятно - спрашиваем smile
PM MAIL ICQ   Вверх
JAPH
Дата 8.12.2007, 00:59 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 312
Регистрация: 8.1.2007
Где: Ленобласть

Репутация: 12
Всего: 23



Считаю нерациональным искать сначала Жорданову форму, так как формула для размера клетки доказывается через схему с башнями.

Ну, допустим, рассмотрю такой пример. Попробую даже вашим способом. Просто он мне кажется интереснее.

Присоединённый файл ( Кол-во скачиваний: 13 )
Присоединённый файл  jordan.rar 30,24 Kb


--------------------
Что непонятно - спрашиваем smile
PM MAIL ICQ   Вверх
Ripper
Дата 9.12.2007, 20:14 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Lonely soul...
**


Профиль
Группа: Участник
Сообщений: 920
Регистрация: 30.6.2004
Где: г. Москва

Репутация: 6
Всего: 23



Блин, все равно ниче не понимаю. Это не та схема как мы решали
1. Находим собственные значения. Ок.
2. находим ЖОРДАНОВУ м-цу. сразу. 
http://www.dep805.ru/education/kk/jmatrix/part4.htm вот как тут все описавается
Только не понятно 1)как решать ур-ия не находя едичную м-цу, сколько и как выбирать переменных которые можно взять за константу, что делать если n-r>1 т.е. когда два лн решения, какой брать и что потом делать, при решении (a-le)h=x. и т.п. :(
 smile 


--------------------
"Он знает: надо смеяться над тем, что тебя мучит, иначе не сохранишь равновесия, иначе мир сведет тебя с ума" - Над кукушкиным гнездом
PM MAIL ICQ   Вверх
cardinal
Дата 9.12.2007, 20:25 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Инженер
****


Профиль
Группа: Экс. модератор
Сообщений: 6003
Регистрация: 26.3.2002
Где: Германия

Репутация: 15
Всего: 99



Совет: подойди к прикольной девушке с твоего курса и убей двух зайцев! smile 


--------------------
Немецкая оппозиция потребовала упростить натурализацию иммигрантов
В моем блоге: Разные истории из жизни в Германии

"Познание бесконечности требует бесконечного времени, а потому работай не работай - все едино".  А. и Б. Стругацкие
PM   Вверх
JAPH
Дата 9.12.2007, 21:14 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 312
Регистрация: 8.1.2007
Где: Ленобласть

Репутация: 12
Всего: 23



Цитата
как решать ур-ия не находя едичную м-цу

Устно smile Подбором smile
Просто метод Гаусса универсальный
Цитата
сколько и как выбирать переменных которые можно взять за константу

Не понял вопроса. Главных переменных столько, каков ранг матрицы, свободных n-r.
Цитата
что делать если n-r>1 т.е. когда два лн решения, какой брать и что потом делать, при решении (a-le)h=x. и т.п.

При нахождении собственных векторов - все решения, при нахождении присоединённых - одно любое.



--------------------
Что непонятно - спрашиваем smile
PM MAIL ICQ   Вверх
Ответ в темуСоздание новой темы Создание опроса
Правила форума "Центр помощи"

ВНИМАНИЕ! Прежде чем создавать темы, или писать сообщения в данный раздел, ознакомьтесь, пожалуйста, с Правилами форума и конкретно этого раздела.
Несоблюдение правил может повлечь за собой самые строгие меры от закрытия/удаления темы до бана пользователя!


  • Название темы должно отражать её суть! (Не следует добавлять туда слова "помогите", "срочно" и т.п.)
  • При создании темы, первым делом в квадратных скобках укажите область, из которой исходит вопрос (язык, дисциплина, диплом). Пример: [C++].
  • В названии темы не нужно указывать происхождение задачи (например "школьная задача", "задача из учебника" и т.п.), не нужно указывать ее сложность ("простая задача", "легкий вопрос" и т.п.). Все это можно писать в тексте самой задачи.
  • Если Вы ошиблись при вводе названия темы, отправьте письмо любому из модераторов раздела (через личные сообщения или report).
  • Для подсветки кода пользуйтесь тегами [code][/code] (выделяйте код и нажимаете на кнопку "Код"). Не забывайте выбирать при этом соответствующий язык.
  • Помните: один топик - один вопрос!
  • В данном разделе запрещено поднимать темы, т.е. при отсутствии ответов на Ваш вопрос добавлять новые ответы к теме, тем самым поднимая тему на верх списка.
  • Если вы хотите, чтобы вашу проблему решили при помощи определенного алгоритма, то не забудьте описать его!
  • Если вопрос решён, то воспользуйтесь ссылкой "Пометить как решённый", которая находится под кнопками создания темы или специальным флажком при ответе.

Более подробно с правилами данного раздела Вы можете ознакомится в этой теме.

Если Вам помогли и атмосфера форума Вам понравилась, то заходите к нам чаще! С уважением, Poseidon, Rodman

 
0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Центр помощи | Следующая тема »


 




[ Время генерации скрипта: 0.0483 ]   [ Использовано запросов: 22 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.