宽带信号DOA估计ism算法 matlab

DOA估计（方向估计）是无线通信领域的一个重要问题，ISM算法（Improved Spatial Smoothing Method）是其中的一种常用算法。

以下是一份MATLAB代码，用于基于ISM算法的宽带信号DOA估计：

clear all; clc;
%% 参数设置
N = 8; % 阵元数目
lamda = 1; % 波长
d = lamda/2; % 阵元间距
theta = 30; % 信源角度
P = 100; % 信号样本数
SNR = 20; % 信噪比

%% 构造信号
S = zeros(N, P);
for k=1:P
    S(:,k) = exp(1j*pi*(0:N-1)'*sin(theta/180*pi)) + sqrt(1/2)*(randn(N,1)+1j*randn(N,1));
end

%% 构造噪声
noise = sqrt(1/2)*(randn(N,P)+1j*randn(N,P));

%% 接收信号
X = S + 10^(-SNR/20)*noise;

%% ISM算法
Rxx = X*X'/P;
Rss = zeros(N,N);
for k=-N+1:N-1
    Rss = Rss + diag(diag(Rxx,k),k);
end
Rss = Rss/(2*N-1);
Rnn = Rxx - Rss;
iRnn = inv(Rnn);
w = iRnn*ones(N,1);
w = w/sqrt(w'*iRnn*w);
theta_est = asin(angle(w)'*(0:N-1)'/(2*pi*d))*180/pi;

%% 结果显示
disp(['实际角度：', num2str(theta)]);
disp(['估计角度：', num2str(theta_est)]);

运行该脚本，可以得到实际角度和估计角度的输出结果。需要注意的是，该代码仅适用于单个信源的情况，若存在多个信源，则需要使用其他多信源DOA估计算法。