fg_09_03.m

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

% fg_09_03   MSE learning curve plot.
m=1; u=.05; N=50; k=[0:N-1];
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(1,N)+(1-2*u).^(2*k)*(c0'*Phiff*c0);
mse20=mse(20)

% Plot the curve.
sp_fig(1,8,5);
set(gca,'fontsize',16);
set(gca,'position',[0.10    0.1365    0.8647    0.78])
line(k,mse,'color','r','marker','o');
title('50 iterations; u=0.05.');
xlabel('k'); ylabel('MSE'); grid on;