📄 dct.c
字号:
/*************************************************************
*Copyright (c) 2005,北京精仪达盛科技有限公司研发部
*All rights reserved
*
*文件名称:fir.c
*文件标示:
*摘 要:本文件内容为tms320lf2407 dsp 有限冲击响应滤波器试验
* 程序
*
*当前版本:1.0
*作 者:王飞
*完成日期:2005年2月17日
*
*取代版本:
*原作者 :
*完成日期:
*************************************************************/
/************************文件预处理***************************/
#include <math.h>
#include "exp3_2407.h"
#define pi 3.1415927
/*************************************************************/
/********************全局变量定义与初始化*********************/
int px[256];
int py[256];
double npass,h[51], x, y, xmid[51];
int m=50;
int n=256;
/*************************************************************/
/*******************函数、子程序声明与定义********************/
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;
/*选择模拟输入通道为4通道*/
* CHSELSEQ1=0x0004;
}
void read_ad() //读AD结果子程序
{
int a;
/*软件启动eoc1--ad转换器*/
* ADCTRL2|=0x2000;
for(a=0;a<256;a++)
{
/*读出ad结果*/
px[a]=(* RESULT0>>6);
}
return;
}
void firdes(int m, double npass)
{
int t;
for (t=0; t<=m; t++)
{
h[t] = sin((t-m/2.0)*npass*pi)/(pi*(t-m/2.0));
}
if (t=m/2) h[t]=npass;
}
/*************************************************************/
/*****************中断服务子程序声明与定义********************/
interrupt void nothing() //哑中断子程序
{
return;
}
/*************************************************************/
/**************************主程序*****************************/
void main(void)
{
int xm,ym;
double fs,fstop,r,rm;
int i,j,p,k,s;
sys_ini();
adc_ini();
k=0;
fs = 16000;
fstop = 2000;
npass = fstop/fs;
for (i=0; i<=m; i++)
{
xmid[i]=0;
}
for(;;)
{
firdes(m, npass);
for (s=0; s<=4; s++)
{
read_ad();
for (i=0; i<=n-1; i++)
{
xm = px[i];
x = xm/1023.0;
for (p=0; p<=m; p++)
{
xmid[m-p] = xmid[m-p-1];
}
xmid[0] = x;
r = 0;
rm= 0;
for (j=0; j<=m; j++)
{
r = xmid[j] * h[j];
rm = rm + r;
}
y = rm;
ym = (int)(1023.0 * y);
py[i] = ym;
}
k=k;
}
}
}
/***************************结束******************************/
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -