循环卷积

1、  题目：循环卷积

2、  代码：

/******************************************************************* * 循环卷积 * V0.1 * * 输入参数：lpData,输入数据；DataLen,输入数据长度；Filter,卷积核； * FilterLengh,卷积核长度；AfterProData,卷积结果；Delay, 卷积核延 * 迟长度。 * * 说明：Delay有输入范围, -DataLen <= Delay < DataLen+2-FilterLen * * 李承宇, lichengyu2345@126.com * * 2010-09-02 *******************************************************************/ #include "stdafx.h" #include "stdio.h" void Convolution(double * lpData, short DataLen, double * Filter, short FilterLengh, double *AfterProData, short Delay) { double * TemFilter; short Temi, Temj, Temk = 0; short temp; double TemVal = 0; TemFilter = new double[FilterLengh]; for( Temi = 0; Temi < FilterLengh; Temi++ ) { TemFilter[Temi] = Filter[FilterLengh - Temi - 1]; } for( Temi = Delay; Temi <= DataLen - 1 + Delay; Temi++ ) { TemVal = 0; for( Temj = 0; Temj < FilterLengh; Temj++ ) { temp = Temj + Temi; if( (temp >= 0) && (temp < DataLen) ) TemVal += TemFilter[Temj] * lpData[temp]; else { if( temp < 0 ) TemVal += TemFilter[Temj] * lpData[temp+DataLen]; else TemVal += TemFilter[Temj] * lpData[temp-DataLen]; } } AfterProData[Temk] = TemVal; Temk ++; } delete [] TemFilter; } void main() { double data[8] = {1, 2, 3, 4, 5, 6, 7, 8 }; double filter[4] = {1, 1, 1, 1}; double AfterProData[8] = {0}; int nDelay; printf("Please type in nDelay:/n"); scanf("%d", &nDelay); Convolution(data, 8, filter, 4, AfterProData, nDelay); printf("/nThe result of convolution is:/n"); for(int i = 0; i < 8; i++) printf("%f/n", AfterProData[i]); }

 

3、  说明：

Delay有输入范围要求：

当Temi > 0时，要满足temp – DataLen < DataLen, 即 Temj(max)+ Temi(max) – DataLen < DataLen, 即 (Filter – 1) + (DataLen – 1 + Delay) – DataLen < DataLen, 所以有Delay < DataLen + 2 – FilterLen;

当Temi < 0 时，要满足 temp + DataLen >= 0, 即 Temj(min) + Temi(min) + DataLen >= 0, 即 0 + Delay + DataLen >= 0, 所以有 Delay >= - DataLen.

综上所述，Delay的取值范围是-DataLen <= Delay < DataLen+2-FilterLen.

如果Delay的值在上述范围之外，将取到lpData[-1]和lpData[DataLen]及之外的值，结果发生错误。

 

4、  改进：

使temp值（Temj + Temi）始终处于0到DataLen-1之间，采用如下方法：

temp = (Temj + Temi) - ( (Temj + Temi) / DataLen ) * DataLen;

       if(temp < 0) temp = DataLen + temp;

 

       改进的代码如下：

/******************************************************************* * 循环卷积 * V0.2 * * 输入参数：lpData,输入数据；DataLen,输入数据长度；Filter,卷积核； * FilterLengh,卷积核长度；AfterProData,卷积结果；Delay, 卷积核延 * 迟长度。 * * 说明：Delay可以是任何长度。 * * 李承宇, lichengyu2345@126.com * * 2010-09-02 *******************************************************************/ void Convolution(double * lpData, short DataLen, double * Filter, short FilterLengh, double *AfterProData, short Delay) { double * TemFilter; short Temi, Temj, Temk = 0; short temp; double TemVal = 0; TemFilter = new double[FilterLengh]; for( Temi = 0; Temi < FilterLengh; Temi++ ) { TemFilter[Temi] = Filter[FilterLengh - Temi - 1]; } for( Temi = Delay; Temi <= DataLen - 1 + Delay; Temi++ ) { TemVal = 0; for( Temj = 0; Temj < FilterLengh; Temj++ ) { temp = Temj + Temi; temp = (Temj + Temi) - ( (Temj + Temi) / DataLen ) * DataLen; if(temp < 0) temp = DataLen + temp; TemVal += TemFilter[Temj] * lpData[temp]; } AfterProData[Temk] = TemVal; Temk ++; } delete [] TemFilter; }

 

 

原代码出：《小波分析及其在图像处理中的应用》，陈武凡主编，科学出版社，2002年4月第一版，程序附录