【算法】LPC共振峰检测

居然在没有线代常识的基础下把这东西搞懂了，虽然只是达到“能用”的程度。

LPC：Linear Prediction Coding，线性预测编码，用已有的信号预测随后的信号。通过最小二乘法获得最小预测误差（意义为声门激励）和脉冲响应（意义为声道滤波器）。（过程相当于解卷？）
主要运算是解矩阵方程，Levinson-Durbin算法可以极大地减少针对这种特殊矩阵的计算开销，其原理是将大规模矩阵运算拆成小一级的矩阵运算来解。

这是一个非常好的介绍：
http://www.emptyloop.com/technotes/A%20tutorial%20on%20linear%20prediction%20and%20Levinson-Durbin.pdf

代码（依赖于CVEDSP）：


对输出的脉冲响应进行FFT可得到频域响应的倒数。通过H(z) = 1/A(z)求得频域响应，相当于频谱包络。    GenerateHamming(Window, 1024);ApplyWindow(InterArray, PFrames[PD.Amount / 2].Data, Window, 0, 1024);   
Boost_FloatSet(LPCCoef, 0, 1024);   
LPC(LPCCoef, InterArray, 1024, 32);   
GenerateBoxCar(Window, 1024);   
GetMagnitudeFromWave(InterArray, LPCCoef, Window, 10);
实测中相比直接取频谱包络能获得更高精度，但仍不能解决共振峰重合问题；而且对低频信号效果较差。




在Rocaloid中应用价值：CVDBStudio自动共振峰标识。

你真是个天才，不夸张地说……

卧槽为啥这么有潜力的项目还没更多大神加入呢QAQ