FPGA数学基础分析及与CORDIC算法计算方式对比

FPGA的最大优势之一是您能够利用其嵌入式DSP模块解决最棘手的数学传递函数。多项式近似法就是此方面的良方。

由于其灵活性与高性能，FPGA已经在众多需要计算复杂数学题或传递函数的工业、科研、军事及其它应用中找到用武之地。苛刻的精度要求与计算时延在更关键的应用中并不少见。

在采用FPGA实现数学函数时，工程师一般选择定点数学（参见：赛灵思中国通讯第80期的《FPGA数学基础》， docs/xcell80/44?e=2232228/2002872 ）。另外，您还可以采用CORDIC等许多算法计算超越函数（参见：赛灵思中国通讯第79期的《如何在FPGA中运用CORDIC算法》， publicaTIons/archives/xcel l/Xcell79.pdf )。

不过，在遇到极为复杂的数学函数时，与在FPGA之中实现精确需求函数相比，还有更高效的方法进行处理。为了理解这些变通方法 – 尤其是其中的多项式近似法，我们首先需要定义相关问题。

FPGA中负责监控铂电阻温度计（PRT）并把PRT电阻转换成温度的复杂数学传递函数就是这样一个例子。这种转换一般采用Callendar-Van Dusen方程实现。通过以下该方程的简化形式，可以确定温度介于0 ℃~660 ℃之间。

式中，R0 为0 ℃时的电阻，a 与b 是PRT的系数，而t 则是温度。

现实中，我们希望从电阻转换到温度。为此，我们需要调整该方程，确保得出的结果是给定电阻下的温度。大多数采用PRT的系统都会设计电子装置、采用电子电路测量PRT的电阻，然后利用FPGA、通过调整后的下式计算温度。

在FPGA中实现此方程即使是经验丰富的FPGA工程师也会望而却步。参考温度绘制所获得的电阻可以获得图1所示图形。从图中可以清晰看出响应的非线性。

直接在FPGA中实现调整后的传递函数会在实际所需设计工作量以及验证方面面临着巨大挑战（确保精度以及跨边界与极端条件函数）。许多工程师会想方设法实现函数，以便减少设计与验证工作量，从而控制项目进度。一个可行的方法是采用查找表保存曲线中的一系列点，同时提供LUT所含点之间的线性插值。

根据相关精度要求和保存在查找表中的元素数量，这种方法有可能满足需求。但是，您仍然需要在设计中包含线性插值函数。此函数在数学上非常复杂，而且往往包含一个非二次幂除法，其会进一步增加复杂性。

相反，您还可以利用另一种方法实现此类传递函数，那就是利用FPGA的本身特性。赛灵思 Spartan-6与7系列ArTIx、Kintex与Virtex等新型FPGA不仅仅只包含传统查找表和触发器，还具有内置DSP Slice、Block RAM、分布式RAM、PCIe®等众多高级IP硬核以及以太网端点、高速串行链路等。

由于其提供的48位累加器，工程师通常把DSP Slice称为DSP48s。不过，这些Slice还提供25 x 18位宽乘法器、加/减功能以及众多其它功能。您可以利用这些内部RAM结构和DSP Slice更轻松实现传递函数。

利用FPGA具有丰富DSP与RAM的结构的一种方法是多项式近似法。为了使用此方法，您必须首先绘出数学函数图，在MATLAB或Excel等数学程序中涵盖输入值范围。然后您可以在相关数据集中添加多项式趋势线，然后可以在FPGA中实现该趋势线的等式，以取代复杂数学函数，只要趋势线等式符合精度要求。

如果一个多项式方程无法针对整个传递函数输入范围提供足够精度，则可以添加更多方程。只要生成一系列在相关输入范围使用的多项式常数，您就能够继续依赖此方法。

能够添加多项式趋势线的数学程序大部分都允许您选择阶次或多项式项的数量。阶次越高，则配合越准确——但是您需要在FPGA中实现更多的项。在针对传递函数示例实施此过程时，我们是采用Microsoft Excel，我们已经获得了图2所示趋势线与等式。本例中采用4次多项式方程。

在获得了适合我们希望实现的传递函数的多项式之后，我们可以采用相同分析工具（在本例中为Excel）针对原始传递函数仔细检查精度。在所述监控温度的例子中，最终测量精度可能会是+/-1C，这并非苛刻的精度要求。尽管如此，根据测量范围和您实现的传递函数，可能证实仍然很难仅用一个多项式方程实现。我们该如何解决这个问题呢？