freqgui.m

trabajos en matlab para verificar son muy buenos

function freqgui(mode)
%FREQGUI Graphical User Interface Program for analog frequency response.


% ADSP Toolbox: Version 2.0 
% For use with "Analog and Digital Signal Processing", 2nd Ed.
% Published by PWS Publishing Co.
%
% Ashok Ambardar, EE Dept. MTU, Houghton, MI 49931, USA
% http://www.ee.mtu/faculty/akambard.html
% e-mail: [e-mail]
% Copyright (c) 1998



%v=matverch;
%if v<5,
%errordlg('This gui requires MATLAB 5.x and above','Version Error');
%return,
%end


if nargin == 0, mode = 0; end

if mode ~= 0,
   f = gcf;
   ui = get(f,'userdata');
   num_edt = ui(3);
   den_edt = ui(5);
   accept_btn = ui(6);
   plot_btn   = ui(17);
   autoplot   = ui(16);
   plot_pop   = ui(13);
   sig_axs    = ui(19);
   opr_axs = ui(29);
   curs_btn = ui(25);
   curs_dsp = ui(21:24);
end

if mode == 0,
   grey = [0.5,0.5,0.5];
   ltgrey = grey*1.5;
   dkgrey = grey*0.5;
   scrn_size = get(0,'screensize');
   fig_w = 620; fig_h = 440;
   flft = (scrn_size(3) - fig_w)/2;
   fbtm = (scrn_size(4) - fig_h)/2;

%   ver_str = version;
ver4=0;   %ver4    = length(findstr('4.',ver_str));
   if ver4,
      int_str = 'yes';
      f = figure('pos',[flft,fbtm,fig_w,fig_h],'menubar','none',...
            'name','Frequency Response','numbertitle','off',...
            'resize','off','vis','off','interruptible',int_str,...
            'color',dkgrey,'defaultaxesfontsize',10,...
            'defaulttextfontsize',10);
   else
      int_str = 'on';
      eval('f = colordef(''new'',''none'');');
      set(f,'pos',[flft,fbtm,fig_w,fig_h],'menubar','none',...
           'name','Frequency Response','numbertitle','off',...
           'resize','off','vis','off','interruptible',int_str,...
           'color',dkgrey,'defaultaxesfontsize',8,...
           'defaultuicontrolback',[0.5,0.5,0.5],...
           'defaulttextfontsize',8);
   end


%%% OPTIONS: %
   uim(1) = uimenu('label','OPTIONS');
   uim(2) = uimenu(uim(1),'label','Help and Info','callback','freqhlp');
   uim(3) = uimenu(uim(1),'label','Export Data','separator','on',...
                          'callback','freqgui(15)');
   uim(4) = uimenu(uim(1),'label','Print Figure','separator','on',...
                          'callback','printdlg(gcf)');
   uim(5) = uimenu(uim(1),'label','Close and Exit','separator','on',...
                          'callback','close(gcf)');
   uim(6) = uimenu('label','   Axis Limits','separator','on',...
                          'callback','freqgui(17)');
   uim(7) = uimenu('label','   Zoom ON','separator','on',...
                          'callback','freqgui(18)');
   uim(8) = uimenu('label','   Zoom OFF','separator','on',...
                          'callback','freqgui(19)');
   uim(9) = uimenu('label','   RESIZE','callback','freqgui(20)');


%%%%%%%


lft =10;btm=fig_h-5;

   ui(1) = uicontrol('style','frame','pos',[lft-5,btm-100,175,85],...
           'back',grey,'fore',ltgrey);
btm = btm-43;
   ui(2) = uicontrol('style','text','pos',[lft,btm,40,17],...
           'string','Num','horiz','right');
   ui(3) = uicontrol('style','edit','pos',[lft+45,btm,120,18],...
           'back','w','string','1');
btm = btm -22;
   ui(4) = uicontrol('style','text','pos',[lft,btm,40,17],...
           'string','Den','horiz','right');
   ui(5) = uicontrol('style','edit','pos',[lft+45,btm,120,18],...
           'back','w','string','1 1');
btm = btm -30;
   ui(6) = uicontrol('style','push','pos',[lft+20,btm,60,20],...
           'string','Accept','callback','freqplt');
   ui(7) = uicontrol('style','push','pos',[lft+90,btm,60,20],...
           'string','Reset','callback','freqgui(2)');

   ui(26) = uicontrol('style','frame','pos',[lft-5,btm-100+10,175,75],...
           'back',grey,'fore',ltgrey);

btm = btm -40;
   ui(18) = uicontrol('style','text','pos',[lft,btm,80,17],...
           'string','Freq Units','horiz','right');
   ui(10) = uicontrol('style','popup','pos',[lft+85,btm,80,17],...
           'string','Hz|rad/s','back','w','callback','freqgui(3)');

