基于MATLAB的GNSS卫星位置的计算及精度分析_GPSmatlab_GNSS_星历文件matlab_

共13个文件

m：9个

16g：1个

16p：1个

5星 · 超过95%的资源 121 浏览量 2021-10-02 04:44:19 上传评论 18 收藏 839KB ZIP 举报

在本项目中，我们主要探讨如何使用MATLAB来处理全球导航卫星系统（Global Navigation Satellite System，简称GNSS）的数据，特别是GPS系统，进行卫星位置的计算以及对其精度进行分析。MATLAB是一种强大的数学计算和数据分析工具，对于处理复杂的数学问题，如 GNSS 数据解析，非常适用。 "广播星历"（Broadcast Ephemeris）是GNSS系统向地面用户广播的一系列参数，用于描述卫星的精确轨道信息。这些文件通常包含关于卫星的轨道参数，如位置、速度、健康状态等，以帮助接收机计算其相对于卫星的位置。在MATLAB中，我们需要编写脚本来读取这些'n'文件，将其中的二进制数据转换为可操作的轨道元素。接下来，计算卫星位置的过程涉及到以下步骤： 1. **数据解析**：从'n'文件中提取出卫星的轨道参数，如偏近点角（Argument of Perigee）、升交点径向距（Right Ascension of Ascending Node）、轨道偏心率（Eccentricity）、平均运动（Mean Motion）等。 2. **坐标转换**：将轨道参数转换为地心坐标系中的位置和速度。这可能涉及地球中心惯性坐标系（ECI）到地心地固坐标系（ECEF）的转换。 3. **时间同步**：确保计算的卫星位置与接收机接收到信号的时间一致。这需要对信号传播时间进行校正，通常通过伪距（Pseudorange）或载波相位测量来实现。 4. **位置解算**：利用卫星的坐标和接收机接收到的信号时间差，应用三角定位原理（多普勒效应或三边测量）来确定接收机的地理位置。在“精度分析”阶段，我们将关注以下几个方面： 1. **误差源分析**：包括信号传播延迟（大气折射、多路径效应）、钟差（卫星和接收机的原子钟不准确）、电离层延迟、对流层延迟等。 2. **定位误差评估**：可以使用重复观测数据来计算位置估计的均方根误差（RMSE），或者与已知参考点进行比较。 3. **精度增强技术**：如差分GPS（DGPS）、精密星历（PPP）和多频段接收等方法，可以显著提高定位精度。 4. **统计分析**：使用统计方法（如协方差分析、假设检验）来评估不同条件下的定位性能。在"基于MATLAB的GNSS卫星位置的计算及精度分析"项目中，你将学习如何实现上述步骤，并理解影响GNSS定位精度的各种因素。通过实际操作，你将能够运用MATLAB的数值计算和可视化功能，对卫星位置计算结果进行深入的分析和解释。压缩包中的文件很可能是包含MATLAB代码和可能的数据文件，用于演示如何进行上述计算和分析。你可以通过运行这些代码，亲身体验GNSS卫星位置计算的整个过程，并理解如何评估和改善定位精度。在实践中，不断提升对MATLAB和GNSS系统的理解，这对于从事相关领域的研究或工作将大有裨益。

资源详情

资源评论

资源推荐

收起资源包目录

基于MATLAB的GNSS卫星位置的计算及精度分析.zip （13个子文件）

基于MATLAB的GNSS卫星位置的计算及精度分析

gpscaculation.m 3KB

readsp3.m 2KB

bdsposition.m 2KB

bdsposition1.m 2KB

brdm0660.16p 4.96MB

glocacuation.m 2KB

f.m 501B

wum18870.sp3 417KB

brdc0660.16n 262KB

gloposition.m 2KB

gpsposition.m 2KB

brdc0660.16g 337KB

bdsgalcaculation.m 9KB

