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

MATLAB解析GPS数据程序文件

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

% 注：本程序可直接在MATLAB 2017a 中运行%该脚本文件用于学习GPS数据的读取，需要做其他用途请自行修改代码%本脚本文件的前面几行代码是要设置的一些参数%默认使用COM3（需视情况修改）%波特率设为GPS模块默认的38400%下面为程序源码clearnum_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';serial3.Terminator = Terminator;%% 打开串口count_opentimes = 1;while contains(serial3.status,'closed') > 0 && count_opentimes < num_MaxTry fopen(serial3); %打开串口count_opentimes = count_opentimes+1;endif contains(serial3.status,'open') < 1disp('open failed!');returnend%% 读取并处理数据% 初始化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); %关闭串口delete(serial3);clear serial3%%%将字符串根据'\r\n'划分成多个子字符串,同时去掉首尾无用的残余字符串function out_strs = split_str2strs(StrData)if contains(class(StrData),'char')uint8(StrData);endrecord = get_pos_enterflag(StrData);if StrData(1) == uint8('$') %开头为'$'的情况flag_start = 1;elseif size(record,2) > 0flag_start = record(1)+2;elseout_strs = cell(0,0);returnendendif StrData(end) == 13flag_end = length(StrData)-1;elseif size(record,2) > 0flag_end = record(end)-1;endendif flag_start >= flag_endout_strs = cell(0,0);returnStrData = 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 > 1out_strs{1,1} = char(StrData(1:record(1)-1));if num_strs == 2out_strs{num_strs,1} = char(StrData(record(1)+2:end));elsefor i = 2 : num_strs-1out_strs{i,1} = char(StrData(record(i-1)+2:record(i)-1));endout_strs{num_strs,1} = char(StrData(record(i)+2:end));endelseout_strs{1,1} = char(StrData);end% 得到字符串中'\r\n'在字符串中的位置(实际为'\r'的位置)function record = get_pos_enterflag(data)record = []; % 记录回车符号位置for ii = 1 : length(data)-1if data(ii) == 13if data(ii+1) == 10record = [record,ii];ii = ii+1;endendendend% 得到具体GPS结构体数据function GPS_Data_tmp = get_GPS_specificData(StrsData)GPS_Data_tmp = [];num_str = size(StrsData,1);for i = 1 : num_strstr_tab = StrsData{i,1};if contains(str_tab,'GGA') > 0GPS_Data_tmp = GNGGA(str_tab);elseif contains(str_tab,'GSA') > 0GPS_Data_tmp = GNGSA(str_tab);elseif contains(str_tab,'GSV') > 0GPS_Data_tmp = GNGSV(str_tab);elseif contains(str_tab,'RMC') > 0GPS_Data_tmp = GNRMC(str_tab);elseif contains(str_tab,'VTG') > 0GPS_Data_tmp = GNVTG(str_tab);elseif contains(str_tab,'GLL') > 0GPS_Data_tmp = GNGLL(str_tab);endendend% 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;GPS_Data_titudeDir = 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_lisec = Time(8:10);GPS_Data_titude.degree = Latitude(1:2); % 纬度GPS_Data_titude.min = Latitude(3:4);tmp = str2double(Latitude(6:9));tmp = tmp*6/1000; % tmp = tmp/10000*60;GPS_Data_titude.sec = num2str(floor(tmp));GPS_Data_lisec = 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_lisec = num2str((tmp-floor(tmp))*10000);% UTC时间转换为时间hour = GPS_Data_tmp.Time.hour;if str2num(hour)+8 >= 24GPS_Data_tmp.Time.hour = num2str(str2num(hour)+8-24);elseGPS_Data_tmp.Time.hour = num2str(str2num(hour)+8);endendfunction GPS_Data_tmp = GNGSA(str_tab)index = strfind(str_tab,',');count = 1;GPS_Data_tmp.LocationMode = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.CurState = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.PRN = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.PDOP = 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.VDOP = str_tab(index(count)+1:index(count+1)-1);count=count+1;% other = str_tab(index(count)+1:end);endfunction GPS_Data_tmp = GNGSV(str_tab)% 此语句为与卫星有关的信息（包括卫星方位，卫星编号）% 暂时用不着，不处理GPS_Data_tmp = [];endfunction GPS_Data_tmp = GNRMC(str_tab)index = strfind(str_tab,',');count = 1;Time = str_tab(index(count)+1:index(count+1)-1);count=count+1; GPS_Data_tmp.LocationState = 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;GPS_Data_titudeDir = 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.speed = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.TrueDir = str_tab(index(count)+1:index(count+1)-1);count=count+1;Date = str_tab(index(count)+1:index(count+1)-1);count=count+1; GPS_Data_tmp.MagneticAngle = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.MagneticDir = 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_lisec = Time(8:10);GPS_Data_titude.degree = Latitude(1:2); % 纬度GPS_Data_titude.min = Latitude(3:4);tmp = str2double(Latitude(6:9));tmp = tmp*6/1000; % tmp = tmp/10000*60;GPS_Data_titude.sec = num2str(floor(tmp));GPS_Data_lisec = 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_lisec = num2str((tmp-floor(tmp))*10000);GPS_Data_tmp.DATE.day = Date(1:2);GPS_Data_tmp.DATE.month = Date(3:4);GPS_Data_tmp.DATE.year = Date(5:6);% UTC时间转换为时间hour = GPS_Data_tmp.Time.hour;if str2num(hour)+8 >= 24GPS_Data_tmp.Time.hour = num2str(str2num(hour)+8-24);elseGPS_Data_tmp.Time.hour = num2str(str2num(hour)+8);endendfunction GPS_Data_tmp = GNVTG(str_tab)index = strfind(str_tab,',');count = 1;GPS_Data_tmp.TrueDir = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.ReferenceTrueDir = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.RelativeDir = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.ReferenceRelativeDir = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.step = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.stepflag = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.velocity = str_tab(index(count)+1:index(count+1)-1);count=count+1;% other = str_tab(index(count)+1:end);endfunction GPS_Data_tmp = GNGLL(str_tab)index = strfind(str_tab,',');count = 1;Latitude = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_titudeDir = 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;Date = str_tab(index(count)+1:index(count+1)-1);count=count+1;GPS_Data_tmp.LocationState = str_tab(index(count)+1:index(count+1)-1);count=count+1;% other = str_tab(index(count)+1:end);% 进一步处理GPS_Data_titude.degree = Latitude(1:2); % 纬度GPS_Data_titude.min = Latitude(3:4);tmp = str2double(Latitude(6:9));tmp = tmp*6/1000; % tmp = tmp/10000*60;GPS_Data_titude.sec = num2str(floor(tmp));GPS_Data_lisec = 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_lisec = num2str((tmp-floor(tmp))*10000); GPS_Data_tmp.DATE.day = Date(1:2);GPS_Data_tmp.DATE.month = Date(3:4);GPS_Data_tmp.DATE.year = Date(5:6);end% 更新获取到的相关数据function GPS_Data = Updata_GPU_Data(GPS_Data,GPS_Data_tmp)% 用不到的数据可以注释掉if isfield(GPS_Data_tmp,'Time') == 1GPS_Data.Time.hour = GPS_Data_tmp.Time.hour;GPS_Data.Time.min = GPS_Data_tmp.Time.min;GPS_Data.Time.sec = GPS_Data_tmp.Time.sec;GPS_lisec = GPS_Data_lisec;endif isfield(GPS_Data_tmp,'DATE') == 1GPS_Data.DATE.day = GPS_Data_tmp.DATE.day;GPS_Data.DATE.month = GPS_Data_tmp.DATE.month;GPS_Data.DATE.year = GPS_Data_tmp.DATE.year;endif isfield(GPS_Data_tmp,'Latitude') == 1GPS_titude.degree = GPS_Data_titude.degree;GPS_titude.min = GPS_Data_titude.min;GPS_titude.sec = GPS_Data_titude.sec;GPS_lisec = GPS_Data_lisec;endif isfield(GPS_Data_tmp,'LatitudeDir') == 1GPS_titudeDir = GPS_Data_titudeDir;endif isfield(GPS_Data_tmp,'Longitude') == 1GPS_Data.Longitude.degree = GPS_Data_tmp.Longitude.degree; GPS_Data.Longitude.min = GPS_Data_tmp.Longitude.min;GPS_Data.Longitude.sec = GPS_Data_tmp.Longitude.sec;GPS_lisec = GPS_Data_lisec; endif isfield(GPS_Data_tmp,'LongitudeDir') == 1GPS_Data.LongitudeDir = GPS_Data_tmp.LongitudeDir;endif isfield(GPS_Data_tmp,'GPSState') == 1GPS_Data.GPSState = GPS_Data_tmp.GPSState;endif isfield(GPS_Data_tmp,'SatelliteNum') == 1GPS_Data.SatelliteNum = GPS_Data_tmp.SatelliteNum;endif isfield(GPS_Data_tmp,'speed') == 1GPS_Data.speed = GPS_Data_tmp.speed;endif isfield(GPS_Data_tmp,'velocity') == 1GPS_Data.velocity = GPS_Data_tmp.velocity;endif isfield(GPS_Data_tmp,'LocationState') == 1GPS_Data.LocationState = GPS_Data_tmp.LocationState;endif isfield(GPS_Data_tmp,'altitude') == 1GPS_Data.altitude = GPS_Data_tmp.altitude;endif isfield(GPS_Data_tmp,'CurState') == 1GPS_Data.CurState = GPS_Data_tmp.CurState;endif isfield(GPS_Data_tmp,'LocationMode') == 1GPS_Data.LocationMode = GPS_Data_tmp.LocationMode;endif isfield(GPS_Data_tmp,'HDOP') == 1GPS_Data.HDOP = GPS_Data_tmp.HDOP;endif isfield(GPS_Data_tmp,'VDOP') == 1GPS_Data.VDOP = GPS_Data_tmp.VDOP;endif isfield(GPS_Data_tmp,'PDOP') == 1GPS_Data.PDOP = GPS_Data_tmp.PDOP;endif isfield(GPS_Data_tmp,'TrueDir') == 1GPS_Data.TrueDir = GPS_Data_tmp.TrueDir;endif isfield(GPS_Data_tmp,'MagneticAngle') == 1GPS_Data.MagneticAngle = GPS_Data_tmp.MagneticAngle;endif isfield(GPS_Data_tmp,'MagneticDir') == 1GPS_Data.MagneticDir = GPS_Data_tmp.MagneticDir;endif isfield(GPS_Data_tmp,'ReferenceTrueDir') == 1GPS_Data.ReferenceTrueDir = GPS_Data_tmp.ReferenceTrueDir;endif isfield(GPS_Data_tmp,'RelativeDir') == 1GPS_Data.RelativeDir = GPS_Data_tmp.RelativeDir;endif isfield(GPS_Data_tmp,'ReferenceRelativeDir') == 1GPS_Data.ReferenceRelativeDir = GPS_Data_tmp.ReferenceRelativeDir; endif isfield(GPS_Data_tmp,'step') == 1GPS_Data.step = GPS_Data_tmp.step;endif isfield(GPS_Data_tmp,'stepflag') == 1GPS_Data.stepflag = GPS_Data_tmp.stepflag;endif isfield(GPS_Data_tmp,'PRN') == 1GPS_Data.PRN = GPS_Data_tmp.PRN;endend% 显示相关GPS数据function show_GPS_Data(GPS_Data)DataAndTime = sprintf('20%02s-%02s-%02s %02s:%02s:%02s:%03s',...GPS_Data.DATE.year,GPS_Data.DATE.month,GPS_Data.DATE.day,...GPS_Data.Time.hour,GPS_Data.Time.min,GPS_Data.Time.sec,GPS_lise c);Location = sprintf(' %s:%02s°%02s′%03s″%04s,%s:%02s°%02s′%03s″%04s',...GPS_titudeDir,GPS_titude.degree,GPS_titude.min,GPS_Data.L atitude.sec,GPS_lisec,...GPS_Data.LongitudeDir,GPS_Data.Longitude.degree,GPS_Data.Longitude.min,GPS_D ata.Longitude.sec,GPS_lisec);% Others = sprintf('GPSState:%s,SatelliteNum:%02s,Speed:%03s,Velocity:%s,LocationState:%s',...%GPS_Data.GPSState,GPS_Data.SatelliteNum,GPS_Data.speed,GPS_Data.velocity,GPS _Data.LocationState);% show_Message_str(strcat(DataAndTime,Location,Others));show_Message_str(strcat(DataAndTime,Location));end% 初始化GPS数据结构体function GPS_Data = GPS_Init()GPS_Data.Time.hour = '0';GPS_Data.Time.min = '0';GPS_Data.Time.sec = '0';GPS_lisec = '0';GPS_Data.DATE.day = '29';GPS_Data.DATE.month = '8';GPS_Data.DATE.year = '18';GPS_titude.degree = '0'; % 纬度GPS_titude.min = '0';GPS_titude.sec = '0';GPS_lisec = '0';GPS_titudeDir = 'N';GPS_Data.Longitude.degree = '0'; % 经度GPS_Data.Longitude.min = '0';GPS_Data.Longitude.sec = '0';GPS_lisec = '0';GPS_Data.LongitudeDir = 'E';GPS_Data.GPSState = '0'; % GPS状态，0：未定位；1：无差分定位；2：带差分定位；3：无效GPS；6：正在估算GPS_Data.SatelliteNum = '0'; % 可用卫星数目GPS_Data.speed = '0';GPS_Data.velocity = '0';GPS_Data.LocationState = 'V';GPS_Data.altitude = '0'; % 海拔高度GPS_Data.CurState = '1'; % 当前状态，1：无定位信息；2：2D；3:3DGPS_Data.LocationMode = 'A'; % 定位模式，'A'：自动，'M'：手动GPS_Data.HDOP = '0.5'; % 水平精度因子GPS_Data.VDOP = '0.5'; % 垂直精度因子GPS_Data.PDOP = '0.5'; % 综合位置精度因子GPS_Data.TrueDir = '0'; % 方位角GPS_Data.MagneticAngle = '0'; % 磁偏角GPS_Data.MagneticDir = 'E'; % 磁偏角参考方向，E/W（东西经）GPS_Data.ReferenceTrueDir = 'T'; % 真实方向的参考方向，T：正北参照系 GPS_Data.RelativeDir = '0'; % 相对方向GPS_Data.ReferenceRelativeDir = 'M'; % M：磁北参照系GPS_Data.step = '0'; % 步长GPS_Data.stepflag = 'N';GPS_Data.PRN = '01'; % 正在使用的卫星PRN码编号end% 可实现静态显示输出容(防止输出容不同换行)function show_Message_str(strData)persistent CurMessage;if isempty(CurMessage)CurMessage = '';endfprintf(1,repmat('\b',1,numel(CurMessage)));fprintf(1,'%s',strData);CurMessage = strData;end。

