nnd10nc.m

《神经网络设计》英文版的配套源代码

function nnd10nc(cmd,arg1,arg2,arg3)
% NND10NC Adaptive noise cancellation demonstration.
 
% Copyright 1994-2005 PWS Publishing Company and The MathWorks, Inc.
% $Revision: 1.6.2.2 $
% First Version, 8-31-95.

%==================================================================

% CONSTANTS
me = 'nnd10nc';
max_t = 0.5;

% DEFAULTS
if nargin == 0, cmd = ''; else cmd = lower(cmd); end

% FIND WINDOW IF IT EXISTS
fig = nndfgflg(me);
if length(get(fig,'children')) == 0, fig = 0; end
  
% GET WINDOW DATA IF IT EXISTS
if fig
  H = get(fig,'userdata');
  fig_axis = H(1);            % window axis
  desc_text = H(2);           % handle to first line of text sequence
  signal_axis = H(3);         % signal axis
  weight_axis = H(4);         % weight axis
  lr_text = H(5);             % text displaying value of learning rate
  lr_bar = H(6);              % learning rate slider
  mc_bar = H(7);              % text displaying value of momentum constant
  mc_text = H(8);             % momentum constant slider
  dc_line = H(9);             % dc line in signal axis
  w_ptr = H(10);              % Pointer to current weights
  w_cont_ptr = H(11);         % Pointer to contour lines in weight axis
  w_pt = H(12);               % Initial weight marker
  P_ptr = H(13);              % Pointer to network input data
  T_ptr = H(14);              % Pointer to network target data
  show_s = H(15);             % "Show signals" radio button
  show_e = H(16);             % "Show error" radio button
  w_line = H(17);             % Weight line in weight axis
  e_line = H(18);             % Estimated signal line in signal axis
  t_line = H(19);             % Target line in signal axis
  d_line = H(20);             % Difference line in signal axis
  S_ptr = H(21);              % Origonal signal
end

%==================================================================
% Activate the window.
%
% ME() or ME('')
%==================================================================

if strcmp(cmd,'')
  if fig
    figure(fig)
    set(fig,'visible','on')
  else
    feval(me,'init')
  end

%==================================================================
% Close the window.
%
% ME() or ME('')
%==================================================================

elseif strcmp(cmd,'close') & (fig)
  delete(fig)

%==================================================================
% Initialize the window.
%
% ME('init')
%==================================================================

