数字信号处理中讲到圆周卷积,试着写了下matlab程序。
根据圆周卷积的公式:
x1(m)不需要变,难点是我们如何将x2((n-m))表示出来,既要沿着y轴反转,还要周期拖延,并且用矩形序列将0-N-1数去读取出来就ok了。
灵活运用了mod函数的功能,就可以达到目的。
这里顺便说下mod与rem的区别,对于整数和0,mod与rem是完全一样的,但对于负数就有区别了,区别在于它们的定义:
rem(X,Y) if Y~= 0, returns X - n.*Y where n= fix(X./Y)
mod(X,Y) if Y~= 0, returns X - n.*Y where n= floor(X./Y)