对多个窄带信号测向，给出matlab代码

以下是一段用于对多个窄带信号进行测向的MATLAB代码示例，其中使用了MUSIC算法：

% 生成多个窄带信号
theta = [20 40 60]; % 信号到达角度
N = 100; % 信号长度
M = 10; % 阵元数
d = 0.5; % 阵元间距
A = zeros(M,length(theta));
for i=1:length(theta)
    A(:,i) = exp(-1j*2*pi*d*(0:M-1)'*sind(theta(i)));
end
S = A*diag([1,1,1])*exp(1j*2*pi*rand(3,N)); % 生成3个信号
X = randn(M,N)+1j*randn(M,N); % 加入噪声

% MUSIC算法测向
R = X*X'/N; % 样本协方差矩阵
[E,D] = eig(R); % R的特征值分解
[~,I] = sort(diag(D),'descend'); % 排序特征值
En = E(:,I(length(theta)+1:end)); % 噪声子空间
theta_range = -90:0.1:90; % 测向范围
Pmusic = zeros(length(theta_range),1);
for i=1:length(theta_range)
    a = exp(-1j*2*pi*d*(0:M-1)'*sind(theta_range(i)));
    Pmusic(i) = 1/(a'*(En*En')*a);
end
Pmusic = 10*log10(abs(Pmusic)/max(abs(Pmusic))); % 转为dB

% 绘制测向图像
figure;
plot(theta_range,Pmusic);
xlabel('Angle (degree)');
ylabel('Spatial Spectrum (dB)');
title('MUSIC Algorithm for DOA Estimation');
grid on;

此代码中首先生成了三个到达角度为20度、40度和60度的窄带信号，并将其合成为一个信号矩阵。然后加入了噪声，并使用MUSIC算法进行信号测向。最后绘制了测向图像。需要注意的是，此代码中的测向范围为-90度到90度，可以根据实际情况进行调整。