合集下载

如何进行MATLAB地图数据处理和可视化

如何进行MATLAB地图数据处理和可视化

如何进行MATLAB地图数据处理和可视化地图数据处理和可视化在很多领域中都是很重要的一项任务，如地理信息系统（GIS）、城市规划、气象学等。

MATLAB作为一种功能强大的科学计算软件，提供了丰富的地图数据处理和可视化工具，能够帮助用户快速、准确地处理和展示地理数据。

本文将介绍如何使用MATLAB进行地图数据处理和可视化。

一、导入地图数据在MATLAB中，我们可以通过导入地理数据文件的方式来获取地图数据。

常见的地理数据文件格式包括Shapefile、GeoTIFF等。

运用MATLAB的Mapping Toolbox中的函数，我们可以轻松地读取这些文件，并将其转换成MATLAB能够处理的数据结构。

例如，如果我们有一个Shapefile文件，可以使用shaperead()函数将其读取为MATLAB中的结构体。

该结构体包含了Shapefile中所有的属性和几何信息。

我们可以通过访问这些属性和几何信息，进行进一步的数据处理和可视化操作。

二、地图数据处理与分析在导入地图数据之后，我们可以进行一系列的数据处理和分析操作。

MATLAB 提供了丰富的地理数据处理函数，如计算地理距离、邻近分析、空间插值等。

