STFT存在频率偏移和幅值错误的问题

旗舰水自然

1. %%短时傅立叶变换
   
2. clear;
   
3. clc;
   
4. fs=1000;%采样频率
   
5. Ts=1/fs;%采样时间间隔
   
6. TT=1;%终止时间
   
7. t=0:Ts:TT;%采样时刻                    
   
8. N=fs*TT;%采样点
   
9. %signal
   
10. x=zeros(size(t));
    
11. x1=sin(2*pi*50*t);x2=sin(2*pi*200*t);
    
12. x=[x1(1:300),x2(301:end)];
    
13. figure(3);
    
14. plot(x);
    
15. title('信号波形图'); %图名
    
16. xlabel('时间'); %x轴
    
17. ylabel('振幅'); %y轴
    
18. %STFFT
    
19. Nw=51;              %窗函数长
    
20. L=Nw/2;             %窗函数每次移动的样点数
    
21. Tn=floor((N-Nw)/L)+1;      %计算把数据x共分成多少段
    
22. nfft=51;
    
23. TF=zeros(Tn,nfft);  %时频矩阵
    
24. for i=1:Tn
    
25. 
    
26. xw=x(25*i-24:25*i+26).*hamming(Nw)';  %加Hamming窗
    
27.     temp=fft(xw,nfft);         %FFT变换
    
28.     temp=fftshift(temp);       %频谱以0频为中心
    
29.     TF(i,:)=temp;              %把谱图存放在TF中
    
30. end
    
31. figure(4);
    
32. fnew=((1:nfft)-nfft/2)*fs/nfft;
    
33. tnew=(1:Tn)*L*Ts;
    
34. [F,T]=meshgrid(fnew,tnew);
    
35. mesh(T,F,abs(TF));
    
36. colorbar;
    
37. title('STFT时频图')
    
38. xlabel('时间  t/ms'); ylabel('频率  f/Hz');
    

复制代码

作出图之后，发现频率不是关于0Hz对称的，而且幅值也不对，不是除以N/2的关系。另外，如何把画只含正频率的图？谢谢！！