ui_mse.m

基于Matlab

% this computes MMSE designs 
%
% function [] = ui_mse();
 
% declare global variables
dfec_global; 
ui_private_global;

if get(h_mse,'Value') 
   newparams;
   figure(fig_mse);

   % Only varying parameter is delay
   if isempty(Nfd) 
      if (delay==-1)
         cursor=find(abs(c)==max(abs(c)));
         semilogy(cursor,MSE); hold on;
         semilogy(cursor,MSE,'o'); hold off;
         grid
         xlabel('delay');
         ylabel('log MSE');
      else
         semilogy(delay,MSE); hold on;
         semilogy(delay,MSE,'o'); hold off;
         grid
         xlabel('delay');
         ylabel('log MSE');
      end;

   % Only varying parameter is equalizer length
   elseif length(delay)==1 
      semilogy(MSE); hold on;
      semilogy(MSE,'o'); hold off;
      grid
      xlabel('Nf');
      ylabel('log MSE');

   % Both length & delay vary
   else
      surfl(log(MSE));
      shading interp
      colormap(bone);
      xlabel('N_f');
      ylabel('delay');
      zlabel('log MSE');
      hold on;
      X=[0 0 Nfd 0];
      if (spacing==1)
         Y=[ceil(Nc-1) max(delay) max(delay) ceil(Nc-1)];
      elseif (spacing==1/2)
         Y=[ceil(Nc/2) max(delay) max(delay) ceil(Nc/2)];
      end;
      Z=log(min(min(MSE)))*[1 1 1 1];
      fill3(X,Y,Z,'black');
      hold off;
      rotate3d on
   end;
   title(['Best: Nf=' num2str(Nf) ', Nd=' num2str(Nd) ...
           ', \delta=' num2str(delta) ', SNR=' num2str(max(SNR))]); 
else
   if ishandle(fig_mse), delete(fig_mse); end;
end;