空间谱测向：利用线阵和面阵实现同频同时多个窄带信号测向。阵元数、信号频率及测向方法任选。给出MATLAB程序代码。使用MATLAB R2020a

以下是一种实现上述要求的MATLAB程序代码，使用了线阵和空间谱测向方法：

%% 参数设置
M = 8; % 阵元数
f = 2e9; % 信号频率
c = 3e8; % 光速
lambda = c/f; % 波长
d = lambda/2; % 阵元间距
theta = [-30 0 15]; % 待测方向角
P = length(theta); % 待测信号数
N = 512; % FFT点数

%% 生成信号数据
rng(1); % 随机数种子
s = randn(N, P); % 生成P个长度为N的随机信号

%% 生成阵列数据
x = zeros(M, N, P); % 频域数据矩阵
for p = 1:P
    for m = 1:M
        x(m,:,p) = exp(-1j*2*pi*f*d*(m-1)*sind(theta(p))/c); % 生成阵列输出数据
    end
end

%% 空间谱测向
Rxx = zeros(M, M); % 信号协方差矩阵
for p = 1:P
    Rxx = Rxx + x(:,:,p)*x(:,:,p)'; % 对所有信号的协方差矩阵求和
end
Rxx = Rxx/P; % 对协方差矩阵求平均
w = beamscan(Rxx, linspace(-90, 90, 181)); % 空间谱测向
w = w/max(abs(w)); % 归一化

%% 绘图
figure;
plot(linspace(-90, 90, 181), abs(w), 'LineWidth', 2);
grid on;
title('空间谱测向');
xlabel('方向角（度）');
ylabel('归一化响应');

上述程序中，先生成了P个长度为N的随机信号，然后通过阵列输出数据的公式生成了阵列数据，其中使用了线阵和待测方向角theta。接着，通过对所有信号的协方差矩阵求和并求平均得到了信号协方差矩阵Rxx，最后使用空间谱测向方法得到了阵列的方向图，并进行了归一化处理。绘图部分使用MATLAB自带的plot函数绘制了方向图。

注意：本程序仅供参考，实际应用可能需要根据具体情况进行修改和调整。