扫频信号的stft变换.docVIP

%% 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’。注意：对于对数扫频，必须有f1f0。 格式：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); figure plot(t,y); axis([0 1 -1 1]); title(线性扫频信号波形图); figure; nfft=;% nfft是fft的长度，对频域分辨率高低有影响 fs=0;%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); figure plot(t,y); axis([0 1 -1 1]); 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); figure plot(t,y); axis([-0.1 0.1 -1.5 1.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，是否有一定规律，还是只能靠试； 2， 是否当时间选的足够长时，后段都将不遵循前段的规律（如指数，二次） 3，t的时间步长如何确定，当选择稍大时程序将无显示。 4，几种扫频信号的波形图都没什么大的区别，为什么？