當(dāng)前位置:首頁 > 工業(yè)控制 > 電子設(shè)計自動化
[導(dǎo)讀]摘要:在學(xué)習(xí)數(shù)字信號處理算法程序中用VC編寫的幾個通用算法程序。 在學(xué)習(xí)信號處理的過程中,看到書上的大部分算法都是用Fortan或者Basic實現(xiàn),于是自己試驗著用VC實現(xiàn)了一下。 1、卷積計算 離散卷積公式的算法實現(xiàn) 圖

摘要:在學(xué)習(xí)數(shù)字信號處理算法程序中用VC編寫的幾個通用算法程序。

在學(xué)習(xí)信號處理的過程中,看到書上的大部分算法都是用Fortan或者Basic實現(xiàn),于是自己試驗著用VC實現(xiàn)了一下。

1、卷積計算

離散卷積公式的算法實現(xiàn)

卷積計算界面


圖1 卷積計算界面

1.1 主程序代碼(省略了部分不關(guān)鍵代碼)

void CInterVolveDlg::CalTheNumByArray() { this->UpdateData(TRUE);
FFuncs funcs[2] = {funch1,funch2};
int n = this->m_ValueN; double* x = new double[2*(n+1)];
//x(n) double* y = new double[2*(n+1)];
//y(n) double* h = new double[2*(n+1)];
//h(n)
//1.init x(n),h(n),y(n) CButton* pbtn = (CButton*) this->GetDlgItem(IDC_RADIO1);
int nChoseItem = 0;
//函數(shù)選擇 if(pbtn->GetCheck()) { nChoseItem = 0; } e
lse { nChoseItem = 1; } for(int i= 0;i<2*(n+1);i++) { if(i< n+1) { x[i] = 1; h[i] = funcs[nChoseItem](i); }
else { x[i] = 0; h[i] = 0; } }
//2.y(i)=SUM(x(m)*h(i-m)) m=0..i for(i=0;i<2*(n+1);i++) { y[i] = Calcy(x,h,i); }
//顯示結(jié)果 delete[] x; delete[] y; delete[] h;}

1.2 各個子函數(shù)實現(xiàn)

typedef double (* FFuncs)(int);
//h1(x) doublefunch1(intn) { doublefbase = (double)4/(double)5; double fr = std::pow(fbase, n); return fr; }
//h2(x)doublefunch2(intn) { doublefpi = 3.1415927; return 0.5*sin((double)0.5*n); }
//y(n)//y(n)=sum(x(m)*y(n-m))m=0..n doubleCalcy(double x[],double h[],int n) {double yvalue = 0;
for(int m= 0;m<=n;m++) { yvalue += x[m]*h[n-m]; }
return yvalue;}

2、DFT與FFT實現(xiàn)

程序界面,具體實現(xiàn)見注釋及代碼:

DFT與FFT實現(xiàn)界面


圖2 DFT與FFT實現(xiàn)界面

2.1 主程序代碼

void CFFTConversionDlg::OnBnClickedBtncal() { this->UpdateData(TRUE);
int nN = this->m_NumN;
float fF = this->m_NumF;
float fT = this->m_NumT;
bool bIsTimesof2 = false;
for(int i= 0;i<100;i++) { if(nN==(2 < < i)) { bIsTimesof2 = true; break; } }
if(!bIsTimesof2) { AfxMessageBox("N請輸入一個以2為底的冪級數(shù)!");
this->GetDlgItem(IDC_EDTN)->SetFocus();
return; } COMP* x = new COMP[nN];
//x(n) COMP* X = new COMP[nN];//X(k) initX(nN,x,fF,fT);
CButton* pRadio = (CButton*)this->GetDlgItem(IDC_RADIODFT);
if(pRadio->GetCheck()) { DFT(nN,x,X); }
else { FFT(nN,x,X); }
char buffer[256];
COMP source = X[nN-1];
sprintf(buffer,"%f+%fi",source.real(),source.imag());
CWnd* pwnd = this->GetDlgItem(IDC_EDTRET);
pwnd->SetWindowText(buffer);
CListCtrl* pList=(CListCtrl*) this->GetDlgItem(IDC_LIST1);
CListOper oper;
oper.FillList(*pList,nN,x,X);
delete[] x;
delete[] X;}

