Gravity Modelling by Using N-Sided Polyhedrons

This repository contains scripts and artificial gravity datasets used and explained in paper "Gravity Modelling of Sedimentary Basins by Using N-Sided Vertical Polyhedrons" by Nedim Gökhan AYDIN and Turgay İŞSEVEN. The scripts and the data are explained below. -intlog5a.m is the MATLAB adaptation of the FORTRAN scripts given by Smith (2000). -gprism.m calculates the gravity contribution of a polyhedron in space. -VoronoiIterations.m applies modelling process for given gravity anomalies, data points and parameters. -The two artificial gravity data (faulty_grav.mat and basin_grav.mat) are 500x500 matrixes for two different underground structures. *The scripts require VoronoiLimit.m function from MATLAB File Exchange to compute Voronoi diagrams for limited area (Sievers, 2023).本仓库包含了Nedim Gökhan AYDIN与Turgay İŞSEVEN在论文《利用N面垂直多面体开展沉积盆地重力模拟》（英文标题："Gravity Modelling of Sedimentary Basins by Using N-Sided Vertical Polyhedrons"）中使用并阐释的脚本与人工重力数据集。下文将对这些脚本与数据进行详细说明。 - `intlog5a.m` 为Smith（2000）所提供FORTRAN脚本的MATLAB适配版本。 - `gprism.m` 用于计算空间中任意多面体的重力贡献值。 - `VoronoiIterations.m` 可针对给定的重力异常、数据点与参数执行建模流程。 - 两份人工重力数据集`faulty_grav.mat`与`basin_grav.mat`均为500×500的矩阵，分别对应两种不同的地下构造。 * 上述脚本需依赖MATLAB文件交换库中的`VoronoiLimit.m`函数，以计算受限区域的沃罗诺伊图（Voronoi diagrams）（Sievers, 2023）。