fg_08_12.m

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

% fg_08_12    TSEmin versus N for different delays.
K=100; k=0:K-1;
sp_fig(1,9,5.5);
set(gca,'fontsize',16);
clr=['r' 'g' 'b'];
for m=1:3,
% Signals.
   s=[1, zeros(1,K-1)];
   f=exp(-.1*k).*sin(.1*k*pi);
   d=[zeros(1,m), s(1:K-m)];
% Least-squares design for each m using covariance.
   TSEmin=ones(1,6);
   for N=1:5,
      Rff=autocovar_mat(f,N);
      rfd=crosscovar(f,d,N);
      b=Rff\rfd;
      TSEmin(N+1)=d*d'+b'*Rff*b-2*b'*rfd;
   end
   line([0:5],TSEmin,'marker','o','color',clr(m),'linewidth',1.2);
end
axis([0 5 0 1]); grid on;
xlabel('Number of weights (N)'); ylabel('TSE_{min}');
text(1.9,.5,'delay=1','fontsize',16);
text(2.9,.5,'delay=2','fontsize',16);
text(3.9,.5,'delay=3','fontsize',16);