下面，我们将介绍两个常用的地图数据处理方法。

1. 计算地理距离MATLAB中使用distance()函数可以计算两点之间的地理距离。

这对于一些需要考虑地球曲率的应用非常重要。

例如，在城市规划中，我们可以使用该函数计算出两个位置之间的实际距离，从而帮助规划交通路线、确定最佳布局等。

2. 空间插值空间插值是指根据有限的采样点数据，通过数学方法推算未知点的数值。

在地图数据处理中，空间插值经常用于填充缺失值、生成连续的等值线图等。

MATLAB中的griddata()函数可以对地理数据进行插值，并生成高质量的连续地图。

三、地图数据可视化地图数据可视化是将处理后的地理数据以图形的形式展示出来。

MATLAB提供了多种可视化方式，如制作等高线图、色彩填充地图、散点图等。

GPS卫星轨道计算及其MATLAB仿真

GPS卫星轨道计算及其MATLAB仿真

GPS卫星轨道计算及其MATLAB仿真黎奇，白征东，李帅，陈波波（清华大学地球空间信息研究所，北京 100084）一、程序设计思路1. 读取RINEX文件（注意：文件路径）2. 计算测量日周积秒（测量日的格里历→GPST）3. 按卫星轨道计算步骤计算WGS-84坐标系坐标（内插）4. 按需要将WGS-84坐标系下坐标转换为所需坐标系坐标5. 画图输出二、n 文件说明及读取程序参考时刻oe t 的RINEX 格式的 “”广播星历文件具体如下：（加粗部分为本次轨道化Ω，率i ，弧度/秒4-22）标svacc ，米）收到的卫星信号解，秒）文件名：RinexNreader.m 输 入：文件地址，卫星编号三、计算测量日的周积秒文件名：GCtoGPS.m （其中调用函数：GCtoJD.m）输入：指定公历的年、月、日、时、分、秒文件名：GCtoJD.m输 入：指定公历的 年、月、日四、GPS 卫星轨道计算步骤及计算程序1. 计算卫星运动的平均角速度n平均角速度()03n =经摄动参数n ∆改正后的平均角速度0n n n =+∆3#61-79），n ∆（2#42-60）；14323.98600510/GM m s =⨯ 2. 计算归化时间k t说明：①广播星历是oe t 时刻的，对应的轨道参数也是oe t 时刻的，而观测时间在t 时刻，显然oe t t ＜。

