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