Matlab实现短时傅里叶变换

实现FFT点数可调、步进可调的短时傅里叶变换；实现加窗功能，窗类型可选。
代码如下
Fs=1000; %采样率1KHz
Ts=1/Fs; %点时间间隔0.001s
L=2048; %2048点数
t=(0:L-1)Ts; %采的时间为20480.001=2.048s
%t1=(0:L/4-13)Ts; %500个点
%t2=(0:L/4+35)Ts; %548个点
t1=(0:L/4-13)Ts; %512-13=499;500个点
t2=(0:L/4+35)Ts; %512+35=547；548个点
%y=sin(502pit)+1; %信号周期:T1=50Hz，T2=120Hz;幅度范围为-2~2;
%y=sin(502pit)+sin(1202pit)+2; %输出无符号十进制,范围0~4
%x = [2cos(2 * pi * 10 * t1+pi/2),2cos(2 * pi * 20 * t1+pi/2),2cos(2 * pi * 30 * t1+pi/2),2cos(2 * pi * 40 * t1+pi/2)]+2;%输出无符号十进制,范围0~4；2048个点
x = [2sin(2 * pi * 10 * t1),2sin(2 * pi * 20 * t1),2sin(2 * pi * 30 * t1),2sin(2 * pi * 40 * t2)]+2;
y = x1023/8192;
sinwave=round((y2^12))