Eka Rinex Gnss

1 clear all;
2 clc;
3 disp(
'| |'
);
4 disp('|Nama: Eka Diah Nur
Safitri |');
5 disp('|NRP:
03311740000006 |');
6 disp('|Matkul: Asistensi GNSS -
A |');
7 disp('|Tugas: Rinex to
Geosentrik |');
8 disp(
'======================================================================================'
); disp(' ');
9 disp(
'======================================================================================'
); disp(' ');
10 disp('= 2.11 N: GPS NAV DATA Reciever Operator 01-SEP-19
00.00.02 =');
11 disp(
'= ='
);
12 disp('= 9 19 9 1 0 0 0.0 -0.200481154025D-0 -0.920863385545D-11
0.000000000000D+00 =');
13 disp('= 0.910000000000D+02 0.312500000000D-01 0.455947563475D-08
0.304179685436D+01 =');
14 disp('= 0.111758708954D-06 0.133976177312D-02 0.137183815241D-04
0.515361626244D+04 =');
15 disp('= 0.000000000000D+00 0.875443220139D-07 -0.265724605959D+01
-0.130385160446D-07 =');
16 disp('= 0.952880411385D+00 0.110343750000D+03 0.168066173056D+01
-0.799176146037D-08 =');
17 disp('= 0.308941440080D-09 0.100000000000D+01 0.206900000000D+04
0.000000000000D+00 =');
18 disp('= 0.240000000000D+01 0.000000000000D+00 0.139698386192D-08
0.910000000000D+02 =');
19 disp('= 0.597618000000D+06 0.400000000000D+01 =');
20 disp(
'======================================================================================'
); disp(' ');
21 disp('Data diatas merupakan dta navigation file GPS format RINEX versi NetR9');
22 disp('pada tanggal 9 Maret 2019 pada jam 00.00.02 oleh satelit nomor 9');
23 disp('Berikut adalah perhitungan pengolahan Data Rinex Navigasi diatas : '); disp(' ');
24 NoSat=9; a0=-0.200481154025D-0; a1=-0.920863385545D-11; a2=0.000000000000D+00;
25 IODE=0.910000000000D+02; Crs=0.312500000000D-01; deltan=0.455947563475D-08; M0=
0.304179685436D+01;
26 Cuc=0.111758708954D-06; e=0.133976177312D-02; Cus=0.137183815241D-04; sqrtA=
0.515364204979D+04;
27 toe=0.000000000000D+00 ; Cic=0.875443220139D-07; KR0=-0.265724605959D+01; Cis=-
0.130385160446D-07;
28 I0=0.952880411385D+00; Crc=0.110343750000D+03; omega=0.168066173056D+01; KR=-
0.799176146037D-08;
29 I=0.308941440080D-09 ; L2=0.100000000000D+01 ; GPSweek=0.206900000000D+04; L2Pcode=
0.000000000000D+00;
30 SatAcc=0.240000000000D+01; SatHea=0.000000000000D+00; TGD=0.139698386192D-08; IODC=
0.910000000000D+02;
31 TransTime=0.597618000000D+06; FitInt=0.400000000000D+01;
32 time=2; toc=(9*3600);
33 format long;
34 disp('1. Untuk menhitung posisi satelit menggunakan broadcast ephemerides,');
35 disp(' diperlukan nilai dari konstanta gravitasi bumi (GM) dan');
36 disp(' kecepatan rotasi bumi (KRe)');
You created this PDF from an application that is not licensed to print to novaPDF printer (http://www.novapdf.com)
37 disp('GM = 3986001.5 x 10^8 m^3/s');
38 disp('KRe = 7.292115147 x 10^-5 rad/s');
39 GM = 3986001.5*10^8; KRe = 7.292115147*10^-5; disp(' ');
40 disp('2. Waktu pengamatan terkoreksi');
41 t=time-(a0+(a1*(time-toc)))+(a2*(time-toc)^2);
42 disp(['t = ',num2str(t),'sekon']); disp(' ');
43 disp('3. Waktu');
44 tk=((5*24*3600)+t)-toe;
45 disp(['tk = ',num2str(tk),'sekon']); disp(' ');
46 disp('4. Sumbu panjang ellipsoid');
47 a=sqrtA^2;
48 disp(['a = ',num2str(a), 'meter']); disp(' ');
49 disp('5. Mean motion terkoreksi;');
50 n=sqrt(GM/(a^3))+deltan;
51 disp(['n = ',num2str(n), 'rad/s']);disp(' ');
52 disp('6. Mean anomaly terkoreksi');
53 Mk=M0+(n*tk);
54 disp(['Mk =',num2str(Mk),'rad']); disp(' ');
55 disp('7. Eccentric anomaly terkoreksi');
56 Ek0=Mk, Ek1=Mk+(e*sin(Ek0)), Ek2=Mk+(e*sin(Ek1)), Ek=Mk+(e*sin(Ek2))
57 disp(['Ek =',num2str(Ek),'rad']); disp(' ');
58 disp('8. True anomaly terkoreksi');
59 vk=2*atan(sqrt((1+e)/(1-e))*(tan(Ek/2)));
60 disp(['vk = ',num2str(vk),'rad']); disp(' ');
61 disp('9. Argument of latitude terkoreksi');
62 Ok=vk+omega; disp(['Ok =', num2str(Ok),'rad']); disp(' ');
63 disp('10. Koreksi argument of latitude');
64 duk=(Cuc*cos(2*Ok))+(Cus*sin(2*Ok));
65 disp(['Koreksi argument of latitude = ',num2str(duk), 'rad']); disp(' ');
66 disp('11. Argument of latitude terkoreksi');
67 uk=Ok+duk;
68 disp(['uk = ',num2str(uk),'rad']); disp(' ');
69 disp('12. Koreksi radius');
70 drk=(Crc*cos(2*Ok))+(Crs*sin(2*Ok));
71 disp(['Koreksi radius = ',num2str(drk), 'meter']); disp(' ');
72 disp('13. Radius terkoreksi');
73 rk=(a*(1-(e*cos(Ek))))+drk;
74 disp(['rk = ',num2str(rk),' meter']); disp(' ');
75 disp('14. Koreksi inklinasi');
76 dik=(Cic*sin(2*Ok))+(Cis*cos(2*Ok));
77 disp(['Koreksi inklinasi = ',num2str(dik), 'rad']); disp(' ');
78 disp('15. Inklinasi Terkoreksi');
79 ik=I0+(I*tk)+dik;
80 disp(['ik = ',num2str(ik),'rad']); disp(' ');
81 disp('16. Posisi pada bidang orbit');
82 xk=rk*cos(uk); yk=rk*sin(uk);
83 disp(['xk = ',num2str(xk), 'meter']);
84 disp(['yk = ',num2str(yk), 'meter']); disp(' ');
85 disp('17. Bujur ascending node terkoreksi');
86 KRk=KR0+((KR-KRe)*tk)-(KRe*toe);
87 disp(['Bujur ascending node terkoreksi = ',num2str(KRk),'rad']); disp(' ');
88 disp('18. Koordinat geosentrik satelit terikat bumi');
89 Xk=(xk*cos(KRk))-(yk*sin(KRk)*cos(ik));
90 Yk=(xk*sin(KRk))+(yk*cos(KRk)*cos(ik));
91 Zk=yk*sin(ik);
92 disp(['Xk = ',num2str(Xk),'meter']);
93 disp(['Yk = ',num2str(Yk),'meter']);
94 disp(['Zk = ',num2str(Zk),'meter']); disp(' ');
95 disp('19. Pembuktian');
96 RK=sqrt(Xk^2+Yk^2+Zk^2)
You created this PDF from an application that is not licensed to print to novaPDF printer (http://www.novapdf.com)

Eka Rinex Gnss

Uploaded by

Copyright:

Available Formats

Eka Rinex Gnss

Uploaded by

Document Information

Original Description:

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Eka Rinex Gnss

Uploaded by

Copyright:

Available Formats

1 clear all;

You might also like