Subversion Repositories pipelined_fft_128

#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

#

# Copyright(c) 2009 by Unicore Systems. All rights reserved

#

# DESIGN           :    FFT128_core

# FILENAME              :        sinerom64_gen.pl

# CREATED               :        10.1:2008

# MODIFIED              :       

# VERSION                       :       1.0

#

# AUTHORS:      Anatolij Sergiyenko.

# HISTORY       :       

#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~            

# DESCRIPTION   :       Test generating PERL file1

# FUNCTION:             Generating sine wave ROM 

#                                        with $n samples  which consists of 4 sine wave sum

#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~            

$aw=7;   ##   - address bit width

$n = 128; ### - period length

$n0 = 0; ### - init address  

$f1 = 1;  ## - first frequency

$f2 = 3;  ## -second frequency

$f3 = 5;  ## - 3-d frequency

$f4 = 7;  ## -4-th frequency

$Pi = 3.14159265358     ;

$ph =$Pi/4; ## initial phase    => 45 degree

    open OutFile, ">Wave_ROM128.v" || die "Cannot open file .V";    # file is open to write

        print OutFile   "//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^\n";

    print OutFile       "//   ROM with ".$n." samples of the sine waves at the frequencies = ".$f1.", ".  $f2.",". $f3." and ".$f4 ."\n";

    print OutFile       "//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^\n";

        print OutFile   "   `timescale 1 ns / 1 ps  \n";

        print OutFile   "module Wave_ROM128 ( ADDR ,DATA_RE,DATA_IM,DATA_REF ); \n";

    print OutFile       "       output [15:0] DATA_RE,DATA_IM,DATA_REF ;     \n";

    print OutFile       "       input [".($aw-1).":0]    ADDR ;     \n";

# # Cosine table generation

    print OutFile       "       reg [15:0] cosi[0:".($n-1)."];    \n";

    print OutFile       "       initial   begin    \n";

        for $j(0..$n/16-1){

          $cos_row=0;

         for $i(0..15){

$CosArr=(32767*cos(2*$Pi*$f1*(($n0+$i+16*$j)/$n)));     #23170

#       $CosArr=(16384*cos(2*$Pi*$f1*(($n0+$i+$j*16)/$n)));

#                        +16383*cos(2*$Pi*$f2*(($n0+$i+$j*16)/$n))); 

#       $CosArr=(8191*cos($ph+2*$Pi*$f1*(($n0+$i+$j*16)/$n))

   #              +8191*cos($ph+2*$Pi*$f2*(($n0+$i+$j*16)/$n))

        #                        +8191*cos($ph+2*$Pi*$f3*(($n0+$i+$j*16)/$n))

           #              +8191*cos($ph+2*$Pi*$f4*(($n0+$i+$j*16)/$n))); 

            $sinv= To_Hex4(int ($CosArr));

        print OutFile   " cosi[".($n0+$i+$j*16)."]=16'h".$sinv.";";

}

print OutFile "\n";

}

print OutFile "     end \n\n";

         print "\n";

# # Sine table generation

    print OutFile       "       reg [15:0] sine[0:".($n-1)."];    \n";

    print OutFile       "       initial   begin    \n";

                for $j(0..$n/16-1){

          $sin_row=0;

         for $i(0..15){

        $SinArr=(32767*sin(2*$Pi*$f1*(($n0+$i+16*$j)/$n)));        #16384          23170

#       $SinArr=(16384*sin(2*$Pi*$f1*(($n0+$i+$j*16)/$n)));        #

#                        +16383*sin(2*$Pi*$f2*(($n0+$i+$j*16)/$n)));    

#       $SinArr=(8191*sin($ph+2*$Pi*$f1*(($n0+$i+$j*16)/$n))

        #                 +8191*sin($ph+2*$Pi*$f2*(($n0+$i+$j*16)/$n))

                #                        +8191*sin($ph+2*$Pi*$f3*(($n0+$i+$j*16)/$n))

               #          +8191*sin($ph+2*$Pi*$f4*(($n0+$i+$j*16)/$n))); 

            $sinv= To_Hex4(int ($SinArr));

        print OutFile   " sine[".($n0+$i+$j*16)."]=16'h".$sinv.";";

}

print OutFile "\n";

}

         print "\n";

print OutFile "      end \n\n";

##Reference table generation

    print OutFile       "       reg [15:0] deltas[0:".($n-1)."];    \n";

    print OutFile       "       initial   begin    \n";

        for $j(0..$n/4-1){

         for $i(0..3){

        print OutFile   " deltas[".($n0+$i+$j*4)."]=16'h0000;";

}

print OutFile "\n";

}

##print OutFile " deltas[".$f1."]=16'h3fff;  deltas[".$f2."]=16'h3fff;\n";

$SinMag=To_Hex4(int (8191*4*sin($ph)));

print OutFile   " deltas[".$f1."]=16'h".$SinMag."; deltas[".$f2."]=16'h".$SinMag."; deltas[".$f3."]=16'h".$SinMag."; deltas[".$f4."]=16'h".$SinMag.";\n";

print OutFile "     end \n\n";

print OutFile " assign DATA_RE=cosi[ADDR];      \n";

print OutFile " assign DATA_IM=sine[ADDR];      \n";

print OutFile " assign DATA_REF=deltas[ADDR];   \n";

print OutFile "endmodule   \n";

   close(OutFile);

   sub To_Hex2{

        my(%Hexnumbers)= (0,0,1,1,2,2,3,3,

                        4,4,5,5,6,6,7,7,

                        8,8,9,9,10,'A',11,'B',

                        12,'C',13,'D',14,'E',15,'F');

        $_[0]=($_[0]>=0) ? $_[0] : 256+$_[0];

        my($h21)=$_[0] % 16;

        my($h22)=int($_[0] / 16);

        return $Hexnumbers{$h22}.$Hexnumbers{$h21};

   };

    sub To_Hex4{

        my(%Hexnumbers)= (0,0,1,1,2,2,3,3,

                        4,4,5,5,6,6,7,7,

                        8,8,9,9,10,'A',11,'B',

                        12,'C',13,'D',14,'E',15,'F');

   $_[0]=($_[0]>=0) ? $_[0] : 65536+$_[0];

        my($h21)=$_[0] % 16;

        my($r21)=int($_[0] / 16);

        my($h22)=$r21 % 16;

        my($r22)=int($r21 / 16);

        my($h23)=$r22 % 16;

        my($h24)=int($r22 / 16);

   return $Hexnumbers{$h24}.$Hexnumbers{$h23}.$Hexnumbers{$h22}.$Hexnumbers{$h21};

  ##return $_[0];

  };

          open (File, "Wave_ROM128.v") || die  "Cannot open r";

          while ($line = <File>)

       {

           print $line;         # Ready file output to the screen

                   }

          close(File);