juanji.m

一个圆周卷积的代码~调用形式为y=juanji(h,x1,M) 其中

function y=juanji(h,x1,M)
x2([1:M])=x1([1:M]);
x2([M+1:2*M])=x1([1:M]);%进行周期延拓
n=length(h);%滤波器的长度
    for i=1:2*M
        sum=0;
        for j=1:n
         if i>j
             sum=sum+h(j)*x2(i-j);%卷积，相乘叠加求和
         end

        end
       
         y(i)=sum; 
        
    end
y=y(81:M+80);%取出需要长度的信号