Board logo

标题: 几个数字信号处理算法程序 [打印本页]

作者: yuyang911220    时间: 2015-10-25 12:55     标题: 几个数字信号处理算法程序

摘要
在学习数字信号处理算法程序中用VC编写的几个通用算法程序。
关键词 离散卷积 FIR
在学习信号处理的过程中,看到书上的大部分算法都是用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; }
  else
  {
    nChoseItem
    = 1; }
  for(int
  i= 0;i<2*(n+1);i++)
  {
    if(i< n+1)
    {
      x = 1;
      h = funcs[nChoseItem](i);
    }
    else
    {
      x = 0;
      h = 0;
    }
  }
  //2.y(i)=SUM(x(m)*h(i-m)) m=0..i
  for(i=0;i<2*(n+1);i++)
  {
    y = Calcy(x,h,i);
  }
  //显示结果
  delete[] x;
  delete[] y;
  delete[] h;
}1.2 各个子函数实现
typedef double
(* FFuncs)(int); //h1(x) double
funch1(int
n) { double
fbase
  = (double)4/(double)5; double fr
  = std::pow(fbase, n); return fr;
  } //h2(x)
double
funch2(int
n) { double
fpi
  = 3.1415927; return 0.5*sin((double)0.5*n);
  } //y(n)
//y(n)=
sum(x(m)*y(n-m))
m=0..n double
Calcy(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实现
程序界面,具体实现见注释及代码:
图2 DFT与FFT实现界面




欢迎光临 电子技术论坛_中国专业的电子工程师学习交流社区-中电网技术论坛 (http://bbs.eccn.com/) Powered by Discuz! 7.0.0