7 应用程序设计.pptVIP

第7章 应用程序设计 7.1.2 卷积算法的MATLAB实现 7.1.3 卷积算法的DSP实现 7.1.3 卷积算法的DSP实现 7.1.3 卷积算法的DSP实现 /* Functions */ extern void load（unsigned int loadValue）; static int processing1（int *output1, int *output2）; static int processing2（int *output2, int *output3）; static int processing3（int *input1,int *output2,int *output4）; static int processing4（int *input2, int *output1）; static void dataIO1（void）; static void dataIO2（void）; 如何优化？ 7.2. FIR滤波器的DSP实现  7.2.1 FIR滤波器的基本结构 7.2.3 FIR滤波器的MATLAB设计 fir1函数 1、线性缓冲区 线性缓冲区的C语言编程实现 线性缓冲区的C语言编程实现 线性缓冲区的C语言编程实现 int fir(int x_in, int *x, int *b, int N) { int y_out; int i; x[0]=x_in; y_out = 0.0; for(i=0;iN;i++) y_out +=b[i]*x[i]; for(i=N-1;i0;i--) x[i] = x[i-1]; return （y_out）; } 其他实现方法： x[0]=x_in; for (i = N-1; i = 0; i--) { y_out += (b[i] * x[i]); x[i] = x[i-1]; } 其他实现方法： 思考题:当FIR滤波器系数为偶对称时，程序实现过程如果修改？ FIR滤波器系数为偶对称时的程序实现 完整的程序 主程序main() 子程序 链接器命令文件 完整的程序 main() { unsigned int i; unsigned int n=0; int index=0; for(i=0;iL;i++) { in[i]=0; out[i]=0; } while(1) { in[index] =input_sample(); out[index] = fir(in, coeff, L, index); index--; if(index==-1) index=L-1; } } 完整的程序 主程序 子函数 链接器命令文件 中断程序结构 Const short N=64; int y_out=0; Short x_buffer[N]; Short buffercount=0; Short y_buffer[N] void main() { while(1); //infinite loop } 中断程序结构 interrupt void c_int11() //ISR { short i; x_buffer[0] = input_sample(); y_out = 0; //initialize filter output for (i = 0; i N; i++) y_out +=(b[i]*x[i]); for (i = N-1; i 0; i--) x_buffer[i] = x_buffer[i-1]; output_sample (y_out); //output filter; y_buffer[buffercount]=y_out; buffercount++; if (buffercount==N) buffercount=0; return; } 完整的程序 主程序 中断服务函数 链接器命令文件 中断向量表 2、循环缓冲区 循环缓冲法 循环缓冲法 循环缓冲法 ARn/CDP---线性/循环寻址