Matlab 短时窗傅里叶变换的计算

有人能提出一种计算给定信号的短时加窗傅里叶变换（也称为加窗傅里叶变换）的方法吗？考虑这个例子数据系列：

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

---

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

这个文件交换解决方案怎么样。。。那你呢？