Matlab音频变速不变调、变调不变速（Phase Vocoder）

其他文章有讲过一些方法，有OLA、WSOLA算法、LSEE-MSTFTM、Phase Vocoder（相位声码器）等等，但都是讲了个大概，没一个能说清楚的，代码就更没有了。找了一个国外的Phase Vocoder变速、变调的方法，原理如下：

http://www.guitarpitchshifter.com/algorithm.html

代码如下：

http://www.guitarpitchshifter.com/matlab.html

其中变调部分（升降n个半音），可以拿来直接运行。

变速可以通过变调代码的前面一大半部分实现，其中变速的倍数用t来表示，t作为函数的输入，函数中的step修改如下：

主函数后面去掉：

% Resample with linearinterpolation
outputTime = interp1((0:(length(outputTimeStretched)-1)),outputTimeStretched,(0:alpha:(length(outputTimeStretched)-1)),'linear');

% Return the result
outputVector = outputTime;

这样即可实现变速不变调和变调不变速。

具体原理，尤其是频率和相位部分真是太™难懂了，不如直接抄代码（逃