qutout.m

附录 光盘说明 本书附赠的光盘包括各章节实例的设计工程与源码

%-------------------------------------------------------------------------------%
%---         程序名称：Rader_Hilber算法的 Quartus 仿真结果     ---%
%---         带奇偶抽样。                                                             ---%
%---             作者：王旭东                                                         ---%
%---             时间：2003年6月25日                                           ---%
%---             版本：v0000.00.02                                                ---%
%---         版权所有违者必究！                                                   ---%
%-------------------------------------------------------------------------------%

%-------------------------------------------------------------------------------%
%--------------- 文件含有三列：ad_out8, In, Qn -----------------------%
%-------------------------------------------------------------------------------%
clear;
clc;

disp('%-----  Radar_Hilbert正交变换算法的 Quartus 仿真结果 -----%');
fid=fopen('..\matlab\rader_hilbert.tbl','r');

%---------------------- 定义常数及变量初始化 ----------------------------%
    OutWidth=16;
    k=0;
    A1='1';
    A2='0';
    A1_next='0';
%-------------------------------------------------------------------------------%
    
while A1_next~='X'
    A1_next=fscanf(fid,'%s',1);
end; 
    A1_next=fscanf(fid,'%s',1);
    A1_next=fscanf(fid,'%s',1);
    A1_next=fscanf(fid,'%s',1);
    A1_next=fscanf(fid,'%s',1);
%    A1_next=fscanf(fid,'%s',1); %Add by me 6.19

while A1_next~='X' %-- 在此处设中断观察A1_next是否取到ad_end的第一个零。
    if (A1_next=='0')
        A1='0';
    else
        A1='1';
        A1_next=fscanf(fid,'%s',1);
        A1_next=fscanf(fid,'%s',1);
        A1_next=fscanf(fid,'%s',1);
        A1_next=fscanf(fid,'%s',1);
        A1_next=fscanf(fid,'%s',1);
    end;
    if(A1=='0')
        A1_next=fscanf(fid,'%s',1);
        A1_next=fscanf(fid,'%s',1);
        A1_next=fscanf(fid,'%s',1);
        A1_next=fscanf(fid,'%s',1);
        A1_next=fscanf(fid,'%s',1);
      if(A1_next=='1') 
        A1='1';
        A1_next1=fscanf(fid,'%s',1);
        A1_next1=fscanf(fid,'%s',1);        
        data_out_r=A1_next1;
        A1_next2=fscanf(fid,'%s',1);
        data_out_i=A1_next2;
        k=k+1;
        A1_next=fscanf(fid,'%s',1); %取数指针归位。
        A1_next=fscanf(fid,'%s',1);
%        A1_next=fscanf(fid,'%s',1);
        
%----------------------  取出实部   -------------------------------------------------%        
        wauk_bin_r=dec2bin(hex2dec(num2str(data_out_r)),OutWidth);
        %wauk_hao(k)=0;
        if (wauk_bin_r(1)==dec2bin(1))
          for j=1:OutWidth
            if wauk_bin_r(j)==dec2bin(1)
               wauk_bin_r(j)=dec2bin(0);
            else
                wauk_bin_r(j)=dec2bin(1);
            end;
          end;
          wauk_hao_r(k)=-(bin2dec(wauk_bin_r)+1);
        else
          wauk_hao_r(k)=hex2dec(num2str(data_out_r));
        end;
%-----------------------  取出虚部  --------------------------------------------------%
        wauk_bin_i=dec2bin(hex2dec(num2str(data_out_i)),OutWidth);
        %wauk_hao(k)=0;
        if (wauk_bin_i(1)==dec2bin(1))
          for j=1:OutWidth
            if wauk_bin_i(j)==dec2bin(1)
               wauk_bin_i(j)=dec2bin(0);
            else
                wauk_bin_i(j)=dec2bin(1);
            end;
          end;
          wauk_hao_i(k)=-(bin2dec(wauk_bin_i)+1);
        else
          wauk_hao_i(k)=hex2dec(num2str(data_out_i));
        end;
        
      end;
    end;
 end;
 fclose(fid);
 
for i=1:k
 outdata(i)=sqrt(wauk_hao_r(i)^2+wauk_hao_i(i)^2);
end;

%-----------------------  Quartus输出数据的图形  ---------------------------------------%
 %figure(3);
 %plot(wauk_hao_r);
 %title 'data-out-r 就是In';
 
 %figure(4);
 %plot(wauk_hao_i,'-r');
 %title 'data-out-i 就是Qn';
 
 i=1:length(wauk_hao_r);
 figure(5);
 plot(i,wauk_hao_r,'-b',i,wauk_hao_i,'-r');
 title 'Radar_Hilbert 的 QuartusII 仿真结果';
%-----------------------------------------------------------------------------------------------%