Navespro.ru

Столкнулся с проблемой, что QUIK использует формат XLTable для передачи данных по DDE.

Данные идут в виде строк или чисел с плавающей запятой - вот с последними основная проблема.
Дело в том, что из пакета данных можно достать массив типа byte (c++), в котором 8 байт отведены под ЧПЗ двойной точности.
Эти 8 байт есть IEEE-754 представление ЧПЗ двойной точности (подробно тут, а тут калькулятор).

Основная проблема, как преобразовать массив байт в число типа double.
Пришлось делать все по битам (здесь красивые картинки) - принимаю комменты по коду, т.к. особо я его не оптимизировал - работает и ладно, пошел дальше.

    //Проверка : 8 байт из посылки DDE: 0x9a99999999a09a40
    //Есть десятичное 1704.15, моя функция дает 1704.1500000000001

    //Массив байт
    byte* data = new byte[8];
   
    data[0] = 0x9a;
    data[1] = 0x99;
    data[2] = 0x99;
    data[3] = 0x99;
    data[4] = 0x99;
    data[5] = 0xa0;
    data[6] = 0x9a;
    data[7] = 0x40;

    //Буфер для перестановки байт
    //Переставляем байты в обратном порядке
    byte buf;
    buf = data[0];
    data[0] = data[7];
    data[7] = buf;
    buf = data[1];
    data[1] = data[6];
    data[6] = buf;
    buf = data[2];
    data[2] = data[5];
    data[5] = buf;
    buf = data[3];
    data[3] = data[4];
    data[4] = buf;

    //Далее представляем 8 байт как сплошные 64 бита для простоты. 63 бит слева, 0 бит справа.
    //63 бит - знак числа
    //если 0 - число положительное, 1 - отрицательное
    int sign = (data[0] & 0x80) >> 7;
    //биты с 62 по 52 - показатель степени, из расшифрованного числа нужно отнять 1023
    int exponent = (int)(((data[0] & 0x7f) << 4) + ((data[1] & 0xf0) >> 4)) - 1023;
    //2, возведенная в нужную степень
    double stepen2 = (int)pow(2.0,exponent);
    //Заранее считаю множитель дробной части, чтобы сократить погрешность
    double footer = stepen2 / 4503599627370496;
    //биты с 51 по 0 - мантисса
   
    //Формула расчета:
    //(-1)^sign * 2^(exponent - 1023) * (1 + mantisse/2^52)
   
    //далее я последнюю скобку умножил на 2^(exponent - 1023),
    //чтобы не пришлось делить на 2^52, а меньше.
    double mantisse =  stepen2  + (
        //51 -48
        ((double)(data[1] & 0x0f) * (double)281474976710656) +
        ((double)data[2] * (double)1099511627776) +
        ((double)data[3] * (double)4294967296) +
        ((double)data[4] * (double)16777216) +
        ((double)data[5] * (double)65536) +
        ((double)data[6] * (double)256) +
        ((double)data[7] )) * footer;

    //Итоговое значение со знаком
    double value = (sign == 0 ? 1 : -1) * mantisse;

Вам это будет интересно!

Последние новости