首页 > 应用 > EDA
[导读]摘要:在学习数字信号处理算法程序中用VC编写的几个通用算法程序。
在学习信号处理的过程中,看到书上的大部分算法都是用Fortan或者Basic实现,于是自己试验着用VC实现了一下。
1、卷积计算
离散卷积公式的算法实现

摘要:在学习数字信号处理算法程序中用VC编写的几个通用算法程序。

在学习信号处理的过程中,看到书上的大部分算法都是用Fortan或者Basic实现,于是自己试验着用VC实现了一下。

1、卷积计算

离散卷积公式的算法实现

卷积计算界面


图1 卷积计算界面

1.1 主程序代码(省略了部分不关键代码)

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;
//函数选择 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); }
//显示结果 delete[] x; delete[] y; delete[] h;}

1.2 各个子函数实现

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实现

程序界面,具体实现见注释及代码:

DFT与FFT实现界面


图2 DFT与FFT实现界面

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为底的幂级数!");
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 子函数代码

说明:其中COMP为复数类型

/*******************************************
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]; } }
//确定级别数
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级数,如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 子函数代码实现

/*********************************************************************
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、结束语

虽然现在DSP算法都有很好C语言实现。但是能够通过自己动手编写代码加深对基础知识的掌握,对自己进行数据采集器件的控制还是有很多益处的。



来源:ks990次

本文引用地址: http://www.21ic.com/app/eda/201806/769781.htm

换一批

延伸阅读

[猎聘集] 西方程序员眼中的东方程序员

西方程序员眼中的东方程序员

你对这些来自东方的程序员有没有一些总结性的看法和观点(比如:他们是否合作,是否能按时提交代码,写出的程序是否有质量?)。依据是什么?......

关键字:程序员

[猎聘集] 重磅!程序猿月薪过7万,可以落户北京了

重磅!程序猿月薪过7万,可以落户北京了

这是北京“近十年最宽松的落户政策”!消息在北京创投圈激起了涟漪。这些年,北京几乎聚集了中国数量最多顶尖的创业公司和创投机构,从IDG资本、红杉中国、经纬创投到百度、京东、美团、今日头条….浩浩荡荡数十万人。这一纸新规有望让他们落地生根。......

关键字:程序猿 月薪 北京

[猎聘集] 面对2018年春节后的跳槽高峰期,程序员该如何更新你的简历

面对2018年春节后的跳槽高峰期,程序员该如何更新你的简历

马上就到一年一度的2018年春节了,金三银四,而春节过后往往是跳槽的高峰期.作为一名程序员,应该如何准备自己的简历呢?我们程序员一般在什么时机需要更新简历呢?......

关键字:程序员 跳槽

[新鲜事] 硅谷的好莱坞化,程序员也有自己的经纪人了?!

硅谷的好莱坞化,程序员也有自己的经纪人了?!

据Businessinsider报道,在硅谷中想当一个出色的程序员,意味着与诸多精英面临直接竞争。公司给你开出了丰厚的薪水,免费的大餐,免费司机接送,硅谷的程序员在某种意义上就像是好莱坞的明星。Altay Guvench,曾经是一名程序员的......

关键字:硅谷 程序员 经纪人

[新鲜事] 8年前用一万个比特币买2个披萨,如今的他又出手了...

8年前用一万个比特币买2个披萨,如今的他又出手了...

2010年5月22日,一名为Laszlo Hanyecz的程序员用一万个比特币购买了两个披萨。当时一枚比特币价值仅为0.003美分,如今则是9669.27美元。......

关键字:比特币 程序员

[疯狂史] 西方技术人眼中的东方技术人

西方技术人眼中的东方技术人

你对这些来自东方的程序员有没有一些总结性的看法和观点(比如:他们是否合作,是否能按时提交代码,写出的程序是否有质量?)。依据是什么?......

关键字:程序员 东方 西方

[猎聘集] 程序员增加收入实用指南

程序员增加收入实用指南

我不算是聪明的人,经历也不算特别成功,但一步一步走来,我认为,我能做到的,你一定也能做到,而且应该还能做得比我更好。......

关键字:程序员 收入

[真心话] 从程序员变成优秀管理者,小米雷军的成长史

从程序员变成优秀管理者,小米雷军的成长史

“什么是成功?成功不是别人觉得你成功就是成功,成功是一种内心深处的自我感受。我只是在追求内心的一些东西,在路上。”学霸很多,能够创办出一流企业的学霸很少;手机公司很多,以软件起家的手机公司很少......

关键字:小米 雷军 程序员
条评论

我 要 评 论

网友评论

技术子站

更多

人才招聘

更多

项目外包

更多

推荐博客