【物理应用】GPS信号捕获跟踪附MATLAB代码

1 内容介绍

​在过去的几年里，紧凑型无线电导航接收器在民用市场上激增，通常以 GPS 或全球定位系统的名义出现。这些商业应用以手持设备的形式出现或内置于汽车、手机和飞机中，起源于冷战时期的军事应用。 GPS 是一个全球无线电导航系统，由 24 颗卫星组成，每颗卫星在地球上空 11,000 海里的轨道上，以及确保卫星正常工作的五个地面站。 GPS 卫星每颗绕地球运行需要 12 个小时。每颗 GPS 卫星以 50 比特/秒的速度连续广播一条导航消息，给出时间、GPS 周数和卫星健康信息（均在消息的第一部分传输）、星历表（在消息的第二部分传输）信息）和年历（信息的后半部分）。 GPS接收器计算位置需要当前时间、卫星位置和接收信号的测量延迟，然后使用三角测量规则确定接收器的位置。 GPS 误差的主要来源是大气效应（电离层、对流层）、卫星时钟误差、多径效应、卫星轨道误差和计算舍入误差。 GPS系统的误差总结在下表中。各个值不是恒定值，而是会发生变化。所有数字均为近似值 [1]：

Table 1. Errors of GPS

2 部分代码

%This Function use Ephemeris Data and Calculate satellite Position 

%CopyRight By Moein Mehrtash

%**************************************************************************                                 *

%**************************************************************************

%**************************************************************************

% Satellite Position By Ephemeris Model 

%Function's Inputs:

    %Pos_SV(m):Satellite Position Matrix

    %Pos_Rcv(m):GPS reciever Position

    %Rho(m):Pseudo Range

    

%Function's Outputs:

    %G:

    %Delta_X:

    %Pos_RCV_N:

    %B:

    

%**************************************************************************

%**************************************************************************

function [G,Delta_X,Pos_Rcv_n,B]=Gen_G_DX_XYZ_B(Pos_SV,Pos_Rcv,Rho);

[m,n]=size(Pos_SV);

d=Distance(Pos_SV,Pos_Rcv);

for i=1:m

    dif=Pos_SV(i,:)-Pos_Rcv;

    unit=dif./d(i);

    for j=1:n

        Unit_Mtrix(i,j)=unit(j);

    end

end

G=[-Unit_Mtrix ones(m,1)];

Delta_Rho=(Rho-d');    

Delta_X=inv(G'*G)*G'*Delta_Rho;

Pos_Rcv_n=(Pos_Rcv'+Delta_X(1:3))';

B=Delta_X(4);

3 运行结果

4 参考文献

[1]刘洋洋, Liu Yang-yang. 基于Matlab的GPS信号捕获与跟踪[C]// 第六届中国卫星导航学术年会. 中国卫星导航系统管理办公室学术交流中心, 2015.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机、雷达通信、无线传感器等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。