DSP中浮點(diǎn)轉(zhuǎn)定點(diǎn)運(yùn)算--舉例及編程中的心得
5.舉例及編程中的心得
5.1舉例
“第3章 DSP芯片的定點(diǎn)運(yùn)算.doc”這篇文章中給了一個(gè)很簡(jiǎn)單有能說明問題的例子,不想動(dòng)大腦了,直接引用過來如下。
這是一個(gè)對(duì)語音信號(hào)(0.3kHz~3.4kHz)進(jìn)行低通濾波的C語言程序,低通濾波的截止頻率為800Hz,濾波器采用19點(diǎn)的有限沖擊響應(yīng)FIR濾波。語音信號(hào)的采樣頻率為8kHz,每個(gè)語音樣值按16位整型數(shù)存放在insp.dat文件中。
例3.7 語音信號(hào)800Hz 19點(diǎn)FIR低通濾波C語言浮點(diǎn)程序
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 語音信號(hào)800Hz 19點(diǎn)FIR低通濾波C語言定點(diǎn)程序
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)換也無非如此。所不同的是,對(duì)于一個(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è)不同嗎。
文章中難免有不妥之處,歡迎大家批評(píng)指正!
以上就是本文的全部內(nèi)容,希望能給大家一個(gè)參考,也希望大家多多支持腳本之家。
- 深入C/C++浮點(diǎn)數(shù)在內(nèi)存中的存儲(chǔ)方式詳解
- 浮點(diǎn)數(shù)在計(jì)算機(jī)中存儲(chǔ)方式是怎樣的
- DSP中浮點(diǎn)轉(zhuǎn)定點(diǎn)運(yùn)算--浮點(diǎn)與定點(diǎn)概述
- DSP中浮點(diǎn)轉(zhuǎn)定點(diǎn)運(yùn)算--浮點(diǎn)數(shù)的存儲(chǔ)格式
- DSP中浮點(diǎn)轉(zhuǎn)定點(diǎn)運(yùn)算--定點(diǎn)數(shù)的加減乘除運(yùn)算
- DSP中浮點(diǎn)轉(zhuǎn)定點(diǎn)運(yùn)算--定點(diǎn)數(shù)模擬浮點(diǎn)數(shù)運(yùn)算及常見的策略
相關(guān)文章
C語言統(tǒng)計(jì)一串字符中空格鍵、Tab鍵、回車鍵、字母、數(shù)字及其他字符的個(gè)數(shù)(Ctrl+Z終止輸入)
這篇文章主要介紹了C語言統(tǒng)計(jì)一串字符中空格鍵、Tab鍵、回車鍵、字母、數(shù)字及其他字符的個(gè)數(shù)(Ctrl+Z終止輸入) ,需要的朋友可以參考下2018-03-03基于C++ list中erase與remove函數(shù)的使用詳解
本篇文章是對(duì)C++ list中erase與remove函數(shù)的使用進(jìn)行了詳細(xì)的分析介紹,需要的朋友參考下2013-05-05c++中的volatile和variant關(guān)鍵字詳解
大家好,本篇文章主要講的是c++中的volatile和variant關(guān)鍵字詳解,感興趣的同學(xué)趕快來看一看吧,對(duì)你有幫助的話記得收藏一下2022-01-01C語言深入分析遞歸函數(shù)的實(shí)現(xiàn)
遞歸(recursive)函數(shù)是“自己調(diào)用自己”的函數(shù),無論是采用直接或間接調(diào)用方式。間接遞歸意味著函數(shù)調(diào)用另一個(gè)函數(shù)(然后可能又調(diào)用第三個(gè)函數(shù)等),最后又調(diào)用第一個(gè)函數(shù)。因?yàn)楹瘮?shù)不可以一直不停地調(diào)用自己,所以遞歸函數(shù)一定具備結(jié)束條件2022-04-04C++:函數(shù)對(duì)象,STL提供的函數(shù)對(duì)象,函數(shù)適配器詳解
這篇文章主要介紹了C++:函數(shù)對(duì)象,STL提供的函數(shù)對(duì)象,函數(shù)適配器的使用詳解,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2021-08-08求斐波那契(Fibonacci)數(shù)列通項(xiàng)的七種實(shí)現(xiàn)方法
本篇文章是對(duì)求斐波那契(Fibonacci)數(shù)列通項(xiàng)的七種實(shí)現(xiàn)方法進(jìn)行了詳細(xì)的分析介紹,需要的朋友參考下2013-05-05