第15章DSP的C语言编程.pptVIP

山东大学控制学院 张东亮 思考题与习题 DSP C语言有哪些特点？ C24x DSP编译器有哪些数据类型？ 如何访问片内外设寄存器的某些位？ 如何直接访问存储器单元？ pragma编译预处理命令有什么用途？ C语言与汇编语言混合编程有哪些方法？ C24x DSP的C编译器扩展了哪几个关键字？ 1个LED指示灯接到DSP的通用I/O引脚GPIOB4。采用通用定时器T1中断方式定时200ms, 用C语言编程使之闪烁, XCLKIN=10MHz, CLKOUT=40MHz。 如何用DSP C语言编程实现常用的FFT与FIR信号处理算法？ //每种蝶形运算在某一级中需要进行N/pow(2,m)次 { //蝶形运算展开，结果的实部和虚部 //分别存储在原实部和虚部位置 X=Xr[k+B]*cos_tab[S]+Xi[k+B]*sin_tab[S]; Y=Xi[k+B]*cos_tab[S]-Xr[k+B]*sin_tab[S]; Xr[k+B]=Xr[k]-X; Xi[k+B]=Xi[k]-Y; Xr[k]=Xr[k]+X; Xi[k]=Xi[k]+Y; } } } for ( m=0;mN/2;m++ ) { w[m]=sqrt(Xr[m]*Xr[m]+Xi[m]*Xi[m]); //计算功率普 } } main() { int i; InitForFFT( ); //FFT初始化函数 MakeWave( ); //波形发生函数 for ( i=0;iN;i++ ) { fWaveR[i]=INPUT[i]; fWaveI[i]=0.0; w[i]=0.0; } Mum=(int)(0.5+log(N)/log(2)); //Mum为蝶形运算的级数, N=2^Mum FFT(fWaveR,fWaveI); for ( i=0;iN;i++ ) DATA[i]=w[i]; while ( 1 ); } void InitForFFT() //FFT初始化函数，建立正余弦函数表 { int i; for ( i=0;iN;i++ ) { sin_tab[i]=sin(PI*2*i/N); cos_tab[i]=cos(PI*2*i/N); } } void MakeWave() //波形发生函数 { int i; for ( i=0;iN;i++ ) { INPUT[i]=sin(PI*2*i/N*3)*1024; //f=3Hz, 正弦函数 } } FFT函数包含的函数finv(N,Xr,Xi)为倒序运算，函数代码如下。 //倒序运算函数finv(N1,Xr,Xi)，对输入序列倒序 //N1为序列长度; Xr[ ], Xi[ ]分别为输入序列的实部和虚部 //倒序原理：倒序数的加1是在最高位加1， //满2向次高位进1，最高位变0，依次往下 //从当前倒序值可求下一倒序值 void finv(int N1, float *xr, float *xi) //倒序运算函数f(N1,Xr,Xi)，对输入序列倒序 { int m, n,N2,k; // m为正序数；n为到序数；k为各个权值；N2为最高位的权值 float T; //临时变量T N2=N1/2; //最高位加1相当于十进制加上最高位的权N1/2 n=N2; //第一个倒序值 for (m=1; m=N1-2;m++) //第0个和最后一个不倒序 { if(mn) //为了避免再次调换，只需对mn的部分调换顺序 { T=xr[m]; xr[m]=xr[n];xr[n]=T; T=xi[m]; xi[m]=xi[n];xi[n]=T; } k=N2; //最高位权值 while (n=k) { n=n-k; //次高位位1，继续上下进位，满2置0 k=(int)(k/2+0.5); //向下权值依次比上级减半 } n=n+k; //得到下一倒序值 } } 13.5.2 FIR数字滤波器 在数字信号处理中，数字滤波占有极其重要的地位。无限冲击响应(Finite Impulse Response, FIR)数字滤波器(Digital Filter)是一种常用数字信号处理算法。利用窗函数法设计FIR滤波器，可以实现线性相位的数字滤波器 。 1. FIR数字滤波器的设计方法 设FIR数字滤波器的单位冲击响应为h(n)， 则传递函数H(z)为 FIR数字滤波器的系数h(n) w(n)为窗函数。理想单位冲击响应h1(n)可以根据给定的理想频率响应求得 FIR数字滤波器的差分方程为 式中，x(i)为输入序列，y(i)为输出序列, N为滤波器阶数