gravity_sphere：球体万有引力场数值模拟工具

资 源 简 介

该项目是一个专门用于计算和模拟球体产生的万有引力场的MATLAB程序，旨在为地球物理、天体物理以及物理教学研究提供精确的数值计算工具。程序的核心功能在于能够根据给定的球体物理参数，高精度地计算并映射球体内部及外部空间中任意点的引力场强度（重力加速度）和引力势能。 通过数值建模方法，该程序可以处理均质球体模型，也支持扩展到具有不同径向密度的分层球体模型。用户可以自由定义球体的半径、总质量或平均密度，并设定三维观测空间。程序会自动应用万有引力定律及其积分形式，生成观测点处的重力矢量数据，并通过MATLAB强大的绘图引擎将抽象的重力场转化为直观的三维等位面、矢量场分布图以及径向重力衰减曲线图。 在应用场景方面，该程序可用于模拟行星重力场、指导地下高密度矿体引发的重力异常分析，或作为基础物理中场论概念的可视化教育辅助工具，帮助用户深入理解重力场随空间距离变化的规律。

详 情 说 明

球体引力场数值计算与模拟项目

项目介绍

本项目是一个基于 MATLAB 开发的专业数值计算工具，专门用于模拟和分析球型天体产生的万有引力场。程序通过物理建模，能够精确计算均质球体内部及外部空间的重力加速度矢量分布和引力势能。该工具不仅适用于地球物理和天体物理的科研计算，也由于其强大的三维可视化功能，可作为物理场论教学的高质量辅助教材。

功能特性

1. 高精度物理建模：基于万有引力定律，实现了均质球体模型的完整引力场解算。
2. 全空间场分布覆盖：程序能够同时处理球体外部（服从平方反比定律）和球体内部（场强与距离成正比）的物理特性。
3. 多物理量解算：同步计算重力加速度（矢量）的三个空间分量以及引力势能（标量）。
4. 四维综合可视化：集成径向剖面图、三维矢量场、等位面以及空间切片图，多维度展示引力场的分布规律。
5. 数值鲁棒性：针对球心奇点进行了数值优化，确保在坐标原点处的计算不会产生空值或程序崩溃。

实现逻辑与算法说明

程序遵循严谨的物理公式进行数值实现，核心逻辑分为以下四个阶段：

1. 物理环境与参数初始化

程序首先定义了万有引力常数（G），并以地球参数为基准设定了球体的半径（约6371km）和平均密度。通过几何关系计算球体的总质量。观测空间被设定为一个 2.5 倍半径的立方体区域，并利用三维网格采样技术生成了高密度的空间坐标阵列。

2. 分段场矢量计算

利用掩码（Masking）技术，程序根据格点到球心的距离（r）将空间划分为两部分进行向量化计算：

• 球体外部 (r ≥ R)：引力加速度强度按照 $g = GM/r^2$ 计算，引力势能 $U = -GM/r$。
• 球体内部 (r < R)：引力加速度强度随半径线性增长，按照 $g = GMr/R^3$ 计算，引力势能遵循抛物线分布 $U = -GM(3R^2 - r^2)/(2R^3)$。
• 矢量方向：通过几何投影将加速度幅度分解到 X、Y、Z 三个坐标轴上，确保重力矢量始终指向球心。

3. 径向特性分析

除了空间网格计算，程序还独立生成了一条从球心向外延伸的径向采样线。通过对该采样线的迭代计算，提取出重力加速度强度和势能在穿过球面边界时的连续性变化数据。

4. 数据可视化实现

程序最后调用 MATLAB 图形引擎，通过以下四个视窗展示计算结果：

• 径向衰减曲线：使用双 Y 轴图表直观展示加速度的线性增长与平方反比衰减过程，以及势能的连续变化。
• 三维重力矢量场分布：利用三维箭头图展示重力场的指向性，并叠加透明球体模型表示物理边界。
• 典型引力势等位面：通过等值面渲染技术展示势能相等的空间包络面。
• 空间剖面强度图：在 X、Y、Z 轴的过中心平面上进行切片投影，通过伪彩色图表征重力场强度的空间衰减。

使用方法

1. 启动 MATLAB 软件。
2. 将程序源代码放置在 MATLAB 的当前工作路径下。
3. 在命令行窗口输入程序主函数名称并回车。
4. 程序将自动执行数值计算并在一个集成窗口中弹出四张分析图表。
5. 用户可以根据需要修改源代码中的半径、密度或网格分辨率（res）等参数来模拟不同的天体模型。

系统要求

• 软件环境：MATLAB R2018b 或更高版本。
• 硬件要求：建议 8GB 以上内存以支持高分辨率的三维网格运算。
• 图形支持：需要支持 OpenGL 加速的图形卡以确保三维等位面和切片图的流畅渲染。