问短时傅立叶变换的问题
%%短时傅立叶变换,解析信号x=Y,Y为0.02s内1024个采样数据。figure(1)
T=0.02;
fs=1024/0.02;
Ts=1/fs;
N=T/Ts;
t=0:N-1;
x=Y;
subplot(221)
plot(t*Ts,x)
X=fft(x);
y=abs(X)
subplot(222)
plot(t(1:200)*fs/N,y(1:200))
Nw=20; %14
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; %27
=meshgrid(fnew,tnew);
mesh(F,T,abs(TF))
subplot(224)
contour(T,F,abs(TF));
最近要用到短时傅立叶变换，但对这方面的知识我还不是很了解，有人能解释第14行到27行的含义吗？谢谢
页:
[1]