shipinzuoye1.m

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

%求(1)的STFT & WVD:
clear all
fs=100;                          %采样率:满足nyquist采样定律
T=1;
N=fs*T;
X=ones(N,1);
%设定STFT参数:
h=boxcar(8);
t_interval=1;        %t_interval是窗函数的移动步长
fft_number=80;      %fft_number定频率轴的采样数目
%调用STFT函数:
st_ft=STFT(X,t_interval,fft_number,h);
%设定WVD参数:
fft_number=80; 
%调用WVD函数:
wvd_trf=WVD(X,fft_number);
%绘图:
subplot(2,1,1);
mesh((1:N),(0:fft_number-1),abs(st_ft));
xlabel('Time');ylabel('Frequence');zlabel('STFT');
subplot(2,1,2);
mesh((1:N),(0:fft_number-1),wvd_trf);
xlabel('Time');ylabel('Frequence');zlabel('WVD');
%wvd_trf=wvd_trf(fft_number/4:3*fft_number/4-1,:);
%subplot(2,1,2);mesh((1:N)',(fft_number/4:3*fft_number/4-1)',abs(wvd_trf));