%
%  This file processes the data from the first trial
%
function data = process_data_1dof;
%
%  enter the data as 
%       frequency (in Hz) input amplitude (in cm/degrees/radians)  output amplitude (in cm/degrees/radians); 
%
%  such as 
%
%  measured = [1.00 1.50 1.740;
%              2.00 1.50 2300;
%              2.50 0.50 1100];
%

  measured = [0.5 5.0 16.5403;
              1.0 3.0 11.62;
              1.5 3.0 18.048;
              2.0 1.0 23.4716;
              2.5 2.0 17.0916;
              3.0 3.0 10.8129;
              3.5 2.0 3.9041;
              4.0 3.0 3.9379;
              4.5 3.0 2.9927;
              5.0 3.0 2.2613;
              5.5 3.0 1.4849;
              6.0 3.0 1.4849;
              6.5 3.0 1.2711;
              7.0 3.0 1.0236;
              7.5 3.0 0.8886;
              2.1 1.0 38.4931;
              2.2 0.5 17.7667;
              2.3 0.5 10.0928;
              1.9 1.0 11.4543;
              2.05 0.5 6.6834];
%
%  Be sure the rows are sorted by frequency
%
   measured = sortrows(measured,1);
%
%  get the magnitude of the transfer function
%
  H = measured(:,3)./measured(:,2);
  f = measured(:,1); 
  data = [f H];
%
  return;