sin.c

CCS环境下

#include<math.h>

#define pi 3.141592654
#define D 30

#pragma DATA_SECTION(OUT,"data_buff1")

float _sinx[256] , _cosx[256];
int I[256] ,  Q[256] , FI[256] ,  FQ[256] ;
int  j_I[511] , j_Q[511]  , OUT[511];


void main(void)
{
  int i,j;
    float t;

  for(i=0;i<256;i++)
  {
    t=i*0.00001;
    _cosx[i]=cos(pi*t*t*D*152588);
    _sinx[i]=sin(pi*t*t*D*152588);

    I[i]=63*_cosx[i];
    Q[i]=63*_sinx[i];
    
  }
  
   /*序列倒序*/ 
   for(i=0;i<256;i++)
   {
       FI[i]=I[255-i];         
       FQ[i]=Q[255-i];   
   }  
   
 /*I路卷积*/       
   for(i=0;i<=510;i++)
     {
       long int t=0;
       for(j=0;j<=255;j++)
         {         
         if(i-j>=0&&i-j<=255)
            t+=I[j]*FI[i-j];        
          }          
       j_I[i]=t*0.01;
     }  
     
 
 /*Q路卷积*/  
   for(i=0;i<=510;i++)
     {
       long int t=0;
       for(j=0;j<=255;j++)
         {         
         if(i-j>=0&&i-j<=255)
            t+=Q[j]*FQ[i-j];        
          }          
       j_Q[i]=t*0.01;
     }  
     
 
     for(j=0;j<511;j++)         /*求模*/
    {
          OUT[j]=j_I[j]*0.03*j_I[j]*0.03+j_Q[j]*0.03*j_Q[j]*0.03;
          OUT[j]=20*sqrt(OUT[j]);
    }   
       
}