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個月之久，現在終于可以了卻自己的愧疚了。從第三部分開始，寫的有寫著急，因為沒有足夠的時間來寫，還想把剩下的寫完，只能這么解決了。我想，文章只是給出學者一個思路，如果自己以后遇到了這樣的問題，還需要自己去探索的，沒有一勞永逸的文章，畢竟問題都各個不同嗎。
文章中難免有不妥之處，歡迎大家批評指正！
以上就是本文的全部內容，希望能給大家一個參考，也希望大家多多支持武林網。
