wvd.m

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

%*********************WVD函数***************************
%  使用说明：
%     输入变量X表示是信号数据（列向量）；
%     fft_number是做FFT使用的离散频率点数，可以用来决定频率采样的精度。
%     生成的wvd数据，行表示时间，列表示频率：即每个时刻的fft放在列中


function  wvd_trf=WVD(X,fft_number)

N=length(X);
wvd_trf=zeros(fft_number,N);

f=(0:fft_number-1);
array1=(1:N)'*f;
w=exp(-j*4*pi/fft_number).^array1;
for i=1:N
    X1=[zeros(N+1,1);X;zeros(N+1,1)];
    X2=[zeros(2*(i),1);X(N:-1:1,1);zeros(2*N+2-2*(i),1)];
       if i<=(N/2)
           Y=X1.*conj(X2);
           Y1=Y((N+2):(N+1+N),1);
       else 
           Y=X1.*conj(X2);
           Y1=Y((1+2*i):(N+2*i),1);
       end
    wvd_trf(:,i)=2*(Y1.'*w).';
    %temp(:,i)=Y1;
    %Y2=fft(Y1,fft_number);
    %Y3=[Y2;Y2];
    %Y4=Y3(2:2:2*fft_number);
    %wvd_trf(:,i)=Y2;
end

wvd_trf=abs(wvd_trf);