yulewalk.asv

DSP 实现基于BURG算法的功率谱估计

fid=fopen('D:\homework\yulewalk\output_next.txt','r');
[A,count]=fscanf(fid,'%f');
%x( );
var2=A(1);
var=sqrt(var2);
for j=1:64,
    x(j)=0;
end
x(1)=1;
for j=1:10
    x(j+1)=A(2*j)+A(2*j+1)*i;
end
figure(1);
[h,w]=freqz(var,x,);
ampl=10*log10(abs(h).*abs(h))
plot(w,ampl);