资源描述:
《基于Matlab的TEQC绘图程序代码》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、%在matlab下新建一个m文件,将以下代码直接拷贝进去,即可执行。%需要一个TEQC生成的plot文件作为参数functionout=teqcplot3(files);%读取TEQC生成的Plot文件,绘制数据图表,支持Copmact、Compact2>Compact3格式%选取一个TEQC的Plot文件%格式说明%*.snl载波L1的信噪比Signaltonoiseratio(S/N)%*.sn2载波L2的信噪比Signaltonoiseratio(S/N)CarrierL2%*.iod*.dl2*.d21电离层延迟观测值变化率
2、(米/秒)Derivativeofionosphericdelayobservable(m/s)%*.ion*.il2*.i21电离层延迟观测值(米)lonosphericdelayobservable(m)%*.mpl*.ml2载波LI的多路径误差MultipathCarrierLI%*.mp2*.m21载波L2的多路径误差MultipathCarrierL2%*.azi卫星方位角(°)Satelliteazimuthaldata(degrees)%*.ele卫星高度角(°)Satelliteelevationdata(degre
3、es)ifnargin==0[filenzpath]=uigetfile('*.snl;*.sn2;*.iod;*.ion;*.mpl;*.mp2;*.azi;*.ele;*.il2;*i21;*.ml2;m21;*.dl2;*d21...'请选择TEQC报告文件:’);else[path,filen,ext]二fileparts(files);path=[pathfilen={[filenext]};end%读取这个文件file=char(filen);%按行读取文件至数组A[A]=importdata([pathfile]J
4、t');%定义SAT,存放卫星数据%GPS有32颗卫星,存放序号1-32,%GLONASS有32颗卫星,存放序号33-64,%BEIDOU有35颗卫星,存放序号65-99SAT(l:length(A),l:99)=NaN;%sats(l:length(A)/l-99)=NaN;%存放采样时间,单位秒tsec(l:length(A))=NaN;%读取文件的第一行filelx=A{l};%判断是哪种格式switchfilelxcase'COMPACT'%读取数据采样I、可隔t_samp=char(A(3));%读取开始时间mjl=cha
5、r(A(4));%读取数据采样间隔T_SAMP=str2num(t_samp(max(find(t_samp=—’)):end));%读取数据采样开始时间MJL_START=str2num(mjl(max(find(mjl=='')):end));%转成时间序列数字,dateserialnumber,从0000年]月1日0时0分0秒开始计算的十进制天数MJD_START=MJL_START+678941.999999741;%i为行号n二l;i=5;case'COMPACT2'%读取数据采用间隔t_samp=char(A(2));%
6、读取开始时间mjl=char(A⑶);%读取数据釆样间隔T_SAMP=str2num(t_samp(max(find(t_samp=—’)):end));%读取数据采样开始时间MJL_START=str2num(mjl(max(find(mjl=='*)):end));%转成时间序列数字,dateserialnumber,从0000年]月1日0时0分0秒开始计算的十进制天数MJDSTART=MJLSTART+678941.999999741;n二l;i=4;case'COMPACT3*%读取开始时间t_start=char(A(2)
7、);t_start=deblank(t_start);s=splitstr(t_start,'**:6);%t_start_time=[char(s{2})'年'char(s{3})'月,char(s{4})'0'char(s{5})'时,char(s{6})分num2str(str2num(char(s{7}))/'%02d,)秒];%获取采样的开始时I、可,2013,12,7,03,05,55t_s_time=[str2num(char(s{2})),str2num(char(s{3})),str2num(char(s{4}))
8、,str2num(char(s{5})),str2num(char(s{6})),str2num(char(s{7}))];n=l;i=3;otherwisedisp(数据格式存在问题J;returnend%n=l;i=3;%sats