clear all;
clc
%以只讀方式打開.tbl文件
fid = fopen('dsb_ask_mod.sim.tbl','r');


%data1 = fscanf(fid,'%s');
%fclose(fid);
%b = find(data == '=');      %找出“＝”的下標
%number = length(b);         %取得文件長度
%j = 0;
%for i = 1:number
 %   if   (data(b(i) - 1) == '0')    %&   (data(b(i) - 2) == '>')  &   (data(b(i) - 3) == '0')   &   (data(b(i) - 4) == '.')  &   (data(b(i) - 5) == '0')  &   ( rem( data(b(i) - 6) ,2 )== '0')  ) %濾除多余的“＝”
 %       j = j+1;
%        %取出幅度值，數據為十六進制數
%        c_s(j,1) = data(b(i) + 1);
%        c_s(j,2) = data(b(i) + 2);
%        c_s(j,3) = data(b(i) + 3);
%        c_s(j,4) = data(b(i) + 4);
%    end
%end
%d_s = hex2dec(c_s);         %將十六進制數轉化為十進制數
%plot(d_s);
%grid on;
j=0;
while(j<1000)  %j<=5000??
    
     data0 = fscanf(fid,'%s',1);
     
    size_data0=size(data0);
%     data0
 if (  data0(size_data0(2))=='>' & size_data0(2)>3 )
        if (data0(size_data0(2)-1)=='0' & data0(size_data0(2)-2)=='.' & data0(size_data0(2)-3)=='0'  )        
        data1 = fscanf(fid,'%s',1);
        data2 = fscanf(fid,'%s',1);
        data2 = fscanf(fid,'%s',1);
        data2 = fscanf(fid,'%s',1);
%         data3 = fscanf(fid,'%s',1);
%         data3 = fscanf(fid,'%s',1);
        data3 = fscanf(fid,'%s',1);
   
        if ( (data1=='1') & (data2=='=') )
            j=j+1;
             data0
             data3
%            temp=str2num(data3);
            temp=hex2dec(data3);
            if temp>=2^13
                temp=temp-2^14;
            end
            Y(j,:)=temp;
        end
    end
 end
end
fclose(fid);


j
%X=0: 20: 20000;

grid on;
% figure(1);
subplot(121)
plot(Y);
% figure(2);
subplot(122)
plot(abs(fft(Y)));