circonv.m

《数字信号处理：基于计算机的方法（第三版）》资源下载第三部分。 著者： （美）米特拉著 作译者： 孙洪等译 ISBN号： 7-121-02652-X 出版日期： 2006-06 出版

function y = circonv(x1, x2)
% Develops a sequence y obtained by the circular
% convolution of two equal-length sequences x1 and x2
L1 = length(x1);
L2 = length(x2);
if L1 ~= L2, error('Sequences of unequal length'), end
y = zeros(1,L1);
x2tr = [x2(1) x2(L2:-1:2)];
for k = 1:L1,
    sh = circshift(x2tr', k-1)';
    h = x1.*sh;
    disp(sh);
    y(k) = sum(h);
end