rcv_tbl.m

baseband解调

function [y1,y2,y3,y4] =  rcv_tbl(points,precision);

fid = fopen('final2_sim.tbl','r');
if fid==-1
    disp('Input file final2_sim.tbl not found');
    return;
end;

br =0;
ts = 1600;

while br == 0
    d = fscanf(fid,'%s',1);
    dd = size(d);
    if (d(dd(2)) == '>')
        br = 1;
    end;
end;

fscanf(fid,'%s',15); % skip until next line
fscanf(fid,'%s',16);
curline=fscanf(fid,'%s',1);
fscanf(fid,'%s',15);
%x = 0;
%read = 0;
%while (read == 0)
 % curline = fscanf(fid,'%s',1); % get timestamp
  %fscanf(fid,'%s',7); % skip until read is next
  %read = fscanf(fid,'%d',1); % get status of read signal
  %fscanf(fid,'%s',6); % skip until next line
  %end;


k = 1;
%strip off anything past decimal point of timestamp
while (curline(k) ~= '.')
  curval(k) = curline(k);
  k = k+1;
end;

curnum = str2num(curval);

next = 1;
while (next < points+1)
  k = 1;
  curline = fscanf(fid,'%s',1); %get timestamp
  fscanf(fid,'%s',7); %skip until next real
  a1 = fscanf(fid,'%x',1); % grab real in HEX format
  a2 = fscanf(fid,'%x',1);% grab imag in HEX format
  b1 = fscanf(fid,'%x',1); % grab real in HEX format
  b2 = fscanf(fid,'%x',1);%grab imag in HEX format
  c1 = fscanf(fid,'%x',1); % grab real in HEX format
  c2 = fscanf(fid,'%x',1);% grab imag in HEX format
  d1 = fscanf(fid,'%x',1); % grab real in HEX format
  d2 = fscanf(fid,'%x',1);%grab imag in HEX format
  %fscanf(fid,'%s',9); % to end of line

  while (curline(k) ~= '.')
    time(k) = curline(k);
    k = k+1; 
  end;

  timenum = str2num(time);

  if (timenum == (curnum + ts))
    curnum = timenum;
      if (a1 > (2^(precision-1)-1))
	    a1 = a1 - 2^precision;
      end;
      if (a2 > (2^(precision-1)-1))
	    a2 = a2 - 2^precision;
      end;
     if (b1 > (2^(precision-1)-1))
	    b1 = b1 - 2^precision;
      end;
      if (b2 > (2^(precision-1)-1))
	    b2 = b2 - 2^precision;
      end; 
      if (c1 > (2^(precision-1)-1))
	    c1 = c1 - 2^precision;
      end;
      if (c2 > (2^(precision-1)-1))
	    c2 = c2 - 2^precision;
      end;
     if (d1 > (2^(precision-1)-1))
	    d1 = d1 - 2^precision;
      end;
      if (d2 > (2^(precision-1)-1))
	    d2 = d2 - 2^precision;
      end; 
    y1(next) = a1+i*a2;
    y2(next) = b1+i*b2;
    y3(next) = c1+i*c2;
    y4(next) = d1+i*d2;
    next = next + 1;
  end;

end;  

%y = y(2:points+1);
%z = z(2:points+1);
fclose(fid);