s_extract.m

spread spectrum communication will be helpful to you!

function [sout,fout] = s_extract(file_name,fmin,fmax,df)

       d=load(file_name,'-mat');
       f=d.data.s(1,1).freq;
       
       
       n = sum(fmin>f);       % low freq index
       m = sum(f<fmax);       % high freq index
       inc = round(df/(f(2)-f(1)));
       if inc >0
          fout = f(n:inc:m);
          v = n:inc:m;
      else
          fout =f(n);
          v = n;
      end;
       for k =1:2
           for l=1:2
                [x,y]      = dbphase2xy(d.data.s(k,l).mag(v),d.data.s(k,l).ph(v));
                sout(k,l).s= x+j*y;
           end;
       end;