%% Чтение RINEX наблюдений из файла в массивы
zimj=readrnx('zimj1500.07o'); % Для станции zimj
onsa=readrnx('onsa1500.07o'); % Для станции onsa
wtzr=readrnx('wtzr1500.07o'); % Для станции wtzr
zimm=readrnx('zimm1500.07o'); % Для станции zimm
%% Доступность спутников для станции zimj C1
imagesc(zimj(:,:,1))
%% Доступность спутников для станции zimj P2
imagesc(zimj(:,:,3))
%% Доступность спутников для станции onsa С1
imagesc(onsa(:,:,1))
%% Доступность спутников для станции onsa P2
imagesc(onsa(:,:,3))
%% Доступность спутников для станции wtzr С1
imagesc(wtzr(:,:,1))
%% Доступность спутников для станции wtzr P2
imagesc(wtzr(:,:,3))
%% Доступность спутников для станции zimj S1
imagesc(zimj(:,:,8))
%% Доступность спутников для станции zimj S2
imagesc(zimj(:,:,9))
%% Доступность спутников для станции onsa S1
imagesc(onsa(:,:,8))
%% Доступность спутников для станции onsa S2
imagesc(onsa(:,:,9))
%% Доступность спутников для станции wtzr S1
imagesc(wtzr(:,:,8))
%% Доступность спутников для станции wtzr S2
imagesc(wtzr(:,:,9))
%%
t=1:2880; % list of epoch index
%% Показания P1 для станции zimj для спутников 1, 11, 21 и 31
plot(t,[zimj(t,1,2),zimj(t,11,2),zimj(t,21,2),zimj(t,31,2)])
%% Показания P2 для станции zimj для спутников 1, 11, 21 и 31
plot(t,[zimj(t,1,3),zimj(t,11,3),zimj(t,21,3),zimj(t,31,3)])
%% Показания L1 для станции zimj для спутников 1, 11, 21 и 31
plot(t,[zimj(t,1,4),zimj(t,11,4),zimj(t,21,4),zimj(t,31,4)])
%% Показания L2 для станции zimj для спутников 1, 11, 21 и 31
plot(t,[zimj(t,1,5),zimj(t,11,5),zimj(t,21,5),zimj(t,31,5)])
%% Показания P1 для станции onsa для спутников 1, 11, 21 и 31
plot(t,[onsa(t,1,2),onsa(t,11,2),onsa(t,21,2),onsa(t,31,2)])
%% Показания P2 для станции onsa для спутников 1, 11, 21 и 31
plot(t,[onsa(t,1,3),onsa(t,11,3),onsa(t,21,3),onsa(t,31,3)])
%% Показания L1 для станции onsa для спутников 1, 11, 21 и 31
plot(t,[onsa(t,1,4),onsa(t,11,4),onsa(t,21,4),onsa(t,31,4)])
%% Показания L2 для станции onsa для спутников 1, 11, 21 и 31
plot(t,[onsa(t,1,5),onsa(t,11,5),onsa(t,21,5),onsa(t,31,5)])
%% Показания P1 для станции wtzr для спутников 1, 11, 21 и 31
plot(t,[wtzr(t,1,2),wtzr(t,11,2),wtzr(t,21,2),wtzr(t,31,2)])
%% Показания P2 для станции wtzr для спутников 1, 11, 21 и 31
plot(t,[wtzr(t,1,3),wtzr(t,11,3),wtzr(t,21,3),wtzr(t,31,3)])
%% Показания L1 для станции wtzr для спутников 1, 11, 21 и 31
plot(t,[wtzr(t,1,4),wtzr(t,11,4),wtzr(t,21,4),wtzr(t,31,4)])
%% Показания L2 для станции wtzr для спутников 1, 11, 21 и 31
plot(t,[wtzr(t,1,5),wtzr(t,11,5),wtzr(t,21,5),wtzr(t,31,5)])
%% Показания D1 для станции zimj для спутников 1, 11, 21 и 31
plot(t,[zimj(t,1,6),zimj(t,11,6),zimj(t,21,6),zimj(t,31,6)])
%% Показания D2 для станции zimj для спутников 1, 11, 21 и 31
plot(t,[zimj(t,1,7),zimj(t,11,7),zimj(t,21,7),zimj(t,31,7)])
%% Показания D1 для станции onsa для спутников 1, 11, 21 и 31
plot(t,[onsa(t,1,6),onsa(t,11,6),onsa(t,21,6),onsa(t,31,6)])
%% Показания D2 для станции onsa для спутников 1, 11, 21 и 31
plot(t,[onsa(t,1,7),onsa(t,11,7),onsa(t,21,7),onsa(t,31,7)])
%% Показания D1 для станции wtzr для спутников 1, 11, 21 и 31
plot(t,[wtzr(t,1,6),wtzr(t,11,6),wtzr(t,21,6),wtzr(t,31,6)])
%% Показания D2 для станции wtzr для спутников 1, 11, 21 и 31
plot(t,[wtzr(t,1,7),wtzr(t,11,7),wtzr(t,21,7),wtzr(t,31,7)])
%% Показания S1 для станции zimj для спутников 1, 11, 21 и 31
plot(t,[zimj(t,1,8),zimj(t,11,8),zimj(t,21,8),zimj(t,31,8)])
%% Показания S2 для станции zimj для спутников 1, 11, 21 и 31
plot(t,[zimj(t,1,9),zimj(t,11,9),zimj(t,21,9),zimj(t,31,9)])
%% Показания S1 для станции onsa для спутников 1, 11, 21 и 31
plot(t,[onsa(t,1,8),onsa(t,11,8),onsa(t,21,8),onsa(t,31,8)])
%% Показания S2 для станции onsa для спутников 1, 11, 21 и 31
plot(t,[onsa(t,1,9),onsa(t,11,9),onsa(t,21,9),onsa(t,31,9)])
%% Показания S1 для станции wtzr для спутников 1, 11, 21 и 31
plot(t,[wtzr(t,1,8),wtzr(t,11,8),wtzr(t,21,8),wtzr(t,31,8)])
%% Показания S2 для станции wtzr для спутников 1, 11, 21 и 31
plot(t,[wtzr(t,1,9),wtzr(t,11,9),wtzr(t,21,9),wtzr(t,31,9)])
%% Зависимость S1 от P1 для станции zimj для спутника 1
plot(zimj(t,1,2),zimj(t,1,8))
%% Зависимость S2 от P2 для станции zimj для спутника 1
plot(zimj(t,1,3),zimj(t,1,9))
%% Зависимость S1 от P1 для станции onsa для спутника 1
plot(onsa(t,1,2),onsa(t,1,8))
%% Зависимость S2 от P2 для станции onsa для спутника 1
plot(onsa(t,1,3),onsa(t,1,9))
%% Зависимость S1 от P1 для станции wtzr для спутника 1
plot(wtzr(t,1,2),wtzr(t,1,8))
%% Зависимость S2 от P2 для станции wtzr для спутника 1
plot(wtzr(t,1,3),wtzr(t,1,9))
%% ЛИНЕЙНЫЕ КОМБИНАЦИИ
%% Ionosphere Linear Combination
f0=10.23*10^6; % Фундаментальная частота
f1=154*f0;
f2=120*f0;
c=299792458;
k1=(f1^2)/(f1^2-f2^2);
k2=-(f2^2)/(f1^2-f2^2);
WaveLenght1=c/f1;
WaveLenght2=c/f2;
%% Независимая от ионосферы линейная комбинация для псевдорасстояний (PI) для станции zimj для спутника 1
plot(t,k1*zimj(t,1,2)+k2*zimj(t,1,3))
%% Независимая от ионосферы линейная комбинация для фазовых измерений (LI) для станции zimj для спутника 1
plot(t,k1*zimj(t,1,4)*WaveLenght1+k2*zimj(t,1,5)*WaveLenght2)
%% Независимая от ионосферы линейная комбинация для псевдорасстояний (PI) для станции onsa для спутника 1
plot(t,k1*onsa(t,1,2)+k2*onsa(t,1,3))
%% Независимая от ионосферы линейная комбинация для фазовых измерений (LI) для станции onsa для спутника 1
plot(t,k1*onsa(t,1,4)*WaveLenght1+k2*onsa(t,1,5)*WaveLenght2)
%% Независимая от ионосферы линейная комбинация для псевдорасстояний (PI) для станции wtzr для спутника 1
plot(t,k1*wtzr(t,1,2)+k2*wtzr(t,1,3))
%% Независимая от ионосферы линейная комбинация для фазовых измерений (LI) для станции wtzr для спутника 1
plot(t,k1*wtzr(t,1,4)*WaveLenght1+k2*wtzr(t,1,5)*WaveLenght2)
%% Независимая от ионосферы линейная комбинация для псевдорасстояний (PI) для станции zimm для спутника 1
plot(t,k1*zimm(t,1,2)+k2*zimm(t,1,3))
%% Независимая от ионосферы линейная комбинация для фазовых измерений (LI) для станции zimm для спутника 1
plot(t,k1*zimm(t,1,4)*WaveLenght1+k2*zimm(t,1,5)*WaveLenght2)
%% Независимая от ионосферы линейная комбинация для псевдорасстояний (PI) для станции zimm для спутника 11
plot(t,k1*zimm(t,11,2)+k2*zimm(t,11,3))
%% Независимая от ионосферы линейная комбинация для фазовых измерений (LI) для станции zimm для спутника 11
plot(t,k1*zimm(t,11,4)*WaveLenght1+k2*zimm(t,11,5)*WaveLenght2)
%% Независимая от ионосферы линейная комбинация для псевдорасстояний (PI) для станции zimm для спутника 21
plot(t,k1*zimm(t,21,2)+k2*zimm(t,21,3))
%% Независимая от ионосферы линейная комбинация для фазовых измерений (LI) для станции zimm для спутника 21
plot(t,k1*zimm(t,21,4)*WaveLenght1+k2*zimm(t,21,5)*WaveLenght2)
%% Независимая от ионосферы линейная комбинация для псевдорасстояний (PI) для станции zimm для спутника 31
plot(t,k1*zimm(t,31,2)+k2*zimm(t,31,3))
%% Независимая от ионосферы линейная комбинация для фазовых измерений (LI) для станции zimm для спутника 31
plot(t,k1*zimm(t,31,4)*WaveLenght1+k2*zimm(t,31,5)*WaveLenght2)
%% MELBOURNE-W AUBBENA LINEAR COMBINATION
%% Widelane linear combination
kmwl1=f1/(f1-f2);
kmwl2=f2/(f1-f2);
kmwp1=f1/(f1+f2);
kmwp2=f2/(f1+f2);
%% Widelane линейная комбинация для псевдорасстояний (PW) для станции zimj для спутника 1
plot(t,kmwp1*zimj(t,1,2)+kmwp2*zimj(t,1,3))
%% Widelane линейная комбинация для фазовых измерений (LW) для станции zimj для спутника 1
plot(t,kmwl1*zimj(t,1,4)*WaveLenght1-kmwl2*zimj(t,1,5)*WaveLenght2)
%% Widelane линейная комбинация для псевдорасстояний (PW) для станции onsa для спутника 1
plot(t,kmwp1*onsa(t,1,2)+kmwp2*onsa(t,1,3))
%% Widelane линейная комбинация для фазовых измерений (LW) для станции onsa для спутника 1
plot(t,kmwl1*onsa(t,1,4)*WaveLenght1-kmwl2*onsa(t,1,5)*WaveLenght2)
%% Widelane линейная комбинация для псевдорасстояний (PW) для станции wtzr для спутника 1
plot(t,kmwp1*wtzr(t,1,2)+kmwp2*wtzr(t,1,3))
%% Widelane линейная комбинация для фазовых измерений (LW) для станции wtzr для спутника 1
plot(t,kmwl1*wtzr(t,1,4)*WaveLenght1-kmwl2*wtzr(t,1,5)*WaveLenght2)
%% Melbourne-W Aubbena линейная комбинация для станции zimj для спутника 1
plot(t,(kmwl1*zimj(t,1,4)*WaveLenght1-kmwl2*zimj(t,1,5)*WaveLenght2)-(kmwp1*zimj(t,1,2)+kmwp2*zimj(t,1,3)))
%% Melbourne-W Aubbena линейная комбинация для станции onsa для спутника 1
plot(t,(kmwl1*onsa(t,1,4)*WaveLenght1-kmwl2*onsa(t,1,5)*WaveLenght2)-(kmwp1*onsa(t,1,2)+kmwp2*onsa(t,1,3)))
%% Melbourne-W Aubbena линейная комбинация для станции wtzr для спутника 1
plot(t,(kmwl1*wtzr(t,1,4)*WaveLenght1-kmwl2*wtzr(t,1,5)*WaveLenght2)-(kmwp1*wtzr(t,1,2)+kmwp2*wtzr(t,1,3)))
%% Melbourne-W Aubbena линейная комбинация для станции zimj для спутника 11
plot(t,(kmwl1*zimj(t,11,4)*WaveLenght1-kmwl2*zimj(t,11,5)*WaveLenght2)-(kmwp1*zimj(t,11,2)+kmwp2*zimj(t,11,3)))
%% Melbourne-W Aubbena линейная комбинация для станции onsa для спутника 11
plot(t,(kmwl1*onsa(t,11,4)*WaveLenght1-kmwl2*onsa(t,11,5)*WaveLenght2)-(kmwp1*onsa(t,11,2)+kmwp2*onsa(t,11,3)))
%% Melbourne-W Aubbena линейная комбинация для станции wtzr для спутника 11
plot(t,(kmwl1*wtzr(t,11,4)*WaveLenght1-kmwl2*wtzr(t,11,5)*WaveLenght2)-(kmwp1*wtzr(t,11,2)+kmwp2*wtzr(t,11,3)))
%% Melbourne-W Aubbena линейная комбинация для станции zimj для спутника 21
plot(t,(kmwl1*zimj(t,21,4)*WaveLenght1-kmwl2*zimj(t,21,5)*WaveLenght2)-(kmwp1*zimj(t,21,2)+kmwp2*zimj(t,21,3)))
%% Melbourne-W Aubbena линейная комбинация для станции onsa для спутника 21
plot(t,(kmwl1*onsa(t,21,4)*WaveLenght1-kmwl2*onsa(t,21,5)*WaveLenght2)-(kmwp1*onsa(t,21,2)+kmwp2*onsa(t,21,3)))
%% Melbourne-W Aubbena линейная комбинация для станции wtzr для спутника 21
plot(t,(kmwl1*wtzr(t,21,4)*WaveLenght1-kmwl2*wtzr(t,21,5)*WaveLenght2)-(kmwp1*wtzr(t,21,2)+kmwp2*wtzr(t,21,3)))