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

15 次查看（过去 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

2 个评论
显示无隐藏无

vimal kumar chawda 2020-5-6

Thank you very much Sir:)

Ameer Hamza 2020-5-6

I am glad to be of help.

请先登录，再进行评论。

类别

MATLAB Graphics 2-D and 3-D Plots Surfaces, Volumes, and Polygons Surface and Mesh Plots

在 Help Center 和 File Exchange 中查找有关 Surface and Mesh Plots 的更多信息

产品

MATLAB

版本

R2020a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by