How to plot the Gravitational Potential vs Radius of earth plot?

4 次查看（过去 30 天）

vimal kumar chawda 2020-5-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/523541-how-to-plot-the-gravitational-potential-vs-radius-of-earth-plot

评论： James Tursa 2020-5-17

clc;
clear all;
close all;
format long g;
G=6.6743015*10^-11;
Density=5.510;
M= (4*pi*G*Density)*10^6/3;
R=[10 20 50 100 6370]; 
%% for the calculation for the gravitational potential
for i =1:length(R)
    
Gravitational_Potential(:,i) = (M*R(:,i));
i=i+1;
end
figure()
plot(R,Gravitational_Potential);
ylabel('Gravitational Potential');
xlabel('Radius of the Sphere');
title('Gravitational potential vs Radius ');
grid on

I have to plot as curve plot as attached

Hint : - consider you have a function y = x^2. And now you have points x=1,2,4,8, from this you can get the corresponding y values: y=1^2,2^2,4^2,8^2. How do you visualize these values?

Now you have the same with a function V=GM/R, and R=10k,… you got V for all these values, you can visualize your results the same way you did above.

5 个评论
显示 3更早的评论隐藏 3更早的评论

vimal kumar chawda 2020-5-17

Thank you. Your links was one of the best explanation ever. Do you have code for all these cases?

James Tursa 2020-5-17

I don't have any code for this, but the equations look farily straightforward so I don't think you should have much trouble writing this.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Ameer Hamza 2020-5-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/523541-how-to-plot-the-gravitational-potential-vs-radius-of-earth-plot#answer_430784

编辑：Ameer Hamza 2020-5-6

在 MATLAB Online 中打开

You are using wrong values of constants and wrong formula for Gravitational potential. Try this code

clc;
clear all;
close all;
format long g;
G=6.6743015e-11;
Density=5510;
R = 6.4e6; % radius of earth
M = 4*pi*R^3*Density/3;
r = linspace(R, 15*R);
%% for the calculation for the gravitational potential
Gravitational_Potential = zeros(size(r));
for i =1:length(r)   
    Gravitational_Potential(i) = -G*M/r(i);
end
figure()
ax = axes();
plot(r, Gravitational_Potential);
ax.XLim(2) = max(r);
ylabel('Gravitational Potential');
xlabel('Radius of the Sphere');
title('Gravitational potential vs Radius ');
grid on