function y=circonvt(x1,x2,N)
%在x1和x2：（時(shí)域）之間的N點(diǎn)循環(huán)卷積
%--------------------------------
%[y]=circonvt(x1,x2,N)
%y=包括循環(huán)卷積的輸出序列
%x1=長(zhǎng)度N1<=N的輸入序列
%x2=長(zhǎng)度N2<=N的輸入序列
%N=循環(huán)緩沖器的大小
%方法y(n)=sum(x1(m)*x2((n-m) mod N))

%檢查x1的長(zhǎng)度
if length(x1)>N
    error('N必須>=x1的長(zhǎng)度')
end
%檢查x2的長(zhǎng)度
if length(x2)>N
    error('N必須>=x2的長(zhǎng)度')
end
x1=[x1 zeros(1,N-length(x1))];
x2=[x2 zeros(1,N-length(x2))];
m=[0:1:N-1];
x2=x2(mod(-m,N)+1);
H=zeros(N,N);
for n=1:1:N
H(n,:)=cirshift(x2,n-1,N);
end
y=x1*H';