matlab数字信号处理

function result = myXcorr(x, y)
    N = length(x);
    M = length(y);
    L = N + M - 1; % 最终结果的长度

    % 补零使得x和y长度相等
    x = [x, zeros(1, M-1)];
    y = [y, zeros(1, N-1)];

    % 线性相关
    linear_corr = conv(x, fliplr(y));

    % 循环相关
    circular_corr = zeros(1, L);
    for i = 1:L
        circular_corr(i) = sum(x .* circshift(y, i-1));
    end

    % 返回结果
    result = struct('linear_corr', linear_corr, 'circular_corr', circular_corr);
end

x = [1, 2, 3, 4, 5];
y = [3, 4, 1, 2, 5];

% 调用线性相关的xcorr函数
result = myXcorr(x, y);
linear_corr = result.linear_corr;

% 调用循环相关的xcorr函数
circular_corr = result.circular_corr;