Таблица поиска в C с диапазонами

У меня есть таблица для расчета температуры на основе диапазона значений АЦП (аналогово-цифрового преобразователя), которые мне нужно реализовать в C. Я не уверен, как это сделать, поскольку это не будет точным значением, это будет находиться в перекрывающемся диапазоне.

Код будет принимать целочисленное значение АЦП (столбцы 2-4 означают низкий уровень АЦП, средний уровень АЦП и высокий уровень АЦП соответственно), считываемый с датчика, и выдавать температуру (столбец 1). Эти значения взяты из документа, который я преобразовал из сопротивления в напряжение в АЦП. Значение ADC будет целым числом, и оно должно наилучшим образом вписываться в диапазон, который, как я полагаю, будет означать самое близкое к центральному значению. Он не должен быть сверхточным, так как это довольно стабильное значение (обычно между 370-350 или 25-26°C), и оно будет использоваться для определения перегрева (41°C).

Вот пример нескольких ячеек таблицы:

Temperature  | ADC Low     |  ADC Center |  ADC High
-------------+-------------+-------------+------------
          25 | 362.1804923 | 372.3636364 | 382.4913871
          26 | 349.9452011 | 359.9395371 | 369.8908548
          27 | 338.1432261 | 347.9502029 | 357.7197732
          28 | 326.7557813 | 336.3737597 | 345.9668118
          29 | 315.7666703 | 325.2012277 | 334.6164426
          30 | 305.1694416 | 314.4195099 | 323.6592884
          31 | 294.9747625 | 304.0429113 | 313.1063265

Любая помощь будет оценена по достоинству.


person smd    schedule 30.11.2016    source источник
comment
Вы действительно имеете в виду, что у вас должна быть таблица поиска? Нет других альтернативных алгоритмов поиска?   -  person kaylum    schedule 30.11.2016
comment
Вы задаетесь вопросом, какую структуру данных использовать для хранения таблицы поиска или как реализовать метод поиска, который принимает текущее значение (в амперах) и дает температуру?   -  person Greg Schmit    schedule 30.11.2016
comment
Итак, если сверхточность не является проблемой, вы можете вычислить среднюю точку центральных значений АЦП между 25 и 26, между 26 и 27, … и округлить это до ближайшего целого числа. Когда входное целое число читается, вы можете найти соответствующую температуру, потому что у вас есть граничное значение, где более высокие числа равны 25ºC, затем 26ºC, затем 27ºC и т. д. Я полагаю, вы могли бы проявить фантазию и поиграть с нечеткой логикой, обрабатывая 3 числа. в виде треугольника, и посмотреть сколько на входе 25ºC и сколько 26ºC и сколько 27ºC и потом придумать составную температуру, но вряд ли оно того стоит.   -  person Jonathan Leffler    schedule 30.11.2016
comment
Таблица поиска была предложена/предпочтительна. Да, меня интересует структура данных (не чтение в файле, а размещение его в виде структуры в коде), а затем использование значения для определения температуры   -  person smd    schedule 30.11.2016


Ответы (2)


Вот некоторый рабочий код. Предполагается, что таблица сопоставления встроена; если вам нужно читать его динамически, это требует некоторого пересмотра кода, но в остальном это не является фундаментальной проблемой.

Как указано в комментарии, код определяет целочисленное значение ADC для каждую температуру таким образом, что если значение больше предела, температура является заданным значением, используя линейную интерполяцию между записями центрального значения АЦП в таблице необработанных данных АЦП. Поиск соответствующей температуры представляет собой простой линейный поиск. Вы можете использовать бинарный поиск, если хотите, или создать таблицу для правдоподобных значений ADC и сопоставить каждое значение ADC с температурой (использует больше места, но дает самый быстрый поиск). Но для диапазона 25-41ºC едва ли стоит беспокоиться о производительности линейного поиска, если только вы не можете продемонстрировать, что это серьезное узкое место, тем более, что при обычном поиске нужно будет просмотреть только 2 или 3 записи в начале. списка (поэтому обычный линейный поиск вполне может превзойти бинарный поиск!).

