Subversion Repositories wdsp

function plotdata(infile,outfile,ibin,fbin,add_g,mult_g,N)

%  plotdata(infile,outfile,ibin,fbin,add_g,mult_g,N)

%   Reads in from infile (test vectors sent into fft_f         iled_tb

%   and outfile (output of fft_filed_tb, usually data.out)

%

%   ibin is integer b its in the input.

%   fbin is fractio na l bits in the input.

%   add_g is the add er growth.

%    mult_g is the mu ltiplier growth.

%    N  is the number of points in the FFT.

%

%   This prog  ram us es:

%       str2fr ac.m

%

clf

fout=fopen(outfile);

ibout=ibin+add_g*log2(  N);

if mult_g >0

    ibout=ibout+floor((  log2(N)-1)/2);

end;

fbout =fbi n+(floor((log2(N)-1)/2)*(mult_g-1));

dp =ibo ut+fbout;

dp2=2*dp;

A=(fscanf(fout,'%  s',[dp2,inf]))';

[nr ,nc]=size(A);

hold on;

offset=N+2  ;

%getting fi rst outputs

if (N+offset)<=nr

    for k=offset+1:offset+N

        out1r(k-offset)= str2frac  (A(k,1:dp),ibout,fbout);

        out1i(k-offset)=str2frac(A(k,     dp+1:dp2),ibout,fbout);

    end

 out1r=out1r';

 out1i=out1i';

  %Bit reversing and r        ecombining outputs

 for k=0:(N-1)

        bitorder=dec2bin(k,log2(N));

        for m=1:(size(bi torder,2)/2)

             temp=bitorder (m);

             bitorder(m)=bitorder(size(bitorder,2)+1-m);

            bitorder(size (bitorder,2)+1-m)=temp;

        end

        bitrevpos=bin2dec(bitorder);

        adata(bitrevpos+1)=out1r(k+1)+i*out1i(k       +1);

    end

    subplot(2,2,1),plot(fftshift(abs(adata)),'k'),axis tight,

    title( 'Magnitude of FFT of Sinc Wave Input'),hold on;

    subplot (2,2,2),plot(fftshift(angle(adata)),'k'),axis tight,

    title('Angle of FFT of  Sinc Wave Input'),hold on;

end

%getting second outputs

offset=offset+N;

if (N +offset)<=nr

     for k=offset+1:offset+N

        out1r(k- offset)=str2frac(A(k,1:dp),ibout,fbout);

         out1i(k-offset)=str2frac(A(k,dp+1:dp2),ibout,fbout);

    end

    out1r=out1r';

 out1i=out1i';

  %Bit reversing and recombining outputs

 for k=0:(N-1)

        bitorder=dec2bin(k,log2(N));

        for m=1:(size(bitorder,2)/2)

            temp=bitorder(m);

            bitorder(m)=bitorder(size(bitorder,2)+1-m);

            bitorder(size(bitorder,2)+1-m)=temp;

        end

          bitrevpos=bin2dec(bitorder);

        bdata(bitrevpos+1)=out1r(k+1)+i*out1i(k+1);

 end

 subplot( 2,2,3),plot(fftshift(real(bdata)),'k'),axis tight,

    title('Real p art of FFT of Square wave input'),hold on;

 subplot(2 ,2 ,4),plot(fftshift(imag(bdata)),'k'),axis tight,

    title('Imaginary part of FF    T of Square wave input'),hold on;

end

%getting input waves from testvec and outputt       ing them on same plots

fin=fopen(infile);

dp=ibin+fbin;

dp2=2*dp+2;

dp=dp+2;

A=(fscanf(fin,'%s',[dp2  ,inf]))';

[nr,nc]=size (A);

offset=2;

%getting first inputs

if (N+offset)<=nr

    for k=offset+1:offse  t+N

        outd ata(k-offset )=str2frac(A(k,3:dp),ibin,fbin) +

i*str2frac(A(k,dp+1:dp 2),ibin,fbin);

    end

    afftdata=fft(outdata );

    subplot(2,2,1),plot(fftshift(abs(afftdata)));

    subpl ot(2,2,2),plot(fftshift(angle(afftdata)));

end

%getting second inputs

offset=offset+N;

if (N+offset)<=nr

    for k=offset+1:offset+ N

         outdata(k-offset)=str2frac(A(k,3:dp),ibin,fbin) +

i*str2frac(A(k,dp+1:dp 2),ibin,fbin);

     end

    bfftda  ta=fft(outdata);

    subplo t(2,2,3),plot(fftshift(real(bfftdata)));

    subplot(2,2,4),plot(fftshift(imag(bfftdata)));

end

fclose('all');