导航技术综合实验报告..doc

导航技术综合实验报告 实验号：实验九、实验十、实验十一、实验十二 2012年4月23日星期一 一、实验目的 理论：基本掌握解码星历表解码星历表解码星历表从帧子帧2和子帧3中得到星历数据转换成十进制格式function [eph, TOW] = ephemeris(bits, D30Star) 1.实验要求： 实验代码见实验9中的ephemeris.m文件，理解实现星历解算的代码，并把函数代码和流程图中的各个环节对应起来。 获得输入的导航数据信息 function [eph, TOW] = ephemeris(bits, D30Star) 利用function形式，调用子函数，输入参数。 检查导航数据的bits长度是否满足导航数据5个子帧的长度要求 if length(bits) 1500 error(The parameter BITS must contain 1500 bits!); end 检查参数bits和D30Star是否是字符串 if ~ischar(bits) error(The parameter BITS must be a character array!); end if ~ischar(D30Star) error(The parameter D30Star must be a char!); end 循环操作，一一对5个子帧进行解码。 for i = 1:5 …end即外层循环。完成 4.1提取第一个子帧的位； subframe = bits(300*(i-1)+1 : 300*i) 4.2检查所提供的30个位的奇偶 for j = 1:10 [subframe(30*(j-1)+1 : 30*j)] = ... checkPhase(subframe(30*(j-1)+1 : 30*j), D30Star); D30Star = subframe(30*j); End 4.3用subframeID = bin2dec(subframe(50:52))提取子帧号 4.4基于子帧的ID号来进行解码，将位信息转换为十进制数。 switch subframeID case 1 eph.weekNumber = bin2dec(subframe(61:70)) + 1024; eph.accuracy = bin2dec(subframe(73:76));%用户测距精度 eph.health = bin2dec(subframe(77:82));%卫星的健康状况 eph.T_GD = twosComp2dec(subframe(197:204)) * 2^(-31);%群的延迟差别评估 该数据为二进制补码，twosComp2dec将补码转换为对应的十进制数 eph.IODC = bin2dec([subframe(83:84) subframe(197:204)]);%数据、时钟的发布号 eph.t_oc = bin2dec(subframe(219:234)) * 2^4;%卫星时钟修正参数 eph.a_f2 = twosComp2dec(subframe(241:248)) * 2^(-55);%卫星时钟修正参数 eph.a_f1 = twosComp2dec(subframe(249:264)) * 2^(-43);%卫星时钟修正参数 eph.a_f0 = twosComp2dec(subframe(271:292)) * 2^(-31);%卫星时钟修正参数 case 2 eph.IODE_sf2 = bin2dec(subframe(61:68));%数据与星历发布号 eph.C_rs = twosComp2dec(subframe(69: 84)) * 2^(-5);%轨道半径的正弦调和修正项的幅度 eph.deltan = ... twosComp2dec(subf