短时傅立叶变换matlab仿真程序Word下载.docx

1、fN = 8000;n = 10000;t = 0:1/fN:（n-1）/fN;Generate signalssig1 = sin（2*pi*500*t）;sig2 = sin（2*pi*1000*t）/2;puls = zeros（size（t）;puls（2100） = 10;puls（2132） = 10;x = sig1+sig2+puls;window widthwin1 = 12.8*fN/1000; %12.8ms 102.4 pointswin2 = 6.4*fN/1000; %6.4ms 51.2 pointswin3 = 3.2*fN/1000; %3.2ms 25.6

2、pointsstepping lengthstepping1 = 16; stepping2 = 32;stepping3 = 64;nfft pointsnfft = min（256,length（x）;painting（Time Frequency Analysis）1.Hamming windows,different window width-figure（1）subplot（311）;specgram（x,nfft,fN,hanning（round（win3）;xlabel（time（ms）ylabel（frequency（kHz）title（STFT Time-Frequency

3、Map,Hamming windows,window width 3.2mssubplot（312）;specgram（x,nfft,fN,hanning（round（win2）;STFT Time-Frequency Map,Hamming windows,window width 6.4mssubplot（313）;specgram（x,nfft,fN,hanning（round（win1）;STFT Time-Frequency Map,Hamming windows,window width 12.8ms%-2.Gauss window ,different window width3

4、.Hamming windows,different stepping length figure（3）specgram（x,nfft,fN,hanning（round（win1）,stepping1）;STFT Time-Frequency Map,Hamming windows,window width12.8ms,stepping length 16specgram（x,nfft,fN,hanning（round（win1）,stepping2）;STFT Time-Frequency Map,Hamming windows,stepping length 32specgram（x,nf

5、ft,fN,hanning（round（win1）,stepping3）;STFT Time-Frequency Map,Hamming windows,window width12.8ms,stepping length 644.Gauss window ,different stepping length figure（4）specgram（x,nfft,fN,gausswin（round（win1）,16）;STFT Time-Frequency Map,Gauss window ,window width12.8ms,stepping length 16specgram（x,nfft,fN,gausswin（round（win1）,32）;STFT Time-Frequency Map,Gauss window ,window width12.8ms,stepping length 32specgram（x,nfft,fN,gausswin（round（win1）,64）;STFT Time-Frequency Map,Gauss window ,window width 12.8ms,stepping length 64