计算VMF1和探空Tm的精度

clc;clear;
name1=['chan','.txt'];
[site,mjd,ah,aw,zhd,zwd,tm]=textread(name1,'%s%s%s%s%s%s%s%*[^\n]','delimiter',' ');%提取前5列数据
n=length(tm);

for i=1:n%在1-n行数据内判断数据是否为空，如果为空则跳出本次循环
    Tm(i,:)=str2num(tm{i});
    s=i/4+1;
    j=floor(s);
    if mod(i,4)==1
        Tm1(j,:)=Tm(i,:);
    elseif mod(i,4)==2
        Tm2(j,:)=Tm(i,:);
    elseif mod(i,4)==3
        Tm3(j,:)=Tm(i,:);
    elseif mod(i,4)==0
        Tm4(j-1,:)=Tm(i,:);
    end
end
GPT3_Tm=xlsread('GPT3_Tm.xlsx');
sou=xlsread('54161_2018.xlsx');
sou_tm=sou(:,1);
nn=length(sou_tm);
for i=1:nn%在1-n行数据内判断数据是否为空，如果为空则跳出本次循环
    Tm(i,:)=str2num(tm{i});
    s=i/2+1;
    j=floor(s);
    if mod(i,2)==1
        sou_Tm1(j,:)=sou_tm(i,:);
    elseif mod(i,2)==0
        sou_Tm2(j-1,:)=sou_tm(i,:);
    end
end

for i=1:365
    vmfTm(i,:)=(Tm1(i,:)+Tm2(i,:)+Tm3(i,:)+Tm4(i,:))/4;
    sou_Tm(i,:)=(sou_Tm1(i,:)+sou_Tm2(i,:))/2;
    errorvmf(i,:)=vmfTm(i,:)-GPT3_Tm(i,5);
    errorsou(i,:)=sou_Tm(i,:)-GPT3_Tm(i,5);
end
mean(errorvmf)
rms(errorvmf)
mean(errorsou)
rms(errorsou)