用c++写的PLC源代码

close all
clear
clc
[xfs]=wavread（‘speech.wav‘）;%读入数据
% 预加重滤波器
xx=double（x）;
y=filter（[1 -0.9495]1xx）;
N=160;
y1=y（1:N）;
w1=hamming（N）;
y2=（y1.*w1）‘;%加窗 取一帧数据
p=10;%预测阶数
%首先求自相关函数
r=zeros（1p+1）;
for k=1:p+1
sum=0;
for m=1:N+1-k
sum=sum+y2（m）.*y2（m-1+k）‘;
end
r（k）=sum;
end
%根据durbin算法求线性预测系数
k=zeros（1p）;
k（1）=r（2）/r（1）;
a=zeros（pp）;
a（11）=k（1）;
e=zeros（1p）;
e（1）=（1-k（1）^2）*r（1）;
%递推过程
for i=2:p
c=zeros（1i）;
sum=0;
for j=1:i-1
sum=sum+（a（i-1j）.*r（i+1-j））;