cirshftf.m

matlab例子

function y=cirshftf(x,m,N)
%长度为N的x序列：（频域）作m个样本的循环移位
%----------------------------------------
%y=cirshftf(x,m,N)
%y=包含循环移位的输出序列
%x=长度<=N的输入序列
%m=移位样本数
%N=循环缓冲器长度
%
%方法：y(n)=idft(dft(x(n))*WN^(mk))
%如果m是个标量，则y是一个序列（行向量）
%如果m是个向量，则y是一个矩阵，每一行是一个循环移位
%x相应于向量m中的各项
%M和x不能为矩阵
if length(x)>N
   error('N必须>=x的长度')
end
x=[x zeros(1,N-length(x))];
k=[0:1:N-1];
WN=exp(-j*2*pi/N);
mk=m*k;
y=idft(dft(x,N).*(WN.^(mk)),N);