Scribd
Unggah
177 tayangan6 halaman

Proyeksi Geodesi: Panjang Busur Meridian

Dokumen tersebut berisi penyelesaian soal hitung proyeksi geodesi yang meliputi perhitungan panjang busur meridian dan paralel dengan menggunakan rumus matematika dan program MATLAB. Hasil perhitungan menunjukkan nilai panjang busur meridian lebih kecil dari panjang busur paralel.

Diunggah oleh

Benita Roseana
Hak Cipta
© All Rights Reserved
Kami menangani hak cipta konten dengan serius. Jika Anda merasa konten ini milik Anda, ajukan klaim di sini.
Format Tersedia
Unduh sebagai DOCX, PDF, TXT atau baca online di Scribd
177 tayangan6 halaman

Proyeksi Geodesi: Panjang Busur Meridian

Dokumen tersebut berisi penyelesaian soal hitung proyeksi geodesi yang meliputi perhitungan panjang busur meridian dan paralel dengan menggunakan rumus matematika dan program MATLAB. Hasil perhitungan menunjukkan nilai panjang busur meridian lebih kecil dari panjang busur paralel.

Diunggah oleh

Benita Roseana
Hak Cipta
© All Rights Reserved
Kami menangani hak cipta konten dengan serius. Jika Anda merasa konten ini milik Anda, ajukan klaim di sini.
Format Tersedia
Unduh sebagai DOCX, PDF, TXT atau baca online di Scribd

Nama : Rika Enjelina Pidu

NIM : 21110115140084
Kelas : B

TUGAS HITUNG PROYEKSI GEODESI

SOAL :

1) Buat gambar irisan normal antara titik p1 dan p2!


2) Hitung panjang busur meridian pada :
A. 0° - 30°
B. 30° - 45°
C. 45° - 60°
D. 60° - 90°
3) Hitung panjang busur paralel pada :
A. 0° - 30°
B. 30° - 45°
C. 45° - 60°
D. 60° - 90°
Tentukan Analisis dari hasil perhitungan tersebut!

JAWABAN :

1) Gambar Irisan Normal antara Titik P1 dan P2 :


(dalam kertas folio)
2) Rumus Panjang Busur Meridian :

B C
Sφ = a(1-e2 )[A(φ2 - φ1) – (sin2φ2 – sin2φ1) + (sin4φ2 – sin4φ1) –
2 4
D E F
(sin6φ2 – sin6φ1) + (sin8φ2 – sin8φ1) – (sin10φ2 – sin10φ1) +
6 8 10
….]

Nilai Panjang Busur Meridian (Sφ) :


A. Untuk 0° - 30° = 191009594.2843
B. Untuk 30° - 45° = 95509571.4908
C. Untuk 45° - 60° = 95513869.0154
D. Untuk 60° - 90° = 191037374.2044
3) Rumus Panjang Busur Paralel :

P = Ncosφ , P = Radius lingkaran parallel

Sλ = (λ2 – λ1) P = (λ2 – λ1) Ncosφ

Nilai Panjang Busur Paralel (Sλ) :


A. Untuk 0° - 30° = 191344110
B. Untuk 30° - 45° = 95672055
C. Untuk 45° - 60° = 95672055
D. Untuk 60° - 90° = 191344110

Analisis :
Dari hasil perhitungan diatas dapat disimpulkan bahwa nilai panjang busur meridian (Sφ)
tidak sama dengan nilai panjang busur parallel(Sλ), akan tetapi panjang busur meridian
lebih kecil dari pada panjang busur parallel (Sφ < Sλ).

Script Matlab :

%Tugas HITUNG PROYEKSI GEODESI


%Nama :Rika Enjelina Pidu
%NIM :21110115140084
%Kelas :B

clear all;
clc;

