esprit算法Matlab实现程序

以下是Esprit算法的Matlab实现程序，仅供参考：

function [theta, A] = esprit(signal, M,)
% signal: 输入信号矩阵，每列为一个传感器的观测值
 M: 估计的信号数
%: 传感器间距离
% theta: 估计出的角度值
% A: 估计出的信号幅值

[N, K] = size(signal);
R = (signal * signal') / K; % 信号协方差矩阵

% 对信号协方差阵进行特征值分解
[V, D] = eig(R);
[~, index] = sort(diag(D), 'descend');
V = V(:, index);

% 构造延迟矩阵
P = zeros(M, N - M + 1);
for i = 1:M
    P(i, :) = diag(V(:, i:N-M+i)' * V(:, 1:N-M+1-i));
end

% 对延迟矩阵进行SVD分解
[U, S, V] = svd(P);
U = U(:, 1:M-1);

% 对U进行特征值分解
[~, index] = sort(abs(fft(U)), 'descend');
theta = angle(index(1:M-1)) / d;

% 计算信号幅值
A = zeros(M, 1);
for i = 1:M
    A(i) = sqrt(sum(abs(signal(i, :)).^2) / K);
end

使用方法：

假设有3个传感器，观测到的信号存储在`signal`矩阵中，每一列为一个传感器的观测值。执行以下命令即可得到估计出的信号角度和幅值：

[theta, A] = esprit(signal, 2, 1);

其中`2`表示估计的信号数，`1`表示传感器间距离。