基于MATLAB实现的对GPS信号多普勒频率和伪码进行二维搜索后，检测峰峰值数量，及判别是否存在欺骗信号+使用说明文档.zip

clear all; close all; %%%%%%%%%% 二维搜索矩阵数据 %%%%%%%%%%%%%%%%%% yym = importdata('data_spoofing.mat');%原始二维搜索矩阵数据 %%%%多普勒频率和伪码二维搜索数据的三维图%%%%% figure; [dop_n,n]= size(yym); [x, y]=meshgrid(1:n, 1:dop_n); mesh(x,y,yym); %%%%%%%%%%%% 峰值及次峰值信息 %%%%%%%%%%%%%%%% data_yym = yym; peak = zeros(3,3); %用于存放峰值和其坐标；及次峰值和其坐标；次次峰值和其坐标 peak_sum = 0; magMean = mean(mean(yym)); %%二维搜索矩阵的均值 % if amp >5.0*magMean % peak_sum = peak_sum + 1; for i = 1:length(peak) [amp crw] = max(max(data_yym')); [amp crn] = max(max(data_yym)); peak(i,:) = [amp,crw,crn]; %存放峰值和其坐标；及次峰值和其坐标 data_yym(crw,crn) = 0; end if amp >5.0*magMean %门限判决 如何设定门限？？？？？？？ gap = abs(peak(1,3)-peak(2,3)); if gap >= 4 %%差值大于两个码片 out = 1; %% 存在欺骗信号 fprintf('存在欺骗信号\n'); else out =0; fprintf('不存在欺骗信号'); end else Acquisition = 0; fprintf('不存在卫星信号'); end % % % % % % % bbb=yym(:,818)'; % % % figure (2); % % % t=1:41; % % % % plot(t,bbb); % % % % % % figure (3); % % t=1:2048; % % % % plot(t,yym); % % % % % % bbb1=yym(26,:)'; % % figure (4); % % t=1:2048; % % % % plot(t,bbb1,'-*'); % % % % % for i = 1: length(peak(:,1)) % % % if peak(i,1)>5.0*magMean % % % peak_sum = peak_sum + 1; % % % end % % % end % % % % % peak_sum % % % % % figure; % % % [dop_n,n]= size(yym); % % % [x, y]=meshgrid(1:n, 1:dop_n); % % % mesh(x,y,yym);