掃頻信號(hào)的stft變換

%%y=chirp（t，f0，t1，f1）%功能：產(chǎn)生一個(gè)線性（頻率隨時(shí)間線性變化）信號(hào)，其時(shí)間軸設(shè)置由數(shù)組t定義。%時(shí)刻0的瞬間頻率為f0，時(shí)刻t1的瞬間頻率為f1.默認(rèn)情況，f0=0Hz，t1=1，f1=100Hz。%格式：y=chirp（t，f0，t1，f1，‘method’）%功能：指定改變掃頻的方法?？捎玫姆椒ㄓ小甽inear’（線性掃頻）、‘quadratic’（二次掃頻）和‘logarithmic’（對(duì)數(shù)掃頻）；%默認(rèn)時(shí)為‘linear’。注意：對(duì)于對(duì)數(shù)掃頻，必須有f1>f0。格式：y=chirp（t，f0，t1，f1，‘method’，phi）%功能：指定信號(hào)的初始相位為phi（單位為度），默認(rèn)時(shí)phi=0%%specgram(x,nfft,fs,window,noverlap);?%[b,f,t]=specgram(data,nfft,Fs,window,numoverlap);??%畫時(shí)頻圖,其中：nfft是fft的長(zhǎng)度，越長(zhǎng)的話，頻域分辨率越高，但是，對(duì)于語(yǔ)音這種時(shí)變信號(hào)，不能過(guò)長(zhǎng)，一般采樣率16k或者8k的取1024或者512，要看采樣率的高低。%fs是采樣率,要比nfft小%window是指窗的長(zhǎng)度，要小于或者等于nff，一般等于nff%numoverlap是指nfft減去步長(zhǎng)，越大越好，但運(yùn)算量越大。一般取nfft的3/4效果就比較好了。overlap是幀重疊的長(zhǎng)度，要<=1/2*nff,取1/4*nff就可以了。%%線性掃頻t1=50;%t1的變化會(huì)影響視頻圖的顯示，t1是參數(shù)變化的關(guān)鍵t=0:0.001:t1;y=chirp(t,0,1,150,'linear');figureplot(t,y);axis([01-11]);title('線性掃頻信號(hào)波形圖');figure;nfft=128;%nfft是fft的長(zhǎng)度，對(duì)頻域分辨率高低有影響fs=100;%fs是采樣率,要比nfft小n=2;%1/n*nfft表示的是幀重疊的長(zhǎng)度spectrogram(y,nfft,fs,nfft,1/2*nfft,'yaxis')title('線性掃頻信號(hào)時(shí)頻圖');%%指數(shù)掃頻t1=10;%t1的變化會(huì)影響視頻圖的顯示，t1是參數(shù)變化的關(guān)鍵t=0:0.001:t1;f0=10;f1=400;y=chirp(t,f0,10,f1,'logarithmic');figureplot(t,y);axis([01-11]);title('指數(shù)掃頻信號(hào)波形圖');figure;nfft=256;%nfft是fft的長(zhǎng)度，對(duì)頻域分辨率高低有影響fs=250;%fs是采樣率,要比nfft小n=2;%1/n*nfft表示的是幀重疊的長(zhǎng)度spectrogram(y,nfft,fs,nfft,1/2*nfft,'yaxis')title('指數(shù)掃頻信號(hào)時(shí)頻圖');%%二次掃頻t1=2;%t1的變化會(huì)影響視頻圖的顯示，t1是參數(shù)變化的關(guān)鍵t=-t1:0.001:t1;y=chirp(t,100,1,200,'quadratic');figureplot(t,y);axis([-0.10.1-1.51.5]);title('二次掃頻信號(hào)波形圖');figure;nfft=256;%nfft是fft的長(zhǎng)度，對(duì)頻域分辨率高低有影響fs=250;%fs是采樣率,要比nfft小n=2;%1/n*nfft表示的是幀重疊的長(zhǎng)度spectrogram(y,nfft,fs,nfft,1/n*nfft,'yaxis')title('二次掃頻信號(hào)時(shí)頻圖');程序運(yùn)行結(jié)果：結(jié)論：1，當(dāng)修改程序中的t1和nfft以及fs時(shí)，運(yùn)行的結(jié)果都有所不同2，對(duì)于線性掃頻，修改參數(shù)對(duì)圖像的分辨率改變3，對(duì)于指數(shù)和二次掃頻，修改參數(shù)會(huì)改變分辨率，以及圖像的形狀，如修改指數(shù)掃頻的t1和nfft，nfft高的頻率分辨就高，其結(jié)果如下圖所示：?jiǎn)栴}：1，如何選擇nfft，是否有一定規(guī)律，還是只能靠試；