扫频信号的stft变换

%%y=chirp（t，f0，t1，f1）%功能：产生一个线性（频率随时间线性变化）信号，其时间轴设置由数组t定义。%时刻0的瞬间频率为f0，时刻t1的瞬间频率为f1.默认情况，f0=0Hz，t1=1，f1=100Hz。%格式：y=chirp（t，f0，t1，f1，‘method’）%功能：指定改变扫频的方法。可用的方法有‘linear’（线性扫频）、‘quadratic’（二次扫频）和‘logarithmic’（对数扫频）；%默认时为‘linear’。注意：对于对数扫频，必须有f1>f0。格式：y=chirp（t，f0，t1，f1，‘method’，phi）%功能：指定信号的初始相位为phi（单位为度），默认时phi=0%%specgram(x,nfft,fs,window,noverlap);?%[b,f,t]=specgram(data,nfft,Fs,window,numoverlap);??%画时频图,其中：nfft是fft的长度，越长的话，频域分辨率越高，但是，对于语音这种时变信号，不能过长，一般采样率16k或者8k的取1024或者512，要看采样率的高低。%fs是采样率,要比nfft小%window是指窗的长度，要小于或者等于nff，一般等于nff%numoverlap是指nfft减去步长，越大越好，但运算量越大。一般取nfft的3/4效果就比较好了。overlap是帧重叠的长度，要<=1/2*nff,取1/4*nff就可以了。%%线性扫频t1=50;%t1的变化会影响视频图的显示，t1是参数变化的关键t=0:0.001:t1;y=chirp(t,0,1,150,'linear');figureplot(t,y);axis([01-11]);title('线性扫频信号波形图');figure;nfft=128;%nfft是fft的长度，对频域分辨率高低有影响fs=100;%fs是采样率,要比nfft小n=2;%1/n*nfft表示的是帧重叠的长度spectrogram(y,nfft,fs,nfft,1/2*nfft,'yaxis')title('线性扫频信号时频图');%%指数扫频t1=10;%t1的变化会影响视频图的显示，t1是参数变化的关键t=0:0.001:t1;f0=10;f1=400;y=chirp(t,f0,10,f1,'logarithmic');figureplot(t,y);axis([01-11]);title('指数扫频信号波形图');figure;nfft=256;%nfft是fft的长度，对频域分辨率高低有影响fs=250;%fs是采样率,要比nfft小n=2;%1/n*nfft表示的是帧重叠的长度spectrogram(y,nfft,fs,nfft,1/2*nfft,'yaxis')title('指数扫频信号时频图');%%二次扫频t1=2;%t1的变化会影响视频图的显示，t1是参数变化的关键t=-t1:0.001:t1;y=chirp(t,100,1,200,'quadratic');figureplot(t,y);axis([-0.10.1-1.51.5]);title('二次扫频信号波形图');figure;nfft=256;%nfft是fft的长度，对频域分辨率高低有影响fs=250;%fs是采样率,要比nfft小n=2;%1/n*nfft表示的是帧重叠的长度spectrogram(y,nfft,fs,nfft,1/n*nfft,'yaxis')title('二次扫频信号时频图');程序运行结果：结论：1，当修改程序中的t1和nfft以及fs时，运行的结果都有所不同2，对于线性扫频，修改参数对图像的分辨率改变3，对于指数和二次扫频，修改参数会改变分辨率，以及图像的形状，如修改指数扫频的t1和nfft，nfft高的频率分辨就高，其结果如下图所示：问题：1，如何选择nfft，是否有一定规律，还是只能靠试；