卫星导航定位算法与程序设计——实验报告
- 格式：doc
- 大小：164.00 KB
- 文档页数：11
fprintf(fid,'纬度幅角的正弦调和项改正的振幅：%d\n轨道半径的余弦调和项改正的振幅：%d\n轨道半径的正弦调和项改正的振幅：%d\n轨道倾角：%d\n轨道倾角变化率：%d\n星历参考时刻：%d\n',cus,crc,crs,i0,idot,toe)
fclose(fid);
、输出结果
2、针对第二部分内容：
fork = 1:NoObs
Obs(u,k) = str2num(lin(2+16*(k-1):16*k-2));
end
end
else% This will typical be Z12 data
Obs = Obs(:,[1 2 3 4 5]);% We cancel the last two columns 6 and 7
GpsSendPosition SdSignPoTemp;
intnTheFitPoint=0;
Position pTemp;
();
if()==0)
{
strErr = ("PRN=%d没有对应星历");
returnfalse;
}
if()==0)
{
returnfalse;
}
cout<<fixed;
for(inti=0;i<();i++)
OC,nGTOC);
ClkBias+
nData[nTheFitPoint].dClkDrift*nData[nTheFitPoint].dClkDriftRate*
pow(ObserveData[j].PRN;
=;
=;
=;
=oData[i].oObserveData[j].dC1;
if(i<10)
cout<<i+1<<" "<<<<" "<<<<" "<<<<" "<<
answer = findstr(line,'ANTENNA: DELTA H/E/N');
if~isempty(answer)
fork = 1:3
[delta, line] = strtok(line);
del = str2num(delta);
ant_delta = [ant_delta del];
<<" "<<" "<<//
(SdSignPoTemp);
}
}
returntrue;
}
3、实验结果
四
计算卫星位置的这部分“C++”代码比较综合，读代码的难度较大，它综合了前面的“C++”O文件读取，N文件读取，然后后面的计算代码较为复杂，需要对卫星轨道计算的公式了解全面，因此对这部分实验，是在老师和同学的帮助下完成的，直接在轨道计算的函数下加入了输出函数，将计算的轨道数据直接输出出来。
一、实验目的
1、进一步熟悉N文件的读入
2、掌握开普勒参数计算卫星轨道的过程
3、编程实现采用广播星历计算卫星轨道
4、掌握MATLAB函数调用步骤
二、
本实验内容包括：
1、调试时间转换函数，熟悉内容，备主函数调用
2、调试广播星历导航文件的读入程序，备主函数调用
3、根据卫星位置计算公式编写主函数，同时调用时间转换、星历读取等的子函数来共同完成卫星位置的计算，最后输出计算结果
Obs_types = [Obs_types ot];
end;
ifound_types = 1;
end;
end;
%fclose(fid);
“grabdata”函数：
functionObs = grabdata(fid, NoSv, NoObs)
%GRABDATA Positioned in a RINEX file at a selected epoch
end;
end
answer = findstr(line,'# / TYPES OF OBSERV');
if~isempty(answer)
[NObs, line] = strtok(line);
NoObs = str2num(NObs);
fork = 1:NoObs
[ot, line] = strtok(line);
、主要函数：
“anheader”函数：
function[Obs_types, ant_delta,ifound_types,approx_XYZ] = anheader(file)
fid = fopen(file,'rt');
eof = 0;
ifound_types = 0;
Obs_types = [];
Obs(u,k) = str2num(lin(1+16*(k-1):16*k-2));
end
% Obs(u,NoObs) = str2num(lin(65:78));
end
end
end
实验结果
四、实验感想
这部分实验是我在之前做的，之前自己有看过“gps_easy”有关的代码，看过相关的“N文件”“O文件”读写函数，并且学会了如何调用这些函数，对里面的输出量有了一点的了解，所以我自己的主要工作就是运用了“fprintf”函数，将读取到matlab中的矩阵写入TXT文档中，这部分工作量不是很大，但较有意义。
{
double*VX=newdouble[4];j++)
{
oGTime);
tr=oGTime;
oData[i].oObserveData[j].PRN,nTheFitPoint);
if(nTheFitPoint==-1)
{
"PRN=%d没?有?D对?应?|星?历¤¨2",oData[i].oObserveData[j].PRN);
通过这部分实验是我加深了对卫星轨道计算公式的了解，以及以后面对复杂的公式应该如何应对。。
2、主要函数
boolCGps::GetGpsPosition()
{
GPSTIME ts;
GPSTIME tr;
GPSTIME ts0;
GPSTIME oGTime;
GPSTIME nGTOC;
vector<GpsSendPosition>SendSignPosition;
GpsPos GpsPTemp;
end;
end
answer = findstr(line,'APPROX POSITION XYZ');
if~isempty(answer)
fork = 1:3
[app_XYZ, line] = strtok(line);
del = str2num(app_XYZ);
approx_XYZ = [approx_XYZ del];
具体的实现函数：
“TIME2JUL”函数：
doubleCTimeTr::TIME2JUL() XT”文件即可。
、主要函数代码:
“get_eph”函数：
functioneph = get_eph(ephemeridesfile)
fide = fopen(ephemeridesfile);
[eph, count] = fread(fide, Inf,'double');
NoObs = 5;
foru = 1:NoSv
lin = fgetl(fid);
lin_doppler = fgetl(fid);
fork = 1:NoObs%%-1
ifisempty(str2num(lin(1+16*(k-1):16*k-2))) == 1, Obs(u,k) = nan;
else%
4、理清程序各模块的功能结构
三
1、实验思路：
在老师提供的“SPP”文件中，直接利用卫星位置计算函数，进行卫星位置的计算，将利用“”文件中的”GetGpsPosition”函数，利用其中的迭代求解卫星位置部分，用“cout”直接将卫星迭代后的位置直接输出，因为星历文件中有较多的星历，所以利用循环语句，将求解出来的卫星位置依次输出出来。
ant_delta = [];
approx_XYZ = [];
while1% Gobbling the header
line = fgetl(fid);
answer = findstr(line,'END OF HEADER');
if~isempty(answer),break;end;
if(line == -1), eof = 1;break;end;
2013 级测绘工程专业
卫星导航定位算法与程序设计
实
验
报
告
实验名称：卫星导航基本程序设计
班 级：
学 号：
姓 名：
实验时间：2016年6月28日~2016年6月30
中 国 矿 业 大 学
实验一 时空基准转换
一、实验目的
1、加深对时空系统及其之间转换关系的理解
2、掌握常用时空基准之间的转换模型与软件实现
noeph = count/22;
eph = reshape(eph, 22, noeph);
“rinexe”函数（部分）：
functionrinexe(ephemerisfile, outputfile)
fide = fopen(ephemerisfile);
head_lines = 0;
fclose(fid);
、输出结果
2、针对第二部分内容：
fork = 1:NoObs
Obs(u,k) = str2num(lin(2+16*(k-1):16*k-2));
end
end
else% This will typical be Z12 data
Obs = Obs(:,[1 2 3 4 5]);% We cancel the last two columns 6 and 7
GpsSendPosition SdSignPoTemp;
intnTheFitPoint=0;
Position pTemp;
();
if()==0)
{
strErr = ("PRN=%d没有对应星历");
returnfalse;
}
if()==0)
{
returnfalse;
}
cout<<fixed;
for(inti=0;i<();i++)
OC,nGTOC);
ClkBias+
nData[nTheFitPoint].dClkDrift*nData[nTheFitPoint].dClkDriftRate*
pow(ObserveData[j].PRN;
=;
=;
=;
=oData[i].oObserveData[j].dC1;
if(i<10)
cout<<i+1<<" "<<<<" "<<<<" "<<<<" "<<
answer = findstr(line,'ANTENNA: DELTA H/E/N');
if~isempty(answer)
fork = 1:3
[delta, line] = strtok(line);
del = str2num(delta);
ant_delta = [ant_delta del];
<<" "<<" "<<//
(SdSignPoTemp);
}
}
returntrue;
}
3、实验结果
四
计算卫星位置的这部分“C++”代码比较综合，读代码的难度较大，它综合了前面的“C++”O文件读取，N文件读取，然后后面的计算代码较为复杂，需要对卫星轨道计算的公式了解全面，因此对这部分实验，是在老师和同学的帮助下完成的，直接在轨道计算的函数下加入了输出函数，将计算的轨道数据直接输出出来。
一、实验目的
1、进一步熟悉N文件的读入
2、掌握开普勒参数计算卫星轨道的过程
3、编程实现采用广播星历计算卫星轨道
4、掌握MATLAB函数调用步骤
二、
本实验内容包括：
1、调试时间转换函数，熟悉内容，备主函数调用
2、调试广播星历导航文件的读入程序，备主函数调用
3、根据卫星位置计算公式编写主函数，同时调用时间转换、星历读取等的子函数来共同完成卫星位置的计算，最后输出计算结果
Obs_types = [Obs_types ot];
end;
ifound_types = 1;
end;
end;
%fclose(fid);
“grabdata”函数：
functionObs = grabdata(fid, NoSv, NoObs)
%GRABDATA Positioned in a RINEX file at a selected epoch
end;
end
answer = findstr(line,'# / TYPES OF OBSERV');
if~isempty(answer)
[NObs, line] = strtok(line);
NoObs = str2num(NObs);
fork = 1:NoObs
[ot, line] = strtok(line);
、主要函数：
“anheader”函数：
function[Obs_types, ant_delta,ifound_types,approx_XYZ] = anheader(file)
fid = fopen(file,'rt');
eof = 0;
ifound_types = 0;
Obs_types = [];
Obs(u,k) = str2num(lin(1+16*(k-1):16*k-2));
end
% Obs(u,NoObs) = str2num(lin(65:78));
end
end
end
实验结果
四、实验感想
这部分实验是我在之前做的，之前自己有看过“gps_easy”有关的代码，看过相关的“N文件”“O文件”读写函数，并且学会了如何调用这些函数，对里面的输出量有了一点的了解，所以我自己的主要工作就是运用了“fprintf”函数，将读取到matlab中的矩阵写入TXT文档中，这部分工作量不是很大，但较有意义。
{
double*VX=newdouble[4];j++)
{
oGTime);
tr=oGTime;
oData[i].oObserveData[j].PRN,nTheFitPoint);
if(nTheFitPoint==-1)
{
"PRN=%d没?有?D对?应?|星?历¤¨2",oData[i].oObserveData[j].PRN);
通过这部分实验是我加深了对卫星轨道计算公式的了解，以及以后面对复杂的公式应该如何应对。。
2、主要函数
boolCGps::GetGpsPosition()
{
GPSTIME ts;
GPSTIME tr;
GPSTIME ts0;
GPSTIME oGTime;
GPSTIME nGTOC;
vector<GpsSendPosition>SendSignPosition;
GpsPos GpsPTemp;
end;
end
answer = findstr(line,'APPROX POSITION XYZ');
if~isempty(answer)
fork = 1:3
[app_XYZ, line] = strtok(line);
del = str2num(app_XYZ);
approx_XYZ = [approx_XYZ del];
具体的实现函数：
“TIME2JUL”函数：
doubleCTimeTr::TIME2JUL() XT”文件即可。
、主要函数代码:
“get_eph”函数：
functioneph = get_eph(ephemeridesfile)
fide = fopen(ephemeridesfile);
[eph, count] = fread(fide, Inf,'double');
NoObs = 5;
foru = 1:NoSv
lin = fgetl(fid);
lin_doppler = fgetl(fid);
fork = 1:NoObs%%-1
ifisempty(str2num(lin(1+16*(k-1):16*k-2))) == 1, Obs(u,k) = nan;
else%
4、理清程序各模块的功能结构
三
1、实验思路：
在老师提供的“SPP”文件中，直接利用卫星位置计算函数，进行卫星位置的计算，将利用“”文件中的”GetGpsPosition”函数，利用其中的迭代求解卫星位置部分，用“cout”直接将卫星迭代后的位置直接输出，因为星历文件中有较多的星历，所以利用循环语句，将求解出来的卫星位置依次输出出来。
ant_delta = [];
approx_XYZ = [];
while1% Gobbling the header
line = fgetl(fid);
answer = findstr(line,'END OF HEADER');
if~isempty(answer),break;end;
if(line == -1), eof = 1;break;end;
2013 级测绘工程专业
卫星导航定位算法与程序设计
实
验
报
告
实验名称：卫星导航基本程序设计
班 级：
学 号：
姓 名：
实验时间：2016年6月28日~2016年6月30
中 国 矿 业 大 学
实验一 时空基准转换
一、实验目的
1、加深对时空系统及其之间转换关系的理解
2、掌握常用时空基准之间的转换模型与软件实现
noeph = count/22;
eph = reshape(eph, 22, noeph);
“rinexe”函数（部分）：
functionrinexe(ephemerisfile, outputfile)
fide = fopen(ephemerisfile);
head_lines = 0;