所以，要想获得t 时刻的轨道参数，需要知道t 与oe t 之间的差值即k t 。

以此，按照oe t 时刻轨道参数，外推t 时刻轨道参数。

②k t 的起算时间是星期六/星期日子夜0点，当302400k t s ＞时，604800k t s -；当302400k t s ＜-时，+604800k t s 。

（604800s=1周） =k oe t t t -，且604800302400604800302400k k k k k k t t t t t t =-⎧⎨=+⎩＞ ＜-已知：oe t （1#4-21）3. 计算观测时刻的平近点角k M0k k M M nt =+已知：0M （2#61-79），n （见1），k t （见2） 4. 计算观测时刻的偏近点角k Esin k k k E M e E =+已知：k M （见3），e （3#23-41）方法：迭代解算，设初值0k k E M =，迭代2次基本收敛。

matlab出租车程序,MATLAB用出租车GPS数据分析深圳道路交通情况源代码

matlab出租车程序,MATLAB用出租车GPS数据分析深圳道路交通情况源代码

matlab出租车程序,MATLAB⽤出租车GPS数据分析深圳道路交通情况源代码%统计满⾜条件的车辆数%tongji.m%求车某时间段车得位置，即位于哪个分区%7点到9点%weizhi79.mformat shortclearclc[filename pathname index]=uigetfile('.txt','请选择结果⽂件','','MultiSelect','on');if indexif iscell(filename)len=length(filename);res=[];col_title=cell(1,len);qishi1=[];l=0;for i=1:lenqishi=[];jieshu=[];fid=fopen([pathname filename{i}]);col_title{i}=strtok(filename{i},'.txt');[val pos]=textscan(fid,'%s%c%c%c%c%c%c%c%c%c%c%c%c%c%c%c%c%c%c%f%f%f%f%f','headerLines',1,'Delimiter',','); if length(find(cell2mat(val(:,22))>2))==0l=l+1;riqi11=cell2mat(val(:,11));riqi1=riqi11-49+1; %⽇期riqi12=cell2mat(val(:,12));riqi2=riqi12-49+1;riqi13=cell2mat(val(:,13));riqi3=riqi13-49+1;hour=10*riqi2+riqi3; %时riqi15=cell2mat(val(:,15));riqi4=riqi15-49+1;riqi16=cell2mat(val(:,16));riqi5=riqi16-49+1;fen=riqi4*10+riqi5; %分n=length(hour);aaa=cell2mat(val(:,20:21)); %第⼀列为经度，第⼆列为纬度% a1=cell2mat(val(:,1));for j=1:n-1if riqi1(j)==3 %对应改riqi1，星期六if hour(j)==7&fen(j)>=0&fen(j)<=10 %对应改上下班% chepai1=[chepai1;a1];qishi=[qishi;aaa(j,:)];endendendendif length(qishi)~=0qishi1=[qishi1;qishi(1,:)];endfclose(fid);iendendend%距离%起始[n1,n2]=size(qishi1);for i=1:n1a1=sqrt((qishi1(i,1)-113.9333)^2+(qishi1(i,2)-22.5338)^2); %与分区⼀的距离a2=sqrt((qishi1(i,1)-114.1164)^2+(qishi1(i,2)-22.5649)^2);a3=sqrt((qishi1(i,1)-114.3222)^2+(qishi1(i,2)-22.6343)^2);a4=sqrt((qishi1(i,1)-114.0418)^2+(qishi1(i,2)-22.6431)^2);a5=sqrt((qishi1(i,1)-114.2457)^2+(qishi1(i,2)-22.7228)^2);a6=sqrt((qishi1(i,1)-114.0522)^2+(qishi1(i,2)-22.5398)^2);a7=sqrt((qishi1(i,1)-113.9091)^2+(qishi1(i,2)-22.7242)^2);a8=sqrt((qishi1(i,1)-113.8698)^2+(qishi1(i,2)-22.5907)^2);aa=[a1 a2 a3 a4 a5 a6 a7 a8];[v,d1]=min(aa);qishi1(i,n2+1)=d1;end%结束bb=qishi1(find(qishi1(:,n2+1)==6),:); %六区车辆信息[n3,n4]=size(bb);R1=0.0234;R2=0.0378;for i=1:n3bbb=sqrt((bb(i,1)-114.0522)^2+(bb(i,2)-22.5398)^2); %距中⼼点的距离if bbb<=R1bb(i,n4+1)=1; %1为最内圆elseif bbb<=R2bb(i,n4+1)=2;elsebb(i,n4+1)=3;endendbbb1=bb(find(bb(:,n4+1)==1),:);。

