IQmath经验.docVIP

在F2812的ROM中，有3K×16位被保留用于存放数学公式表以及未来的开发。主要应用于高速度和高精度的实时计算，比同等程度的ANSIC C语言效率更高，同时可以节省用户更多的设计和调试时间。在此我对IQmath的应用做一下总结，纯属个人见解，仅供大家参考。 1???????? IQmath库 为了应用IQmath，首先要从TI官方网站下载IQmath库，文档名称为SPRC087。我们主要应用库里面的：IQmath.cmd，IQmathLib.h，IQmath.lib。新建一个工程，将IQmath.lib，IQmath.cmd添加到工程，同时在main（）函数之前增加语句：#include “IQmathLib.h”。注意：rts2800.lib和DSP281x_Headers_nonBIOS.cmd也要加到工程里面。 当然也可以不用IQmath.cmd文件，而用自己的CMD文件，只要在你的CMD里面添加以下代码即可： MEMORY { ?? PAGE 0: ?? BOOTROM (RW) : origin = 0x3ff000 , length = 0x000fc0 ?? RAMH0 (RW)?? : origin = 0x3f8000 , length = 0x002000 ? } SECTIONS { ?? IQmathTables?? : load = BOOTROM , type = NOLOAD , PAGE=0 ?? IQmath??????? : load = RAMH0 , PAGE=0??????????? } 以上代码的红色部分可以适当修改。IQmathTables段必须设置为NOLOAD型。 2???????? IQmath应用 完成以上几步之后就可以在你的主函数里应用IQmath提供的函数进行计算了。建议在应用之前把IQmathLib.h浏览一下，了解各个函数是怎样实现的。下面举一个简单的例子： #include DSP281x_Device.h#include IQmathLib.h#define ?PI? 3.14159_iq ??sinout_iq;float ?sinout_flt;void main(void){? InitSysCtrl();? InitXintf();? DINT;? IER=0X0000;? IFR=0X0000;? ? sinout_iq=_IQ29sin(_IQ29mpy(_IQ29(0.25),_IQ29(PI)));? sinout_flt=_IQ29toF(sinout_iq);? for(;;){}} ? 上述代码的功能是计算sin(π/4)的值,然后赋给sinout_flt，计算结果如图2.1所示。 图2.1? 实例计算结果 ? sinout_iq值的格式为_iq29类型（也就是long型，参见IQmathLib.h），所以要通过函数_IQ29toF(sinout_iq)转化为float类型，才是我们需要的最终结果。 ? 除个别函数外，一般情况下，计算公式里的所有变量都为一个iq类型，如上述主函数的第6行语句，全为_iq29类型。函数的具体讲解请参考IQmath手册，在此不再啰嗦。 ? 3???????? Q格式的选择 ? IQmath一共提供了30种Q格式，具体选择哪种格式要兼顾精度和值的大小依据下表而定： 例如将数5.0转为Q格式，只能从_iq1~_iq28里面选择，而不能转化为_iq29和_iq30表示，因为_iq29能转化的最大值为3.999999998，否则会发生溢出。所以在定Q格式时要对数的范围做一下估计。也正是由于这个原因，有些三角函数不能采用_iq30格式。 4???????? 计算arcsin与arccos 很多人有这样的疑问，函数表里面为什么没有提供arccos()和arcsin()函数呢？怎样才能计算这两个函数呢？其实只要你的高等数学过关的话，就很容易理解。因为arccos()和arcsin()可以通过反正切函数atan()间接求得，而函数表里面恰好提供了反正切函数，参见以下两个公式： arcsin(X) = atan(X / sqr(-X * X + 1))?? ??????????????????????反正弦 arccos(X) = atan(-X / sqr(-X * X + 1)) + 2 * atan(1)?????????????? 反余弦 arcsec(X) = atan(X / sqr(X * X - 1)) + sgn((X) - 1) * (2 * atan(1))??? 反正割 arccosec(X) = atan(X / sqr(X * X - 1