当前位置：文档之家 › MATLAB解析GPS数据程序

MATLAB解析GPS数据程序

  • 格式：doc
  • 大小：58.50 KB
  • 文档页数：16

页眉内容 页脚内容16 % 注：本程序可直接在MATLAB 2017a 中运行 %该脚本文件用于学习GPS数据的读取，需要做其他用途请自行修改代码 %本脚本文件的前面几行代码是要设置的一些参数 %默认使用COM3（需视情况修改） %波特率设为GPS模块默认的38400 %下面为程序源码 clear num_execute = 100; % 执行次数 num_SingleRead = 150; %单次从串口读取的字节数(最好设置足够大（最低大概设为80），保证单次读取的数据包含一条完整的GPS数据)

Timedelay = 0.2; % 用于延时读取串口数据 BaudRate = 38400; % 读取数据的波特率 Terminator = 'CR'; num_MaxTry = 5; %打开串口的最多尝试次数 BytesAvailableFcnCount = 1000; %% 设置参数 % delete(instrfindall); % 串口打开失败时使用此句 % delete(s);clear s % 串口打开失败时使用此句 serial3 = serial('COM3'); % 串口设置 serial3.BytesAvailableFcnMode = 'byte'; % serial3.InputBufferSize = 38400; % 输出波特率 serial3.BaudRate = BaudRate; % 读入波特率 % serial3.OutputBufferSize = 1024; serial3.BytesAvailableFcnCount = BytesAvailableFcnCount; serial3.ReadAsyncMode = 'continuous'; 页眉内容 页脚内容16 serial3.Terminator = Terminator; %% 打开串口 count_opentimes = 1; while contains(serial3.status,'closed') > 0 && count_opentimes < num_MaxTry fopen(serial3); %打开串口 count_opentimes = count_opentimes+1; end if contains(serial3.status,'open') < 1 disp('open com failed!'); return end %% 读取并处理数据 % 初始化 GPS_Data = GPS_Init(); while(num_execute > 0) GPS_DataStrs = fread(serial3,num_SingleRead,'char'); %一次读出10个字符 GPS_DataStrs = reshape(GPS_DataStrs,1,[]); GPS_DataStrs = split_str2strs(GPS_DataStrs); GPS_Data_tmp = get_GPS_specificData(GPS_DataStrs); GPS_Data = Updata_GPU_Data(GPS_Data,GPS_Data_tmp); show_GPS_Data(GPS_Data); pause(Timedelay); % 延时 num_execute = num_execute-1; end % fprintf(s,'abcd'); %给串口的发送数据 % fscanf(s); %从串口的接收缓存读数据 %% 关闭串口并删除相关数据 fclose(serial3); %关闭串口 页眉内容 页脚内容16 delete(serial3); clear serial3 %% %将字符串根据'\r\n'划分成多个子字符串,同时去掉首尾无用的残余字符串 function out_strs = split_str2strs(StrData) if contains(class(StrData),'char') uint8(StrData); end record = get_pos_enterflag(StrData); if StrData(1) == uint8('$') %开头为'$'的情况 flag_start = 1; else if size(record,2) > 0 flag_start = record(1)+2; else out_strs = cell(0,0); return end end if StrData(end) == 13 flag_end = length(StrData)-1; else if size(record,2) > 0 flag_end = record(end)-1; end end if flag_start >= flag_end out_strs = cell(0,0); return 页眉内容 页脚内容16 end StrData = StrData(flag_start:flag_end); % 截取有效数据，方便下面划分子字符串 record = get_pos_enterflag(StrData); num_strs = size(record,2)+1; out_strs = cell(num_strs,1); if num_strs > 1 out_strs{1,1} = char(StrData(1:record(1)-1)); if num_strs == 2 out_strs{num_strs,1} = char(StrData(record(1)+2:end)); else for i = 2 : num_strs-1 out_strs{i,1} = char(StrData(record(i-1)+2:record(i)-1)); end out_strs{num_strs,1} = char(StrData(record(i)+2:end)); end else out_strs{1,1} = char(StrData); end % 得到字符串中'\r\n'在字符串中的位置(实际为'\r'的位置) function record = get_pos_enterflag(data) record = []; % 记录回车符号位置 for ii = 1 : length(data)-1 if data(ii) == 13 if data(ii+1) == 10 record = [record,ii]; ii = ii+1; end end end 页眉内容 页脚内容16 end end % 得到具体GPS结构体数据 function GPS_Data_tmp = get_GPS_specificData(StrsData) GPS_Data_tmp = []; num_str = size(StrsData,1); for i = 1 : num_str str_tab = StrsData{i,1}; if contains(str_tab,'GGA') > 0 GPS_Data_tmp = GNGGA(str_tab); elseif contains(str_tab,'GSA') > 0 GPS_Data_tmp = GNGSA(str_tab); elseif contains(str_tab,'GSV') > 0 GPS_Data_tmp = GNGSV(str_tab); elseif contains(str_tab,'RMC') > 0 GPS_Data_tmp = GNRMC(str_tab); elseif contains(str_tab,'VTG') > 0 GPS_Data_tmp = GNVTG(str_tab); elseif contains(str_tab,'GLL') > 0 GPS_Data_tmp = GNGLL(str_tab); end end end % GPS字符串解析 function GPS_Data_tmp = GNGGA(str_tab) index = strfind(str_tab,','); count = 1; Time = str_tab(index(count)+1:index(count+1)-1);count=count+1; Latitude = str_tab(index(count)+1:index(count+1)-1);count=count+1; 页眉内容 页脚内容16 GPS_Data_tmp.LatitudeDir = str_tab(index(count)+1:index(count+1)-1);count=count+1;

