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

%在 matlab下新建一个m文件，以下代码直接拷贝，即可执行。 一个TEQC生成的plot文件作为参数 function out=teqcplot3(files); 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 %读取开始时间 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) 秒] ; %获取采样的开始时间，2013,12，7,03,05,55 t_s_time=[str2num(char(s{2})), str2num(char(s{3})),str2num(char(s{4})),str2num(char(s{5})),str2num(char(s{6})),str2num(char(s{7}))]; n=1;i=3; otherwise disp(数据格式存在问题); return end %n=1;i=3; %sats=str2num(A{i}); snyggfilen=strrep(filen,_,-); %生成进度条 h=waitbar(0,[正在读取数据，请