elseif strcmp(cmd,'init') & (~fig)

  % CHECK FOR FILES
  % ODJ 12/23/07 Check for NN Toolbox removed
  %if ~nnfexist(me),return,end

  % CONSTANTS
  N = 3;                      % Samples per cycle
  f = 60;                     % Frequency of noise
  s = N*f;                    % Samples per second
  ts = s*max_t+1;             % Number of samples
  A = 0.1;                    % Gain on noise
  theta = pi/2;               % Phase shift of noise

  k = 0.2;                    % Signal amplitude
  %signal = k*rands(1,ts);     % Signal
  signal = k*(2*rand(1,ts)-1);     % Signal

  i = 1:ts;
  noise = 1.20*sin(2*pi*(i-1)/N);
  filtered_noise = A*1.20*sin (2*pi*(i-1)/N + theta);
  delayed_noise = [noise; 0 noise(1:length(noise)-1)];

  noisy_signal = signal + filtered_noise;

  w = [0 -2];
  time = [1:ts]/ts*max_t;     % Simulation time points

  % THE FUNCTION
  P = delayed_noise;
  T = noisy_signal(1:ts);
  A=2*P*P';
  d=-2*P*T';   
  c=T*T';

  % CONTOUR DATA
  xx = -2.1:0.15:2.1;
  yy = xx;
  [XX,YY] = meshgrid(xx,yy);
  F = (A(1,1)*XX.^2+(A(1,2)+A(2,1))*XX.*YY+A(2,2)*YY.^2)/2 +...
    d(1)*XX + d(2)*YY + c;

  % NEW DEMO FIGURE
  fig = nndemof2(me,'DESIGN','Adaptive Noise Cancellation','Chapter 10','');
  set(fig, ...
    'windowbuttondownfcn',nncallbk(me,'down'), ...
    'BackingStore','off',...
    'nextplot','add');
  H = get(fig,'userdata');
  fig_axis = H(1);
  desc_text = H(2);

  % ICON
  nndicon(10,458,363,'shadow')

  % SLIDE BARS
  lr = 0.2;
  text(180,140,'Learning Rate:',...
    'color',nndkblue,...
    'fontw','bold',...
    'fontsize',12,...
    'horizontalalignment','left')
  lr_text = text(340,140,num2str(lr),...
    'color',nndkblue,...
    'fontw','bold',...
    'fontsize',12,...
    'horizontalalignment','right');
  text(180,95,'0.1',...
    'color',nndkblue,...
    'fontw','bold',...
    'fontsize',12,...
    'horizontalalignment','left')
  text(340,95,'1.5',...
    'color',nndkblue,...
    'fontw','bold',...
    'fontsize',12,...
    'horizontalalignment','right');
  lr_bar = uicontrol(...
    'units','points',...
    'position',[180 110 160 16],...
    'style','slider',...
    'backg',nnltgray,...
    'callback',[me '(''lr'')'],...
    'min',0.1,...
    'max',1.5,...
    'value',lr);

  mc = 0.0;
  text(180,70,'Momentum:',...
    'color',nndkblue,...
    'fontw','bold',...
    'fontsize',12,...
    'horizontalalignment','left')
  mc_text = text(340,70,num2str(mc),...
    'color',nndkblue,...
    'fontw','bold',...
    'fontsize',12,...
    'horizontalalignment','right');
  text(180,25,'0.0',...
    'color',nndkblue,...
    'fontw','bold',...
    'fontsize',12,...
    'horizontalalignment','left')
  text(340,25,'1.0',...
    'color',nndkblue,...
    'fontw','bold',...
    'fontsize',12,...
    'horizontalalignment','right');
  mc_bar = uicontrol(...
    'units','points',...
    'position',[180 40 160 16],...
    'style','slider',...
    'backg',nnltgray,...
    'callback',[me '(''mc'')'],...
    'value',mc);

  % SIGNAL AXIS
  color_order = [nnred; nndkblue; nngreen; [1 0 1]];
  signal_axis = axes(...
    'units','points',...
    'position',[40 210 300 120],...
    'box','on', ...
    'color',nnltyell, ...
    'xcolor',nndkblue, ...
    'ycolor',nndkblue, ...
    'zcolor',nndkblue, ...
    'fontsize',10,...
    'nextplot','add',...
    'colororder',color_order,...
    'xlim',[0 max_t],...
    'ylim',[-2.1 2.1]);
  xlabel('Time');
  set(get(signal_axis,'xlabel'),'fontw','bold')
  ylabel('Amplitude');
  set(get(signal_axis,'ylabel'),'fontw','bold')
  title('Original (blue) and Estimated (red) Signals')
  set(get(signal_axis,'title'),...
    'color',nndkblue,...
    'fontw','bold',...
    'fontsize',12)
  d_line = line(time,time*0,...
    'visible','off',...
    'erasemode','none',...
    'color',nnred);
  t_line = line(time,signal,...
    'erasemode','none',...
    'color',nnltyell);
  e_line = line(time,time*0,...
    'erasemode','none',...
    'color',nnltyell);
  dc_line = plot3([0 max_t],[0 0],[1 1],'--',...
    'color',nndkblue,...
    'erasemode','none');
  set(get(signal_axis,'xlabel'),'fontw','bold','color',nndkblue)
  set(get(signal_axis,'ylabel'),'fontw','bold','color',nndkblue)
  set(get(signal_axis,'zlabel'),'fontw','bold','color',nndkblue)
  view(2)

  % WEIGHT AXIS
  weight_axis = axes(...
    'units','points',...
    'position',[40 40 120 120],...
    'box','on', ...
    'color',nnltyell, ...
    'xcolor',nndkblue, ...
    'ycolor',nndkblue, ...
    'zcolor',nndkblue, ...
    'fontsize',10,...
    'nextplot','add',...
    'colororder',color_order,...
    'xlim',[-2.1 2.1],...
    'ylim',[-2.1 2.1]);
  xlabel('W(1,1)');
  set(get(weight_axis,'xlabel'),...
    'fontw','bold')
  ylabel('W(1,2)');
  set(get(weight_axis,'ylabel'),...
    'fontw','bold')
  title('Adaptive Weights')
  set(get(weight_axis,'title'),...
    'color',nndkblue,...
    'fontw','bold',...
    'fontsize',12)
  [dummy,w_cont] = contour(xx,yy,F,10);
  for i=(w_cont')
    set(i,'erasemode','none');
  end
  w_pt = plot3(w(1),w(2),1,'.',...
    'color',nnred,...
    'erasemode','none',...
    'markersize',20);
  w_line = line(w(1),w(2),...
    'erasemode','none',...
    'color',nnred);
  set(get(weight_axis,'xlabel'),'fontw','bold','color',nndkblue)
  set(get(weight_axis,'ylabel'),'fontw','bold','color',nndkblue)
  set(get(weight_axis,'zlabel'),'fontw','bold','color',nndkblue)
  view(2)

  % RADIO BUTTONS
  show_s = uicontrol(...
    'units','points',...
    'position',[180 160 70 20],...
    'style','radiobutton',...
    'string','Signals',...
    'callback',[me '(''show_s'')'],...
    'background',nnltgray,...
    'value',1);
  show_e = uicontrol(...
    'units','points',...
    'position',[250 160 100 20],...