How to plot and create a Gravity-Altitude graph?
4 次查看(过去 30 天)
显示 更早的评论
I have been tasked to create or make a gravitational force-altitude graph from the surface to the top of the atmosphere. However, I am confused and I do not know how to start it because I don't have anything yet aside from the essential variables to be used such as the equation to get the gravity at any latitude, the equatorial radius, and other needed constant variables and a code (acquired from Google) that seemingly fails. Can anyone help me start?
By the way, this is the code:
function [ g ] = Gravity( z )
% this function used to get g at any altitude
%% inputs:
% z : the altitude above the earth (km)
%% outputs:
% g : the gravity at altitude z
g0 = 9.81; % m/s
Re = 6378; % km
g = g0*Re^2/(Re+z)^2; % m/s^2
end
0 个评论
回答(1 个)
KSSV
2021-3-19
You need to input z values....Do you have those values along with lat, lon?
You may check the below example with dummy data.
clc; clear all ;
x = -10:1:10 ; nx = length(x) ; % lat
y = -10:1:10 ; ny = length(y) ; % lon
[X,Y] = meshgrid(x,y) ; % mesh
z = randi(1000,ny,nx) ; % random altitudes within 1000 units
% this function used to get g at any altitude
%% inputs:
% z : the altitude above the earth (km)
%% outputs:
% g : the gravity at altitude z
g0 = 9.81; % m/s
Re = 6378; % km
g = g0*Re^2./(Re+z).^2; % m/s^2
surf(X,Y,g)
4 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Earth and Planetary Science 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!