定点仿真.pdf

C 语言中的定点计算仿真 1 定点计算原理 见王灵光的《定点数运算原理》。 2 定点仿真 2.1 算法原理 2.1.1 程序流程 在 C 语言中，没有直接可以调用的函数来进行定点计算，所以需要自己编 写函数来模拟定点计算。一般来说，我们都是用整数运算来模拟定点运算的，即 C 语言中的int 类型。对于任何一个定点程序，都有下图所示的三个过程： double/float 类型转变为int 类型 主程序，里面所进行的所有关于定点的计算都通过调用函数来实现 int 类型转变为double/float 定点计算和浮点计算的差别很大，定点数的数值范围与数值精度都有严格的 要求。对于一个I +B +1位的定点系统，它的整数位宽为I ，小数位宽为B ，剩下 的一位为符号位，所表示的定点数范围从max 2I =−2−B 到min −2I 。使用定点 计算的系统中所有的数值都不能超过此范围，超过此范围的数将被限幅。而在此 范围之间的数值，其数值的最小精度也被限定在2−B ，如果小数位大于B 的话， 也要四舍五入为B 位小数，精度控制在2−B 。 2.1.2 基本运算 在C 语言中，当输入一个double/float 变量x 的时候，首先我们把它转为int 变量X ： ⎧⎢2B ⋅x +0.5⎥ if x ≥0 ⎪⎣ ⎦ X ⎨ ⎪⎩⎡⎢2B ⋅x −0.5⎤⎥ if x 0 此处采用了四舍五入的取整方法对其化整，化整之后，再对其限幅。此时有： I: 输入的整数二进制位宽 B ：输入的小数二进制位宽 max 2I +B =−1 ：系统所能表示的最大值； min =−2B +I ：系统所能表示的最小值； max 和min 都是int 类型变量。所以限幅过程如下： max if X max ⎧ ⎪ X ⎨min if X =min ⎪ ⎩X if other 此时的X 代表了x 在定点系统中的二进制表示值。 在 C 语言中，必须对所有的运算重新定义才能使用，定义的方法主要通过 调用函数来实现。有两个定点数X 和Y ，它们的和为定点数Z ，那么有 Z=X+Y max if Z max ⎧ ⎪ 且Z ⎨min if Z min ⎪ ⎩ Z if other 我们可以看到，当两个数的和Z 大于最大值max 的时候，它就被限幅到等于max ； 而当它小于最小值min 的时候，它被限幅到min 。X 、Y 本身就是定点数，所以 它们的精度为2−B ，因此其和值z 的精度也就维持在2−B 了。减法有类似的定义。 而乘除法的定义则难度稍微大些。我们仍然设两个定点数 X 和 Y ，它们的 积为定点数Z ，那么有： ⎧ × ⎢X Y ⎥ if X Y