dct.c

DSP2407实现的离散余弦变换和逆变换,5个通道内置定时AD采集,在语音和图象处理方面广泛使用,算法成熟,有助于在运用中理论学习

/*
;***********************************************************************
;*				  						 北京达盛科技有限公司
;*                                     	    研    发    部
;*
;*                                 	http://www.techshine.com
;*
 *文件名称：dct.c
 *文件标示：
 *摘    要：本文件内容为tms320lf2407 dsp 离散余弦变换试验程序     
 *************************************************************/

/************************文件预处理***************************/
#include <math.h>
#include "exp4_2407.h"
#define pi 3.1415927
/*************************************************************/

/********************全局变量定义与初始化*********************/
  int px[128];
  int pz[128];		
  double x[128],y[128],z[128];
  int n=128;
/*
  x为输入序列，y为DCT变换输出序列，z为重构序列，序列长度为128
  n为DCT变换序列长度
  dct1c2为DCT变换子程序，idct1c2为DCT逆变换子程序
*/ 
/*************************************************************/

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

void adc_ini()                  //ADC初始化子程序
{    
    /*仿真器挂起时,完成最后一次ad,
      高中断优先级,双排序模式,禁用其他模式*/
     * ADCTRL1=0x25C0;          
    /*禁用中断模式，清中断标志*/
     * ADCTRL2=0x0202;          
    /*启用排序器1，最大转换数为1*/
     * MAXCONV=0x0000;          
    /*选择模拟输入通道为5通道*/
     * CHSELSEQ1=0x0005;        
}

void read_ad()                  //读AD结果子程序
{
     int a;
    /*软件启动eoc1--ad转换器*/
     * ADCTRL2|=0x2000;         
     for(a=0;a<256;a++)
     {       
      /*读出ad结果*/
       px[a]=(* RESULT0>>6);    
     }
     return;
}

void dct1c2 (double x[], double y[], int n)
  {
  		double s,t;
  		int i,j;
  		s=0.0;
   	for (i=0; i<=n-1; i++)
   	  {
   		s=s+x[i];
   	  }
   	
   	for (i=1; i<=n-1; i++)
      {
   		t=0.0;
   		for (j=0; j<=n-1; j++)
   		  {
   			t=t+x[j]*cos((2.0*j+1)*i*pi/(2.0*n));
   		  }
   		y[i]=sqrt(2.0/n)*t;
   	   }
   	y[0]=s/(sqrt(n/1.0));
  }
  
  void idct1c2(double y[], double z[], int n)
  {
  		double t;
  		int i,j;
  		for (i=0; i<=n-1; i++)
   	 	  {
		    t=0.0;
		    for (j=0; j<=n-1; j++)
		     {
		  	   t=t+y[j]*cos((2.0*i+1)*j*pi/(2.0*n));
		     }
      		z[i]=y[0]/(sqrt(n/1.0))+sqrt(2.0/n)*t;
      	  }
  		
  }
/*************************************************************/

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

/**************************主程序*****************************/
void main(void)
{
		int k=0;
		int i,j;		
		int xm,zm; 
		double xmid=0;
  sys_ini();
  adc_ini();  
  for(;;)  
  {
		  xmid=0;
		  for (j=0; j<=2; j++)
		  {
			read_ad();
			for (i=0; i<=127; i++)
			{
				xm = px[i];
				xmid=xmid+xm;
				x[i] = xm;
			}		
		    
		    xmid=xmid/128.0;
		    for(i=0;i<128;i++)
		    {
		      x[i]=x[i]-xmid;
		    }
			dct1c2 (x,y,n);		
			idct1c2 (y,z,n);
				
			for (i=0; i<=127; i++)
			{
				zm=(int)z[i];
				pz[i] = zm;
			}	
		   k=k;
	     }		
  }
}
/***************************结束******************************/