btm = btm -22;
   ui(8) = uicontrol('style','text','pos',[lft,btm,80,17],...
           'string','Freq Range','horiz','right');
   ui(9) = uicontrol('style','edit','pos',[lft+85,btm,80,18],...
           'back','w','string','0,0.5',...
           'callback','freqgui(3)');

btm = btm -22;
   ui(15) = uicontrol('style','text','pos',[lft,btm,80,17],...
           'string','Freq Scale','horiz','right','userdata',0);

   ui(11) = uicontrol('style','popup','pos',[lft+85,btm,80,17],...
           'string','Linear|Log','back','w','callback','freqgui(3)');

   ui(27) = uicontrol('style','frame','pos',[lft-5,btm-100+30+5,175,50],...
           'back',grey,'fore',ltgrey);


btm = btm -38;
   ui(13) = uicontrol('style','check','pos',[lft,btm,165,17],...
           'string','Magnitude in dB','value',0,...
           'callback','freqplt');

btm = btm-22;
   ui(12) = uicontrol('style','check','pos',[lft,btm,165,17],...
           'string','Unwrap Phase','value',0,...
           'callback','freqplt');

btm = btm -30;
   ui(16) = uicontrol('style','check','pos',[lft,btm,95,17],...
           'string','Insta-Plot','fore','c','value',0,'callback','freqgui(5)');

   ui(17) = uicontrol('style','push','pos',[lft+100,btm,65,17],...
           'string','Plot','callback','freqplt');

btm = btm -100;
   ui(25) = uicontrol('style','check','pos',[lft+90,btm,70,20],...
           'string','Cursor','back','y','callback','freqgui(7)');


   pos = get(f,'pos');
   lef = pos(3)-166;
   bot = pos(4)-21;

   ui(21) = uicontrol('style','text','pos',[lef,bot,20,20],...
           'string','X ','horiz','right');
   ui(22) = uicontrol('style','text','pos',[lef+20,bot,60,20],...
           'string','','back','w');
   ui(23) = uicontrol('style','text','pos',[lef+85,bot,20,20],...
           'string','Y ','horiz','right');
   ui(24) = uicontrol('style','text','pos',[lef+105,bot,60,20],...
           'string','','back','w');
   set(ui(21:24),'vis','off');


   ui(19) = axes('units','pixel','pos',[210,250,390,160],'box','on',...
            'xlim',[0,0.5],'color','k',...
            'xcolor','w','ycolor','w','xgrid','on','ygrid','on');
   title('Magnitude Spectrum');
   xlabel('Frequency  [Hz]');
   ui(20) = line('xdata',0,'ydata',0,'color','c');
   set(ui(19),'userdata',ui(20));

   ui(29) = axes('units','pixel','pos',[210,40,390,160],'box','on',...
            'xlim',[0,0.5],'color','k',...
            'xcolor','w','ycolor','w','xgrid','on','ygrid','on');
   title('Phase Spectrum (deg)');
   xlabel('Frequency  [Hz]');
   ui(30) = line('xdata',0,'ydata',0,'color','c');
   set(ui(29),'userdata',ui(30));



   set(f,'userdata',ui);
   freqplt;
   drawnow;
   set(f,'vis','on');

elseif mode == 2, % clear
   set(num_edt,'string','1');
   set(den_edt,'string','1 1');
   freqplt;

elseif mode == 3, % edit fields
   if get(autoplot,'value'),freqplt,end


elseif mode == 5, % autoplot
      if get(autoplot,'value'),
      set(plot_btn,'enable','off');
      freqplt,
   else
      set(plot_btn,'enable','on');
   end
   

elseif mode == 7, % cursor
    if get(curs_btn,'value')
    set(curs_dsp,'vis','on');

    set(f,'WindowButtonMotionFcn','freqgui(8)');

    else

  set(f,'Pointer','arrow','WindowButtonMotionFcn','');


   set(curs_dsp,'vis','off')
   set(curs_btn,'value',0),

    end
elseif mode == 8
       h=gca;
       set(f,'currentaxes',h);
       lim=get(h,'pos');  % get current axis position
       c_p = [lim(1) lim(1)+lim(3) lim(2) lim(2)+lim(4)];
       f_p=get(f,'CurrentPoint'); % get co-ordinates of current point

     if f_p(1) < c_p(1) | f_p(1) > c_p(2) | f_p(2) < c_p(3) | f_p(2) > c_p(4)
        set(f,'Pointer','arrow');

        set(curs_dsp(2),'string','');
        set(curs_dsp(4),'string','');

     else
        curr_pt=get(h,'CurrentPoint');  % get current mouse position
        set(f,'Pointer','crosshair');
        set(curs_dsp(2),'string',num2str(curr_pt(1,1)));
        set(curs_dsp(4),'string',num2str(curr_pt(1,2)));
     end


