Matlab分析NMEA0183报文源代码(含高斯6度带投影计算)

2024-05-13 06:08

本文主要是介绍Matlab分析NMEA0183报文源代码(含高斯6度带投影计算),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

clc
clear all
close all

format long g

FILE_NAME = ‘D:/data6/RxCOM6_20190528gps.dat’;

fid = fopen(FILE_NAME,‘rb’);
if(fid == -1)
error(‘open file failure!’);
end

cntL_GGA = 0;
%cntL_RMC = 0;

Failure_cnt = 0;
line_cnt = 1;
GGA_line_cnt = 0; %文件中GGA所在行
RMC_line_cnt = 0;

while(1)
Line = fgetl(fid);
if(Line == -1)
break;
end

if(strncmp(Line,'$GPGGA',6))cntL_GGA = cntL_GGA + 1;            %GGA信息计数Line = strrep(Line,',',' ');GGA_temp = sscanf(Line, '%*s %f %f %*s %f %*s %d %d %f %f %*s %f', [1, inf]);GGA_array(cntL_GGA,:) = GGA_temp;%高斯6度投影iPI =pi/180;ZoneWide = 6; %6度带宽   a = 6378245; f = 1/298.257223563; %WGS84坐标系参数  ProjNo = floor((floor(GGA_temp(3)/100)+rem(GGA_temp(3),100)/60)/ZoneWide);%投影带号计算longitude0 = ProjNo * ZoneWide + ZoneWide / 2;longitude0 = longitude0 * iPI;latitude0 = 0;longitude1 = (floor(GGA_temp(3)/100)+rem(GGA_temp(3),100)/60) * iPI; %经度转换为弧度  latitude1 = (floor(GGA_temp(2)/100)+rem(GGA_temp(2),100)/60) * iPI; %纬度转换为弧度  e2 = 2 * f - f * f;ee = e2 * (1 - e2);NN = a /sqrt(1 -e2*sin(latitude1)*sin(latitude1));T = tan(latitude1) *tan(latitude1);C = ee * cos(latitude1) * cos(latitude1);A = (longitude1 - longitude0) * cos(latitude1);%参数转换M = a * ((1 - e2 / 4 - 3 * e2 * e2 / 64 - 5 * e2 * e2 * e2 / 256) * latitude1 - ...(3 * e2 / 8 + 3 * e2 * e2 / 32 + 45 * e2 * e2 * e2 / 1024) *...sin(2 * latitude1) + (15 * e2 * e2 / 256 + 45 * e2 * e2 * e2 / 1024) * ...sin(4 * latitude1) - (35 * e2 * e2 * e2 / 3072) * sin(6 * latitude1));xval = NN * (A + (1 - T + C) * A * A * A / 6 + (5 - 18 * T + T * T + 72 * C - 58 * ee) * A * A * A * A * A / 120);yval = M + NN *tan(latitude1) * (A * A / 2 + (5 - T + 9 * C + 4 * C * C) * A * A * A * A / 24 + (61 - 58 * T + T * T + 600 * C - 330 * ee) * A * A * A * A * A * A / 720);GGA_array(cntL_GGA,2)= xval+500000;GGA_array(cntL_GGA,3)= yval;

% %经纬度转空间直角坐标
% a=6378137;
% f=1/298.257223563;
% b = a * (1 - f);
% e = sqrt(aa-bb)/a;
% N = a / sqrt(1 - e * e * sin((floor(GGA_temp(2)/100)+…
% rem(GGA_temp(2),100)/60) *pi/ 180) * sin((floor(GGA_temp(2)/100)+rem(GGA_temp(2),100)/60) * pi / 180));
%
% %GGA语句中的以度/分表示的经纬度转换成度表示
% GGA_array(cntL_GGA,2) = (N +GGA_temp(7))*cos((floor(GGA_temp(2)/100)+rem(GGA_temp(2),100)/60)pi/180)
% cos((floor(GGA_temp(3)/100)+rem(GGA_temp(3),100)/60)*pi/180);
% GGA_array(cntL_GGA,3) = (N +GGA_temp(7))*cos((floor(GGA_temp(2)/100)+rem(GGA_temp(2),100)/60)pi/180)
% sin((floor(GGA_temp(3)/100)+rem(GGA_temp(3),100)/60)*pi/180);

