我在仿真混合源通过四阶累积量的DOA估计的过程中，出现了当远场源大于20度，就测不出现的现象。这篇论文给的信源设置中用到了归一化频率，具体如下：
混合源信号是等能量，相对独立的。并且形式为非高斯Sk(t) = e ^（ j * 2pi * fk * t +j * Φk ）
其中fk是归一化频率，Φk是相互独立的随机相位，均匀分布在[0，2π]上。
假设噪声是加性空间白复高斯噪声，第k个信号的SNR定义为SNR＝10log10（σ^2(k)/σ^2(n)）, 其中σ^2(k)和σ^2(n)是第k个源的功率和噪声。
下面是我编写的，烦请大家看看哪里出问题了
clear；
degToRad = pi/180; % 角度->弧度
radToDeg = 180/pi; % 弧度->角度
ERM1 = 11; % UCA阵元数量 2*5+1
Sensor = - (ERM1 - 1)/2 : 1 : (ERM1 - 1)/2;
snr = 15; % 信噪比
c = 3e8; %声速
d = 7.5e-2;%阵元间距
dd = 0.25;
% D = d * ERM1; %整个阵列的孔径
% dlambda = -(ERM1 - 1)/2 * dd :dd: (ERM1 - 1)/2 * dd;
N = 200; %快拍数
Kn = 1; % 近场源数 Kn<=M
Kf = 2; % 远场源数
K = Kn + Kf ; %总信源数
theta=[ -15 5 20]; %DOA
f = c/(4*d); %f=1e8
lambda = c / f;
% r = [3 * lambda 8 * lambda Inf]; %2
r = [3 * lambda Inf Inf ]; %1
fs = 3e9; % 采样频率
Ts = 1/fs;
t = Ts * [1:1:N];
%% Date model
Fk = zeros(K,1);
% Fk(:,:) = [0.21 0.22 0.23 ]; %归一化
Fk(:,:) = [2e8 1.5e8 2.5e8 ]; %非归一化
S = exp(j * 2 * pi * Fk * t +j * rand(K,1)*2*pi );
