DSP中浮點(diǎn)轉(zhuǎn)定點(diǎn)運(yùn)算--舉例及編程中的心得

更新時(shí)間：2016年06月17日 10:42:44 作者：ganxingming

本文主要講解DSP浮點(diǎn)轉(zhuǎn)定點(diǎn)運(yùn)算舉例及編程中的心得，具有參考價(jià)值，需要的朋友可以參考一下。

5．舉例及編程中的心得

5．1舉例

　　“第3章 DSP芯片的定點(diǎn)運(yùn)算.doc”這篇文章中給了一個(gè)很簡單有能說明問題的例子，不想動(dòng)大腦了，直接引用過來如下。

　　這是一個(gè)對語音信號(0.3kHz~3.4kHz)進(jìn)行低通濾波的C語言程序，低通濾波的截止頻率為800Hz，濾波器采用19點(diǎn)的有限沖擊響應(yīng)FIR濾波。語音信號的采樣頻率為8kHz，每個(gè)語音樣值按16位整型數(shù)存放在insp.dat文件中。

例3.7 語音信號800Hz 19點(diǎn)FIR低通濾波C語言浮點(diǎn)程序

復(fù)制代碼代碼如下:

#include <stdio.h>

const  int length = 180 /*語音幀長為180點(diǎn)＝22.5ms@8kHz采樣*/

void  filter(int xin[ ],int xout[ ],int n,float h[ ]); /*濾波子程序說明*/

/*19點(diǎn)濾波器系數(shù)*/

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];

/*低通濾波浮點(diǎn)子程序*/

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);        /*取一幀語音數(shù)據(jù)*/

        filter(indata,outdata,19,h);            /*調(diào)用低通濾波子程序*/

        for(i=0;i<length;i++)  putw(outdata[i],fp2);     /*將濾波后的樣值寫入文件*/

    }

    fcloseall( );                   /*關(guān)閉文件*/

    return(0);

}

例3.8 語音信號800Hz 19點(diǎn)FIR低通濾波C語言定點(diǎn)程序

復(fù)制代碼代碼如下:

#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];

/*低通濾波定點(diǎn)子程序*/

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];

}

主程序與浮點(diǎn)的完全一樣。

5．2編程中的心得

　　通過上邊的例子，相信大家都大致明白了轉(zhuǎn)換的過程。一個(gè)函數(shù)模塊如此，一個(gè)大的軟件項(xiàng)目的轉(zhuǎn)換也無非如此。所不同的是，對于一個(gè)項(xiàng)目，很少有一個(gè)個(gè)單獨(dú)的模塊，一般都是幾個(gè)或數(shù)個(gè)模塊都有聯(lián)系，千一發(fā)而動(dòng)全身，所以，再進(jìn)行轉(zhuǎn)換之前，最好理清各個(gè)軟件模塊之間的關(guān)系。

　　協(xié)調(diào)，協(xié)調(diào)，再協(xié)調(diào)。

　　寫到這里，與afreez一起學(xué)習(xí)DSP中浮點(diǎn)轉(zhuǎn)定點(diǎn)運(yùn)算這個(gè)問題的討論終于可以收筆了。拖拖拉拉，拉了有3個(gè)月之久，現(xiàn)在終于可以了卻自己的愧疚了。從第三部分開始，寫的有寫著急，因?yàn)闆]有足夠的時(shí)間來寫，還想把剩下的寫完，只能這么解決了。我想，文章只是給出學(xué)者一個(gè)思路，如果自己以后遇到了這樣的問題，還需要自己去探索的，沒有一勞永逸的文章，畢竟問題都各個(gè)不同嗎。

　　文章中難免有不妥之處，歡迎大家批評指正！

以上就是本文的全部內(nèi)容，希望能給大家一個(gè)參考，也希望大家多多支持腳本之家。

您可能感興趣的文章: