function fn=circonvt(x1,x2,N)
%circonvt函數實現輸入序列x1和x2的循環卷積,fn為輸出序列
%N 為循環卷積長度
%實現方法:fn=sum(x1(m)*x2((n-m) mod N))
if (length(x1)>N|length(x2)>N)%判斷輸入信號的長度
error('N的長度必須大于輸入數據的長度');
end
x1=[x1,zeros(1,N-length(x1))];
x2=[x2,zeros(1,N-length(x2))];
fn=ifft(fft(x1,N).*fft(x2,N),N);