function y=cirshift(x,m,N)
%長度為N的x序列:(時域)作m個樣本的循環移位
%-------------------------------------
%[y]=cirshift(x,m,N)
%y=包含循環移位的輸出序列
%x=長度<=N的輸入序列
%m=移位樣本數
%N=循環緩沖器的長度
%方法：y(n)=x((n-m) mod N)

%檢查x的長度
if length(x)>N
    error('N必須>=x的長度')
end
x=[x zeros(1,N-length(x))];
n=[0:1:N-1];
n=mod(n-m,N);
y=x(n+1);