将IEEE float转换为TI TMS320C30 32bits浮动在python中

时间:2013-07-24 12:02:45

标签: python floating-point signal-processing ieee-754 texas-instruments

我需要按照惯例将python float转换为TI DSP TMS320C30浮点表示:

http://www.ti.com/lit/an/spra400/spra400.pdf#page=13

我已经尝试了一些东西,但我似乎无法绕过提出的算法。 我还发现了C version算法,但它看起来像是在TI DSP中运行的版本,因此有些操作我无法弄清楚(例如,反转符号)。

我在下面有一个非常天真的实现,但它不起作用......

# see http://www.ti.com/lit/an/spra400/spra400.pdf

import math

def twos_comp(val, bits):
    """compute the 2's complement of int value val"""
    if( (val&(1<<(bits-1))) != 0 ):
        val = val - (1<<bits)
    return val

def float_to_ti( f ):
    m,e = math.frexp( f )
    # print m, e,
    mantissa = int(str(m)[2:])
    exponent = twos_comp((e - 127), 8) if e != 0 else (-128)
    sign = 1 if f < 0 else 0
    # print sign, mantissa, exponent
    return ((sign << 30) + (exponent << 24) + mantissa) & 0xffffffff

预期值的一些示例:

# Float     TI Decimal value of the resulting 32bits
#################################################
# 0.0         2147483648
# 1.0         0
# 100         105381888
# 0.000021    4029688104
# 10          52428800
# -1.0        4286578688
# -10.0       65011712
# -0.0021     4160118745
# -84.5487    114747153

我认为归结为python返回尾数/有效数的方式,但我不确定。
你会如何在这里开始做事?

注意:我发现this question可能是相关的,我会查看struct pack并解压缩..

仅供参考:我使用C程序检索理论值,我在DSP中加载如下:

{
  float f = -1.0; printf("F: %f -> %u", f, *(unsigned int*)&f);
  f = -10.0;      printf("F: %f -> %u", f, *(unsigned int*)&f);
  f = -0.0021;    printf("F: %f -> %u", f, *(unsigned int*)&f);
  f = -84.5487;   printf("F: %f -> %u", f, *(unsigned int*)&f);
}

工作实施

在Armin的回答之后,我让它为负数工作:

def float_to_ti(f):
    m, e = math.frexp(f)
    ma = m
    m = abs(m * 2)
    e -= 1
    sign = (math.copysign(1.0, f) < 0)
    man = int((m - 1.0) * (1 << 23))
    if sign:
        man *= -1
    if e == 0 and sign == 1:
        e = 255
    if f == 0.0:
        return (128 << 24) | (sign << 23)
    return ((e & 0xff) << 24) | (sign << 23) | man & 0x7fffff

1 个答案:

答案 0 :(得分:1)

以下代码通过了您的测试:

def float_to_ti( f ):
    m, e = math.frexp(f)
    m *= 2
    e -= 1
    sign = (math.copysign(1.0, f) < 0)
    if f == 0.0:
        return (128 << 24) | (sign << 23)
    if sign:
        m += 2.0
        if m == 1.0:
            m = 0.0
            e -= 1
    else:
        m -= 1.0
    assert 0.0 <= m < 1.0
    return ((e & 0xff) << 24) | (sign << 23) | int(m * (1 << 23) + 0.5)

注意不同的顺序(指数,符号,尾数)。另请注意,math.frexp()不返回IEEE格式的任何内容,因此此代码不会担心任何IEEE详细信息:(e & 0xff)将指数(作为有符号字符)转换为无符号数。最后,请注意C30格式不支持非正规格,这意味着它隐含了尾数的顶部位(因此m - 1.0)。