%Samp8_1
randn('state',0);   %設(shè)置隨機函數(shù)的狀態(tài)
a=[1 0.1 0.1 0.1 0.1];   %濾波器分母多項式系數(shù)
x=impz(1,a,10)+randn(10,1)/20;  %求得帶噪聲的濾波器脈沖響應(yīng)
%采用第一種方法
r=xcorr(x);            %求信號x的自相關(guān)r，參考第9章的該函數(shù)用法。
r(1:length(x)-1)=[];%該序列的前面部分(相當(dāng)于x的長度-1)受到邊界的影響，扣除
aa=levinson(r,4)    %用Levison-Durbin遞推算法采用自相關(guān)序列求出系數(shù)aa
%采用第二種方法
a1=lpc(x,4)    %直接調(diào)用lpc（Linear Predictor Coefficients）函數(shù)求得結(jié)果