correlation.c

本文件内容为tms320lf2407 dsp的中断向量表 将对应中断的_nothing改成眺转到指定服务子程序的指令 即可以使用此中断

/*************************************************************
 *Copyright (c) 2005,北京精仪达盛科技有限公司研发部
 *All rights reserved
 *
 *文件名称：correlation.c
 *文件标示：
 *摘    要：本文件内容为tms320lf2407 dsp相关运算实验程序     
 *          
 *
 *当前版本：1.0
 *作    者：王飞
 *完成日期：2005年2月17日
 *
 *取代版本：
 *原作者  ：
 *完成日期：
 *************************************************************/

/************************文件预处理***************************/
#include <math.h>
#include "exp3_2407.h"
#define PI  3.1415927
#define Length  128            /*input array x,y length*/
#define Lengthcor  2*Length-1  /*ouput array r length*/
/*************************************************************/

/********************全局变量定义与初始化*********************/
int i,k,j,mode;
double sum_real, t_real, temp_real;
double x_real[Length],y_real[Length];
double r_real[Lengthcor]; 
/*************************************************************/

/*******************函数、子程序声明与定义********************/
void sys_ini()               //系统初始化子程序
{
 /*关总中断*/
  asm(" setc INTM");         
 /*抑制符号位扩展*/
  asm(" clrc SXM");          
 /*累加器中结果正常溢出*/
  asm(" clrc OVM");          
 /*禁止看门狗*/
  * WDCR=0x00E8;             
 /*配置时钟锁相为4倍频CLKOUT=4*10M=40M,使能ADC时钟*/
  * SCSR1=0x0081;    
 /*io、ram、program都设为0等待读写*/
  WSGR=0x0000;                         
 /*清除所有中断标志,"写1清0"*/ 
  * IFR=0xFFFF; 
 /*使能cpu中断1*/
  * IMR=0x0001;                          
}   
/*************************************************************/

/*****************中断服务子程序声明与定义********************/
interrupt void nothing()     //哑中断子程序
{
  return;         
}
/*************************************************************/

/**************************主程序*****************************/
void main(void)
{
   sys_ini();

   k=0;
   sum_real=0;
   t_real=0;
   temp_real=0;
   
/* 
	when mode=1, result is biased estimate;
   	when mode=0, result is unbiased estimate
*/                               
   
   mode=1;
   
   for(i=0; i<Length-1; i++)   /*Initialize*/
   {
      x_real[i]=0;
      y_real[i]=0;
      r_real[i]=0;
   }

   for(i=0; i<=Length-1; i++)  /*Input x,y*/
   {
      x_real[i]=sin(2*PI*i/(Length-1));
      y_real[i]=cos(2*PI*i/(Length-1));
   }

   for(k=0; k<=Length-1; k++)  /* Caculation */
   {
   	  sum_real=0;
   	  
   	  for(j=0; j<=Length-1-k; j++)
      {
         t_real=x_real[j+k]*y_real[j];
         sum_real=sum_real+t_real;
      }

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

      else
      {
         r_real[Length-1-k]=sum_real/Length;
      }
   }
   for(k=0; k<=Length-1; k++)  
   {
      sum_real=0;

      for(j=0; j<=Length-1-k; j++)
      {
         t_real=x_real[j]*y_real[j+k];
         sum_real=sum_real+t_real;
      }

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

      else
      {
         r_real[Length-1+k]=sum_real/Length;
      }
   }
   for (;;)
   {
     i=0;
   }
}
/***************************结束******************************/