基于Matlab地TEQC绘图程序代码.docxVIP

%在 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;*.d1 2;*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-99 SAT(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.999999741; %i为行号 n=1;i=5; case COMPACT2 %读取数据采用间隔 t_samp=char(A(2)); %读取开始时间 mjl=char(A(3)); %读取数据采样间隔 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.999999741; n=1;i=4; case COMPACT3 %读取开始时间 t_start=char(A(2)); t_start = deblank(t_start); s = splitstr(t_start,**,6); %t_start_time=[char(s{2}) 年 char(s{3}) 月 char(s{4}) 日 char(s{5}) 时 char(s{6}) 分 num2str(str2num(char(s{7})),%02d) 秒] ; %获