get_hk_from_hk_abs.m

（本文是在论坛已有一篇文章《SIRP法相干相关K分布雷达杂波的建模与仿真》的基础上修改而来

function [Hk] = Get_Hk_From_Hk_Abs(Hk_Abs)

N = length(Hk_Abs);
Hk = zeros(1, N);

% N 为偶数时
if(mod(N, 2) == 0)
    for w1 = 0 : 1 : N/2 - 1
        Hk(1, w1+1) = Hk_Abs(1, w1+1) * exp((-j*2*pi/N) * w1 * (N-1)/2);
    end;
    Hk(N/2) = 0;
    for w1 = N/2+1 : 1 : N-1
        Hk(1, w1+1) = Hk_Abs(1, w1+1) * exp((j*2*pi/N) * (N-w1) * (N-1)/2);
    end;
else
     %  如果是奇数
     for w1 = 0 : 1 : (N-1)/2
        Hk(1, w1+1) = Hk_Abs(1, w1+1) * exp((-j*2*pi/N) * w1 * (N-1)/2);
    end;
    for w1 = (N+1)/2 : 1 : N-1
        Hk(1, w1+1) = Hk_Abs(1, w1+1) * exp((j*2*pi/N) * (N-w1) * (N-1)/2);
    end;
end;