#include <assert.h>
#include <stdio.h>
#include <stdlib.h>
#include <time.h>

struct raw_adc
{
    int temp;
    double adc_lo;
    double adc_mid;
    double adc_hi;
};

/* Assumed sorted in order of increasing temperature */
/* Assumed monotonic in order of adc_mid values in sorted order */
/* Assumed lower temperature has higher adc_mid value */
/* Assumed ADC readings always positive (0 used as dummy) */
/* Assume contiguous series of temperatures */
static struct raw_adc raw_adc_data[] =
{
    { 25, 362.1804923, 372.3636364, 382.4913871 },
    { 26, 349.9452011, 359.9395371, 369.8908548 },
    { 27, 338.1432261, 347.9502029, 357.7197732 },
    { 28, 326.7557813, 336.3737597, 345.9668118 },
    { 29, 315.7666703, 325.2012277, 334.6164426 },
    { 30, 305.1694416, 314.4195099, 323.6592884 },
    { 31, 294.9747625, 304.0429113, 313.1063265 },
};
enum { NUM_RAW_ADC = sizeof(raw_adc_data) / sizeof(raw_adc_data[0]) };

struct map_adc
{
    int temp;
    int adc_value;
};

static struct map_adc adc_map[NUM_RAW_ADC];

static void map_raw_adc_values(void)
{
    int i;
    for (i = 0; i < NUM_RAW_ADC - 1; i++)
    {
        adc_map[i].temp = raw_adc_data[i].temp;
        /* Optionally add 0.5 before assigning */
        adc_map[i].adc_value = (raw_adc_data[i].adc_mid + raw_adc_data[i+1].adc_mid) / 2;
    }
    /* Last value is deemed to be hotter than the last recorded value */
    adc_map[i].temp = adc_map[i-1].temp + 1;
    adc_map[i].adc_value = 0;
}

static int temp_from_adc(int adc_value)
{
    int i;
    for (i = 0; i < NUM_RAW_ADC; i++)
    {
        /* Use of > determined by inspection of data - colder temps have higher ADC value */
        if (adc_value > adc_map[i].adc_value)
            return adc_map[i].temp;
    }
    assert(0);      /* Can't happen! */
    return 300;     /* If it gets here, the machine is too hot! */
}

static void dump_adc_map(void)
{
    for (int i = 0; i < NUM_RAW_ADC; i++)
    {
        assert(raw_adc_data[i].temp == adc_map[i].temp);
        printf("T = %.3dºC ADC = (%6.2f:%6.2f:%6.2f) Limit = %d\n",
                raw_adc_data[i].temp, raw_adc_data[i].adc_lo,
                raw_adc_data[i].adc_mid, raw_adc_data[i].adc_hi,
                adc_map[i].adc_value);
    }
}

int main(void)
{
    map_raw_adc_values();
    dump_adc_map();
    srand(time(0));
    for (int i = 0; i < 20; i++)
    {
        /* Range of ADC values in table is 294-382 */
        /* Generate random value in that range */
        int adc_rdg = rand() % (382 - 294) + 294;
        printf("ADC: %d = %d ºC\n", adc_rdg, temp_from_adc(adc_rdg));
    }
    return 0;
}

Пример запуска:

