Noted, where cos f =cos (L)x cos (lat), and value of L equal with absolute value from Les-Ls. So from that information we can find value of E using MATLAB with this source code:
function Elv = satelev(lat,lgtd) Re=6378; %Radius Bumi Ro=35786; %Ketinggian Altitude satelit Pada Orbit Geostasioner a=lat*(pi/180); %latitude stasiun bumi (lintang) b=lgtd*(pi/180); %Selisih longitude satelit (bujur) Elv=(atan(((cos(a)*cos(b))-(Re/(Re+Ro)))/((1-(cos(a)*cos(b))^2)^(1/2))))*(180/pi); end
save this function to satelev.m (this matlab file). You have to know, This program is a function program, so when you want to use it. you have to declare the latitude of your ES and the value of L, using this source code
lgtd=value of lgtd; lat=value of lat;
after that you can run the function and find the value with this source code
E= satelev(lat,lgtd);
0 Comments
Post a Comment