基于Matlab的TEQC绘图程序代码.docVIP

%在 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-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