Realizing the Calculation of a Fully Normalized Associated Legendre Function Based on an FPGA

A large number of fully normalized associated Legendre function (fnALF) calculations are required to compute Earth’s gravity field elements using ultra high-order gravity field coefficient models. In the surveying and mapping industry, researchers typically rely on CPU-based systems for these calcul...

Full description

Saved in:
Bibliographic Details
Main Authors: Yuxiang Fang, Qingbin Wang, Yichao Yang
Format: Article
Language:English
Published: MDPI AG 2024-11-01
Series:Sensors
Subjects:
Online Access:https://www.mdpi.com/1424-8220/24/22/7262
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:A large number of fully normalized associated Legendre function (fnALF) calculations are required to compute Earth’s gravity field elements using ultra high-order gravity field coefficient models. In the surveying and mapping industry, researchers typically rely on CPU-based systems for these calculations, which leads to limitations in execution speed and power efficiency. Although modern CPUs improve instruction execution efficiency through instruction-level parallelism, the constraints of a shared memory architecture impose further limitations on the execution speed and power efficiency. This results in exponential increases in computation time as demand rises alongside high power consumption. In this article, we present a new computational implementation of an fnALF based on the ZYNQ platform. We design a task-parallel “pipeline” architecture which converts the original serial logic into a more efficient hardware implementation, and we utilize a redundant calculation layer to handle repetitive coefficient computations separately. The experimental results demonstrate that our system achieved accurate and rapid calculations. Under the only one geocentric residual latitude condition, we measured the computation times for spherical harmonic coefficient degrees of 360, 720, and 1080 to be 0.155922 s, 0.520950 s, and 1.401609 s, respectively. In the case of the multiple geocentric residual latitudes condition, our design generally yielded efficiency gains of over three times those of MATLAB R2020b implementation. Additionally, our calculated results were used to determine the geoid height in the field with an error of less than <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mo>±</mo><mn>0.1</mn><mspace width="0.166667em"></mspace><mi mathvariant="normal">m</mi></mrow></semantics></math></inline-formula>, confirming the reliability of our computations.
ISSN:1424-8220