本科毕设做雷达微动的时频分析时,最常用的就是STFT了。matlab已经有了时频分析工具箱,调用里面的tfrstft函数即可轻松实现。今天看了源代码感觉写的比较复杂,于是结合微动,尝试自己实现一下。
STFT是在信号时域加窗函数,并通过移动窗函数进行FFT以获得信号的时变频率特性。可见,窗函数长度越短,对信号信号的时域划分越细,时域分辨性越好,但是频率分辨性越差。因此要合理选择窗函数长度。
代码实现的思路:先求
这个认为姑且认为是关于x与g的瞬时相关函数吧。其中要注意,和都是时间范围,但却是两个变量。
第二步就是求瞬时相关函数关于