MATLAB技术GPS数据处理教程

MATLAB技术GPS数据处理教程

MATLAB技术GPS数据处理教程导言GPS（全球定位系统）以其高精度的定位功能在各个领域得到广泛应用。

在我们的日常生活中，GPS定位已经成为了司空见惯的技术，无论是导航系统、运输管理还是环境监测等。

本文旨在介绍使用MATLAB技术处理GPS数据的方法和技巧，帮助读者更好地利用现有的GPS数据进行分析和应用。

一、GPS数据处理的基本概念在开始介绍MATLAB技术处理GPS数据之前，我们首先需要了解一些GPS数据处理的基本概念。

GPS接收器通过接收卫星发射的信号来计算接收器与卫星之间的距离，进而确定接收器的位置。

GPS数据通常包括接收器的经纬度坐标、时间戳和位置精度等信息。

二、数据准备与读取要处理GPS数据，首先需要准备好相应的数据文件。

通常，GPS数据会以文本文件（如CSV或TXT）的形式进行存储。

在MATLAB中，可以使用load或readtable函数来读取存储在文本文件中的GPS数据，并将其转换为MATLAB中的矩阵或表格形式，方便进行后续的处理。

三、数据预处理与清洗在进行GPS数据处理之前，常常需要对数据进行预处理和清洗，以去除异常值和噪音，确保数据的准确性。

MATLAB提供了一系列的函数和工具箱，例如数据滤波技术、异常值检测算法等，可以用于数据的预处理和清洗。

此外，还可以利用MATLAB内置的绘图工具对数据进行可视化，快速发现异常数据。

四、数据分析与可视化数据处理的核心在于数据分析和可视化。

MATLAB提供了丰富的函数和工具箱，可以对GPS数据进行多维度的分析和可视化。

例如，可以利用MATLAB的统计工具箱对数据进行统计分析，找出数据中的趋势和规律。

同时，也可以使用MATLAB的绘图函数绘制位置轨迹图、速度图以及加速度图等，直观地展示数据的空间和时序特征。

五、轨迹重建与预测GPS数据处理的另一个重要应用是轨迹重建和预测。

基于历史GPS数据，我们可以通过建立合适的模型，使用MATLAB的预测函数对未来的轨迹进行预测。