#include "math.h"

#define  n    128             /*input array x,y length*/
#define  m    2*n-1
#define  pi  4.0*atan(1.0)

int i,k,j,mode;
float sum_real,t_real,temp_real;
float x_real[n],y_real[n];
float r_real[m];   

extern  c54init();
main()
{    
   unsigned int t;
   
   k=0;
   sum_real=0;   
   t_real=0;   
   temp_real=0; 
                            
   mode=1;                     /*0-無偏相關(guān)估計;1-有偏相關(guān)估計;*/

   c54init();
   
   for(i=0;i<n;i++)
   {
      x_real[i]=0;     
      y_real[i]=0;     
      r_real[i]=0;
      
   }

   for(i=0;i<n;i++)
   {
      x_real[i]=1.0*(sin(2*pi*i/n));
      y_real[i]=1.0*(cos(2*pi*i/n));      
   }

    for(k=0;k<n;k++)
   {
      sum_real=0;  
      
      for(j=0;j<n-k;j++)
      {
         temp_real=y_real[j];         

         t_real=x_real[j+k]*temp_real;         

         sum_real=sum_real+t_real;             
      }

      if(mode==0)
      {
         r_real[n-1-k]=sum_real/(float)(n-k);        
      }

      else
      {
         r_real[n-1-k]=sum_real/(float)n;        
      }
   }  
   
   for(k=0;k<n;k++)
   {
      sum_real=0;  
      
      for(j=0;j<n-k;j++)
      {
         temp_real=y_real[j+k];         

         t_real=x_real[j]*temp_real;         

         sum_real=sum_real+t_real;             
      }

      if(mode==0)
      {
         r_real[n-1+k]=sum_real/(float)(n-k);        
      }

      else
      {
         r_real[n-1+k]=sum_real/(float)n;        
      }
   }    
   
   for(;;) 
   {
     t++;
   }
}