Longitude = str_tab(index(count)+1:index(count+1)-1);count=count+1;

GPS_Data_tmp.LongitudeDir = str_tab(index(count)+1:index(count+1)-1);count=count+1;

GPS_Data_tmp.GPSState = str_tab(index(count)+1:index(count+1)-1);count=count+1;

GPS_Data_tmp.SatelliteNum = str_tab(index(count)+1:index(count+1)-1);count=count+1;

GPS_Data_tmp.HDOP = str_tab(index(count)+1:index(count+1)-1);count=count+1;

GPS_Data_tmp.altitude = str_tab(index(count)+1:index(count+1)-1);count=count+1;

% other = str_tab(index(count)+1:end); % 进一步处理 GPS_Data_tmp.Time.hour = Time(1:2); GPS_Data_tmp.Time.min = Time(3:4); GPS_Data_tmp.Time.sec = Time(5:6); GPS_Data_tmp.Time.millisec = Time(8:10); GPS_Data_tmp.Latitude.degree = Latitude(1:2); % 纬度 GPS_Data_tmp.Latitude.min = Latitude(3:4); tmp = str2double(Latitude(6:9)); tmp = tmp*6/1000; % tmp = tmp/10000*60; GPS_Data_tmp.Latitude.sec = num2str(floor(tmp)); GPS_Data_tmp.Latitude.millisec = num2str((tmp-floor(tmp))*10000); GPS_Data_tmp.Longitude.degree = Longitude(1:3); % 经度 GPS_Data_tmp.Longitude.min = Longitude(4:5); tmp = str2double(Longitude(7:10)); tmp = tmp*6/1000; % tmp = tmp/10000*60; GPS_Data_tmp.Longitude.sec = num2str(floor(tmp)); GPS_Data_tmp.Longitude.millisec = num2str((tmp-floor(tmp))*10000);

合集下载

matlab 解析格式

matlab 解析格式

matlab 解析格式MATLAB 是一种高级编程语言和交互式环境，主要用于数值计算、可视化和算法开发。

在MATLAB 中，解析格式通常指的是解析数据文件的格式，如文本文件、CSV 文件、Excel 文件、JSON 文件等。

以下是MATLAB 中解析几种常见文件格式的基本方法：1.文本文件（.txt）:使用textscan或fscanf函数来解析文本文件。

matlab复制代码fid = fopen('filename.txt', 'r'); % 打开文件formatSpec = '%f %f %f'; % 格式化字符串，根据文件内容设定data = textscan(fid, formatSpec); % 解析数据fclose(fid); % 关闭文件2.逗号分隔值文件（.csv）:使用csvread、readtable或readmatrix函数来读取CSV 文件。

