DSP中浮点转定点运算--举例及编程中的心得

DSP中浮点转定点运算--举例及编程中的⼼得
5．举例及编程中的⼼得
5．1举例
“第3章 DSP芯⽚的定点运算.doc”这篇⽂章中给了⼀个很简单有能说明问题的例⼦，不想动⼤脑了，直接引⽤过来如下。

这是⼀个对语⾳信号(0.3kHz~3.4kHz)进⾏低通滤波的C语⾔程序，低通滤波的截⽌频率为800Hz，滤波器采⽤19点的有限冲击响应FIR滤波。

语⾳信号的采样频率为8kHz，每个语⾳样值按16位整型数存放在insp.dat⽂件中。

例3.7 语⾳信号800Hz 19点FIR低通滤波C语⾔浮点程序
复制代码代码如下:
#include <stdio.h>
const int length = 180 /*语⾳帧长为180点＝22.5ms@8kHz采样*/
void filter(int xin[ ],int xout[ ],int n,float h[ ]); /*滤波⼦程序说明*/
/*19点滤波器系数*/
static float h[19]=
{0.01218354,-0.009012882,-0.02881839,-0.04743239,-0.04584568,
-0.008692503,0.06446265,0.1544655,0.2289794,0.257883,
0.2289794,0.1544655,0.06446265,-0.008692503,-0.04584568,
-0.04743239,-0.02881839,-0.009012882,0.01218354};
static int x1[length+20];
/*低通滤波浮点⼦程序*/
void filter(int xin[ ],int xout[ ],int n,float h[ ])
{
int i,j;
float sum;
for(i=0;i<length;i++) x1[n+i-1]=xin[i];
for (i=0;i<length;i++)
{
sum=0.0;
for(j=0;j<n;j++) sum+=h[j]*x1[i-j+n-1];
xout[i]=(int)sum;
}
for(i=0;i<(n-1);i++) x1[n-i-2]=xin[length-1-i];
}
/*主程序*/
void main( )
{
FILE *fp1,*fp2;
int frame,indata[length],outdata[length];
fp1=fopen(insp.dat,"rb"); /*输⼊语⾳⽂件*/
fp2=fopen(outsp.dat,"wb"); /*滤波后语⾳⽂件*/
frame=0;
while(feof(fp1)==0)
{
frame++;
printf("frame=%d/n",frame);
for(i=0;i<length;i++) indata[i]=getw(fp1); /*取⼀帧语⾳数据*/
filter(indata,outdata,19,h); /*调⽤低通滤波⼦程序*/
for(i=0;i<length;i++) putw(outdata[i],fp2); /*将滤波后的样值写⼊⽂件*/
}
fcloseall( ); /*关闭⽂件*/
return(0);
}
例3.8 语⾳信号800Hz 19点FIR低通滤波C语⾔定点程序
复制代码代码如下:
#include <stdio.h>
const int length=180;
void filter(int xin[ ],int xout[ ],int n,int h[ ]);
static int h[19]={399,-296,-945,-1555,-1503,-285,2112,5061,7503,8450,
7503,5061,2112,-285,-1503,-1555,-945,-296,399}; /*Q15*/
static int x1[length+20];
/*低通滤波定点⼦程序*/
void filter(int xin[ ],int xout[ ],int n,int h[ ])
{
int i,j;
long sum;
for(i=0;i<length;i++) x1[n+i-1]=xin[i];
for (i=0;i<length;i++)
{
sum=0;
for(j=0;j<n;j++) sum+=(long)h[j]*x1[i-j+n-1];
xout[i]=sum>>15;
}
for(i=0;i<(n-1);i++) x1[n-i-2]=xin[length-i-1];
}
主程序与浮点的完全⼀样。

5．2编程中的⼼得
通过上边的例⼦，相信⼤家都⼤致明⽩了转换的过程。

⼀个函数模块如此，⼀个⼤的软件项⽬的转换也⽆⾮如此。

所不同的是，对于⼀个项⽬，很少有⼀个个单独的模块，⼀般都是⼏个或数个模块都有联系，千⼀发⽽动全⾝，所以，再进⾏转换之前，最好理清各个软件模块之间的关系。

协调，协调，再协调。

写到这⾥，与afreez⼀起学习DSP中浮点转定点运算这个问题的讨论终于可以收笔了。

拖拖拉拉，拉了有3个⽉之久，现在终于可以了却⾃⼰的愧疚了。

从第三部分开始，写的有写着急，因为没有⾜够的时间来写，还想把剩下的写完，只能这么解决了。

我想，⽂章只是给出学者⼀个思路，如果⾃⼰以后遇到了这样的问题，还需要⾃⼰去探索的，没有⼀劳永逸的⽂章，毕竟问题都各个不同吗。

⽂章中难免有不妥之处，欢迎⼤家批评指正！
以上就是本⽂的全部内容，希望能给⼤家⼀个参考，也希望⼤家多多⽀持。