% 提取北斗与伽利略卫星参数 clear,clc navdata=struct('prn',nan,'year',nan,'month',nan,'day',nan,'hour',nan,'min'... ,nan,'sec',nan,'a0',nan,'a1',nan,'a2',nan,'aode',NaN,'crs',NaN,'dn',NaN,... 'm0',NaN,'cuc',NaN,'e',NaN,'cus',NaN,'a',NaN,'toe',... NaN,'cic',NaN,'om0',NaN,'cis',NaN,'i0',NaN,'crc',NaN,'w',NaN,... 'dom',NaN,'di', NaN,'cflg12',NaN,'weekno',NaN,'pflg12',NaN,... 'svaccuracy',NaN,'svhealth',NaN,'tgd',NaN,'aodc',NaN,'transmission',... NaN,'s1',NaN,'s2',NaN,'s3',NaN); fid = fopen('brdm0660.16p'); k=1; position=[]; bdsagal=[]; format long e while 1 line = fgetl(fid); answer =findstr(line,'END OF HEADER'); if ~isempty(answer) break; end end while 1 line = fgetl(fid); if ~isempty( findstr(line,'E11')) navdata(k).prn = str2num(line(2:4)); navdata(k).year= str2num(line(5:9)); navdata(k).month = str2num(line(10:12)); navdata(k).day = str2num(line(13:15)); navdata(k).hour = str2num(line(16:18)); navdata(k).min = str2num(line(19:21)); navdata(k).sec = str2num(line(22:23)); navdata(k).a0 = str2num(line(24:42)); navdata(k).a1 = str2num(line(44:61)); navdata(k).a2 = str2num(line(62:80)); line = fgetl(fid); % second line navdata(k).aode = str2num(line(5:23)); navdata(k).crs = str2num(line(24:42)); navdata(k).dn = str2num(line(43:61)); navdata(k).m0 = str2num(line(62:80)); line = fgetl(fid); navdata(k).cuc = str2num(line(5:23)); navdata(k).e = str2num(line(24:42)); navdata(k).cus = str2num(line(43:61)); navdata(k).a = str2num(line(62:80)) line = fgetl(fid); navdata(k).toe = str2num(line(5:23)); navdata(k).cic = str2num(line(24:42)); navdata(k).om0 = str2num(line(43:61)); navdata(k).cis = str2num(line(62:80)); line = fgetl(fid); navdata(k).i0 = str2num(line(5:23)); navdata(k).crc = str2num(line(24:42)); navdata(k).w = str2num(line(43:61)); navdata(k).dom = str2num(line(62:80)); line = fgetl(fid); navdata(k).di = str2num(line(5:23)); navdata(k).cflg12 = str2num(line(24:42)); navdata(k).weekno = str2num(line(43:61)); navdata(k).pflg12 = str2num(line(62:80)); line = fgetl(fid); navdata(k).svaccuracy = str2num(line(5:23)); navdata(k).svhealth= str2num(line(24:42)); navdata(k).tgd = str2num(line(43:61)); navdata(k).aodc = str2num(line(62:80)); line = fgetl(fid); navdata(k).transmission = str2num(line(5:23)); navdata(k).s1 = str2num(line(24:42)); navdata(k).s2 = str2num(line(43:61)); navdata(k).s3 = str2num(line(62:80)); if navdata(k).min==0 position=gpsposition( navdata(k).year,navdata(k).month,navdata(k).day,navdata(k).hour,navdata(k).min,navdata(k).sec,navdata(k).a0,navdata(k).a1,navdata(k).a2,navdata(k).crs,navdata(k).dn,navdata(k).m0,navdata(k).cuc,navdata(k).e,navdata(k).cus, navdata(k).a,navdata(k).toe,navdata(k).cic,navdata(k).om0,navdata(k).cis,navdata(k).i0,navdata(k).crc,navdata(k).w,navdata(k).dom,navdata(k).di); bdsagal(k,1)=position(1); bdsagal(k,2)=position(2); bdsagal(k,3)=position(3); k=k+1; end elseif ~isempty( findstr(line,'C01')) navdata(k).prn = str2num(line(2:4)); navdata(k).year= str2num(line(5:9)); navdata(k).month = str2num(line(10:12)); navdata(k).day = str2num(line(13:15)); navdata(k).hour = str2num(line(16:18)); navdata(k).min = str2num(line(19:21)); navdata(k).sec = str2num(line(22:23)); navdata(k).a0 = str2num(line(24:42)); navdata(k).a1 = str2num(line(43:61)); navdata(k).a2 = str2num(line(62:80)); line = fgetl(fid); % second line navdata(k).aode = str2num(line(5:23)); navdata(k).crs = str2num(line(24:42)); navdata(k).dn = str2num(line(43:61)); navdata(k).m0 = str2num(line(62:80)) line = fgetl(fid); navdata(k).cuc = str2num(line(5:23)); navdata(k).e = str2num(line(24:42)); navdata(k).cus = str2num(line(43:61)); navdata(k).a = str2num(line(62:80)) line = fgetl(fid); navdata(k).toe = str2num(line(5:23)); navdata(k).cic = str2num(line(24:42)); navdata(k).om0 = str2num(line(43:61)); navdata(k).cis = str2num(line(62:80)); line = fgetl(fid); navdata(k).i0 = str2num(line(5:23)); navdata(k).crc = str2num(line(24:42)); navdata(k).w = str2num(line(43:61)); navdata(k).dom = str2num(line(62:80)); line = fgetl(fid); navdata(k).di = str2num(line(5:23)); navdata(k).cflg12 = str2num(line(24:42)); navdata(k).weekno = str2num(line(43:61)); navdata(k).pflg12 = str2num(line(62:80)); line = fgetl(fid); navdata(k).svaccuracy = str2num(line(5:23)); navdata(k).svhealth= str2num(line(24:42)); navdata(k).tgd = str2num(line(43:61)); navdata(k).aodc = str2num(line(62:80)); line = fgetl(fid); navdata(k).transmission = str2num(line(5:23)); navdata(k).s1 = str2num(line(24:42)); navdata(k).s2 = str2num(line(43:61)); navdata(k).s3 = str2num(line(62:80)); position=bdsposition( navdata(k).year,navdata(k).month,navdata(k).day,navdata(k).hour,navdata(k).min,navdata(k).sec,navdata(k).a0,navdata(k).a1,navdata(k).a2,navdata(k).crs,navdata(k).dn,navdata(k).m0,navdata(k).cuc,navdata(k).e,navdata(k).cus, navdata(k).a,navdata(k).toe,navdata(k).cic,navdata(k).om0,navdata(k).cis,navdata(k).i0,navdata(k).crc,navdata(k).w,navdata(k).dom,navdata(k).di); bdsagal(k,1)=position(1); bdsagal(k,2)=position(2); bdsagal(k,3)=position(3); k=k+1; elseif ~isempty( findstr(line,'C06')) navdata(k).prn = str2num(line(2:4)); navdata(k).year= str2num(line(5:9)); navdata(k).month = str2num(line(10:12)); navdata(k).day = str2num(line(13:15)); navdata(k).hour = str2num(line(16:18)); navdata(k).min = str2num(line(19:21)); navdata(k).sec = str2num(line(22:23)); navdata(k).a0 = str2num(line(24:42)); navdata(k).a1 = str2num(line(43:61)); navdata(k).a2 = str2num(line(62:80)); line = fgetl(fid); % second line navdata(k).aode = str2num(line(5:23)); navdata(k).crs = str2num(line(24:42)); navdata(k).dn = str2num(line(43:61)); navdata(k).m0 = str2num(line(62:80)); line = fgetl(fid); navdata(k).cuc = str2num(line(5:23)); navdata(k).e = str2num(line(24:42)); navdata(k).cus = str2num(line(43:61)); navdata(k).a = str2num(line(62:80)); line = fgetl(fid); navdata(k).toe = str2num(line(5:23)); navdata(k).cic = str2num(line(24:42)); navdata(k).om0 = str2num(line(43:61)); navdata(k).cis = str2num(line(62:80)); line = fgetl(fid); navdata(k).i0 = str2num(line(5:23)); navdata(k).crc = str2num(line(24:42)); navdata(k).w = str2num(line(43:61)); navdata(k).dom = str2num(line(62:80)); line = fgetl(fid); navdata(k).di = str2num(line(5:23)); navdata(k).cflg12 = str2num(line(24:42)); navdata(k).weekno = str2num(line(43:61)); navdata(k).pflg12 = str2num(line(62:80)); line = fgetl(fid); navdata(k).svaccuracy = str2num(line(5:23)); navdata(k).svhealth= str2num(line(24:42)); navdata(k).tgd = str2num(line(43:61)); navdata(k).aodc = str2num(line(62:80)); line = fgetl(fid); navdata(k).transmission = str2num(line(5:23)); navdata(k).s1 = str2num(line(24:42)); navdata(k).s2 = str2num(line(43:61)); navdata(k).s3 = str2num(line(62:80)); position=bdsposition1( navdata(k).year,navdata(k).month,navdata(k).day,navdata(k).hour,navdata(k).min,navdata(k).sec,navdata(k).a0,navdata(k).a1,navdata(k).a2,navdata(k).crs,navdata(k).dn,navdata(k).m0,navdata(k).cuc,navdata(k).e,navdata(k).cus, navdata(k).a,navdata(k).toe,navdata(k).cic,navdata(k).om0,navdata(k).cis,navdata(k).i0,navdata(k).crc,navdata(k).w,navdata(k).dom,navdata(k).di); bdsagal(k,1)=position(1); bdsagal(k,2)=position(2); bdsagal(k,3)=position(3); k=k+1; elseif ~isempty( findstr(line,'C11')) navdata(k).prn = str2num(line(2:4)); navdata(k).year= str2num(line(5:9)); navdata(k).month = str2num(line(10:12)); navdata(k).day = str2num(line(13:15)); navdata(k).hour = str2num(line(16:18)); navdata(k).min = str2num(line(19:21)); navdata(k).sec = str2num(line(22:23)); navdata(k).a0 = str2num(line(24:42)); navdata(k).a1 = str2num(line(43:61)); navdata(k).a2 = str2num(line(62:80)); line = fgetl(fid); % second line navdata(k).aode = str2num(line(5:23)); navdata(k).crs = str2num(line(24:42)); navdata(k).dn = str2num(line(43:61)); navdata(k