为了克服MATLAB时频工具箱在计算大容量的地震数据的分频处理时，存在计算速度慢、内存不足以及ＦＦＴ的计算长度严重影响时频分辨率等问题，使其适应于实际海量地震数据的时频分析。对MATLAB时频工具箱中基于离散短时傅里叶变换的时频算法进行了改进，采用Fortran语言在微机上编程实现了该算法。理论实算表明，改进后的Fortran程序能够高效地计算大容量、短时窗地震数据的分频处理。