ui_traj.m

基于Matlab

% keep all variables local in scope
function [] = ui_traj();
 
% declare global variables
dfec_global; 
ui_private_global;

if get(h_traj,'Value')
   figure(fig_traj);

   % Real plots
   if (isreal(fahist) & isreal(fbhist) & (isreal(dhist)))
      if (spacing==1)
         if (Nd>0)
            m=2;
         else
            m=1;
         end;

         subplot(m,1,1);
         plot(fhist');
         zoom on;
         title('Forward Filter Trajectories');
         xlabel('reps');
         ylabel('f(k)');

         if (Nd>0)
            subplot(m,1,2);
            plot(dhist');
            zoom on;
            title('Feedback Filter Trajectories');
            xlabel('reps');
            ylabel('d(k)');
         end;

      else
         if (Nd>0)
            m=3;
         else
            m=2;
         end;

         subplot(m,1,1);
         plot(fahist');
         zoom on;
         title('Forward Filter (a) Trajectories');
         xlabel('reps');
         ylabel('fa(k)');
         
         subplot(m,1,2);
         plot(fbhist');
         zoom on;
         title('Forward Filter (b) Trajectories');
         xlabel('reps');
         ylabel('fb(k)');

         if (Nd>0)
            subplot(m,1,3);
            plot(dhist');
            zoom on;
            title('Feedback Filter Trajectories');
            xlabel('reps');
            ylabel('d(k)');
         end;
      end;
 
   % Complex plots 
   else
      if (spacing==1)
         if (Nd>0)
            m=2;
         else
            m=1;
         end;

         subplot(m,2,1);
         plot(real(fhist.'));
         zoom on;
         title('Real part: Forward Filter Trajectories');
         xlabel('reps');
         ylabel('Re f(k)');

         subplot(m,2,2);
         plot(imag(fhist.'));
         zoom on;
         title('Imag part: Forward Filter Trajectories');
         xlabel('reps');
         ylabel('Im f(k)');

         if (Nd>0)
            subplot(m,2,3);
            plot(real(dhist.'));
            zoom on;
            title('Real part: Feedback Filter Trajectories');
            xlabel('reps');
            ylabel('Re d(k)');

            subplot(224);
            plot(imag(dhist.'));
            zoom on;
            title('Imag part: Feedback Filter Trajectories');
            xlabel('reps');
            ylabel('Im d(k)');
         end;
         
      else
         if (Nd>0)
            m=3;
         else
            m=2;
         end;

         subplot(m,2,1);
         plot(real(fahist.'));
         zoom on;
         title('Real part: Forward Filter (a) Trajectories');
         xlabel('reps');
         ylabel('Re fa(k)');
 
         subplot(m,2,2);
         plot(imag(fahist.'));
         zoom on;
         title('Imag part: Forward Filter (a) Trajectories');
         xlabel('reps');
         ylabel('Im fa(k)');

         subplot(m,2,3);
         plot(real(fbhist.'));
         zoom on;
         title('Real part: Forward Filter (b) Trajectories');
         xlabel('reps');
         ylabel('Re fb(k)');

         subplot(m,2,4);
         plot(imag(fbhist.'));
         zoom on;
         title('Imag part: Forward Filter (b) Trajectories');
         xlabel('reps');
         ylabel('Im fb(k)');


         if (Nd>0)
            subplot(m,2,5);
            plot(real(dhist.'));
            zoom on;
            title('Real part: Feedback Filter Trajectories');
            xlabel('reps');
            ylabel('Re d(k)');

            subplot(m,2,6);
            plot(imag(dhist.'));
            zoom on;
            title('Imag part: Feedback Filter Trajectories');
            xlabel('reps');
            ylabel('Im d(k)');
         end;
      end;
   end;
else
   if ishandle(fig_traj), delete(fig_traj); end;
end;