%在 matlab 下新建一个 m 文件,将以下代码直接拷贝进去,即可执行。％需要一个 TEQC 生成的 plot 文件作为参数function out=teqcplot3(files）;％读取 TEQC 生成的 Plot 文件，绘制数据图表，支持 Copmact、Compact2、Compact3 格式%选取一个 TEQC 的 Plot 文件％ 格式说明％ ＊。sn1 载波 L1 的信噪比 Signal to noise ratio （S/N） ％ *。sn2 载波 L2 的信噪比 Signal to noise ratio （S/N） Carrier L2％ ＊。iod ＊。d12 ＊。d21 电离层延迟观测值变化率(米/秒） Derivative of ionospheric delay observable （m/s）％ ＊。ion ＊。i12 *.i21 电离层延迟观测值(米）Ionospheric delay observable (m）％ *。mp1 *.m12 载波 L1 的多路径误差 Multipath Carrier L1% ＊。mp2 ＊。m21 载波 L2 的多路径误差 Multipath Carrier L2％ ＊.azi 卫星方位角(°） Satellite azimuthal data (degrees）% ＊.ele 卫星高度角（°) Satellite elevation data （degrees）if nargin==0 [filen,path ］=uigetfile（’＊。 sn1；＊ .sn2；＊ .iod；＊。 ion；*.mp1; ＊。 mp2 ；＊.azi;*.ele;*。i12；＊i21;＊.m12；m21；＊。d12；＊d21’,.。. ’请选择 TEQC 报告文件：')；else ［path，filen，ext]=fileparts（files)； path=[path ’\’］; filen=｛［filen ext］};end%读取这个文件file=char（filen）；%按行读取文件至数组 A［A]=importdata([path file］,’\t’)；%定义 SAT，存放卫星数据%GPS 有 32 颗卫星，存放序号 1-32，%GLONASS 有 32 颗卫星，存放序号 33—64,％BEIDOU 有 35 颗卫星，存放序号 65—99SAT（1:length（A），1:99)=NaN;％sats（1:length（A）,1：99)=NaN；％存放采样时间，单位秒tsec（1：length(A))=NaN；％读取文件的第一行filelx=A{1}；%推断是哪种格式switch filelx case ’COMPACT' %读取数据采样间隔 t_samp=char(A（3)）; ％读取开始时间 mjl=char（A（4)）； %读取数据采样间隔 T_SAMP=str2num（t_samp（max（find(t_samp==’ '））：end）)； ％读取数据采样开始时间 MJL_START=str2num（mjl(max（find(mjl==’ ’)）：end）); ％转成时间序列数字,date serial number，从 0000 年 1 月 1 日 0 时 0 分 0 秒开始计算的十进制天数 MJD_START=MJL_START+678941。9...
