VPF::[Matlab] Собственные числа матрицы

Rudolf

Дата 2.12.2007, 18:33 (ссылка)

(нет голосов)

Загрузка ...

Новичок

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

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

Вобщем мне надо написать функции нахождения максимального и минимального собственных чисел матрицы степенным и градиентным методом. Максимальное собственное число степенным методом я нашел:

Код


function [lambda,X,k]=eig_powermax(A,epsilon,max_k) 
% Вход     - A - матрица n x n 
%          - epsilon - допустимая погрешность 
%          - max_k - максимальное число итераций 
% Выход    - lambda - максимальное по модулю собственное значение 
%          - X - соответствующий собственный вектор 

% инициализация параметров 

k=0; % номер итерации 
[m,n]=size(A); % размерность матрицы 
if m~=n 
    error('Матрица не квадратная') 
end 

X=ones(n,1); % единичное начальное приближение 

while 1 
    k=k+1;% счетчик итераций 
    if k>max_k 
       error('Число итераций больше заданного') 
    end 
    Y=A*X; 
    if k>=2 
        lambda_s=lambda; 
    end 
    lambda=(Y'*X)/(X'*X); 
    % нормализация Y 
    X=(1/norm(Y,inf))*Y; 
    % проверка условия окончания 
    if k>=2 
        if abs(lambda-lambda_s)/abs(lambda)<=epsilon 
            break 
        end 
    end 
end % конец цикла 
return

А минимальное не получается....как я понял там надо решить систему
A*Y^(k+1)=X^(k), после чего получится максимальное собственное число м-цы A^-1, а его обратное значение будет минимальным у м-цы A. Я пишу:

Код


function [lambda,X,k]=eig_powermin(A,epsilon,max_k) 
% Вход     - A - матрица n x n 
%          - epsilon - допустимая погрешность 
%          - max_k - максимальное число итераций 
% Выход    - lambda - минимальное по модулю собственное значение 
%          - X - соответствующий собственный вектор 

% инициализация параметров 

k=0; % номер итерации 
[m,n]=size(A); % размерность матрицы 
if m~=n 
    error('Матрица не квадратная') 
end 

X=ones(n,1); % единичное начальное приближение 

while 1 
    k=k+1;% счетчик итераций 
    if k>max_k 
       error('Число итераций больше заданного') 
    end 
    Y=X/A; 
    if k>=2 
        lambda_s=lambda; 
    end 
    lambda=(Y'*X)/(X'*X); 
    % нормализация Y 
    X=(1/norm(Y,inf))*Y; 
    % проверка условия окончания 
    if k>=2 
        if abs(lambda-lambda_s)/abs(lambda)<=epsilon 
            break 
        end 
    end 
end % конец цикла 
lambda=1/lambda; 
return

Матлаб ругается...Подскажите пожалуйста где ошибка.

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