% Чтение RINEX
gope=readrnx('gope1500.07o');
onsa=readrnx('onsa1500.07o');
kiru=readrnx('kiru1500.07o');
zimm=readrnx('zimm1500.07o');
%% Дост-ть сп-в gope C1
imagesc(gope(:,:,1))
%% gope P2
imagesc(gope(:,:,3))
%% gope S1
imagesc(gope(:,:,8))
%% gope S2
imagesc(gope(:,:,9))
%% onsa С1
imagesc(onsa(:,:,1))
%% onsa P2
imagesc(onsa(:,:,3))
%% onsa S1
imagesc(onsa(:,:,8))
%% onsa S2
imagesc(onsa(:,:,9))
%% kiru С1
imagesc(kiru(:,:,1))
%% kiru P2
imagesc(kiru(:,:,3))
%% kiru S1
imagesc(kiru(:,:,8))
%% kiru S2
imagesc(kiru(:,:,9))
%%
t=1:2880; % list of epoch index
%% Пок-я P1 ст. gope сп-в 10, 20, 30 и 40
plot(t,[gope(t,10,2),gope(t,20,2),gope(t,30,2),gope(t,40,2)])
%% P2 для gope
plot(t,[gope(t,10,3),gope(t,20,3),gope(t,30,3),gope(t,40,3)])
%% L1 для gope
plot(t,[gope(t,10,4),gope(t,20,4),gope(t,30,4),gope(t,40,4)])
%% L2 для gope
plot(t,[gope(t,10,5),gope(t,20,5),gope(t,30,5),gope(t,40,5)])
%% P1 для onsa
plot(t,[onsa(t,10,2),onsa(t,20,2),onsa(t,30,2),onsa(t,40,2)])
%% P2 для onsa
plot(t,[onsa(t,10,3),onsa(t,20,3),onsa(t,30,3),onsa(t,40,3)])
%% L1 для onsa
plot(t,[onsa(t,10,4),onsa(t,20,4),onsa(t,30,4),onsa(t,40,4)])
%% L2 для onsa
plot(t,[onsa(t,10,5),onsa(t,20,5),onsa(t,30,5),onsa(t,40,5)])
%% P1 для kiru
plot(t,[kiru(t,10,2),kiru(t,20,2),kiru(t,30,2),kiru(t,40,2)])
%% P2 для kiru
plot(t,[kiru(t,10,3),kiru(t,20,3),kiru(t,30,3),kiru(t,40,3)])
%% L1 для kiru
plot(t,[kiru(t,10,4),kiru(t,20,4),kiru(t,30,4),kiru(t,40,4)])
%% L2 для kiru
plot(t,[kiru(t,10,5),kiru(t,20,5),kiru(t,30,5),kiru(t,40,5)])
%% D1 для gope
plot(t,[gope(t,10,6),gope(t,20,6),gope(t,30,6),gope(t,40,6)])
%% D2 для gope
plot(t,[gope(t,10,7),gope(t,20,7),gope(t,30,7),gope(t,40,7)])
%% D1 для onsa
plot(t,[onsa(t,10,6),onsa(t,20,6),onsa(t,30,6),onsa(t,40,6)])
%% D2 для onsa
plot(t,[onsa(t,10,7),onsa(t,20,7),onsa(t,30,7),onsa(t,40,7)])
%% D1 для kiru
plot(t,[kiru(t,10,6),kiru(t,20,6),kiru(t,30,6),kiru(t,40,6)])
%% D2 для kiru
plot(t,[kiru(t,10,7),kiru(t,20,7),kiru(t,30,7),kiru(t,40,7)])
%% S1 для gope
plot(t,[gope(t,10,8),gope(t,20,8),gope(t,30,8),gope(t,40,8)])
%% S2 для gope
plot(t,[gope(t,10,9),gope(t,20,9),gope(t,30,9),gope(t,40,9)])
%% S1 для onsa
plot(t,[onsa(t,10,8),onsa(t,20,8),onsa(t,30,8),onsa(t,40,8)])
%% S2 для onsa
plot(t,[onsa(t,10,9),onsa(t,20,9),onsa(t,30,9),onsa(t,40,9)])
%% S1 для kiru
plot(t,[kiru(t,10,8),kiru(t,20,8),kiru(t,30,8),kiru(t,40,8)])
%% S2 для kiru
plot(t,[kiru(t,10,9),kiru(t,20,9),kiru(t,30,9),kiru(t,40,9)])
%% Зав-ть S1 от P1 ст. gope сп-ка 10
plot(gope(t,10,2),gope(t,10,8))
%% S2 от P2 для gope
plot(gope(t,10,3),gope(t,10,9))
%% S1 от P1 для onsa
plot(onsa(t,10,2),onsa(t,10,8))
%% S2 от P2 для onsa
plot(onsa(t,10,3),onsa(t,10,9))
%% S1 от P1 для kiru
plot(kiru(t,10,2),kiru(t,10,8))
%% S2 от P2 для kiru
plot(kiru(t,10,3),kiru(t,10,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 ст. gope сп. 10
plot(t,k1*gope(t,10,2)+k2*gope(t,10,3))
%% Незав-я от ионосферы лин. комб. фазовых измерений LI ст. gope сп. 10
plot(t,k1*gope(t,10,4)*WaveLenght1+k2*gope(t,10,5)*WaveLenght2)
%% PI для onsa
plot(t,k1*onsa(t,10,2)+k2*onsa(t,10,3))
%% LI для onsa
plot(t,k1*onsa(t,10,4)*WaveLenght1+k2*onsa(t,10,5)*WaveLenght2)
%% PI для kiru
plot(t,k1*kiru(t,10,2)+k2*kiru(t,10,3))
%% LI для kiru
plot(t,k1*kiru(t,10,4)*WaveLenght1+k2*kiru(t,10,5)*WaveLenght2)
%% PI для zimm
plot(t,k1*zimm(t,10,2)+k2*zimm(t,10,3))
%% LI для zimm
plot(t,k1*zimm(t,10,4)*WaveLenght1+k2*zimm(t,10,5)*WaveLenght2)
%% PI для zimm сп. 20
plot(t,k1*zimm(t,20,2)+k2*zimm(t,20,3))
%% LI для zimm
plot(t,k1*zimm(t,20,4)*WaveLenght1+k2*zimm(t,20,5)*WaveLenght2)
%% PI для zimm сп. 30
plot(t,k1*zimm(t,30,2)+k2*zimm(t,30,3))
%% LI для zimm
plot(t,k1*zimm(t,30,4)*WaveLenght1+k2*zimm(t,30,5)*WaveLenght2)
%% PI для zimm сп. 40
plot(t,k1*zimm(t,40,2)+k2*zimm(t,40,3))
%% LI для zimm сп. 40
plot(t,k1*zimm(t,40,4)*WaveLenght1+k2*zimm(t,40,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 ст. gope сп. 10
plot(t,kmwp1*gope(t,10,2)+kmwp2*gope(t,10,3))
%% Widelane лин. комб. фазовых измерений LW ст. gope сп. 10
plot(t,kmwl1*gope(t,10,4)*WaveLenght1-kmwl2*gope(t,10,5)*WaveLenght2)
%% PW onsa
plot(t,kmwp1*onsa(t,10,2)+kmwp2*onsa(t,10,3))
%% LW onsa
plot(t,kmwl1*onsa(t,10,4)*WaveLenght1-kmwl2*onsa(t,10,5)*WaveLenght2)
%% PW kiru
plot(t,kmwp1*kiru(t,10,2)+kmwp2*kiru(t,10,3))
%% LW kiru
plot(t,kmwl1*kiru(t,10,4)*WaveLenght1-kmwl2*kiru(t,10,5)*WaveLenght2)
%% Melbourne-W Aubbena лин. комб. ст. gope сп. 10
plot(t,(kmwl1*gope(t,10,4)*WaveLenght1-kmwl2*gope(t,10,5)*WaveLenght2)-(kmwp1*gope(t,10,2)+kmwp2*gope(t,10,3)))
%% onsa
plot(t,(kmwl1*onsa(t,10,4)*WaveLenght1-kmwl2*onsa(t,10,5)*WaveLenght2)-(kmwp1*onsa(t,10,2)+kmwp2*onsa(t,10,3)))
%% kiru
plot(t,(kmwl1*kiru(t,10,4)*WaveLenght1-kmwl2*kiru(t,10,5)*WaveLenght2)-(kmwp1*kiru(t,10,2)+kmwp2*kiru(t,10,3)))
%% Melbourne-W Aubbena лин. комб. ст. gope сп. 20
plot(t,(kmwl1*gope(t,20,4)*WaveLenght1-kmwl2*gope(t,20,5)*WaveLenght2)-(kmwp1*gope(t,20,2)+kmwp2*gope(t,20,3)))
%% onsa
plot(t,(kmwl1*onsa(t,20,4)*WaveLenght1-kmwl2*onsa(t,20,5)*WaveLenght2)-(kmwp1*onsa(t,20,2)+kmwp2*onsa(t,20,3)))
%% kiru
plot(t,(kmwl1*kiru(t,20,4)*WaveLenght1-kmwl2*kiru(t,20,5)*WaveLenght2)-(kmwp1*kiru(t,20,2)+kmwp2*kiru(t,20,3)))
%% Melbourne-W Aubbena лин. комб. ст. gope сп. 30
plot(t,(kmwl1*gope(t,30,4)*WaveLenght1-kmwl2*gope(t,30,5)*WaveLenght2)-(kmwp1*gope(t,30,2)+kmwp2*gope(t,30,3)))
%% onsa
plot(t,(kmwl1*onsa(t,30,4)*WaveLenght1-kmwl2*onsa(t,30,5)*WaveLenght2)-(kmwp1*onsa(t,30,2)+kmwp2*onsa(t,30,3)))
%% kiru
plot(t,(kmwl1*kiru(t,30,4)*WaveLenght1-kmwl2*kiru(t,30,5)*WaveLenght2)-(kmwp1*kiru(t,30,2)+kmwp2*kiru(t,30,3)))