matlab反幂法迭代矩阵

反幂法迭代矩阵是一种求解特征值和特征向量的方法，具体步骤如下：

1. 输入矩阵A和初始向量x0，设A的特征值均为实数且按绝对值大小递减排列，且特征值λ1 > λ2 > … > λn。

2. 计算逆矩阵B = (A - μI)^(-1)，其中μ为猜测的特征值。

3. 迭代计算x(k+1) = Bx(k)，直到满足收敛条件，即|x(k+1) - x(k)| < ε，其中ε为预设的误差。

4. 计算特征值λ = 1/α + μ，其中α为x(n+1)与x(n)的内积。

5. 计算特征向量v = x(n+1) / ||x(n+1)||。

以下是MATLAB代码示例：

function [lambda, v] = inverse_power_method(A, x0, mu, eps, max_iter)
% A: 待求特征值和特征向量的矩阵
% x0: 初始向量
% mu: 猜测的特征值
% eps: 收敛误差
% max_iter: 最大迭代次数
% lambda: 特征值
% v: 特征向量

n = size(A, 1);
x = x0 / norm(x0);
B = inv(A - mu * eye(n));

for k = 1:max_iter
    x_old = x;
    x = B * x;
    x = x / norm(x);
    if norm(x - x_old) < eps
        alpha = x' * x_old;
        lambda = 1 / alpha + mu;
        v = x;
        return;
    end
end

error('无法收敛');

注意：反幂法只能求解矩阵A的逆矩阵存在的情况下的特征值和特征向量。如果A的特征值中有重复值，则反幂法可能会失效。如果A的特征值中有复数，则需要使用复数版本的反幂法。