matlab复制代码data = csvread('filename.csv'); % 读取 CSV 文件到矩阵% 或者T = readtable('filename.csv'); % 读取 CSV 文件到表格3.Microsoft Excel 文件（.xls 或.xlsx）:使用xlsread（针对旧版.xls 文件）或readtable / readmatrix（针对新版 .xlsx 文件）函数来读取Excel 文件。

matlab复制代码data = xlsread('filename.xls'); % 读取旧版 Excel 文件% 或者T = readtable('filename.xlsx'); % 读取新版 Excel 文件到表格4.JSON 文件:MATLAB R2016b 之后的版本支持JSON 文件的读写。

可以使用jsonread和jsonwrite函数。

基于RENIX3.02格式读取卫星观测数据O文件MATLAB源代码

基于RENIX3.02格式读取卫星观测数据O文件MATLAB源代码

function [APXYZ DeltaHEN signal_type obsdata] = read_obs(filename)% %%function [signal_type obsdata] = read_obs(filename)%程序描述：读取观测数据文件，rinex格式%---------------------------------------------------------------%-- Inputs:%--filename%-- Outputs:%--obsdata %-- observation message% ---modified by Zhang 2017-12-1%---------------------------------------------------------------%-----Open the file-------%fid = fopen(filename);disp(['读取卫星观测数据文件...... : ',filename]);% [filename,pathname]=uigetfile('*.**o','读取卫星观测数据文件');% fid=fopen(strcat(pathname,filename),'rt');% if(fid==-1)% error('error to open the '+ filename);% end%% -----读取头部-----%hang = zeros(1,80);head_lines = 0;while 1head_lines = head_lines + 1;hang = fgetl(fid);head(head_lines).str = str2mat(hang);if(~isempty(findstr(hang,'APPROX POSITION XYZ'))) %obtain approximate marker position (WGS-84)，获得观测点大概位置（WGS-84）APXYZ=[0,0,0];APXYZ(1) = str2num(hang(1:14));APXYZ(2) = str2num(hang(15:28));APXYZ(3) = str2num(hang(29:42));endif(~isempty(findstr(hang,'ANTENNA: DELTA H/E/N'))) %obtain Antenna height、Eccentricities of antenna center relative to marker to the east and north (all units in meters) %天线高度：天线底部相交于标记点的高度，天线中心相对标记点东向和北向距离（单位：米）DeltaHEN=[0,0,0];DeltaHEN(1) = str2num(hang(1:14));DeltaHEN(2) = str2num(hang(15:28));DeltaHEN(3) = str2num(hang(29:42));endif(~isempty(findstr(hang,'WAVELENGTH FACT L1/2'))) %获取L1和L2载波的缺省系数，1:?整周；2:半周；0(L2载波):单频接收机L1_Wavalength = str2num(hang(1:6));L2_Wavalength = str2num(hang(7:12));endif(~isempty(findstr(hang,'# / TYPES OF OBSERV'))) %文件中观测类型数量%格式例如：21 L1 P1 C1 L2 P2 D1 D2 S1 S2# / TYPES OF OBSERV% L5 C5 D5 S5 L7 C7 D7 S7 L8# / TYPES OF OBSERV% C8 D8 S8 # / TYPES OF OBSERVtype = str2num(hang(1:6));if(type <= 9)%观测类型数量少于9个for t = 1:typesignal_type(t).types = hang(6 * t +5 : 6 * t + 6);endelse%观测类型数量多于9个，分多行显示for t = 1:9%由于第1行已经读取，先存储第1行signal_type(t).types = hang(6 * t +5 : 6 * t + 6);endfor m=1:floor(type/9)%再循环存储其他行head_lines = head_lines + 1;hang = fgetl(fid);head(head_lines).str = str2mat(hang);for t = 9*m+1:9*m+9if(t<=type)signal_type(t).types = hang(6 * (t - 9*m) +5 : 6 * (t - 9*m) + 6);endendendendend%文件中观测类型数量读取完毕，相应说明：% TYPES OF OBSERV：在本数据文件中所存储不同观测值类型的数量和观测值类型列表，如果超过9种观测值类型, 则使用续行。

利用MATLAB神经网络实现GPS高程转换设计

利用MATLAB神经网络实现GPS高程转换设计


