当你编辑了你的问题,志愿者通常停止回答你。这是一个共享资源:我们回答每个人的学习。如果你想私下问你的问题你可以雇个人来回答。
写一个matlab代码计算和情节的变化潜在的球形地球内部和外部。地球和它的半径已知的质量。考虑地球内部的密度均匀
27日视图(30天)
显示旧的评论
好的,我试着在代码,这是我迄今为止。当我运行我的代码,它情节图没有曲线。我应该得到一个这样的图。请帮助。
%这个代码将计算和情节的变化重力pontential里面
%和外部的球形地球
% U代表引力势
% G =引力常数(m ^ 3公斤^ 1 s ^ 2)
G = 6.67 * 10 e-11;
% M =地球的质量(千克)
M = 5.97 * 10 e24;
% r =距离地球中心质点(m)
%假设一个质点位于一个半径r = 1000000米
%绘制结果
r = 10000000;
%半径R =地球表面(m)
R = 6371000;
% p =密度(公斤/米^ 3)
%利用泊松方程我们可以得出一个公式,可以给我们
%在球形地球引力势。让U1代表
%潜在的范围内。
为x = 0:1000000: r
而r <
U1 = - (G *米)* (3 * R ^ 2 R ^ 2) / (R ^ 2 * 3);
打破
结束
结束
%,我们可以找到一个方程推导出拉普拉斯方程,这将给我们
%引力势在地球表面上的一个点。让U2
%代表地球上的潜力。
为x = 0:1000000: r
而r >
U2 = - * M (G) / r;
打破
结束;
结束;
%的潜在领域的中心是通过把r = 0:让加州大学
%代表了潜在的地球的中心。
加州大学= - (3 G * *米)/ 2 * r;
v = r / r;
情节(v v, U1, U2, v,加州大学);
包含(' r / r ');
ylabel (“U (m ^ 2 s ^ 2)”);
标题(重力势的情节作为半径的函数);
集(gca),“XAxisLocation”,“高级”);
接受的答案
沃尔特·罗伯森
2015年9月22日
看看你的代码
为x = 0:1000000: r
而r <
U1 = - (G *米)* (3 * R ^ 2 R ^ 2) / (R ^ 2 * 3);
打破
结束
结束
你的循环是一个for循环/ x,但你不使用x体内。
你最初的r >所以每次while循环立即退出。
如果你最初的r <你就写每个迭代的U1值但价值并不取决于前面的值,和你不保存U1值。效果只似你浪费时间一段时间然后跑过去迭代循环。
你为你的第二个“for”也有类似的问题。
当你结束计算v r和r, r和r只有一个值,所以你v是标量。