T = 025ºC ADC = (362.18:372.36:382.49) Limit = 366
T = 026ºC ADC = (349.95:359.94:369.89) Limit = 353
T = 027ºC ADC = (338.14:347.95:357.72) Limit = 342
T = 028ºC ADC = (326.76:336.37:345.97) Limit = 330
T = 029ºC ADC = (315.77:325.20:334.62) Limit = 319
T = 030ºC ADC = (305.17:314.42:323.66) Limit = 309
T = 031ºC ADC = (294.97:304.04:313.11) Limit = 0
ADC: 298 = 31 ºC
ADC: 358 = 26 ºC
ADC: 343 = 27 ºC
ADC: 315 = 30 ºC
ADC: 358 = 26 ºC
ADC: 352 = 27 ºC
ADC: 374 = 25 ºC
ADC: 322 = 29 ºC
ADC: 372 = 25 ºC
ADC: 376 = 25 ºC
ADC: 333 = 28 ºC
ADC: 334 = 28 ºC
ADC: 356 = 26 ºC
ADC: 307 = 31 ºC
ADC: 304 = 31 ºC
ADC: 305 = 31 ºC
ADC: 324 = 29 ºC
ADC: 358 = 26 ºC
ADC: 324 = 29 ºC
ADC: 322 = 29 ºC
person Jonathan Leffler    schedule 30.11.2016
comment
если бы таблица изменилась с -40 до 125, вы бы порекомендовали другой метод поиска скорости или вы думаете, что этого тоже будет достаточно? - person smd; 30.11.2016
comment
Это будет зависеть от того, испытываете ли вы нагрузку на ЦП или нет. В худшем случае при линейном поиске вы найдете менее 170 сравнений для температуры 125 и выше; это будет 9 сравнений для бинарного поиска (но сравните 16 и 4 для линейного и бинарного поиска в диапазоне 25-41). Так что да, я бы, вероятно, использовал бинарный поиск в гораздо большем диапазоне температур, но это сопряжено с некоторой дополнительной сложностью. Однако я сомневаюсь, что это вызовет нагрузку на любую систему. Сказав это, еще одной проблемой будет доступ к памяти — ошибки страниц и т. д. Это зависит от контекста! - person Jonathan Leffler; 30.11.2016
comment
Я бы, вероятно, не стал хранить в памяти таблицу «необработанных данных АЦП» — я бы предварительно вычислил ее. Даже если бы вам нужно было перекалибровать АЦП, я бы, вероятно, использовал рабочую программу для чтения таблицы карты АЦП из файла при запуске и по запросу, а также программу для предварительного вычисления (повторного предварительного вычисления) карты АЦП. Предположительно, такие повторные калибровки будут редкостью. - person Jonathan Leffler; 30.11.2016
comment
Также обратите внимание, что вам понадобится пользовательская функция двоичного поиска. Стандартный (bsearch()) ищет точное совпадение, но вы необходимо иметь возможность возвращать значение, когда показания АЦП не отображаются в таблице. (Я также упомяну, что последние 4 из 6 дробных цифр в необработанной таблице ADC вряд ли будут иметь значение — отсюда и выбранный формат печати. ​​У меня есть оговорки по поводу любых дробных цифр, но определенно последние 4 из них — тем более, что АЦП сообщает только целые значения.) - person Jonathan Leffler; 30.11.2016

Я думаю: используйте двоичный поиск для обхода массива структур, которые объявляют диапазон:

struct Range {
    float low;
    float center;
    float high;
    int   temperature;
};

struct RangeTable {
    struct Range* list;
    size_t length;
}

// TODO: Read the file here

struct RangeTable* table = ...
// NOTE: Ensure table.list is sorted in ascending `center` value, either in the file, or use a sorting function on the array after loading.

/////////////

/// <summary>Looks up the <paramref="adcReading" /> value in the provided table and returns a pointer to the best-matching Range struct value.</summary>
struct Range* lookup( struct RangeTable * table, float adcReading ) {

    // do a binary-search on the `center` value, which will return the Range object with the closest `center` value.
    struct Range* closestCenter = binarySearch( table->list, table->length, adcReading );

    if( closestCenter == NULL ) {
        // error condition.
        exit( 1 );
    }
    if( adcReading < closestCenter->low || adcReading > closestCenter->high ) {
        // out of range
    }
    else {
        // inside the range
    }
}

Реализация функций сортировки и бинарного поиска, конечно же, зависит от вас.

person Dai    schedule 30.11.2016