2.2 子函數(shù)代碼

說明:其中COMP為復(fù)數(shù)類型

/*******************************************
Name :DFT* Function
:Disperse Fuliye Transformation* Params
:N -- Total count of sampling points* X -- Input sequence* Return
:XN(k)=sum[x(n)*Pow(e,j2*Pi/N)]
* k,n
:0..N-1
*******************************************
/void DFT(int N,COMP x[],COMP XK[]){ double C = (2*pi)/N;
COMP t(0,0),ret(0,0);
for(int k=0;k < N;k++)
{ ret = COMP(0,0);
for(int i=0;i< N;i++) { t = COMP(cos(C*k*i),-sin(C*k*i));
ret += x[i]*t; } XK[k] = ret; } }/
*******************************************
Name
:FFT* Function
:Fast Fuliye Transformation* Params
:N -- Total count of sampling points* X -- Input sequence* Return
:XN(k)=sum[x(n)*Pow(e,j2*Pi/N)] * k,n
:0..N-1
*******************************************
/void FFT(int N,COMP X[],COMP XK[]){ int j=0; COMP U=0,W=0;
COMP* A = XK;
//Adjust sequence for(int i=0;i< N;i++) { if(i==0) { A[0] = X[0]; }
else { j=GetInverse(N,j);
A[i] = X[j]; } }
//確定級別數(shù)
for(int M=0;M< N;M++) { if((1<< M)==N) break; }
for(int L=1;L<=M;L++)//1-M級依次確定 { int LE = (int)pow(2,L);
//間隔 int LE1 = LE/2;
//W級數(shù),如W0,W1,W2... W=COMP(cos(pi/LE1),-sin(pi/LE1));
U=COMP(1,0); for(j=0;j< LE1;j++)
// { i=j; while(i< N) { int IP = i+LE1; COMP T=A[IP]*U;
A[IP]=A[i]-T;//蝶形計算 A[i]=A[i]+T; i+=LE; } U=U*W;
//不同的W次冪 } }}void initX(int N,COMP x[],float F,float T){
for(int i=0;i< N;i++) { x[i] = COMP(cos(2*pi*F*T*i),0);
}}

3.2 子函數(shù)代碼實現(xiàn)

/*********************************************************************
Name :
FuncHd* Function:
Hd()--Required frequency response function
***********************************************************************
/COMP FuncHd(double LowLimit,double UpperLimit,COMP x){
if(x.real()>UpperLimit||x.real() < LowLimit) return 0;
else return 1;
}void FIR(double LowLimit,double UpperLimit,int N,COMP Hn[])
{ int M = 2*N;
for(int i=0;i < N;i++) { Hn[i] = COMP(0,0);
for(int k=0;k < M;k++) { COMP C = COMP(cos(2*pi*i*k/(double)M),sin(2*pi*i*k/(double)M));
Hn[i] += C*FuncHd(LowLimit,UpperLimit,COMP(cos(2*pi*k/(double)M),sin(2*pi*k/(double)M)));
} Hn[i] = Hn[i]*COMP(1/(double)M,0);
}}

4、結(jié)束語

雖然現(xiàn)在DSP算法都有很好C語言實現(xiàn)。但是能夠通過自己動手編寫代碼加深對基礎(chǔ)知識的掌握,對自己進行數(shù)據(jù)采集器件的控制還是有很多益處的。



來源:ks990次

本站聲明: 本文章由作者或相關(guān)機構(gòu)授權(quán)發(fā)布,目的在于傳遞更多信息,并不代表本站贊同其觀點,本站亦不保證或承諾內(nèi)容真實性等。需要轉(zhuǎn)載請聯(lián)系該專欄作者,如若文章內(nèi)容侵犯您的權(quán)益,請及時聯(lián)系本站刪除。
換一批
延伸閱讀

9月2日消息,不造車的華為或?qū)⒋呱龈蟮莫毥谦F公司,隨著阿維塔和賽力斯的入局,華為引望愈發(fā)顯得引人矚目。

關(guān)鍵字: 阿維塔 塞力斯 華為

加利福尼亞州圣克拉拉縣2024年8月30日 /美通社/ -- 數(shù)字化轉(zhuǎn)型技術(shù)解決方案公司Trianz今天宣布,該公司與Amazon Web Services (AWS)簽訂了...

關(guān)鍵字: AWS AN BSP 數(shù)字化

倫敦2024年8月29日 /美通社/ -- 英國汽車技術(shù)公司SODA.Auto推出其旗艦產(chǎn)品SODA V,這是全球首款涵蓋汽車工程師從創(chuàng)意到認(rèn)證的所有需求的工具,可用于創(chuàng)建軟件定義汽車。 SODA V工具的開發(fā)耗時1.5...

關(guān)鍵字: 汽車 人工智能 智能驅(qū)動 BSP

北京2024年8月28日 /美通社/ -- 越來越多用戶希望企業(yè)業(yè)務(wù)能7×24不間斷運行,同時企業(yè)卻面臨越來越多業(yè)務(wù)中斷的風(fēng)險,如企業(yè)系統(tǒng)復(fù)雜性的增加,頻繁的功能更新和發(fā)布等。如何確保業(yè)務(wù)連續(xù)性,提升韌性,成...

關(guān)鍵字: 亞馬遜 解密 控制平面 BSP

8月30日消息,據(jù)媒體報道,騰訊和網(wǎng)易近期正在縮減他們對日本游戲市場的投資。

關(guān)鍵字: 騰訊 編碼器 CPU

8月28日消息,今天上午,2024中國國際大數(shù)據(jù)產(chǎn)業(yè)博覽會開幕式在貴陽舉行,華為董事、質(zhì)量流程IT總裁陶景文發(fā)表了演講。

關(guān)鍵字: 華為 12nm EDA 半導(dǎo)體

8月28日消息,在2024中國國際大數(shù)據(jù)產(chǎn)業(yè)博覽會上,華為常務(wù)董事、華為云CEO張平安發(fā)表演講稱,數(shù)字世界的話語權(quán)最終是由生態(tài)的繁榮決定的。

關(guān)鍵字: 華為 12nm 手機 衛(wèi)星通信

要點: 有效應(yīng)對環(huán)境變化,經(jīng)營業(yè)績穩(wěn)中有升 落實提質(zhì)增效舉措,毛利潤率延續(xù)升勢 戰(zhàn)略布局成效顯著,戰(zhàn)新業(yè)務(wù)引領(lǐng)增長 以科技創(chuàng)新為引領(lǐng),提升企業(yè)核心競爭力 堅持高質(zhì)量發(fā)展策略,塑強核心競爭優(yōu)勢...

關(guān)鍵字: 通信 BSP 電信運營商 數(shù)字經(jīng)濟

北京2024年8月27日 /美通社/ -- 8月21日,由中央廣播電視總臺與中國電影電視技術(shù)學(xué)會聯(lián)合牽頭組建的NVI技術(shù)創(chuàng)新聯(lián)盟在BIRTV2024超高清全產(chǎn)業(yè)鏈發(fā)展研討會上宣布正式成立。 活動現(xiàn)場 NVI技術(shù)創(chuàng)新聯(lián)...

關(guān)鍵字: VI 傳輸協(xié)議 音頻 BSP

北京2024年8月27日 /美通社/ -- 在8月23日舉辦的2024年長三角生態(tài)綠色一體化發(fā)展示范區(qū)聯(lián)合招商會上,軟通動力信息技術(shù)(集團)股份有限公司(以下簡稱"軟通動力")與長三角投資(上海)有限...

關(guān)鍵字: BSP 信息技術(shù)
關(guān)閉
關(guān)閉