Matlab GPS信号存在转发式欺骗干扰时，做出欺骗信号和真实卫星信号的相关峰图，其中伪码生成采用weil码

clear;clc; PRN1=Weil(1,1); %产生一个周期的伪码,10230长度,1ms周期 PRN1(PRN1==1)=-1;%计算自相关、互相关时必须将二进制转换为电平值 PRN1(PRN1==0)=1; fc=10.23e6; %码速率 fs=102.3e6; %采样率 T=1; %积分时间1ms PRN=repmat(PRN1,1,T+1);%复制扩展PRN1至1行、T+1列，每行列的元素都是PRN1 N=0.001*fs*T; %采样点数 n=[0.0000001,1:N-1]; t=n/fs; %采样点时间 A=1.5; %欺骗信号相对真实信号功率比 BB_real=PRN(ceil(t*fc)); %真实采样信号,ceil(z)取大于z的最小整数 %傅里叶变换，时域卷积=频域相乘 no_spoofing=ifft(fft(BB_real).*conj(fft(BB_real)));%conj(z)=real(Z)-i*imag(Z),即z的复共轭。.*表对应元素一一相乘 no_spoofing=fftshift(no_spoofing)/N;%fftshift可以让0频率点回到坐标原点，让图像正确显示。而fft的坐标会有问题 delay1=1/fc*3; %欺骗信号的延时时间 BB_spoofing1=sqrt(A)*PRN(ceil((t+delay1)*fc)); %欺骗信号没有重叠 spoofing1=ifft(fft(BB_real).*conj(fft(BB_spoofing1+BB_real))); spoofing1=fftshift(spoofing1)/N; delay2=1/fc*0.5; %欺骗信号的延时时间 BB_spoofing2=sqrt(A)*PRN(ceil((t+delay2)*fc)); %欺骗信号部分重叠 spoofing2=ifft(fft(BB_real).*conj(fft(BB_spoofing2+BB_real))); spoofing2=fftshift(spoofing2)/N; delay3=0; %欺骗信号的延时时间 BB_spoofing3=sqrt(A)*PRN(ceil((t+delay3)*fc)); %欺骗完全重叠 spoofing3=ifft(fft(BB_real).*conj(fft(BB_spoofing3+BB_real))); spoofing3=fftshift(spoofing3)/N; m=80; figure(1) subplot(4,1,1); plot(abs(no_spoofing(N/2-m:N/2+m))); title('(a)没有加欺骗信号时的相关峰','fontsize',14); axis([0 160 0 3]); hold on; subplot(4,1,2); plot(abs(spoofing1(N/2-m:N/2+m))); title('(b)欺骗信号与真实信号没有重叠','fontsize',14); axis([0 160 0 3]); hold on; subplot(4,1,3); plot(abs(spoofing2(N/2-m:N/2+m))); title('(c)欺骗信号与真实信号部分重叠','fontsize',14); axis([0 160 0 3]); hold on; subplot(4,1,4); plot(abs(spoofing3(N/2-m:N/2+m))); title('(d)欺骗信号与真实信号完全重叠','fontsize',14); axis([0 160 0 3]); hold on;