culc_coor.m

来自「跳频图案是跳频通信成功的基础」· M 代码 · 共 11 行

function [max_coor] = culc_coor(a, b)
% 计算等长序列a,b的循环互相关并绘图
L = length(a);
K = floor(L/2);
a_span = [a(L-K+1:L),a,a(1:K)];
b_span = [b(L-K+1:L),b,b(1:K)];
for k = -K : K % 周期互相关函数PACF
    PACFab(k+K+1) = 1/L * a(1:L)*b_span(k+K+1:k+K+L)';
end
PACFab_abs = abs(PACFab);
max_coor = max(PACFab_abs);