cirshift.m

计DFT

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);