一维卷积的实现

ouble convfft(vector<double> &a, vector<double> &b, vector<double> &cx) { unsigned int sz = a.size() + b.size() - 1; cx.resize(sz); int b_len = b.size(); int a_len = a.size(); for (int n = 0; n < sz; n++) { int kmin = (n >= b_len - 1) ? n - (b_len - 1) : 0; int kmax = (n < a_len - 1) ? n : a_len - 1; for (int k = kmin; k <= kmax; k++) { cx[n] += a[k] * b[n - k]; } }

return 0;

这种是根据卷积的特点进行计算的，还有就是通过傅里叶进行求卷积；先将原来的数据扩展，长度为信号的长度+滤波器的长度-1;之后求各自的傅里叶变换，之后求乘积，在求逆傅里叶变换。这种还没有实现。