esprit算法matlab程序

以下是Esprit算法的Matlab程序示例：

% 定义信号和参数
N = 100; % 信号长度
M = 4; % 阵元数
d = 0.5; % 阵元间距
theta = [20, 40]; % 信号入射角
w = 2 * pi * [0.1, 0.15]; % 信号频率
snr = 10; % 信噪比

% 生成信号
A = exp(1i * (0:M-1)' * d * cosd(theta) * w);
s = A * (randn(length(w), N) + 1i * randn(length(w), N)) * 10^(-snr/20);

% Esprit算法
Rxx = s * s' / N;
[U, S, V] = svd(Rxx);
U1 = U(:, 1:M-1);
U2 = U(:, 2:M);
P = U1 * U2' / (U1 * U1');
phi = eig(P);
theta_hat = acosd(phi) / pi * 180;

% 显示结果
disp(['真实入射角：', num2str(theta)]);
disp(['估计入射角：', num2str(theta_hat)]);

该程序实现了一个4阵元线性阵列接收到两个信号的Esprit算法，其中`N`为信号长度，`M`为阵元数，`d`为阵元间距，`theta`为信号入射角，`w`为信号频率，`snr`为信噪比。程序首先生成信号，然后使用Esprit算法估计信号入射角，并输出真实入射角和估计入射角。