m=[1 2;3 4];
[q r]=qr(m);
eig_vec=q;
eig_val=r*q;
flag=true;
iter=0;
while(flag==1)
    prev_eig_val=eig_val;
    [q r]=qr(eig_val);
    eig_val=r*q;
    eig_vec=eig_vec*q;
    err=abs(diag(eig_val)-diag(prev_eig_val));
    if(mean(err)<0.00001)
        flag=false;
    end;
end;
eig_val=diag(eig_val);
 
disp('QR-алгоритм');
eig_val
eig_vec
disp('Стандартная функция eig()');
[eig_vec eig_val]=eig(m)

Программа выводит следующие данные:

QR-алгоритм

eig_val =

5.3723
-0.3723

eig_vec =

0.4160 -0.9094
0.9094 0.4160

Стандартная функция eig()

eig_vec =

-0.8246 -0.4160
0.5658 -0.9094

eig_val =

-0.3723 0
0 5.3723

Собственные значения посчиталась правильно, но вектора посчитались не верно. Не могу понять почему. Что может быть не так?

И второй вопрос по поводу критерия остановки, насколько верно то как я оцениваю точность решения, чтобы выйти из цикла?

Заказать работу у эксперта

Programming

Эксперт

94731 / 64177 / 26122

Регистрация: 12.04.2006

Сообщений: 116,782

24.09.2015, 12:01

Ответы с готовыми решениями:

Разработать м-файл для нахождения собственных чисел и собственных векторов матрицы
1 1 3 -1 2 0 2 1 0

Нахождение собственных значений и собственных векторов, построение графика
Доброго времени суток, форумчане! Помогите, пожалуйста, решить задачу в Матлаб. Условие следующее:...

Функция hmatrixevd из пакета alglib для поиска собственных значений и собственных векторов
Добрый день. Нужно найти собственные значения и вектора комплексной эрмитовой матрицы. Под решение...

Критерий остановки QR-алгоритма для поиска собственных чисел или векторов
Добрый день. Я пишу функцию, которая ищет собственные вектора через QR-алгоритм. Но нигде не могу...

	24.09.2015, 12:01

IT_Exp

Эксперт

87844 / 49110 / 22898

Регистрация: 17.06.2006

Сообщений: 92,604

24.09.2015, 12:01

Помогаю со студенческими работами здесь

Вычисление собственных значений и собственных векторов матрицы. Метод скалярных произведений
Курсовая работа!! Задание: "Вычисление собственных значений и собственных векторов матрицы. Метод...

Метод итерации нахождения собственных чисел и собственных векторов матрицы
Помогите пожалуйста

Поиск собственных чисел и собственных векторов матрицы 3 на 3 с фиксированным шагом
Здраввствуйте, дорогие форумчане. Задача: собственных чисел и собственных векторов матрицы 3 на...

Решение СДУ методом определения собственных значений и собственных векторов
Дана система уравнений: \begin{cases} & \dot{x}= 2x-y \\ & \dot{y}=x+2{e}^{t} \end{cases}...

Метод итераций: нахождение собственных векторов и собственных значений матрицы
Доброго времени суток. метод итерации нахождение собственных векторов и собственных значений...

Нахождение собственных чисел и собственных векторов степенным методом
Здравствуйте. Не могу понять что не так, вроде бы по теории все как и должно быть, но результаты не...

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:

Наверх