目录
1.算法描述
最早的GPS包含L1和L2两个频段，其中L1上调制CA码，P码以及导航电文，L2上调制P码和导航电文。在实际接收到的GPS信号中，我们除了能够接受到CA码和P码外，还能检测到L1和L2两种载波信号。GPS双频发送器的基本构架如下所示：
使用CA码和P码两种码来区分双频道中的两种不同的频道。但是P码周期非常长，美国用P码周期是140多天，而简化后的民用版本也要7天多，我们无法在仿真或者实际测试的时候花那么多时间去验证P码的捕获，所以这里，我们将P码部分做了下简化，使用伪随机序列周期为2048bit，来代替P码部分。
一般情况下，导航电文的频率为50hz，CA码的频率为1M，P码的频率为10M。这里，为了测试的需要，我们需要降低频率来进行测试。整个发送端，按如下的结构设计：
通过捕获模块获得初始的频偏值，然后进行载波同步。与此同时，通过码同步，完成相位的捕获。最后进入跟踪阶段。
先介绍一下捕获的基本原理：
本地码生成器以C/A码标称频率产生C/A码与接收到的采样信号相关累加，一个积分周期（通常1个码周期）后，相关峰与检测门限比较，如果相关峰大于门限，则认为捕获成功，得到对应的码相位估计；如果相关峰小于门限，码发生器自动将本地码码相位向前或向后跳动1/2或1/4个码片，然后继续相关累加检测，最多在 或 个伪码周期后找到与本地伪码同步的输入伪码的相位状态（ 即为一个码周期内码片的数目），以实现伪码的捕获。下面对步进相关法进行简单介绍，其原理图见图1。
跟踪部分：
其内部详细结构如下所示：
2.仿真效果预览
算法仿真:MATLAB 2010b
FPGA设计：ISE12.2
FPGA仿真:Modelsim6.5SE
导航电文和CA码及P码异或之后的信号。
导航电文和CA码及P码异或之后的信号，通过成型滤波器之后的效果。
这个最后发送出去的QPSK,BPSK以及相加之后的射频信号。
最后捕获跟踪之后的信号，放大看如下所示：
一开始的逐渐变大的过程就是锁定过程
3.MATLAB核心程序
-
`
timescale
1
ns /
1
ps
-
module
GPS_Rec(
-
i_clk,
-
i_rst,
-
i_QPSK,
-
i_BPSK,
-
//Capture
-
o_Ca_index,
-
o_CA,
-
o_abs_addCA,
-
o_P_index,
-
o_P,
-
o_abs_addP,
-
o_fre_est_Ca,
-
o_fre_est_P,
-
//Tracking
-
o_I_L1,
-
o_Q_L1,
-
o_I_L2,
-
o_Q_L2,
-
o_Dwen_rec_L1,
-
o_Dwen_rec_L2
-
);
-
-
input
i_clk;
-
input
i_rst;
-
input
signed
[15:0]
i_QPSK;
-
input
signed
[15:0]
i_BPSK;
-
//Capture
-
output
[9:0]
o_Ca_index;
-
output
signed
[1:0]
o_CA;
-
output
signed
[21:0]
o_abs_addCA;
-
output
[10:0]
o_P_index;
-
output
signed
[1:0]
o_P;
-
output
signed
[23:0]
o_abs_addP;
-
output
signed
[23:0]
o_fre_est_Ca;
-
output
signed
[23:0]
o_fre_est_P;
-
-
//Tracking
-
output
signed
[15:0]
o_I_L1;
-
output
signed
[15:0]
o_Q_L1;
-
output
signed
[15:0]
o_I_L2;
-
output
signed
[15:0]
o_Q_L2;
-
output
signed
[1:0]
o_Dwen_rec_L1;
-
output
signed
[1:0]
o_Dwen_rec_L2;
-
-
wire
clk_ca;
-
wire
clk_ca_2code;
-
wire
clk_p;
-
wire
clk_p_2code;
-
-
//2 time ca clock
-
CLOCK_DCM2
CLOCK_DCM2_u(
-
.i_clk (i_clk),
-
.i_rst (i_rst),
-
.o_clk_dwen (),
-
.o_clk_ca (clk_ca),
-
.o_clk_ca_2code (clk_ca_2code),
-
.o_clk_p (clk_p),
-
.o_clk_p_2code (clk_p_2code)
-
);
-
-
-
//CAPTURE
-
//CAPTURE
-
//frequency capture
-
wire
signed
[9:0]
o_Ca_index;
-
wire
signed
[1:0]
o_CA;
-
wire
signed
[21:0]
o_abs_addCA;
-
-
frequency_capture_channel1
frequency_capture_channel1_u(
-
.i_clk (i_clk),
-
.i_clk_ca (clk_ca),
-
.i_clk_ca2times (clk_ca_2code),
-
.i_rst (i_rst),
-
.i_QPSK (i_QPSK),
-
.o_fre_est (o_fre_est_Ca),
-
.o_I_filter (),
-
.o_Q_filter (),
-
.o_Ca_index (o_Ca_index),
-
.o_CA (o_CA),
-
.o_abs_addIQ (o_abs_addCA)
-
);
-
-
-
wire
signed
[10:0]
o_p_index;
-
wire
signed
[1:0]
o_p;
-
wire
signed
[23:0]
o_abs_addP;
-
frequency_capture_channel2
frequency_capture_channel2_u (
-
.i_clk (i_clk),
-
.i_clk_p (clk_p),
-
.i_clk_p2times (clk_p_2code),
-
.i_rst (i_rst),
-
.i_BPSK (i_BPSK),
-
.o_fre_est (o_fre_est_P),
-
.o_I_filter (),
-
.o_Q_filter (),
-
.o_p_index (o_P_index),
-
.o_p (o_P),
-
.o_abs_addIQ (o_abs_addP)
-
);
-
-
-
-
//TRACKING
-
//TRACKING
-
-
-
wire
signed
[15:0]
o_I_filter1;
-
wire
signed
[15:0]
o_Q_filter1;
-
wire
signed
[15:0]
o_I_filter2;
-
wire
signed
[15:0]
o_Q_filter2;
-
-
Frequency_track_tops
Frequency_track_tops_u(
-
.i_clk (i_clk),
-
.i_rst (i_rst),
-
.i_QPSK (i_QPSK),
-
.i_BPSK (i_BPSK),
-
.
i_FRE_index1(o_fre_est_Ca),
-
.
i_FRE_index2(o_fre_est_P),
-
.o_I_filter1 (o_I_L1),
-
.o_Q_filter1 (o_Q_L1),
-
.o_I_filter2 (o_I_L2),
-
.o_Q_filter2 (o_Q_L2)
-
);
-
-
//CA capture
-
wire
signed
[1:0]
Dwen_rec_L1;
-
CA_early_late_track_tops
CA_early_late_track_tops_u(
-
.i_clk (clk_ca_2code),
-
.i_clk_ca (clk_ca),
-
.i_rst (i_rst),
-
.i_CA_index (o_Ca_index+
1),
-
.i_Idin (o_I_L1),
-
.i_Qdin (o_Q_L1),
-
.o_dout (o_Dwen_rec_L1)
-
);
-
-
//P capture
-
wire
signed
[1:0]
Dwen_rec_L2;
-
P_early_late_track_tops
P_early_late_track_tops_u (
-
.i_clk (clk_p_2code),
-
.i_clk_p (clk_p),
-
.i_rst (i_rst),
-
.i_P_index (o_P_index+
1),
-
.i_Idin (o_I_L2),
-
.i_Qdin (o_Q_L2),
-
.o_dout (o_Dwen_rec_L2)
-
);
-
endmodule
-
01
_118m
4.完整MATLAB
V
转载：https://blog.csdn.net/hlayumi1234567/article/details/128258392