shipinzuoye2.m

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

%求信号(2)的STFT分布
f1=10;f2=25;f3=40;               %定义3个分量的模拟频率fa(且fd=fa/fft_number)stft中的频率轴是fft的k，即有：
                                 %频率轴是fft的k，即有：(2*pi/fft_number)*k=2*pi*fa/fft_number,=>k=fa
fs=100;                          %采样率:满足nyquist采样定律
T=1;
N=fs*T;
t1=(0:1/N:T/4-1/N)';
t2=(T/4:1/N:T/2-1/N)';
t3=(T/2:1/N:T-1/N)';
X1=exp(j*2*pi*f1*t1);
X2=exp(j*2*pi*f2*t2);
X3=exp(j*2*pi*f3*t3);
%设定STFT参数:
X=[X1;X2;X3];
h=boxcar(10);
t_interval=1;        %t_interval是窗函数的移动步长
fft_number=N;      %fft_number定频率轴的采样数目
%调用STFT函数:
st_ft=STFT(X,t_interval,fft_number,h);
%设定WVD参数:
fft_number=N; 
%调用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),abs(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));