% %GGA语句中的以度/分表示的经纬度转换成度表示
% GGA_array(cntL_GGA,2) = floor(GGA_temp(2)/100)+… %度(纬度)
% rem(GGA_temp(2),100)/60 ; %以度表示的分(纬度)
% GGA_array(cntL_GGA,3) = floor(GGA_temp(3)/100)+… %度(经度)
% rem(GGA_temp(3),100)/60 ; %以度表示的分(经度)
GGA_line_cnt(line_cnt) = 1;
end

%if(strncmp(Line,'$BDRMC',6))%cntL_RMC = cntL_RMC + 1;            %RMC信息计数%Line = strrep(Line,',',' ');%RMC_temp = sscanf(Line, '%*s %f %s %f %*s %f %*s %f %f %f %f %*s %f', [1, inf]);%RMC_array(cntL_RMC,:) = RMC_temp;%RMC_line_cnt(line_cnt) = 1;
%end   line_cnt = line_cnt + 1;

end
GGA_line_cnt = GGA_line_cnt’;
%RMC_line_cnt = RMC_line_cnt’;

%% ----------------------平面轨迹----------------------------
h=figure(1);set(h,‘name’,['平面轨迹 ’ FILE_NAME]);
plot(GGA_array(:,3),GGA_array(:,2),‘b.-’)
xlabel(‘Longitude [deg]’,‘fontsize’,12,‘fontweight’,‘bold’);
ylabel(‘Latitude [deg]’,‘fontsize’,12,‘fontweight’,‘bold’);
title({‘平面轨迹’,FILE_NAME},‘fontsize’,12,‘fontweight’,‘bold’);
grid;
%% -----------------------三维轨迹--------------------------
% h=figure(2);set(h,‘name’,['三纬轨迹 ’ FILE_NAME]);
% plot3(GGA_array(:,3),GGA_array(:,2),GGA_array(:,7),‘b.-’)
% xlabel(‘Longitude [deg]’,‘fontsize’,12,‘fontweight’,‘bold’);
% ylabel(‘Latitude [deg]’,‘fontsize’,12,‘fontweight’,‘bold’);
% zlabel(‘Height [m]’,‘fontsize’,12,‘fontweight’,‘bold’);
% title({‘三维轨迹’,FILE_NAME},‘fontsize’,12,‘fontweight’,‘bold’);
% grid;
%% -------------------------地速---------------------------------
%h=figure(3);set(h,‘name’,[‘地面速度’ FILE_NAME]);
%plot(RMC_array(:,5),‘b.-’);
% xlabel(‘Epoch [s]’,‘fontsize’,12,‘fontweight’,‘bold’);
% ylabel(‘SOG [knot]’,‘fontsize’,12,‘fontweight’,‘bold’);
% title({‘Speed of Ground’,FILE_NAME},‘fontsize’,12,‘fontweight’,‘bold’);
% grid;
%% --------------------------地面航向------------------------------
% h=figure(4);set(h,‘name’,[‘地面航向’ FILE_NAME]);
% plot(RMC_array(:,6),‘b.-’);
% xlabel(‘Epoch [s]’,‘fontsize’,12,‘fontweight’,‘bold’);
% ylabel(‘COG [knot]’,‘fontsize’,12,‘fontweight’,‘bold’);
% title({‘Course of Ground’,FILE_NAME},‘fontsize’,12,‘fontweight’,‘bold’);
% grid;
%% ----------------------------高度----------------------------
h=figure(5);set(h,‘name’,[‘高程’ FILE_NAME]);
plot(GGA_array(:,7),‘b.-’);
hold on;
plot(GGA_array(:,4),‘b.-’);
hold on;
plot(GGA_array(:,5),‘b.-’);
hold on;
xlabel(‘Epoch [s]’,‘fontsize’,12,‘fontweight’,‘bold’);
ylabel(‘Height[m]’,‘fontsize’,12,‘fontweight’,‘bold’);
title({‘height’,FILE_NAME},‘fontsize’,12,‘fontweight’,‘bold’);
grid;

fclose(fid);

这篇关于Matlab分析NMEA0183报文源代码(含高斯6度带投影计算)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/984883

相关文章

kotlin中const 和val的区别及使用场景分析

《kotlin中const和val的区别及使用场景分析》在Kotlin中,const和val都是用来声明常量的,但它们的使用场景和功能有所不同,下面给大家介绍kotlin中const和val的区别,... 目录kotlin中const 和val的区别1. val:2. const:二 代码示例1 Java

Go标准库常见错误分析和解决办法

《Go标准库常见错误分析和解决办法》Go语言的标准库为开发者提供了丰富且高效的工具,涵盖了从网络编程到文件操作等各个方面,然而,标准库虽好,使用不当却可能适得其反,正所谓工欲善其事,必先利其器,本文将... 目录1. 使用了错误的time.Duration2. time.After导致的内存泄漏3. jsO

Spring事务中@Transactional注解不生效的原因分析与解决

《Spring事务中@Transactional注解不生效的原因分析与解决》在Spring框架中,@Transactional注解是管理数据库事务的核心方式,本文将深入分析事务自调用的底层原理,解释为... 目录1. 引言2. 事务自调用问题重现2.1 示例代码2.2 问题现象3. 为什么事务自调用会失效3

找不到Anaconda prompt终端的原因分析及解决方案

《找不到Anacondaprompt终端的原因分析及解决方案》因为anaconda还没有初始化,在安装anaconda的过程中,有一行是否要添加anaconda到菜单目录中,由于没有勾选,导致没有菜... 目录问题原因问http://www.chinasem.cn题解决安装了 Anaconda 却找不到 An

Spring定时任务只执行一次的原因分析与解决方案

《Spring定时任务只执行一次的原因分析与解决方案》在使用Spring的@Scheduled定时任务时,你是否遇到过任务只执行一次,后续不再触发的情况?这种情况可能由多种原因导致,如未启用调度、线程... 目录1. 问题背景2. Spring定时任务的基本用法3. 为什么定时任务只执行一次?3.1 未启用

C++ 各种map特点对比分析

《C++各种map特点对比分析》文章比较了C++中不同类型的map(如std::map,std::unordered_map,std::multimap,std::unordered_multima... 目录特点比较C++ 示例代码 ​​​​​​代码解释特点比较1. std::map底层实现:基于红黑

Spring、Spring Boot、Spring Cloud 的区别与联系分析

《Spring、SpringBoot、SpringCloud的区别与联系分析》Spring、SpringBoot和SpringCloud是Java开发中常用的框架,分别针对企业级应用开发、快速开... 目录1. Spring 框架2. Spring Boot3. Spring Cloud总结1. Sprin

Spring 中 BeanFactoryPostProcessor 的作用和示例源码分析

《Spring中BeanFactoryPostProcessor的作用和示例源码分析》Spring的BeanFactoryPostProcessor是容器初始化的扩展接口,允许在Bean实例化前... 目录一、概览1. 核心定位2. 核心功能详解3. 关键特性二、Spring 内置的 BeanFactory

MyBatis-Plus中Service接口的lambdaUpdate用法及实例分析

《MyBatis-Plus中Service接口的lambdaUpdate用法及实例分析》本文将详细讲解MyBatis-Plus中的lambdaUpdate用法,并提供丰富的案例来帮助读者更好地理解和应... 目录深入探索MyBATis-Plus中Service接口的lambdaUpdate用法及示例案例背景

MyBatis-Plus中静态工具Db的多种用法及实例分析

《MyBatis-Plus中静态工具Db的多种用法及实例分析》本文将详细讲解MyBatis-Plus中静态工具Db的各种用法,并结合具体案例进行演示和说明,具有很好的参考价值,希望对大家有所帮助,如有... 目录MyBATis-Plus中静态工具Db的多种用法及实例案例背景使用静态工具Db进行数据库操作插入