hmain.m

gps/galileo信号源频谱分析

global time_unit; % 一个完整扩频码周期
global time; % 数据跳变时间单位
global fs; % 采样率
global nn; % 一个完整扩频周期采样点数
global kk; % 数据总采样点
global F_if;   % 载波中频
global F_if_before; 
global CA_freq; % PN码速率 BOC(m,n)中n对应的速率
global tc; 
global CA ;
global F_Carrier;
global NOISE_POWER;
global CA_length;
% %==================== boc(2,2)========================%

time = 20*(10^(-3));
time_unit = 4*(10^(-3));

F_if_before = 4.309*(10^6);
F_if =1.405*(10^6);
F_Carrier = 1575.42*(10^6);
CA_freq = 1.023*(10^6);
fc = 2*CA_freq;
CA_length = 1023/CA_freq;

fs = 5*(10^6);
nn = CA_length*fs;
kk = (time/CA_length)*nn;

bz = 1; 
PN = CreateBOC(bz);
CA = [];

k = 5;
for n = 1:length(PN)
CA((1+k*(n-1)):k*n) = PN(n)*ones(1,k);
end
CA;
tc = 1/(k*fc); 

NOISE_POWER = 100;

% %==================== boc(2,1)========================%
% time_unit = 1*(10^(-3));
% time = 4*(10^(-3));
% fs = 1.023*4*3*(10^6);
% nn = 1*(10^(-3))*fs;
% kk = (time/time_unit)*nn;
% F_if = 1.023*4*(10^6); % carrier_if 
% CA_freq = 1.023*(10^6);
% fc = 1.023*4*(10^6);
% bz = 2;
% PN = CreateBOC(bz);
% CA = [];
% k = 6;
% for n = 1:length(PN)
% CA((1+k*(n-1)):k*n) = PN(n)*ones(1,k);
% end
% CA;
% tc = 1/(k*fc)