fg_10_03.m

英文书《Digital Signal Processing with Examples in MATLAB》附带的MATLAB实例

% fg_10_03    Predictive coding of v1.
v1=[1:7,6:-1:0,1:4,3:-1:1,2:7,6:-1:2];
N=2;
R=autocovar_mat(v1,N);
r=crosscovar([0,v1(1:31)],v1,N);
b=R\r;
e=round(v1-filter(b,1,v1));
Entropy_of_e=entropy(e)
Count_of_e=32*freq(e)
sp_fig(1);
axes('position',[.05,.2,.4,.4]);
set(gca,'fontname','times','fontsize',16);
line([0:31],v1,'marker','o','color','k'); grid on;
axis([0 32 0 7]);
xlabel('Sample (k)'); ylabel('x_k');
set(gca,'xtick',[0:31]);
set(gca,'xticklabel','0|||||5|||||10|||||15|||||20|||||25|||||30||');
set(gca,'ytick',[0:7]);
title('Signal');

axes('position',[.55,.2,.4,.4]);
set(gca,'fontname','times','fontsize',16);
line([0:31],e,'marker','o','color','k'); grid on;
axis([0 32 -2 2]);
xlabel('Sample (k)'); ylabel('e_k');
set(gca,'xtick',[0:31]);
set(gca,'xticklabel','0|||||5|||||10|||||15|||||20|||||25|||||30||');
title('Prediction error');