English
- function e=exlpc(x,p);
- %线性预测编码LPC的实现，x为要编码的原始数列，P为x内的周期。没有再对x分块，
- %x就是一块
- %出自书《小波与傅里叶分析基础》的第0章最后一题
- %v=[a1,a2,...,ap];
- %y=[x(P+1),....,x(N)];
- %%%%%%%%%%%%%%%%%%%%%%%%注意：当p大于等于真实周期时，准确率为100%。
- %%%%%%%%%%%%%%%%%%%当n小于真实周期时，准确率急剧下降，越小越不准确
- %%%%%%%%%%%%%%%%%%%真实周期处是个临界点
- n=numel(x);
- y=x((p+1):n)';
- for i=1:(n-p)
- z(i,:)=x((i+p-1):-1:i);
- end
- v=(z'*z)^(-1)*z'*y;
- e=y-z*v;
