fg_09_05.m

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

% fg_09_05    Two MSE learning curve plots.
m=1; u=.05; K1=51; k1=[0:K1-1]';

% Ideal mse convergence.
phi=cos(2*pi*[0:12]/12);
Phiff=[phi(1) phi(2); phi(2) phi(1)];
phifd=[phi(m+1); phi(m+2)];
bopt=Phiff\phifd;
mse_min=phi(1)-phifd'*bopt;
c0=[1; -2.75]-bopt;
mse=mse_min*ones(K1,1)+(1-2*u).^(2*k1)*(c0'*Phiff*c0);

% LMS mse convergence.
b0=[1; -2.75]; u=.05;
K2=51; k2=[0:K2-1]';
d=sqrt(2)*sin(2*pi*k2/12);
f=[0; d(1:K2-1)];
e=d-lms_filter(b0,u,d,f);

% Plot the curves.
sp_fig(1,9,5.5);
set(gca,'fontsize',16);
line(k1,mse,'color','b','linestyle','--');
line(k2,e.^2,'color','r','marker','o');
title('50 iterations; u=0.05.');
xlabel('k'); ylabel('e_k^2'); grid on;