wht.m

给定入射信号角度分别为 -22度

%************ WVD_Hough变换 ***************

%使用说明:
%      wvd为要做Hough变换的wvd矩阵;
%      M为斜率采样点数; 
%      N为频率轴截距采样点数,一般取值范围为频率轴一个周期(m/2); 
%      得到的Hough变换的'截距'轴的对应值为初始频率;'斜率'轴的K值
%      对应的实际斜率值为k_interval*(-M/2+(K-1));

function  WVD_Hough=WHT(wvd,M,N)

[m,n]=size(wvd);
WVD_Hough=zeros(M,N);

k_interval=m/n/M ;
array1=ones(M,1)*(1:N);
array2=k_interval*(-M/2:M/2-1)'*ones(1,N);

rem_array=ones(M,N)*m;


for i=1:n
    arrayi=array1+i*array2;
    arrayi=abs(rem(arrayi,rem_array));
    arrayi=round(arrayi)+ones(M,N);

    array3=wvd(arrayi(:),i);
    WVD_Houghi=reshape(array3,M,N);
    WVD_Hough=WVD_Hough+ WVD_Houghi;
end