对样 本集 Ｐ进行学 习 ， 立 映射关 系 ： 建
＝
厂 ，） （ Ｙ
（） ３
式 中 ， Ｙ为平 面坐标 ； 、 为高 程异 常 。
２ ２ Ｂ 网络 设 计 技 巧 ． Ｐ
Ｇ Ｓ大 地 高 向正 常 高 的转 换 。如 何 才 能 提 高 Ｇ Ｓ Ｐ Ｐ 大地 高 向正常 高 转换 的实 际 精 度 问题 ， 直是 Ｇ Ｓ 一 Ｐ 测绘 应用 领域 在探索 与讨 论 的问题 。
Ｐ ＝ ｛ ｌＰ … Ｐ ｝ Ｐ ，２ 式 中 ， ， ） ｉ ， … ， 。 Ｐ ＝（ Ｙ ， ， ＝１２， ｎ （） ２
式中， 日为大地高 ； ｈ为正常高 ； 为高程异常。 为充 分发 挥 Ｇ Ｓ高精 度 三 维 坐标 的优 越性 ， Ｐ 尤
其是 Ｇ Ｓ高程测 量 在工程 测 量 中的应 用 ， 须 实现 Ｐ 必
数学 者将 Ｒ Ｆ神经 网络应 用 于 Ｇ Ｓ高程 转换 ， Ｂ Ｐ 取得
了一 定 的成果 。但是 ， 何 运用 ＭＡ Ｌ Ｂ神 经 网络 如 ＴＡ 工具 箱实现 这两 种神 经 网络 的 Ｇ Ｓ高 程 转换 ， Ｐ 以及 在实 现过 程 中需注 意 的 问 题 ， 却还 没 有 文 章 作 专 门
２ １ 问题 描 述 ．
大地 高 ， 我 国所采 用 的 高 程是 相 对 于 似 大地 水 准 而
面 的正常 高 。在忽 略 垂线 偏 差 的情 况下 ， 如下 关 有
系式 ¨： Ｈ ＝ｈ＋ （） １
用 已知点 的 （ ， 和 高 程 异 常值 建立 神 经 Ｙ）
网络 的已知样 本集 ：
维普资讯
中 图 分 类 号 ： ２８ ３ Ｐ ２ ． 文 献 标 识 码 ：Ｂ 文 章 编 号 ：１７ —０４ ２ ０ ） ３０４ —４ ６ １３４ （ ０ ８ ０ —０ ８０
matlab卫星轨迹

matlab卫星轨迹

MATLAB卫星轨迹引言卫星轨迹是描述卫星在地球或其他天体之间运动的路径。

通过了解卫星轨迹，我们可以预测和控制卫星的运行，以及计划卫星的任务。

在本文档中，我们将使用MATLAB来分析和绘制卫星轨迹。

确定卫星轨迹方程卫星在空间中的运动可以通过多种方式描述，其中一种常用的方法是使用开普勒问题的解析解。

开普勒问题是描述两个质点间引力相互作用的运动方程。

对于一个质点沿着椭圆轨道运动的情况，其运动方程可以表示为：轨道方程轨道方程其中，r是卫星与中心天体（例如地球）之间的距离，a和b是椭圆的半长轴和半短轴，e是离心率，θ是卫星相对于半长轴的偏移角。

在MATLAB中，我们可以使用以下代码来计算卫星的轨道方程：function [x, y] = compute_orbit(a, e, theta)r = a * (1 - e^2) ./ (1 - e * cos(theta));x = r .* cos(theta);y = r .* sin(theta);end绘制卫星轨迹在计算了卫星的轨道方程之后，我们可以使用MATLAB的绘图工具将卫星轨迹可视化。

以下是绘制卫星轨迹的代码示例：a = 6378; % 半长轴e = 0.1; % 离心率theta = linspace(0, 2*pi, 1000); % 角度范围[x, y] = compute_orbit(a, e, theta);figure;plot(x, y);axis equal;title('卫星轨迹');xlabel('X轴');ylabel('Y轴');在上述代码中，我们假设半长轴为6378千米，离心率为0.1，并生成1000个角度点。

然后，我们使用compute_orbit 函数计算卫星的轨道坐标，并使用plot函数绘制这些坐标。

最后，我们使用axis equal命令来确保图形的横纵比例相等，以保持轨道的形状准确。