my_read_tbl_5m_psk.asv

pie edcode编码 程序设计

clear all;
clc
%以只读方式打开.tbl文件
fid = fopen('pie_code.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<160)  %j<=5000??

     data0 = fscanf(fid,'%s',1);
    size_data0=size(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);
        data_in = fscanf(fid,'%s',1);
        data2 = fscanf(fid,'%s',1);
        data_out = fscanf(fid,'%s',1);
%         cos_out2 = fscanf(fid,'%s',1);
%         pre_fir_q = fscanf(fid,'%s',1);
%         pre_fir_i = fscanf(fid,'%s',1);
%         out_Q = fscanf(fid,'%s',1);
%         out_I = fscanf(fid,'%s',1);        
%         pre_mult_Q = fscanf(fid,'%s',1);
%         pre_mult_I = fscanf(fid,'%s',1);
%         after_mult_Q3 = fscanf(fid,'%s',1);
%         after_mult_Q4 = fscanf(fid,'%s',1);
%         pre_compare_out = fscanf(fid,'%s',1);
%         wave_Q = fscanf(fid,'%s',1);
%         wave_I = fscanf(fid,'%s',1);

        if ( (data1=='0') & (data2=='=') )
            j=j+1;
%            temp=str2num(data3);
             temp=hex2dec(data_in);
%             if temp>=2^15                               %bu ma!
%                 temp=temp-2^16;
%             end
             Y(j,:)=temp;
            
            %            temp=str2num(data3);
             temp1=hex2dec(data_out);
             if temp1>=2^16                               %bu ma!
                 temp1=temp1-2^17;
             end
             Y1(j,:)=temp1;
%             
%             %            temp=str2num(data3);
%             temp2=hex2dec(cos_out2);
%             if temp2>=2^13                               %bu ma!
%                 temp2=temp2-2^14;
%             end
%             Y2(j,:)=temp2;
%             
%             %            temp=str2num(data3);
%             temp3=hex2dec(pre_fir_q);
%             if temp3>=2^13                               %bu ma!
%                 temp3=temp3-2^14;
%             end
%             Y3(j,:)=temp3;
%             
%             %            temp=str2num(data3);
%             temp4=hex2dec(pre_fir_i);
%             if temp4>=2^13                               %bu ma!
%                 temp4=temp4-2^14;
%             end
%             Y4(j,:)=temp4;
%             
%             %            temp=str2num(data3);
%             temp5=hex2dec(out_Q);
%             if temp5>=2^25                               %bu ma!
%                 temp5=temp5-2^26;
%             end
%             Y5(j,:)=temp5;
%             
%             %            temp=str2num(data3);
%             temp6=hex2dec(out_I);
%             if temp6>=2^25                               %bu ma!
%                 temp6=temp6-2^26;
%             end
%             Y6(j,:)=temp6;
%             
%             %            temp=str2num(data3);
%             temp7=hex2dec(pre_mult_Q);
%             if temp7>=2^15                               %bu ma!
%                 temp7=temp7-2^16;
%             end
%             Y7(j,:)=temp7;
%             
%             %            temp=str2num(data3);
%             temp8=hex2dec(pre_mult_I);
%             if temp8>=2^15                               %bu ma!
%                 temp8=temp8-2^16;
%             end
%             Y8(j,:)=temp8;
            
            %            temp=str2num(data3);
%             temp9=hex2dec(wave_Q);
% %             if temp>=2^15                               %bu ma!
% %                 temp=temp-2^16;
% %             end
%             Y9(j,:)=temp9;
%             
%             %            temp=str2num(data3);
%             temp10=hex2dec(wave_I);
% %             if temp>=2^15                               %bu ma!
% %                 temp=temp-2^16;
% %             end
%             Y10(j,:)=temp10;
%             
%             Y11(j,:)=temp9+temp10;
%             
%             %            temp=str2num(data3);
%             temp11=hex2dec(pre_compare_out);
% %             if temp>=2^15                               %bu ma!
% %                 temp=temp-2^16;
% %             end
%             Y11(j,:)=temp11;
        end
        end
        end
 end
fclose(fid);


j
%X=0: 20: 20000;

grid on;
% figure(1);
% subplot(221)
% plot(Y1);
% % figure(2);
% subplot(222)
% plot(abs(fft(Y1)));
% 
% % figure(3);
% subplot(223)
% plot(Y2);
% % figure(4);
% subplot(224)
% plot(abs(fft(Y2)));
% 
% 
% figure(5);
% subplot(221)
% plot(Y3);
% % figure(6);
% subplot(222)
% plot(abs(fft(Y3)));
% 
% % figure(7);
% subplot(223)
% plot(Y4);
% % figure(8);
% subplot(224)
% plot(abs(fft(Y4)));
% 
% figure(9);
% subplot(221)
% plot(Y5);
% % figure(10);
% subplot(222)
% plot(abs(fft(Y5)));
% 
% % figure(11);
% subplot(223)
% plot(Y6);
% % figure(12);
% subplot(224)
% plot(abs(fft(Y6)));
% 
% figure(13);
% subplot(221)
% plot(Y7);
% % figure(14);
% subplot(222)
% plot(abs(fft(Y7)));
% 
% % figure(15);
% subplot(223)
% plot(Y8);
% % figure(16);
% subplot(224)
% plot(abs(fft(Y8)));
% 

% % figure(18);
% subplot(222)
% plot(abs(fft(Y9)));
% 
%  figure(19);
 figure(1);
subplot(211)
 plot(Y);
 subplot(212)
plot(Y1);
figure(2);
subplot(211) 
plot(abs(fft(Y)));
subplot(212)
 plot(abs(fft(Y1)));
figure(3);
subplot(211) 
plot(10*log10(abs(fft(Y))));
subplot(212)
plot(10*log10(abs(fft(Y1))));
 
% % figure(20);
%  subplot(122)
% plot(abs(fft(Y10)));

%  figure(21);
%  plot(Y11);
% figure(22);
% plot(abs(fft(Y11)));