f_start= str2num(get(handles.fstart_edit,'string'));  % start freq in MHz must be < fb !!!        f_stop = str2num(get(handles.fstop_edit,'string'));   % stop freq        range= 1+ (f_stop>13) + (f_stop>130);  % clever those hp dudes         switch range           case 1              % < 13 MHz              fprintf(RF,['R1W1FA',strrep(sprintf('%4.0f',f_start*100),' ','0'),'FB',strrep(sprintf('%4.0f',f_stop*100),' ','0'),'E']);             case 2              fprintf(RF,['R2W1FA',strrep(sprintf('%4.0f',f_start*10),' ','0'),'FB',strrep(sprintf('%4.0f',f_stop*10),' ','0'),'E']);             case 3              fprintf(RF,['R3W1FA',strrep(sprintf('%4.0f',f_start),' ','0'),'FB',strrep(sprintf('%4.0f',f_stop),' ','0'),'E']);          end;    % note the "E" = Execute     end;      fprintf(DISPS,'C3I1M1E');              %  turn electrical line compensation off     fprintf(DISPS,'C3I2M1E');       fprintf(DISPS,'L');                    % learn DISPlay module setup   disp_str = fscanf(DISPS);   ref_c1=str2num(disp_str(15:20))/100;   % channel 1 reference level  (div by 100 for dB)   ref_c2=str2num(disp_str(31:36))/10;    % channel 2 reference level  (div by 10 for phase)   per_div.ph    =[180 90 45 20 10  5  2  1];   per_div.mag   = [20 10  5  2  1 .5 .2 .1];   mscale=2;         % start with 10 dB per div   pscale=3;         % and 45 deg/div for phase    autorange=get(handles.autorange,'value');   d=0.5;   % need to pause or ref_c1 setting gets hosed in autoscale mode   for k=1:1+autorange       fprintf(SN,'AB');   % stop current measurement in storage normalizer        switch param          case 's11'                fprintf(SPAR,'1');           % 1 is forward, 2 is reverse                pause(d)                fprintf(DISPS,['C1I4M2S',int2str(mscale),'E']);  % A/R, Mag                fprintf(DISPS,['C2I4M3S',int2str(pscale),'E']);  % A/R, Phase          case 's12'                fprintf(SPAR,'2');           % 1 is forward, 2 is reverse                pause(d)                fprintf(DISPS,['C1I4M2S',int2str(mscale),'E']);  % A/R, Mag                fprintf(DISPS,['C2I4M3S',int2str(pscale),'E']);  % A/R, Phase                       case 's21'                fprintf(SPAR,'1');           % 1 is forward, 2 is reverse                pause(d)                fprintf(DISPS,['C1I5M2S',int2str(mscale),'E']);  % B/R, Mag                fprintf(DISPS,['C2I5M3S',int2str(pscale),'E']);  % B/R, Phase          case 's22'                fprintf(SPAR,'2');           % 1 is forward, 2 is reverse                pause(d)                fprintf(DISPS,['C1I5M2S',int2str(mscale),'E']);  % B/R, Mag                fprintf(DISPS,['C2I5M3S',int2str(pscale),'E']);  % B/R, Phase       end;       fprintf(SN,'S2C1I2C2I2TS2'); %turn c1 & c2 on, 2 sweeps (ts2) to get mag & phase       fprintf(SN,'EX');       fprintf(SN,'P3');    % select Page 3 , implicit AB        ch1str=fscanf(SN);   % fetch mag in dB result       fprintf(SN,'P4');    % select Page 4 , implicit AB        ch2str=fscanf(SN);   % fetch phase result            % results should be the same size, but sometimes thay are not       %     lch1=length(ch1str)      %     lch2=length(ch2str)      %     [ch1str(1:20);ch2str(1:20)]      %     ch1str(2960:end)      %     ch2str(2960:end)      % end;              x_end=2980;                                         % empirically determined ending point        x_beg=5;                                            % ditto        dB_scale =per_div.mag(mscale)*10/500;               % if mag in dB and 10 db per div 100/500=0.2          x1=dB_scale*str2num(ch1str(x_beg:x_end))+ref_c1;    % trim bogus data, add ref offset       deg_scale=per_div.ph(pscale)*8/400; %          x2=deg_scale*str2num(ch2str(x_beg:x_end))+ref_c2;   % trim bogus data, add ref offset             N= min(length(x2),length(x1));       result.mag= x1(1:N);       result.ph = x2(1:N);       result.f  = f_start+(0:(N-1))'*(f_stop-f_start)/N;          % need this if autoranging       % mag and phase spread        if autorange & (k==1)          min_mag = min(result.mag);             max_mag= max(result.mag);          ref_c1 = round((min_mag+max_mag)/2);        % new reference level          if ref_c1>=0 s='+'; else  s='-'; end;          sm=['C1M2R',s,sprintf('%04i',abs(round(ref_c1*100))),'E'];          fprintf(DISPS,sm);          usedmag_per_div =  (max_mag-min_mag)/10;   % 10 dB per div initial measurement          mscale=sum(usedmag_per_div<=per_div.mag);  % new mscale                  % do same for phase          [min_ph,i1] = min(result.ph);          [max_ph,i2] = max(result.ph);                   % ph_min_max=[min_ph,max_ph]          ref_c2 = round((min_ph+max_ph)/2);                  if ref_c2>=0 s='+'; else  s='-'; end;          sp=['C2M3R',s,sprintf('%04i',abs(round(ref_c2*10))),'E'];          fprintf(DISPS,sp);          usedph_per_div=(max_ph-min_ph)/8;         % 360 phase should map to 45 deg / div          pscale=sum(usedph_per_div <= per_div.ph); % new pscale      end;   end;      % restore condictions that were read back from panel setup   fprintf(DISPS,[disp_str,'E']);   fprintf(RF,[rf_str,'E']);   fprintf(SN,'EX');    % resume normal measurement process   close_network_analyzer(SN,DISPS,RF,SPAR);   set(handles.status,'string',[param,' complete']);% --------------------------------------------------------------------% --------------------------------------------------------------------function varargout = open_Callback(h, eventdata, handles, varargin)      if isempty(eventdata)          [file_n,path_n]=uigetfile([pwd,'\*.rfa'],'Open File',0.5,0.5);           path_file=[path_n,file_n];          if isempty(file_n)             return          end;       else           path_file=eventdata;       end;      s=load(path_file,'-mat');      data=s.data;      for k=1:2          for p=1:2              set(handles.s(k,p).mag,'xdata',data.s(k,p).freq, 'ydata',data.s(k,p).mag);              set(get(handles.s(k,p).mag,'parent'),'xlim',data.s(k,p).freq([1,end]));              set(handles.s(k,p).ph,'xdata',data.s(k,p).freq, 'ydata',data.s(k,p).ph);              set(get(handles.s(k,p).ph,'parent'),'xlim',data.s(k,p).freq([1,end]));          end;      end% --------------------------------------------------------------------function varargout = save_Callback(h, eventdata, handles, varargin)% --------------------------------------------------------------------function varargout = save_as_Callback(h, eventdata, handles, varargin)                  for k=1:2                      for p=1:2                          data.s(k,p).freq = get(handles.s(k,p).mag,'xdata');  % this has cal factor                          data.s(k,p).mag  = get(handles.s(k,p).mag,'ydata');  % if cal factor was applied                          data.s(k,p).ph   = get(handles.s(k,p).ph,'ydata');   % which may be the way to go                      end;                  end                  [file_n,path_n]=uiputfile([pwd,'\*.rfa'],'Save File',0.5,0.5);                   if isempty(file_n)                      return                  end;                                    % check for rfa extension                  p=findstr(file_n,'.');                  if ~isempty(p)                     % force it                     temp=file_n(1:(p-1));                     file_n=temp;                  end;                  file_n=[file_n,'.rfa'];                  file_path=[path_n,file_n];                  save(file_path,'data','-mat');                  % --------------------------------------------------------------------function varargout = ps_Callback(h, eventdata, handles, varargin)                     name=eventdata;                     tag=[eventdata,'_plot'];                     h_fig=findobj('tag',tag);                     if isempty(h_fig)                        h_fig=rfna_plot;                        set(h_fig,'tag',tag,'name',name);                     else                        figure(h_fig);                      end;                     k = str2num(name(2)); p = str2num(name(3));                     data.freq = get(handles.s(k,p).mag,'xdata');                     data.mag  = get(handles.s(k,p).mag,'ydata');                     data.ph   = get(handles.s(k,p).ph,'ydata');                     data.name = name;                     rfna_plot('plot_data',h_fig,data,guidata(h_fig));% --------------------------------------------------------------------function varargout = fstart_edit_Callback(h, eventdata, handles, varargin)% --------------------------------------------------------------------function varargout = fstop_edit_Callback(h, eventdata, handles, varargin)% --------------------------------------------------------------------function varargout = inst_pannel_Callback(h, eventdata, handles, varargin)% --------------------------------------------------------------------function varargout = autorange_Callback(h, eventdata, handles, varargin)% --------------------------------------------------------------------function varargout = correct_Callback(h, eventdata, handles, varargin)