画三维图形的时频分析的程序怎么解释？-面包板社区

http://www.chinavib.com/forum/thread-29541-1-9.html

这段程序怎么解释？

各位高手，这是本论坛里面的一个帖子里面其中的一段程序，是用来画三维图形的时频分析。程序如下：
%%短时傅立叶变换
%解析信号x=exp(j*pi*k*t.^2)
clear,clc,close all
figure(1)
k=4;T=5;
fc=k*T;
fs=3*fc;
Ts=1/fs;
N=T/Ts;
x=zeros(1,N);
t=0:N-1;
x=exp(j*k*pi*(t*Ts).^2);
% x=awgn(x,-3,'measured');
subplot(221)
plot(t*Ts,real(x))
X=fft(x);
X=fftshift(X);
subplot(222)
plot((t-N/2)*fs/N,abs(X))
Nw=20;
L=Nw/2;
Tn=(N-Nw)/L+1;
nfft=32;
TF=zeros(Tn,nfft);
for i=1:Tn
xw=x((i-1)*10+1:i*10+10);
temp=fft(xw,nfft);
temp=fftshift(temp);
TF(i,

=temp;
end
subplot(223)
fnew=((1:nfft)-nfft/2)*fs/nfft;
tnew=(1:Tn)*L*Ts;
[F,T]=meshgrid(fnew,tnew);
mesh(F,T,abs(TF))
subplot(224)
contour(F,T,abs(TF))
其中我对后面画三维图那里的程序不是很明白，比如：NW是什么意思？L又是什么意思，为什么要NW/2？Tn=（N-NW）/L+1又如何解释？nfft的取值有什么说法没？

-------------------------------------

songzy41：

这一段程序中的关键是进行短时傅里叶变换，然而构成一个三维谱图。在取数时每
段数据的窗函数是移动10个数据样点，而窗函数每次取20个数据样点。
Nw=20;             %窗函数长
L=Nw/2;          %窗函数每次移动的样点数
Tn=(N-Nw)/L+1;    %计算把数据x共分成多少段
nfft=32;          %FFT的长度
TF=zeros(Tn,nfft);  %将存放三维谱图，先清零
for i=1:Tn
xw=x((i-1)*10+1:i*10+10);  %取一段数据
temp=fft(xw,nfft);       %FFT变换
temp=fftshift(temp);    %频谱以0频为中心
TF(i,

=temp; %把谱图存放在TF中
end

---------------------------------------

楼上的高手，先谢谢你给出的解释！现在我要解决的就是我有一段采集的数据，采样频率为500Hz，我想取其中的5秒来分析，不知道我的程序改按照上面的程序怎么改？能帮帮我吗？万分感激！

我试着用上面的程序编了一下，出来个这样的图（附图1），是双边频的，怎么把他改成单边频的呢？还有别人的图是这样的（附图2），这个程序怎么改才能变成图2那样呢？
程序如下：
Nw=20;fs=500;t=5;
Ts=1/fs;N=t/Ts;
L=Nw/2;
Tn=(N-Nw)/L+1;
nfft=1024;
TF=zeros(Tn,nfft);
for i=1:Tn
xw=x((i-1)*10+1:i*10+10);
temp=fft(xw,nfft);
temp=fftshift(temp);
TF(i,

=temp;
end
figure(2);
fnew=((1:nfft)-nfft/2)*fs/nfft;
tnew=(1:Tn)*L*Ts;
[F,T]=meshgrid(fnew,tnew);
mesh(T,F,abs(TF));