lama juga ga buka blog ini.
hha
iseng-iseng
ngerjain tugas dan sebenarnya cuman butuh kalkulator buat ngerjain aja sih...
tp berhubung iseng..
script "keisengan" :
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
clear
clc
format long g
%% Input
tet = [6 53 29.49737];
alp = [107 36 43.21995];
h= 816.5963;
tet = dms2degrees(tet);
alp = dms2degrees(alp);
a = 6378137;
b = 6356752;
%soal 2.b.1
as1 = 6378165;
bs1 = 6356783.3;
%soal 2.b.2
as2 = 6356760;
bs2 = 6378145;
%soal 2.b.3
Es3 = 0.082136016;
%soal 2.b.4
fs4 = 1/151;
%% program
%soal 2.b.1 &b.2 nilai tet,alp,h -> X,Y,Z (WGS'84) -> te,alp,h (WGS'2.b.1&2)
% tet,alp,h -> X,Y,Z
E = ((a^2 - b^2)/(b^2))
N = ((a)/(sqrt(1-(E*sind(tet)))))
Xwgs84 =(N + h)*cosd(tet)*cosd(alp)
Ywgs84 =(N + h)*cosd(tet)*sind(alp)
Zwgs84 =(N*(1-E)+h)*sind(tet)
%soal 2.b.1
%X,Y,Z -> tet,alp,h (wgs's1)
Es1 = ((as1^2 - bs1^2)/(bs1^2))
Ns1 = ((as1)/(sqrt(1-(E*sind(tet)))))
Ms1 = ((as1*(1-Es1))/((1-(E*sind(tet)^(3/2)))))
p = sqrt((Xwgs84*Xwgs84)+(Ywgs84*Ywgs84))
alp_s1 = atand (Ywgs84/Xwgs84)+180
tet01 = atand (Zwgs84/(p*(1-Es1)))
h01 = (sqrt((p^2)+(Zwgs84^2)))-((as1*sqrt(1-Es1))/(sqrt(1-Es1*((cosd(tet))^2))))
Z = (Ns1*(1-Es1)+h01)*sind(tet01)
P = (Ns1+h01)*cosd(tet01)
dZ = Zwgs84 - Z
dP = p - P
dtet = (dZ*cosd(tet01)-dP*sind(tet01))/(Ms1+h01)
dh = dZ*sind(tet01)+dP*cosd(tet01)
tets1= tet01 + dtet
hs1 = h01 +dh
% Tranformasi ke geosentrik
tetGeosentrik = atand ((sqrt(1-Es1))*tand(tets1));
tetGS = degrees2dms (tetGeosentrik)
% soal 2.b.2
%X,Y,Z -> tet,alp,h (wgs's1)
Es2 = ((as2^2 - bs2^2)/(bs2^2))
Ns2 = ((as2)/(sqrt(1-(E*sind(tet)))))
Ms2 = ((as1*(1-Es2))/((1-(E*sind(tet)^(3/2)))))
p2 = sqrt((Xwgs84*Xwgs84)+(Ywgs84*Ywgs84))
alp_s2 = atand (Ywgs84/Xwgs84)+180
tet02 = atand (Zwgs84/(p*(1-Es2)))
h02 = (sqrt((p^2)+(Zwgs84^2)))-((as1*sqrt(1-Es1))/(sqrt(1-Es1*((cosd(tet))^2))))
Z = (Ns1*(1-Es2)+h02)*sind(tet02)
P2 = (Ns1+h02)*cosd(tet02)
dZ2 = Zwgs84 - Z
dP2 = p2 - P2
dtet2 = (dZ*cosd(tet02)-dP2*sind(tet02))/(Ms2+h02)
dh = dZ2*sind(tet02)+dP2*cosd(tet02)
tets2= tet02 + dtet2
hs1 = h02 +dh
% Tranformasi ke geosentrik
tetGeosentrik2 = atand ((sqrt(1-Es2))*tand(tets2))
tetGSen2 = degrees2dms(tetGeosentrik2)
%soal 2.b.3
Es3 = 0.082136016
tetGeosentrik3 = atand ((sqrt(1-Es3))*tand(tet))
tetGS3 = degrees2dms (tetGeosentrik3)
%soal 2.b.4
Es4 = (1-(1-fs4)^2)
tetGeosentrik4 = atand ((sqrt(1-Es4))*tand(tet))
tetGS3 = degrees2dms (tetGeosentrik4)
%soal 2.c.
% koordinat XYZ dalam geosentrik
tetGS = atand ((sqrt(1-E))*tand(tet))
% tetGS = degrees2dms (tetGS)
N = ((a)/(sqrt(1-(E*sind(tetGS)))))
Xwgs84 =(N + h)*cosd(tetGS)*cosd(alp)
Ywgs84 =(N + h)*cosd(tetGS)*sind(alp)
Zwgs84 =(N*(1-E)+h)*sind(tetGS)
%soal 2.d
Esd = ((10^2 - 8^2)/(8^2))
N = ((10)/(sqrt(1-(Esd *sind(tet)))))
Msd = ((10*(1-Esd))/((1-(Esd *sind(tet)^(3/2)))))
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
sebenarnya ini ada yang salah sih.. kira-kira apa ya?
