GNSS电离层与对流层大气误差分析MATLAB代码

clc;clear;

A=importdata('GMGT-ZQGT.ddatm');

B=importdata('GMGT-ZQGT.ddatm343');

dzwd=zeros(length(A),1);dion=zeros(length(A),1);

count=1;count1=1;

for i=1:length(B)

    if(strcmp(B{i}(2:4),'ZWD'))

            C=B{i}(45:49);

        for j=1:length(A)

                if(strcmp(A{j}(2:4),'ZWD'))

                D=A{j}(45:49);

                    if(strcmp(C,D))

                        dzwd(count,1)=str2double(A{j}(85:91))-str2double(B{i}(85:91));

                        count=count+1;

                    end

                end

        end

    end

            if(strcmp(B{i}(2:4),'ION')&&strcmp(B{i}(12:14),'C02'))

                E=B{i}(45:49);

                if(strcmp(A{j}(2:4),'ION')&&strcmp(A{j}(12:14),'C02'))

                    F=A{j}(45:49);

                    if(strcmp(E,F))

                        dion(count1)=str2double(A{j}(85:91))-str2double(B{i}(85:91));

                        count1=count1+1;

                    end

                end

            end

end

f1=figure;

plot(dion(1:count1),'r');

title('C02的ION之差');

f2=figure;

plot(dzwd(1:count),'y');

title('ZWD之差');