samp8_4.m

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

%Samp8_4
randn('seed',0);   %设置随机数种子
x=impz(1,[1 0.1 0.1 0.1],10)+randn(10,1)/20;  %
a1=lpc(x,3);    %运用lpc函数建立传递函数系数
[b2,a2]=prony(x,3,3); %运用prony函数建立传递函数系数
[b3,a3]=stmcb(x,3,3); %运用stmcb函数建立传递函数系数
%计算误差
disp('函数LPC输出:')
err1=x-impz(1,a1,10); %x与lpc函数确定模型脉冲响应的每个数的误差
sumerr1=sum(err1.^2)  %x与lpc函数确定模型脉冲响应的总误差
disp('函数PRONY输出:')
err2=x-impz(b2,a2,10); %x与prony函数确定模型脉冲响应的每个数误差
sumerr2=sum(err2.^2) %x与prony函数确定模型脉冲响应的总误差
disp('函数STCMB输出:')
err3=x-impz(b3,a3,10); %x与stmcb函数确定模型脉冲响应的每个数误差
sumerr3=sum(err3.^2) %x与stmcbc函数确定模型脉冲响应的总误差