function y = cirshftt(x,m,N)
% 長度為 N 的x序列: (時域)作m采樣點圓周移位
% -------------------------------------------------------------------
% [y] = cirshftt(x,m,N)
% y = 包含圓周移位的輸出序列
% x = 長度 <= N的輸入序列
% m = 移位采樣數 
% N = 圓周緩沖器長度 
% 方法: y(n) = x((n-m) mod N)

% Check for length of 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);