如何从区间获得对数分布

时间:2016-02-23 18:59:18

标签: scale distribution intervals logarithm

我目前正在尝试将间隔切割成不等宽切片。实际上我希望每个切片的宽度遵循对数规则。例如,第一个间隔应该大于第二个间隔,等等。

我很难记住我的数学讲座。假设我知道 a b 分别是我的间隔 I n 的下边界和上边界是切片的数量: 如何找到每个切片的下边界和上边界(遵循对数刻度)?

换句话说,这就是我为了获得等宽间隔所做的工作:

for (i = 1; i< p; i++) {
    start = lower + i -1 + ((i-1) * size_piece);

    if (i == p-1 ) {
      end = upper;
    } else {  
      end = start + size_piece;
    }
    //function(start, end)
  }

其中: p-1 =切片数量, size_piece = | b-a |

我现在想要的是 start end 值,但遵循对数刻度而不是算术刻度(将在某些函数中调用 for 循环)。

提前感谢您的帮助。

1 个答案:

答案 0 :(得分:1)

如果我理解了您的问题,这个C ++程序将向您展示可以使用的算法的实际示例:

#include <iostream>
#include <cmath>

void my_function( double a, double b ) {
    // print out the lower and upper bounds of the slice
    std::cout << a << " -- " << b << '\n';
}

int main() {

    double start = 0.0, end = 1.0;
    int n_slices = 7;

    // I want to create 7 slices in a segment of length = end - start
    // whose extremes are logarithmically distributed:
    //     |         1       |     2    |   3  |  4 | 5 |6 |7|
    //     +-----------------+----------+------+----+---+--+-+
    //   start                                              end

    double scale = (end - start) / log(1.0 + n_slices);
    double lower_bound = start;
    for ( int i = 0; i < n_slices; ++i ) {
        // transform to the interval (1,n_slices+1):
        //     1                 2          3      4    5   6  7 8
        //     +-----------------+----------+------+----+---+--+-+
        //   start                                              end

        double upper_bound = start + log(2.0 + i) * scale;

        // use the extremes in your function
        my_function(lower_bound,upper_bound);

        // update
        lower_bound = upper_bound;
    }

    return 0;
}

输出(切片的极值)是:

0 -- 0.333333
0.333333 -- 0.528321
0.528321 -- 0.666667
0.666667 -- 0.773976
0.773976 -- 0.861654
0.861654 -- 0.935785
0.935785 -- 1