elseif mode == 15, %Export data
freqplt;
fstr=get(ui(10),'string');fval=get(ui(10),'value');fstr=fstr(fval,:);
mval=get(ui(13),'value');
if mval==0, mstr='Linear';else,mstr='Decibel';end
pval=get(ui(12),'value');
if pval==0, pstr='';else,pstr='Unwrapped';end
fll=get(ui(11),'string');lval=get(ui(11),'value');fll=fll(lval,:);%lin/log

str1='Export transfer function numerator as';
str2='Export transfer function denominator as';
str3=['Export ' fll ' freq in ' fstr ' as'];
if lval==2,str3=['(Use semilogx to plot) ' str3];end
str4=['Export ' mstr ' magnitude as'];
str5=['Export ' pstr ' phase (deg) as'];
exp_str={str1,str2,str3,str4,str5};
exp_var={'ctfnum','ctfden','cfreq','Mctf','Pctf'};

expect=inputdlg(exp_str,'Save Variables as',1,exp_var);
if ~isempty(expect)
hn=get(ui(21),'userdata');%%% FIX FOR UIs 
hd=get(ui(22),'userdata');%%% FIX FOR UIs 
mtf=get(ui(23),'userdata');
ptf=get(ui(24),'userdata');

if ~isempty(expect{1}),assignin('base',expect{1},hn);end
if ~isempty(expect{2}),assignin('base',expect{2},hd);end
if ~isempty(expect{3}),assignin('base',expect{3},mtf(:,1));end
if ~isempty(expect{4}),assignin('base',expect{4},mtf(:,2));end
if ~isempty(expect{5}),assignin('base',expect{5},ptf);end
end


elseif mode == 16, % Axis Control Help
msgbox([' First click in a plot window.                             ',... 
'                  Then, come back and execute'],...
'Axis Control Info','help')


elseif mode == 17, % Axis Control
cur_axs=[];   cur_obj = gco;
   cur_par = get(cur_obj,'parent');
if ~isempty(cur_par)
   if cur_par == f,
      cur_axs = cur_obj;
   else
      cur_axs = cur_par;
   end
end

if isempty(cur_axs),freqgui(16);
else
set(curs_btn,'value',0);
freqgui(7); %%% REPLACE number BY CURSOR MODE 
if any(cur_axs==ui([19,29])),%%%% CHANGE UI numbers
if cur_axs==ui(19),txt='Magnitude';%%%% CHANGE UI numbers
elseif cur_axs==ui(29),txt='Phase';%%%% CHANGE UI numbers
end
axlimdlg(txt,[1 1]); %[1 1]=auto +linear/log,  [1 0]=auto, no linear/log
else
freqgui(16);
end
end


elseif mode == 18, % zoom on
set(curs_btn,'value',0);
freqgui(7);                    % Use CURSOR MODE 
if get(ui(15),'userdata')==0  %%Change ui(15) to whatever
zoom on
set(ui(15),'userdata',1);
else
msgbox('Zoom is already ON.','Zoom Info','help')
end

elseif mode == 19, % zoom off
set(curs_btn,'value',0);
freqgui(7);                    % Use CURSOR MODE 
if get(ui(15),'userdata')==1

%%%%%Insert old zoom off here
set(f,'currentaxes',ui(19));
zoom out,zoom reset,zoom off
set(f,'currentaxes',ui(29));
zoom out,zoom reset,zoom off
set(ui(15),'userdata',0);
else
msgbox('Zoom is already OFF.','Zoom Info','help')
end



elseif mode == 20, % mouse functionality
cur_axs = [];   cur_obj = gco;
   cur_par = get(cur_obj,'parent');

if ~isempty(cur_par)
   if cur_par == f,
      cur_axs = cur_obj;
   else
      cur_axs = cur_par;
   end
end


if isempty(cur_axs),freqgui(16);return,end

   if any(cur_axs == [sig_axs,opr_axs]),
      axes(cur_axs);
         axs_pos = get(cur_axs,'pos');
         if axs_pos(4) == 160, % need to maximize
               set([sig_axs,opr_axs],'vis','off');
            set(cur_axs,'pos',[210,40,390,375],'vis','on');

         else % minimization
              set(sig_axs,'pos',[210,250,390,160],'vis','on');
              set(opr_axs,'pos',[210,40,390,160],'vis','on');
         end
   else
         freqgui(16);
   end


end