fg_08_08.m

来自「英文书《Digital Signal Processing with Examp」· M 代码 · 共 27 行

% fg_08_08    3-d MSE plot
N=20; m=1;
phi=cos(2*pi*[0:12]/12);
Phiff=[phi(1) phi(2); phi(2) phi(1)]
phifd=[phi(m+1); phi(m+2)]
b=Phiff\phifd
min_mse=phi(1)-phifd'*b

db=1.99;
b0=linspace(b(1)-db,b(1)+db,N);
b1=linspace(b(2)-db,b(2)+db,N);
mse=zeros(N,N);
for i=1:N,
   for j=1:N,
      mse(i,j)=1+[b0(i) b1(j)]*Phiff*[b0(i); b1(j)]-2*[b0(i) b1(j)]*phifd;
   end
end
sp_fig(1,9,5.5);
set(gca,'position',[.2 .2 .75 .75]);
set(gca,'fontsize',16); 
surf(b0,b1,mse,'CDataMapping','direct'); grid on;
h=line([4 1.732],[-1 -5],[0 0],'color','r');
set(h,'linestyle','none','marker','*','markersize',20);
xlabel('b_0'); ylabel('b_1'); zlabel('MSE');
view(30,15);
colormap hsv;
caxis([-2 4]);