samp8_1.m

来自「数字信号处理的Matlab实现的随书光盘」· M 代码 · 共 11 行

%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）函数求得结果