disp('-------------------------------------------------------------------');
disp('_________Progam Hitung Panjang Busur Meridian dan Pararel__________');
disp('____________________Hitung Proyeksi Geodesi________________________');
disp('_____________Oleh Rika Enjelina Pidu (21110115140084)______________');
disp('-------------------------------------------------------------------');
disp(' ');
disp('PILIH SOAL:');
disp('1. Hitung Panjang Busur Meridian (Sl)');
disp('2. Hitung Panjang Busur Paralel (Sb)');
soal=input('Pilih soal (1/2) ===> ');
switch soal
case 1
disp('# Panjang Busur Meridian');
disp('Pilih data lintang');
disp('1. 0 - 30');
disp('2. 30 - 45');
disp('3. 45 - 60');
disp('4. 60 - 90');
Lin=input('pilih sudut lintang (1/2/3/4)==> ');
%Datum WGS 84
disp(' ');
a=6378137;
b=6356752;
e=((a^2-b^2)/(a^2))^(1/2);
%Perhitungan Busur Meridian
switch Lin
case 1
l1=0;
l2=30;

G=1+((3/4)*(e^2))+((45/64)*(e^4))+((175/256)*(e^6))+((11025/16384)*(e^8))+((4
3659/65536)*(e^10));

H=(3/4)*(e^2)+(15/16)*(e^4)+(525/512)*(e^6)+(2205/2048)*(e^8)+(72765/65536)*(
e^10);

I=(15/64)*(e^4)+(105/256)*(e^6)+(2205/2048)*(e^8)+(10395/16384)*(e^10);
J=(35/512)*(e^6)+(315/2048)*(e^8)+(31185/131072)*(e^10);
K=(315/16384)*(e^8)+(3465/65536)*(e^10);
L=(693/131072)*(e^10);
Sl=a*(1-(e^2))*(G*(l2-l1)-(H/2)*(sind(2*l2)-
sind(2*l1))+(I/4)*(sind(4*l2)-sind(4*l1))-(J/6)*(sind(6*l2)-
sind(6*l1))+(K/8)*(sind(8*l2)-sind(8*l1))-(L/10)*(sind(10*l2)-sind(10*l1)));
disp(['Nilai Panjang Busur Meridian (S) = ' num2str(Sl)]);
case 2
l1=30;
l2=45;

G=1+(3/4)*(e^2)+(45/64)*(e^4)+(175/256)*(e^6)+(11025/16384)*(e^8)+(43659/6553
6)*(e^10);

H=(3/4)*(e^2)+(15/16)*(e^4)+(525/512)*(e^6)+(2205/2048)*(e^8)+(72765/65536)*(
e^10);

I=(15/64)*(e^4)+(105/256)*(e^6)+(2205/2048)*(e^8)+(10395/16384)*(e^10);
J=(35/512)*(e^6)+(315/2048)*(e^8)+(31185/131072)*(e^10);
K=(315/16384)*(e^8)+(3465/65536)*(e^10);
L=(693/131072)*(e^10);
Sl=a*(1-(e^2))*(G*(l2-l1)-(H/2)*(sind(2*l2)-
sind(2*l1))+(I/4)*(sind(4*l2)-sind(4*l1))-(J/6)*(sind(6*l2)-
sind(6*l1))+(K/8)*(sind(8*l2)-sind(8*l1))-(L/10)*(sind(10*l2)-sind(10*l1)));
disp(['Nilai Panjang Busur Meridian (S) = ' num2str(Sl)]);
case 3
l1=45;
l2=60;

G=1+(3/4)*(e^2)+(45/64)*(e^4)+(175/256)*(e^6)+(11025/16384)*(e^8)+(43659/6553
6)*(e^10);

H=(3/4)*(e^2)+(15/16)*(e^4)+(525/512)*(e^6)+(2205/2048)*(e^8)+(72765/65536)*(
e^10);

I=(15/64)*(e^4)+(105/256)*(e^6)+(2205/2048)*(e^8)+(10395/16384)*(e^10);
J=(35/512)*(e^6)+(315/2048)*(e^8)+(31185/131072)*(e^10);
K=(315/16384)*(e^8)+(3465/65536)*(e^10);
L=(693/131072)*(e^10);
Sl=a*(1-(e^2))*(G*(l2-l1)-(H/2)*(sind(2*l2)-
sind(2*l1))+(I/4)*(sind(4*l2)-sind(4*l1))-(J/6)*(sind(6*l2)-
sind(6*l1))+(K/8)*(sind(8*l2)-sind(8*l1))-(L/10)*(sind(10*l2)-sind(10*l1)));
disp(['Nilai Panjang Busur Meridian (S) = ' num2str(Sl)]);
case 4
l1=60;
l2=90;

