Question

就像标题所说的那样，我正在尝试使用二进制搜索方法来搜索排序后的向量，以获得最接近的给定值并返回其索引。我试图使用lower / upper_bound（），但是返回的值是第一个或最后一个矢量值，或者是“ 0”。下面是我已将温度和电压读入矢量的txt文件。

1.4 1.644290    -12.5
1.5 1.642990    -13.6
1.6 1.641570    -14.8
1.7 1.640030    -16.0
1.8 1.638370    -17.1

这是我当前有效的线性搜索

double Convert::convertmVtoK(double value) const
{
    assert(!mV.empty());
    auto it = std::min_element(mV.begin(), mV.end(), [value] (double a, double b) {
        return std::abs(value - a) < std::abs(value - b);
    });
    assert(it != mV.end());
    int index = std::distance(mV.begin(), it);
    std::cout<<kelvin[index];
    return kelvin[index];
}

这是我目前正在努力提高性能的算法。

double Convert::convertmVtoK(double value)
{
    auto it = lower_bound(mV.begin(), mV.end(), value);

    if (it == mV.begin())
    {
        it = mV.begin();
    }
    else
    {
        --it;
    }

    auto jt = upper_bound(mV.begin(), mV.end(), value), out = it;

    if (it == mV.end() || jt != mV.end() && value - *it > *jt - value)
    {
        out = jt;
    }

     cout<<"This is conversion mV to K"<<" "<< *out;

任何建议将不胜感激。我认为问题可能在于向量按降序排序，但是我需要顺序保持不变以便比较值。

已解决，感谢@John。对于将来需要此功能的任何人，这里都是有效的。

double Convert::convertmVtoK(double value) const
{

    auto it = lower_bound(mV.begin(), mV.end(), value, [](double a, double b){ return a > b; });
    int index = std::distance(mV.begin(), it);
    std::cout<<kelvin[index];
}

Answer 1

由于范围不是递增的（按降序排列），因此可以使用std :: lower_bound并使用大于运算符，如注释中所述。但是，这只会使您获得的第一个结果超过或等于您的数字。这并不意味着它是“最接近的”，这就是您要的。

相反，我将使用std :: upper_bound，因此您不必检查相等性（将其加倍只是为了使其更糟），然后放回其中一个以获得另一个边界数据点，并计算哪个是实际上更近了。连同一些边界检查：

#include <vector>
#include <algorithm>
#include <iostream>
#include <functional>
#include <iterator>

// for nonincreasing range of double, find closest to value, return its index
int index_closest(std::vector<double>::iterator begin, std::vector<double>::iterator end, double value) {
    if (begin == end){
        // we're boned
        raise std::exception("index_closest has no valid index to return");
    }

    auto it = std::upper_bound(begin, end, value, std::greater<double>());

    // first member is closest
    if (begin == it)
        return 0;

    // last member is closest. end is one past that.
    if (end == it)
        return std::distance(begin, end) - 1;

    // between two, need to see which is closer
    double diff1 = abs(value - *it);
    double diff2 = abs(value - *(it-1));
    if (diff2 < diff1)
        --it;
    return std::distance(begin, it);
}

int main()
{
    std::vector<double> data{ -12.5, -13.6, -14.8, -16.0, -17.1 };
    for (double value = -12.0; value > -18.99; value = value - 1.0) {
        int index = index_closest(data.begin(), data.end(), value);
        std::cout << value << " is closest to " << data[index] << " at index " << index << std::endl;
    }
}

输出

-12 is closest to -12.5 at index 0
-13 is closest to -12.5 at index 0
-14 is closest to -13.6 at index 1
-15 is closest to -14.8 at index 2
-16 is closest to -16 at index 3
-17 is closest to -17.1 at index 4
-18 is closest to -17.1 at index 4

请注意，例如-14比-14.8更接近-13.6，这是您当前工作点的具体反例。还要注意两个端点输入的重要性。

从这里开始，您可以服用开尔文[i]。当您不需要使用外部数据数组作为函数的返回值时，我不满意，所以我只返回了索引。

Answer 2

此方法有效，但不能100％确保始终提供最接近的值。合并了@KennyOstrom方法的一部分。

    double Convert::convertmVtoK(double value) const
{

    auto it = lower_bound(mV.begin(), mV.end(), value, [](double a, double b){ return a > b; });
    int index = std::distance(mV.begin(), it);

    if(value>mV[0] || value < mV.back())
    {
        std::cout<<"Warning: Voltage Out of Range"<<"\n";
    }

    else if(value==mV[0] || value==mV.back()
            ||fabs(value - mV[index]) <= 0.0001 * fabs(value))
    {
        std::cout<<kelvin[index];
        return kelvin[index];
    }

    else
    {
        double diff1 = std::abs(value - mV[index]);
        double diff2 = std::abs(value - mV[index-1]);

        if (diff2 < diff1)
        {
            std::cout<<kelvin[index-1];
            return kelvin[index-1];
        }

        else
        {
            std::cout<<kelvin[index];
            return kelvin[index];
        }
    }

}

Answer 3

您可以使用以下命令获取具有最接近值的迭代器：

auto FindClosest(const std::vector<double>& v, double value)
{
    // assert(std::is_sorted(v.begin(), v.end(), std::greater<>{}));
    auto it = std::lower_bound(v.begin(), v.end(), value, std::greater<>{});

    if (it == v.begin()) {
        return it;
    } else if (it == v.end()) {
        return it - 1;
    } else {
        return std::abs(value - *it) < std::abs(value - *(it - 1)) ?
               it : it - 1;
    }
}

最接近值C ++的二进制搜索向量

3 个答案: