计算短时窗口傅里叶变换

Question

有人能建议一种计算给定信号的短时窗口傅里叶变换的方法，也称为加窗周期变换信号的窗口傅里叶变换？请考虑此示例数据系列：

t = 1:365;
x = cos((2*pi)/12*t)+randn(size(t));

我理解该方法背后的原理，但尝试在matlab中实现它是另一回事。我这个例子不足以提供答案，对具有不同周期性的信号的任何解决方案都将非常感激。