P4码仿真(含自相关结果)
P1-P3码的原理介绍见之前的博客:https://blog.csdn.net/m0_47074531/article/details/113621351
https://blog.csdn.net/m0_47074531/article/details/120544775
https://blog.csdn.net/m0_47074531/article/details/120545006
clc;
clear all;
close all;
%————————————code of P4—————————————%
N=64;
for k = 1:N
phase(k)=((pi/N)*(k-1)^2)-(pi*(k-1)); % P4 Phase
end
figure(1);
plot(phase);
title(['P4 Code Phase Shift']);
grid on;
%————————————画出相位图—————————————%
xx=0:length(phase)-1;
figure(2);
stairs(xx,phase);grid
title([ 'P4码时相特性曲线']);
xlabel('k');
ylabel('P4 相位 (rad)');
%——————————求 P4 自相关特性——————————%
un=rem(phase, 2*pi);
P4_signal=exp(j*un(1:1:N)); %P4 signal
[a,b]=xcorr(P4_signal);
d=abs(a);
d=20*log10(d+1e-6);
figure(3)
plot(b,d);
title('P4码自相关函数');
grid on;