G=1+(3/4)*(e^2)+(45/64)*(e^4)+(175/256)*(e^6)+(11025/16384)*(e^8)+(43659/6553
6)*(e^10);

H=(3/4)*(e^2)+(15/16)*(e^4)+(525/512)*(e^6)+(2205/2048)*(e^8)+(72765/65536)*(
e^10);

I=(15/64)*(e^4)+(105/256)*(e^6)+(2205/2048)*(e^8)+(10395/16384)*(e^10);
J=(35/512)*(e^6)+(315/2048)*(e^8)+(31185/131072)*(e^10);
K=(315/16384)*(e^8)+(3465/65536)*(e^10);
L=(693/131072)*(e^10);
Sl=a*(1-(e^2))*(G*(l2-l1)-(H/2)*(sind(2*l2)-
sind(2*l1))+(I/4)*(sind(4*l2)-sind(4*l1))-(J/6)*(sind(6*l2)-
sind(6*l1))+(K/8)*(sind(8*l2)-sind(8*l1))-(L/10)*(sind(10*l2)-sind(10*l1)));
disp(['Nilai Panjang Busur Meridian (S) = ' num2str(Sl)]);
otherwise
disp('ERROR!!!');
end
case 2
disp('# Panjang Busur Pararel');
disp('Pilih data busur');
disp('1. 0 - 30');
disp('2. 30 - 45');
disp('3. 45 - 60');
disp('4. 60 - 90');
Buj=input('pilih sudut busur (1/2/3/4)==> ');
L=input('Besar Lintang = ');
disp(' ');
%Datum WGS 84
disp(' ');
a=6378137;
b=6356752;
e=((a^2-b^2)/(a^2))^(1/2);
N=a/((1-(e^2)*((sind(L))^2))^0.5);
%Perhitungan Busur Pararel
switch Buj
case 1
b1=0;
b2=30;
p=N*cosd(L);
Sb=(b2-b1)*p;
disp(['Nilai Panjang Busur Paralel = ' num2str(Sb)]);
case 2
b1=30;
b2=45;
p=N*cosd(L);
Sb=(b2-b1)*p;
disp(['Nilai Panjang Busur Paralel = ' num2str(Sb)]);
case 3
b1=45;
b2=60;
p=N*cosd(L);
Sb=(b2-b1)*p;
disp(['Nilai Panjang Busur Paralel = ' num2str(Sb)]);
case 4
b1=60;
b2=90;
p=N*cosd(L);
Sb=(b2-b1)*p;
disp(['Nilai Panjang Busur Paralel = ' num2str(Sb)]);
otherwise
disp('ERROR!');
end
otherwise
disp('ERROR!');
end
disp('-----------------------------------------------------------');
disp('-----------------------------------------------------------');
disp('Terima kasih sudah memakai program ini, Pakai lagi?');
disp('YA (1)');
disp('TIDAK (2)');
pilih= input('Pilih =======> ');
switch pilih
case 1
run C:\Users\OPTION\Documents\RIKA.m;
case 2
disp('BYE');
otherwise
disp('Tidak dapat diproses');
end

Pembuktian
------------------------------------------------------------------
_________Progam Hitung Panjang Busur Meridian dan Pararel__________
____________________Hitung Proyeksi Geodesi________________________
_____________Oleh Rika Enjelina Pidu (21110115140084)______________
-------------------------------------------------------------------

PILIH SOAL:
1. Hitung Panjang Busur Meridian (Sl)
2. Hitung Panjang Busur Paralel (Sb)
Pilih soal (1/2) ===> 2
# Panjang Busur Pararel
Pilih data busur
1. 0 - 30
2. 30 - 45
3. 45 - 60
4. 60 - 90
pilih sudut busur (1/2/3/4)==> 2
Besar Lintang = 20

Nilai Panjang Busur Paralel = 89937546.2755


-----------------------------------------------------------
-----------------------------------------------------------

Anda mungkin juga menyukai